fortran77 自宅ではコンパイルできないので

質問

課題が出ているのですが、
自宅のパソコンではコンパイルできないので、
このプログラムで動くか教えてほしいです

dx/dt=1+x^2
初期条件x=0,t=0
を０＜ｔ＜１．５まで
１０等分で解く

オイラー法
f(x,t)=1+x**2
x=0,t=0
a=1.5
b=10
h=a/b

do i=0,a,h
x=x+f(x,t)*h
write(*,*) x,t
end do

２次ルンゲクッタ法
do i=0,a,h
x=x+h{(1-m)*f(x,t)+mf(x+h/(2m),t++h/(2m)*f(x,y))}
write(*,*)x,t
end do

４次ルンゲクッタ法
do i=0, a,h
k1=h*f(x,t)
k2=h*f(x+h/2,t+k1/2)
k3=h*f(x+h/2,t+k2/2)
k4=h*f(x+h/2,t+k3/2)
x=x+(k1+2*k2+2*k3+k4)/6
write(*,*) x,t
end do

これでx(t)=tantであることが求められるでしょうか？？
間違っている箇所を指摘していただけると嬉しいです。

あと２階常微分でのルンゲクッタ法というのが全く理解できません。
どなたか詳しく教えてください。

お願いします。

回答No.1

noname#129118 · ベストアンサー

Vectorに無料コンパイラがありませんか。なければ、basicインタープリタでもOKですから、実行してみると良いのですが、、、
(fortranとbasicは言語仕様が似ているので便利かと思います)

●doループは整数刻みで、期待する値とは異なるはず…即ちバグ
●その他は文法エラーはなさそうです

回答No.2

Tacosan

そりゃぁもともと BASIC は「TSS で動いていた FORTRAN の簡略版」ではあるけど, さすがに今「fortranとbasicは言語仕様が似ているので便利かと思います」はないと思うよ＞#1. フリーなコンパイラだって, GCC を含めいくつか存在するんだし.

とりあえず「間違っている箇所を指摘していただけると嬉しいです」というなら完全なプログラムを出してほしい. まさか, これを「FORTRAN のプログラム」というつもりはないよね?

あと, 最後の「２階常微分でのルンゲクッタ法」は基本的に「高階常微分方程式は連立 1階常微分方程式に直せる」ことを使っただけ. もちろん式の形によっては (ちょっとだけ) 精度がいい形にもなるけど.

お礼

ありがとうございます