; = (float *) malloc (M * sizeof (float)); = fopen ("d: rezultat.txt", "w" );
// обчислення і запис у файл нульового слоя_sloya (0); (j = 0; j <= M; j + +)
{[j] = fi (j * h); _v_file (u0 [j]);
}
// обчислення і запис у файл першого слоя_sloya (1); (j = 0; j <= M; j + +)
{[j] = psi (j * h) * l + u0 [j]; _v_file (u1 [j]);
}
/* обчислення і запис у файл все з інших верств i = 2, .. N
Схема тришарова, для обчислення одного шару необхідні значення функції у вузлах на попередніх 2-х шарах.
Для обчислення, наприклад 2го шару використовуються значення на нульовому і першому, які зберігаються в масивах u0 і u1 відповідно. Значення ж третього шару записуються в масив u2. p align="justify"> Для обчислення на третьому шарі використовуються значення на першому і другому, які зберігаються в масивах u1 і u2 відповідно. Значення ж третього шару записуються в масив u0. p align="justify"> Для обчислення на четвертому шарі використовуються значення на другому і третьому шарах, які зберігаються в масивах u2 і u0 відповідно. Значення ж третього шару записуються в масив u1. p align="justify"> І так далі до останнього шару. */(i = 2; i <= N; i + +)
{
octatok = fmod (i, 3); (octatok)
{//обчислюються значення на шарах з i = 3,6,9, ... 0:
// обчислюються значення на шарах з i = 4,7,11, ... 1:
// обчислюються значення на шарах з i = 2,5,8, ... 2: zapis_sloya (i); (j = 1; j }
} (u0); (u1); (u2); (); 0;
}