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

}

void FillF(void) // calc function F(x,y) at area Dh_

{

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

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

if (arrayP[i][j].BelongsToDh_==1)

arrayP[i][j].F=arrayP[i][j].xx*pow(arrayP[i][j].yy,2);

else arrayP[i][j].F=0;

}

int IsInner(int i,int j) //Is point inner?

{

if ((arrayP[i-1][j].BelongsToDh_==1) &&

(arrayP[i+1][j].BelongsToDh_==1) &&

(arrayP[i][j+1].BelongsToDh_==1) &&

(arrayP[i][j-1].BelongsToDh_==1)) return 1;

else return 0;

}

void CreateDh(void) //Create area Dh

{

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

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

if ((arrayP[i][j].BelongsToDh_==1) &&

IsInner(i,j))

arrayP[i][j].BelongsToDh=1;

}

void FillF_() //calc new appr. values of F

{

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

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