수다닷컴

소스코드

솔직히 어떤 오류인지는 모르겠어요(이제 입문단계입니다..)


오류가 떴는데 그 다음에 어떻게 해야할지를 모르겠습니다.#include stdio.h
#include math.hvoid main()
{

file *out;int i, j, n;
double u[22][6];
double ru[22][6];
double r[22][6];

u[0][0]=4.00; u[0][1]=4.00;u[0][2]=4.00; u[0][3]=4.00;u[0][4]=4.00;u[0][5]=4.00;u[0][6]=4.00;u[0][7]=4.00;u[0][8]=4.00;u[0][9]=4.00;u[0][10]=4.00;
u[0][11]=0.00;u[0][12]=0.00;u[0][13]=0.00;u[0][14]=0.00;u[0][15]=0.00;u[0][16]=0.00;u[0][17]=0.00;u[0][18]=0.00;u[0][19]=0.00;u[0][20]=0.00;u[0][21]=0.00;
u[1][0]=3.90; u[1][1]=3.88;u[1][2]=3.86; u[1][3]=3.84;u[1][4]=3.80;u[1][5]=3.75;u[1][6]=3.70;u[1][7]=3.62;u[1][8]=3.54;u[1][9]=4.00;u[1][10]=3.40;
u[1][11]=0.80;u[1][12]=0.70;u[1][13]=0.60;u[1][14]=0.50;u[1][15]=0.45;u[1][16]=0.36;u[1][17]=0.27;u[1][18]=0.22;u[1][19]=0.14;u[1][20]=0.08;u[1][21]=0.06;
u[2][0]=3.85; u[2][1]=3.79;u[2][2]=3.66; u[2][3]=3.56;u[2][4]=3.46;u[2][5]=3.39;u[2][6]=3.24;u[2][7]=3.01;u[2][8]=2.76;u[2][9]=2.20;u[2][10]=2.00;
u[2][11]=2.00;u[2][12]=1.69;u[2][13]=1.60;u[2][14]=1.30;u[2][15]=1.10;u[2][16]=0.89;u[2][17]=0.68;u[2][18]=0.49;u[2][19]=0.38;u[2][20]=0.20;u[2][21]=0.15;
u[3][0]=3.80; u[3][1]=3.68;u[3][2]=3.45; u[3][3]=3.38;u[3][4]=3.24;u[3][5]=3.02;u[3][6]=2.89;u[3][7]=2.76;u[3][8]=2.45;u[3][9]=2.22;u[3][10]=2.00;
u[3][11]=2.00;u[3][12]=1.72;u[3][13]=1.66;u[3][14]=1.34;u[3][15]=1.20;u[3][16]=1.01;u[3][17]=0.86;u[3][18]=0.65;u[3][19]=0.43;u[3][20]=0.22;u[3][21]=0.18;
u[4][0]=3.75; u[4][1]=3.59;u[4][2]=3.40; u[4][3]=3.28;u[4][4]=3.16;u[4][5]=2.98;u[4][6]=2.87;u[4][7]=2.72;u[4][8]=2.41;u[4][9]=2.16;u[4][10]=2.00;
u[4][11]=2.00;u[4][12]=1.78;u[4][13]=1.70;u[4][14]=1.46;u[4][15]=1.26;u[4][16]=1.06;u[4][17]=0.88;u[4][18]=0.66;u[4][19]=0.46;u[4][20]=0.24;u[4][21]=0.22;
u[5][0]=3.70; u[5][1]=3.56;u[5][2]=3.38; u[5][3]=3.20;u[5][4]=3.12;u[5][5]=2.96;u[5][6]=2.80;u[5][7]=2.66;u[5][8]=2.36;u[5][9]=2.10;u[5][10]=2.00;
u[5][11]=2.00;u[5][12]=1.84;u[5][13]=1.77;u[5][14]=1.52;u[5][15]=1.30;u[5][16]=1.10;u[5][17]=0.90;u[5][18]=0.69;u[5][19]=0.48;u[5][20]=0.30;u[5][21]=0.26;out = fopen(output.dat , w);/* boundary condition */for(i=0; i11; i++){
ru[0][i]=4.00;
}
for(i=11; i21; i++){
ru[0][i]=0.00;
}
for(j=2; j6; j++){
ru[j][10]=2.00;
ru[j][11]=2.00;
}
/* main loof */r[5][0]=0.5*(u[4][0]+u[5][1]);
r[5][21]=0.5*(u[4][21]+u[5][20];for(j=1; j5; j++){
for(i=1; i10; i++){
ru[j][i]=0.25*(u[j][i-1]+u[j][i+1]+u[j-1][i]+u[j+1][i]);
r[j][i]=0.25*(u[j][i-1]+u[j][i+1]+u[j-1][i]+u[j+1][i])-ru[j][i];
}
}for(j=1; j5; j++){
for(i=12; i21; i++){
ru[j][i]=0.25*(u[j][i-1]+u[j][i+1]+u[j-1][i]+u[j+1][i]);
r[j][i]=0.25*(u[j][i-1]+u[j][i+1]+u[j-1][i]+u[j+1][i])-ru[j][i];
}
}
for(j=1; j5; j++){
ru[j][0]=0.25*(2.0*u[j][1]+u[j-1][0]+u[j+1][0]);
r[j][0]=0.25*(2.0*u[j][1]+u[j-1][0]+u[j+1][0])-ru[j][0];
}
for(j=1; j5; j++){
ru[j][21]=0.25*(2.0*u[j][20]+u[j-1][21]+u[j+1][21]);
r[j][21]=0.25*(2.0*u[j][20]+u[j-1][21]+u[j+1][21])-ru[j][21];
}
for(j=1; j2; j++){
ru[j][10]=0.25*(2.0*u[j][9]+u[j-1][10]+u[j+1][10]);
r[j][10]=0.25*(2.0*u[j][9]+u[j-1][10]+u[j+1][10])-ru[j][10];
}
for(i=1; i10; i++){
ru[5][i]=0.25*(2.0*u[4][i]+u[5][i-1]+u[5][i+1]);
r[5][i]=0.25*(2.0*u[4][i]+u[5][i-1]+u[5][i+1])-ru[5][i];
}for(i=12; i21; i++){
ru[5][i]=0.25*(2.0*u[4][i]+u[5][i-1]+u[5][i+1]);
r[5][i]=0.25*(2.0*u[4][i]+u[5][i-1]+u[5][i+1])-ru[5][i];
}for(j=0; j6; j++){
for(i=0; i22; i++){
ru[j][i]=u[j][i];
}
}
if(r[j][i]0.001)
break;for(j=0; j6; j++){
for(i=0; i14; i++){
fprintf(out,%f ,ru[j][i]);}
}
}

신청하기