Visual C++での数値計算のプログラミング
質問初めてになります。
プログラミングにあまり詳しくない大学院の数理科の学生です。
学校で熱方程式の陽解法のプログラミングのレポートが出されたのですが、困っています。
レポートで詰まっている点は
windows Visual C++ Expressionでの陽解法のプログラミングです。
indows Visual C++ Expression2008が良くわからないのでコマンドプロンプトで実行しています。
熱方程式の初期値問題の陽解法のプログラムをなんとか組んでいます。
私が組んだプログラムではLinuxでは通るのですがwindows Visual C++ Expressionでは通りません。
このプログラムをwindows Visual C++ Expression 2008で通すにはどのように直したがいいのか教えて頂きたいと思います。よろしくお願いします。
以下は私なりに組んだプログラミングです。
πの値、u1[N+1]、u2[N+2]
の3箇所でエラーがでてると思われます。
#include<stdio.h>
#include<stdio.h>
#include<math.h>
int main(){
int i,k,kmax;
int N=10;
double dx
double dt=0.01;
double u1[N+1];
double u2[N+1];
double r=dt/(dx*dx);
double T=1.0;
kmax=T/dt;
u1[0]=0;
u1[N]=0;
for(i=1;i<=N;i++){
u1[i]=sin(M_PI*i*dx);
}
for(k=1;k<=kmax;k++){
for(i=1;i<=N-1;i++){
u2[i]=r*u1[i-1]+(1-2*r)*u1[i]+r*u1[i+1];
}
for(i=1;i<=N-1;i++){
u1[i]=u2[i];
}
for(i=0;i<=N;i++){
printf("%f %f %f\n",dt*k,dx*i,u1[i]);
}
printf("\n");
}
return 0;
}
お礼
購読していませんでした。 やってみます。 情報ありがとうございます。