Гармонические колебания пружин

>	restart:with(plots):with(plottools): N:=20:

>	pru:=proc(x1,x2,y)local K,d,r,l: r:=6: K:=30: l:=(x2-x1)/3: d:=l/K:  PLOT(CURVES([[x1,y],seq([x1+l+j*d,r*sin(Pi*j/2)+y],j= 0..K),[x2,y]])):

>

end:

Warning, the name changecoords has been redefined

Warning, the name arrow has been redefined

>	for t to N do

>	m:=20*evalf(2+sin(4*Pi*t/N));

>	H[t]:=pru(0,m,0):od:

>	for t to N do

>	m1:=30*evalf(2+sin(2*Pi*t/N));

>	H2[t]:=pru(0,m1,15):od:

>	a:=display(seq(H[t],t=1..N),insequence=true):

>	b:=display(seq(H2[t],t=1..N),insequence=true):

>	display(a,b,axes=none,scaling=constrained);

>