function dYdt = spring_eq(t , Y)
global N1 N2 N3
dYdt = [Y(2) ; -N1^2*Y(1) ; Y(4) ; -N2^2*Y(3) ; Y(6) ; -N3^2*Y(5)];