! 線形加速度法 real a(101), v(101),x(101),as(101) dt=0.1 time=0. a(1)=0. v(1)=0. x(1)=0. gal=1 DO 100 i=1,100 x(i+1)=x(i)+v(i)*dt+(2*a(i)+1)*dt*dt/6. v(i+1)=v(i)+(a(i)+1.)*dt/2. a(i+1)=-1*gal-2*v(i+1)-x(i+1) DO 200 J=1,100 as(j)=a(i+1) x(i+1)=x(i)+v(i)*dt+(2.*a(i)+as(j))*dt*dt/6. v(i+1)=v(i)+(a(i)+as(j))*dt/2. a(i+1)=-1*gal-2.*v(i+1)-x(i+1) 200 CONTINUE time=time+dt write(6,*) a(i+1), time 100 CONTINUE END