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

arrayY[i+1]=tmp;}}

LastUnreg=LastUnreg-1; }

for (i=0;i<33;i++) { printf("%g ",arrayY[i]); }

printf("\n");} // End of Regulation

void CrMtrD(void) //Create general Matrix

{

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

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

arrayP[i][j].BelongsToDh=0;}

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

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

arrayP[i][j].xx=arrayX[i];

arrayP[i][j].yy=arrayY[j];

}

// printf("%g %g",arrayP[12][0].xx,arrayP[12][0].yy);

// printf("\n");

}

int IsFit(point Par) //does point belong to area D?

{

if ((Par.xx<=4) && (Par.xx>=1.99) && (Par.yy>=Par.xx)

&& (Par.yy<=Par.xx+4)) return 1;

else return 0;

}

void CreateDh_(void) //Create area Dh_

{

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

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

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

cout << arrayP[1][1].BelongsToDh_<< "\n";