Реферат: Решение задачи Дирихле для уравнения Лапласа методом сеток

arrayP[i][j].F_=(arrayP[i-1][j].F+arrayP[i+1][j].F+

arrayP[i][j-1].F+arrayP[i][j+1].F)/4;

else arrayP[i][j].F_=0; }

}

void CountDif() // find maximal difference abs(F-F_)

{

k=0;

for(i=0;i<13;i++)

for(j=0;j<33;j++)

{ if (arrayP[i][j].BelongsToDh==1) {

diff[k]=fabs(arrayP[i][j].F_-arrayP[i][j].F);

k++;}}

E1=diff[0];

for (k=1;k<500;k++) {

if (diff[k]>E1) E1=diff[k];}

}

void MakeFile()

{

ofstream f;

FILE *f1=fopen("surf.dat","w1");

fclose(f1);

f.open("surf.dat",ios::out,0);

for(i=0;i<13;i++)

for(j=0;j<33;j++) { if (arrayP[i][j].BelongsToDh==1) {

f<<arrayP[i][j].xx<<" "<<arrayP[i][j].yy<<

" "<<arrayP[i][j].F_<<"\n";}}

f.close() ;

}

4. ГРАФИКИ РЕШЕНИЙ