09/12/03 21:10:12
fortranの初心者です
彗星の軌道をプログラミングしたいのですが上手くいきません
どこがおかしいか教えてくださいませんか？
program main
integer i,j
real ms,me,x,y,vx,vy,t,dt,G
ms=1.989*(10**30)
dt=100
G=6.673*(10**(-11))
K=G*ms/R**2
x=0.8826*(10**(11))
y=0.0
R=sqrt(x**2+y**2)
vx=0.0
vy=55000
do i=0,1000
x=x+vx*dt*i
y=y+vy*dt*i
R=sqrt(x**2+y**2)
vx=vx-K*(x/R)*dt*i
vy=vy-K*(y/R)*dt*i
write(6,*) x,y
enddo
end