[shikha5]

Hi, I have written code for 2d diffusion problem and want to get its isotherm contour and temperature at mid y plane.

/*----------------------output---------------------------*/
fprintf(b,"I=%d\tJ=%d \n",nx+1,ny+1);
for(j=0;j<=ny;j++)
{
for(i=0;i<=nx;i++)
{
fprintf(b,"%f\t%f\t%f\n",x[i],y[j],th[i][j]);
}
}
for(i=0;i<=nx;i++)
{
fprintf(c,"\n%f\t%f",x[i],th[i][ny/2]);
}

Is it a right way to write this part of code to get the output, because I am not getting desired results. Any suggestions plz.