Re: Превод с фортрана на СИ
Добавлено: 03 ноя 2009, 16:18
void main()
{
double e[50][1000], f[50][1000], t[50][50][1000],b1,b2,b,b3,b4,b5,sumtau,a1,a2,t1,l122,t2,l111,l112,l121,l211,l222,l221,l212,c,tsr,tw,rrx,rry,rr,tau,r,dy,dx,df,deltax,deltaf1,deltay,deltaf2;
double x1,y1,tl,tlfx,tlfy,nx,mx,ny,my,nnx,nny,ny1,nx1,kfs,j,k,i,bn,rf,mws,kj,mks,ts;
//ввод данных
//выбор шага сетки по Х
a1=4.95E-6;
a2=5.78E-6;
deltax=x1/1000;
deltaf1=tlfx/1000;
dx=sqrt(6*a1*tau);
df=sqrt(4*a2*tau);
nx=int((deltax/dx)+0.5)+1;
mx=int((deltaf1/df)+0.5)+nx+1;
nnx=mx-nx;
printf ("Число узлов в металле по Х=",nx);
printf ("Число узлов в форме по Х=",nnx);
//выбор шага сетки по У
deltaf2=tlfy/1000;
deltay=y1/1000;
dy=sqrt(4*a1*tau);
ny=int((deltay/dy)+0.5)+1;
my=int((deltaf2/df)+0.5)+ny+1;
nny=my-ny;
printf ("Число узлов в металле по Y=",ny);
printf ("Число узлов в форме по Y=",nny);
//задание начальных условий в отливке
k=1;
do
{ j=1;
do
{ i=1;
t[j][k]=t1;
}
while(i=nx);
}
while(i=ny);
{
double e[50][1000], f[50][1000], t[50][50][1000],b1,b2,b,b3,b4,b5,sumtau,a1,a2,t1,l122,t2,l111,l112,l121,l211,l222,l221,l212,c,tsr,tw,rrx,rry,rr,tau,r,dy,dx,df,deltax,deltaf1,deltay,deltaf2;
double x1,y1,tl,tlfx,tlfy,nx,mx,ny,my,nnx,nny,ny1,nx1,kfs,j,k,i,bn,rf,mws,kj,mks,ts;
//ввод данных
//выбор шага сетки по Х
a1=4.95E-6;
a2=5.78E-6;
deltax=x1/1000;
deltaf1=tlfx/1000;
dx=sqrt(6*a1*tau);
df=sqrt(4*a2*tau);
nx=int((deltax/dx)+0.5)+1;
mx=int((deltaf1/df)+0.5)+nx+1;
nnx=mx-nx;
printf ("Число узлов в металле по Х=",nx);
printf ("Число узлов в форме по Х=",nnx);
//выбор шага сетки по У
deltaf2=tlfy/1000;
deltay=y1/1000;
dy=sqrt(4*a1*tau);
ny=int((deltay/dy)+0.5)+1;
my=int((deltaf2/df)+0.5)+ny+1;
nny=my-ny;
printf ("Число узлов в металле по Y=",ny);
printf ("Число узлов в форме по Y=",nny);
//задание начальных условий в отливке
k=1;
do
{ j=1;
do
{ i=1;
t[j][k]=t1;
}
while(i=nx);
}
while(i=ny);