フォートランで正しく表示されません。間違いを教えてください原子(He)を飛ばして、原点にある原子(Si)に影響されて軌道を変えて飛んでいくグラフを書きたいのですが、エラーが表示されてしまいます。 implicit real*8 (a-h,o-z) real*8 M1,M2,x(1001),y(1001),vx(1001),vy(1001) open(unit=10,file 'a.dat',status='unknown') Z1=2.0 M1=4.0 Z2=14.0 M2=28.0 u=M1*M2/(M1+M2) x0=-5E-8 y0=0.1E-8 vx=3.0E10*0.07 vy=0.0 dt=0.0001 do i=1,1000 x(i+1)=x(i)+vx(i)*dt y(i+1)=y(i)+vy(i)*dt vx(i+1)=vx(i)*fx(x(i),y(i))*dt/u (*) vy(i+1)=vy(i)*fy(x(i),y(i))*dt/u enddo close(10) stop endfunction fx(x,y) (*)function fy(x,y) (*) e2=14.4*1.6E-12 fx=Z1*Z2*e2/x(i) fy=Z1*Z2*e2/y(i) return end----------------------------------------------(*)の部分エラーがでます。わかる方教えてくださいm(_ _)m
![]()