fortranのプログラム

質問

fortranのプログラム

現在、fortranの勉強をしております。
そこで、質問があるのですが、
ある関数ｆ(x,y,z)の座標（x,y,z）の値がデータとして与えられているとき、

S=10+f

をfortranで計算したいと考えております。
ただ、関数ｆは複数（ｆ１、ｆ２、ｆ３）あり、次々とｆに代入してＳを計算したいのですが、どのようにプログラムしたらいいか思いつきません。
どなたか、ヒントだけでもいいので、教えてください。

ちなみに、私が考えたプログラムは(下のプログラムはポイントだけ書いてあります。endやその他関係ないと思われるところは省いております。)

do 100 k=1,3
S=S+fk(x,y,z)
continue

function f1
f1(s,t,u)=・・・
return
end

f2(s,t,u)=・・・
return
end

使用しているバージョンは、fortran77(本当は90を使っているのですが、77だけで書いています)です。

回答No.1

f272 · ベストアンサー

私が77でやるなら

program main
S=0
x=1
y=2
z=3
do 100 k=1,3
S=S+f(k,x,y,z)
100 continue
print*,S
end

function f(k,s,t,u)
if (k.eq.1) then
f=s+t+u
else if (k.eq.2) then
f=s+t+u*2
else if (k.eq.3) then
f=s+t+u*3
endif
return
end

こんな感じ。

お礼

早速のご回答ありがとうございます。
このプログラミングを参考に書いてみます。