Реферат Рішення завдання Неймана для рівняння Пуассона в прямокутній області

0, h, tau,

eps = 0.00000001,

** A_u, ** A_u_accur,

** A_u_u_accur, ** r,

** z, ** z_1d, * f_1d;

char fname [] = "out.txt"; // ім'я вихідного файлу

setlocale (LC_CTYPE, "Russian"); // ставимо російську мову консолі

printf ("Введіть кількість вузлів сітки:") ;//

scanf_s ("% d", & n); //

n -; //

h = 1./n; // обчислюємо крок сітки

printf ("Виведення даних у файл% s ... n", fname);

file = fopen (fname, "w");

// ініціалізація змінних початок

mask = new int * [n +1];

u = new double * [n +1];

u_accur = new double * [n +1];

f = new double * [n +1];

res_S = new double * [n +1];

res_A = new double * [n +1];

A_u = new double * [n +1];

A_u_accur = new double * [n +1];

A_u_u_accur = new double * [n +1];

r = new double * [n +1];

z = new double * [n +1];

for (int i = 0; i <= n; i + +)

{

mask [i] = new int [n +1];

u [i] = new double [n +1];

u_accur [i] = new double [n +1];

f [i] = new double [n +1];

res_S [i] = new double [n +1];

res_A [i] = new double [n +1];

A_u [i] = new double [n +1];

A_u_accur [i] = new double [n +1];

A_u_u_accur [i] = new double [n +1];

r [i] = new double [n +1];

z [i] = new double [n +1];

}

for (int i = 0; i <= n; i + +)

{

for (int j = 0; j <= n; j + +)

{

mask [i] [j] = 0;

u [i] [j] = 0;

u_accur [i] [j] = 0;

f [i] [j] = 0;

res_S [i] [j] = 0;

res_A [i] [j] = 0;

z [i] [j] = 0;

}

}

// ініціалізація змінних кінець

// створення маски квадратної області з вирізом і умовою Діріхле на кордонах початок

int iStep, jStep;

iStep = n/2;