pproxim_correct;// Значення аппроксім. з поправкою [3] [4];/​​/ Система рівнянь для пошуку коеф. квіт. функції
double Parab [3] = {0,0,0};// Коефіцієнти апроксимуючої функції
// ---------------------------------------- -----------------------------------
// Прототипи функцій
doubleMy (double, double);// Функція у правій частині діф.уравненія
double ApproxFunc (double);// апроксимується функція
// ---------------------------------------- -----------------------------------
__fastcall TForm1 :: TForm1 (TComponent * Owner)
: TForm (Owner)
{
}
{
solutions = newdouble [n];// Виділення пам'яті під масиви значень
approxim = new double [n]; = new double [n]; = new double [n]; _half = new double [n * 2]; = new double [n]; _correct = new double [n]; _half = new double [n * 2]; (int i = 0; i
// Обчислення вузлів сітки з півкроком (int i = 0; i <2 * n; i + +) knots_half [i] = x0 + i * h_half;
}
{i; k1, k2, k3, k4; [0] = y0;// Початкове условіе_half [0] = y0; -> Enabled = true;// Дозвіл раніше недоступних об'єктів управління
Solution-> Enabled = true;
// Обчислення методом Рунге-Кутта 4 порядку
for (i = 0; i
}
// Обчислення методом Рунге-Кутта 4 порядку з півкроком
for (i = 0; i <2 * n-1; i + +)
}
// Висновок рішення на екран
Screen-> Lines-> Add ("Рішення диф. рівняння");
Lines-> Add ("x =" + AnsiString (knots [i]) +
" t f (x, y) =" +
// Висновок рішення з півкроком на екран
Screen-> Lines-> Add ("Рішення диф. рівняння c півкроком");
Lines-> Add ("x =" + AnsiString (knots_half [i]) +
" t f (x, y) =" +
}
{cur = 0; i, j, k, p, q;
// Обчислення коефіцієнтів системи рівнянь для обчислення квіт. функції
// Відповідно до методу найменших квадратів
// a11 a12 a13 | b1
// a21 a22 a23 | b2
// a31 a32 a33 | b3tmp = 0; [3] [4];
CalcDy-> Enabled = true;// Дозвіл раніше недоступних органів управління
