微分方程式のプログラムについて

質問

次の問題を解こうとしましたが、プログラムが動かなかったので実行できませんでした。今、非常に困っています。分かりにくいプログラムですみません。

問　温度Taの空気中に置かれている物体の温度Tは、次の常微分方程式に従う。
　　　dT/dt=-k(T-Ta)
t=0でT=100[℃]、Ta=20[℃]、k=1×10^-2[s^-1]として、物体の温度が60℃にまで冷却される時間を求めよ。厳密解とも比較すること。

自分が作ったプログラムは
#include <stdio.h>
#include <math.h>

double k = 1e-2, tmpa=20.0;
double f(double t, double tmp)
{
return -k*(tmp-tmpa);
}

int main(void)
{
double t0 = 0.0, tmp0 = 100, t1, tmp1, ttmp1;
double k1, k2, k3, k4, h;

printf("刻み幅=");
scanf("%lf", h);
t1 = t0+h;
ttmp1=tmpa+(tmp0-tmpa)+exp(-k*(t0-t0));
printf("%8.4lf %8.4lf %8.4lf\n", t0, tmp0, ttmp1);
while(t0 <= 100.0){
k1 = f(t0, tmp0);
k2 = f(t0+0.5*h, tmp0+0.5*h*k1);
k3 = f(t0+0.5*h, tmp0+0.5*h*k2);
k4 = f(t0+h, tmp0+h*k3);
tmp1 = tmp0+h*(k1+2.0*k2+2.0*k3+k4)/6.0;
printf("%8.4lf %8.4lf %8.4lf\n", t1, tmp1, ttmp1);
if(tmp1 <= 60.0) break;
t0 = t1;
tmp0 = tmp1;
t1 = t0+h;
}
return 0;
}

あと厳密解はT=e^-kt(T0-Ta)+Ta求められるそうです。
変な間違いが多いと思いますが、このプログラムのどこが間違っているか教えてください。厳密解だけでもいいです。お願いします。

回答No.1

kochory · ベストアンサー

「プログラムが動かなかった」というのがどういう意味なのかはっきりさせなさい。
コンパイルできないのか、できても実行できないのか、実行できても動作がおかしいのか。
エラーメッセージが出る場合はそれも質問で示すこと。
以上はプログラムに関する質問をするときの常識。

プログラム上おかしいところを言うと、scanfの引数をhでなく&hにすること。
あと関数f(t,tmp)内で変数tはまったく使っていないので、渡す必要がない。

お礼

ありがとうございました。何とかプログラムは動きました。

補足

すみません。実行はできました。そのあと忘れましたが英語のメッセージが出てきました。書き忘れましたがこのプログラムはルンゲクッタ法です。