tify"> 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] == LeftBottomPoint)
{
u_accur [i] [j] =
}
}
}
// обчислення u_accur, f залежно від маски кінець
int matrixSize = 0;
for (int i = 0; i <= n; i + +)
{
for (int j = 0; j <= n; j + +)
{
if (mask [i] [j]! = FictivePoint)
{
matrixSize + +;
}
}
}
z_1d = new double * [matrixSize];
for (int i = 0; i <= matrixSize; i + +)
{
z_1d [i] = new double [matrixSize];
}
// Перевірка спряженості - початок
int p = 0, k = 0;
for (int i = 0; i <= n; i + +)
{
for (int j = 0; j <= n; j + +)
{
if (i> 0)
{
z [i-1] [n] = 0;
z [i-1] [n-1] = 0;
z [2] [1] = 0;
z [2] [2] = 0;
}
if (mask [i] [j]! = FictivePoint)
{
z [i] [j-1] = 0;
z [i] [j] = 1;
A (n, mask, z, res_A);
for (int l = 0; l <= n; l + +)
{
for (int m = 0; m <= n; m + +)
{
if (mask [l] [m]! = FictivePoint)
{
z_1d [p] [k] = res_A [l] [m];
p + +;
}
}
}
k + +;
p = 0;
}
}
}
f_1d = new double [matrixSize];
p = 0;
for (int l = 0; l <= n; l + +)
{
for (int m = 0; m <= n; m + +)
{
if (mask [l] [m]! = FictivePoint)
{
f_1d [p] = f [l] [m];
p + +;
}
}
}
print (matrixSize, z_1d, f_1d);
// Перевірка ...
