gn="justify"> jStep = n/2;
for (int i = 0; i <= n; i + +)
{
for (int j = 0; j <= n; j + +)
{
if ((i> iStep) | | ((i <= iStep) && (j> jStep)))
{
mask [i] [j] = ActualPoint;
}
if (i == 0) mask [i] [j] = FictivePoint; // ліва межа
if (i == n) mask [i] [j] = FictivePoint; // права межа
if (j == 0) mask [i] [j] = FictivePoint; // нижня межа
if (j == n) mask [i] [j] = FictivePoint; // верхня межа
}
}
// створення маски квадратної області з вирізом і умовою Діріхле на кордонах кінець
// вивід маски початок
fprintf (file, "МАСКА ОБЛАСТІ! n");
for (int j = n; j> = 0; j -)
{
for (int i = 0; i <= n; i + +)
{
fprintf (file, "% d", mask [i] [j]);
}
fprintf (file, " n");
}
// вивід маски кінець
// обчислення u_accur, f залежно від маски початок
for (int i = 0; i <= n; i + +)
{
for (int j = 0; j <= n; j + +)
{
if (mask [i] [j] == ActualPoint)
{
u_accur [i] [j] = sin (M_PI * i * h) * sin (M_PI * j * h);
f [i] [j] = - (-2 * M_PI * M_PI * sin (M_PI * i * h) * sin (M_PI * j * h));
}
if (mask [i] [j] == TopPoint)
{
u_accur [i] [j] = sin (M_PI * i * h) * sin (M_PI * j * h);
f [i] [j] = - (M_PI * sin (M_PI * i * h) * cos (M_PI * j * h)/h);
}
if (mask [i] [j] == RightPoint)
{
u_accur [i] [j] = sin (M_PI * i * h) * sin (M_PI * j * h);
f [i] [j] = - (M_PI * cos (M_PI * i * h) * sin (M_PI * j * h)/h);
}
if (mask [i] [j] == BottomPoint)
{
u_accur [i] [j] = sin (M_PI * i * h) * sin (M_PI * j * h);
f [i] [j] = (M_PI * sin (M_PI * i * h) * cos (M_PI * j * h)/h);
}
if (mask [i] [j] == LeftPoint)
{
