编了个程序，编译链接都没问题，最后运行怎么就出现这种问题了，哪位高手有空给俺看看辣!!多谢了
程序如下：
#include <stdio.h>
#include <math.h>
double abstr(double);
int sign(double);
void main(void)
{
double m[500],c[500],p[500],d[500],ecn[500];
double mom[100],coc[100];
int i;
for(i=0;i<500;i++)
{
m[i]=0.0;
c[i]=0.0;
p[i]=0.0;
d[i]=0.0;
ecn[i]=0.0;
}
for(i=0;i<100;i++)
{
mom[i]=0.0;
coc[i]=0.0;
}
//**************Enter Data To Store In Input.dat******************
FILE *file1,*file2,*file3;
double fy,es,esh;
double fc,fct;
double as1,as;
double l,a,h,b;
double ao[2];
int sn,ln;
file1=fopen("input.dat","r");
fscanf(file1,"%f%f%f",&fy,&es,&esh);
fscanf(file1,"%f%f",&fc,&fct);
fscanf(file1,"%f%f",&as1,&as);
fscanf(file1,"%f%f%f%f",&l,&a,&h,&b);
fscanf(file1,"%f%f",&ao[0],&ao[1]);
fscanf(file1,"%d%d",&sn,&ln);
fclose(file1);
//***************end of inputing data**********************
double decn=0.000002,de=0.00005,ee,em;
double sf1=0.0,sf2=0.0,dsf;
double ffc=0.0,fs=0.0,fps=0.0;
double mi=0.0,mic,mis,mip;
double z,e,s,r;
double eo,eu,cc=0.0;
double lp,hh,hn,aas,etop;
double esy,da;
double eps,ep,sps,epo,ep1,fpu,aps;
double dd,dsn,dl;
double mm,mo,co,dp,tan;
int j,k,n,ii,jj,t;
int jmax1=500,jmax2=0,jmax3=0;
esy=fy/es;
ep1=195000;epo=214000;fpu=1860;aps=1668;ep=275;t=4;
hn=h/sn;
eu=-0.0035;
eo=-0.002;
ee=-0.0001;
j=0;ii=0;
ecn[0]=-0.00013;
do
{
j++;
ecn[j]=ecn[j-1]+decn;
do
{
ii++;
mic=0.0;
mis=0.0;
c[j]=2.0*(ecn[j]-ee)/h;
//**************calculate the force of concrete***********
ffc=0.0;
for(i=0;i<sn;i++)
{
z=h/2-i*hn-hn/2;
e=ee-z*c[j];
if(e>0.0)s=0.0;
else if(e>eo)s=-fc*(2*e/eo-(e/eo)*(e/eo));
else if(e>eu)s=-fc;
ffc=ffc+s*b*hn;
mic=mic+s*b*z*hn;
}
//*****************calculate the force of steel**************
fs=0.0;
for(k=0;k<2;k++)
{
z=ao[k]-h/2.0;
e=ee-z*c[j];
if(abstr(e)<=esy)s=e*es;
else s=sign(e)*fy;
if(z<0.0)aas=as1;
else aas=as;
fs=fs+s*aas;
mis=mis+s*aas*z;
}
//*****************calculate the force of precost steel**********
fps=0.0;
eps=ee+ep*c[j]+0.00519;
if(eps<0.00668)sps=ep1*eps;
else sps=(epo*eps)/pow(1+pow(epo*eps/fpu,t),1.0/t);
fps=sps*aps;
mip=fps*ep;
fs=fs+fps;
mis=mis+mip;
if(ii==1)
{
sf1=ffc+fs;
ee=ee+de;
em=de;
}
else
{
sf2=ffc+fs;
dsf=sf2-sf1;
em=-sf2*em/dsf;
if(em==0.0)em=de;
ee=ee+em;
sf1=sf2;
}
mm=mic+mis;
}while(abstr(sf2)>=0.1);
if(e>esy)
{
if(j<=jmax1)jmax1=j;
}
etop=ee-(h/2.0-hn/2.0)*c[j];
m[j]=mm;
if(m[j]>=mi)
{
mi=m[j];
jmax2=j;
}
}while(etop<=1.2*etop||fps<=fpu);
jmax3=j;
//*************put out the answer of moment and curvature**********
file2=fopen("out1.dat","w");
for(j=0;j<jmax3;j++)
{
fprintf(file2,"%18.8g,%18.8g \n",c[j],m[j]);
}
fclose(file2);
}
int sign(double num)
{
int q;
if(num<0.0)q=-1;
if(num>0.0)q=1;
return(q);
}
double abstr(double x)
{
if(x>=0)x=x;
else x=-x;
return x;
}
