求助:这到底是什么错误
新手上路,向DX请教
#include"stdio.h"
#include"math.h"
#include"iostream.h"
#include"fstream.h"
#define M 5
#define N 5
#define x 0.01
#define y 0.01
#define p 370
#define cp 2.5
#define k 6.5e-4
#define hlv 2500
#define r 0.3
#define s 0.02
#define e 0.3
#define dm 6e-9
#define tin 383
#define tout 383
#define uin 0.04
#define uout 0.04
#define tb 283
#define ub 0.86
void main()
{
double t[M+1][N+1],u[M+1][N+1],errx=0.0,erry=0.0,a,b,c,d,A,B,C,D;
double ain=2.25,aout=2.25,bin=2.5e-6,bout=2.5e-6;
int i,j;
int step=0; //迭代次数
ofstream o_file("result.dat"); //输出dat文件
for(j=0;j<=N;j++)
for(i=0;i<=M;i++)
{ t[i][j]=0; //初始化
u[i][j]=0;
}
for(i=0;i<=M;i++) //初始条件
{
t[i][0]=tb;
u[i][0]=ub;
}
for(j=0;j<N;j++) //边界条件
{
t[0][j+1]=2*k*y*(t[1][j]-t[0][j])/(p*cp*x*x)+t[0][j]+r*hlv*(u[0][j+1]-u[0][j])/cp-2*y*(ain*(t[0][j]-tin)+(1-e)*hlv*bin*(u[0][j]-uin))/(p*cp*x);
u[0][j+1]=2*dm*y*(u[1][j]-u[0][j]+s*(t[1][j]-t[0][j]))/(x*x)+u[0][j]-2*y*bin*(u[0][j]-uin)/x;
}
for(j=0;j<=N;j++)
{
t[M][j+1]=2*k*y*(t[M-1][j]-t[M][j])/(p*cp*x*x)+t[M][j]+r*hlv*(u[M][j+1]-u[M][j])/cp-2*y*(aout*(t[M][j]-tout)+(1-e)*hlv*bout*(u[M][j]-uout))/(p*cp*x);
u[M][j+1]=2*dm*y*(u[M-1][j]-u[M][j]+s*(t[M-1][j]-t[M][j]))+u[M][j]-2*bout*y*(u[M][j]-uout)/x;
do
{ errx=0;
erry=0;
for(j=0;j<N;j++)
for(i=1;i<M;i++)
{ a=t[i][j+1];
b=u[i][j+1];
A=k/(p*cp);
B=r*hlv/cp;
C=k*dm/(k+dm*s*p*r*hlv);
D=dm*s*p*cp/(k+dm*s*p*r*hlv);
t[i][j+1]=y*A*(t[i+1][j]+t[i-1][j])/(x*x)+(x*x-2*y*A)*t[i][j]/(x*x)+B*(u[i][j+1]-u[i][j]);
u[i][j+1]=C*y*(u[i+1][j]+u[i-1][j])/(x*x)+(x*x-2*y*C)*u[i][j]/(x*x)+D*(t[i][j+1]-t[i][j]);
c=fabs(t[i][j+1]-a);
d=fabs(u[i][j+1]-b);
errx=errx+c;
erry=erry+d;
}
step++;
}while(errx+erry>=0.0001); //误差值设为0.0001
o_file<<"VARIABLES="<<"\"t\","<<"\"u\","<<endl;
o_file<<"ZONE "<<"I="<<M+1<<", J="<<N+1<<", F=POINT"<<endl;
for(j=0;j<=N;j++)
for(i=0;i<=M;i++)
{
cout<<"t"<<i<<j<<"="<<t[i][j]<<" u"<<i<<j<<"="<<u[i][j]<<endl;
o_file<<t[i][j]<<" "<<u[i][j]<<endl;
}
printf("step=%d\n",step);
}
Compiling...
1.cpp
d:\program files\microsoft visual studio\vc98\include\math.h(92) : error C2059: syntax error : 'constant'
d:\program files\microsoft visual studio\vc98\include\math.h(92) : error C2238: unexpected token(s) preceding ';'
f:\论文相关\tdma\1.cpp(100) : fatal error C1004: unexpected end of file found
Error executing cl.exe.
1.obj - 3 error(s), 0 warning(s)