偏微分のプログラムについてです.

質問

(∂^２ u)/(∂x^2 )+(∂^2 u)/(∂y^2 )=2(x^2+y^2 )
境界上でu(0,y)=u(x,0)=0,u(1,y)=y^2,u(x,1)=x^2である.点(0.5,0.5)の値を求めよ.その際きざみ幅を幾つか変化させることで,反復回数がどのように変化するか確かめよ.
という問題です.自分で以下のプログラムを作ったのですが,コンパイルされません.ぜひどこがおかしいか教えてください.

#include<stdio.h>
#define EPS 10e-6
#define N 100
#define M 100

main(void){
int i,j;
int n=19,m=19,r=0;
double h,x0,S;
double u1[N][M],u2[N][M];

h=2.0/(n+1);
for(i=0;i<n+2;i++){
for(j=0;j<m+2;j++){
u1[i][j]=0.0;
}
}

do{
for(i=0;i<n+2;i++){
u2[i][0]=0.0;
u2[i][m+1]=0.0;
}

for(j=0;j<m+2;j++){
u2[0][j]=0.0;
u2[n+1][j]=0.0;
}

for(i=1;i<n+1;i++){
for(j=1;j<m+1;j++){
u2[i][j]=(u1[i-1][j]+u1[i][j+1]+u1[i+1][j]+u1[i][j-1]+2.0*h*h)/4.0;
}
}

r++;

S=fabs(u1[(n+1)/2][(m+1)/2]-u2[(m+1)/2][(m+1)/2]);
for(i=0;i<n+2;i++){
for(j=1;j<m+1;j++){
u1[i][j]=u2[i][j];
}
}
}
while(S>EPS);
x0=u2[(n+2)/2][(m+2)/2];

printf("きざみ幅=%lf\n",h);
printf("反復回数=%d\n",r);
printf("点(0.5,0.5)=%lf\n",x0);
}

回答No.1

asuncion · ベストアンサー

> コンパイルされません.

コンパイル時のエラーメッセージを提示してください。
ちなみに、当方では警告は出ましたがエラーは出ませんでした。

> while(S>EPS);

SがEPSより大きい間、単にループするだけとなっていますが、
意図通りですか？

補足

ご返答ありがとうございます．
　このプログラムはc080725.cで保存しているのですが，まず cl c080725.c と入力した時点ではエラーはでないのですが，いざ次に c080725で実行すると，普通printfのところが出てきて入力できるとおもうのですが，出てきません．Crt c で画面を戻しています．
　また，while(S>EPS)のところは意図的にではありません．どうしたらいいのでしょうか？

回答No.2

asuncion

> while(S>EPS)のところは意図的にではありません．

失礼いたしました。当方の誤りです。
do ～ while 構文ですので、
当然
while(S>EPS);
です。申し訳ありません。

fabs() のために math.h をインクルードしたり、
main関数の戻り値が int 型なので最後に return 0; を加えたりして
コンパイル時に警告もエラーも出ないようにすると、当方では

きざみ幅=0.100000
反復回数=542
点(0.5,0.5)=0.587418

という結果を得ました。正しいかどうかはわかりません。

お礼

言って頂いた通りにすると,きちんとコンパイルできました.知識のない私に丁寧に教えてくださり,本当に感謝します.ありがとうございました.