Мои наработки программы, но она неправильно работает. Подскажите, пожалуйста, что не так в алгоритме и вообще?
void __fastcall TForm1::Button2Click(TObject* Sender)
{
x = StrToFloat(Edit1->Text);
y = StrToFloat(Edit2->Text);
z = StrToFloat(Edit3->Text);
esp = StrToFloat(Edit4->Text);
double a[n][n], a1[n][n], b[n], det, invdet, f[n],
K[n]; //a[3][3]-матрица Якоби;
//a1[3][3]-обратная матрица; b[3]-произведение а1[3] на f[3]; f[3]-массив функций; K[3]-массив переменных
double xk, yk, zk, xn, yn, zn,
max; //изменение значения переменных при k=0,1,2..
int k = 0; //кол-во итераций
//Матрица Якоби
{
a[0][0] = (-2 * x - 1);
a[0][1] = (2 * z);
a[0][2] = (2 * y);
a[1][0] = (-3 * z);
a[1][1] = (2 * y - 1);
a[1][2] = (-3 * x);
a[2][0] = (-2 * y);
a[2][1] = (-2 * x);
a[2][2] = (-2 * z - 1);
//Определитель
det = a[0][0] * a[1][1] * a[2][2] + a[0][1] * a[1][2] * a[2][0] +
a[0][2] * a[1][0] * a[2][1] -
a[0][2] * a[1][1] * a[2][0] - a[0][0] * a[1][2] * a[2][1] - a[0][1] *
a[1][0] * a[2][2];
invdet = 1 / det;
//Обратная матрица
for (i = 0; i < n; i++)
{
for (j = 0; j < n; j++)
{
a1[0][0] = a[1][1] * a[2][2] - a[1][2] * a[2][1] * invdet;
a1[0][1] = -a[1][0] * a[2][2] + a[2][0] * a[1][2] * invdet;
a1[0][2] = a[1][0] * a[2][1] - a[1][1] * a[2][0] * invdet;
a1[1][0] = -a[0][1] * a[2][2] + a[0][2] * a[2][1] * invdet;
a1[1][1] = a[0][0] * a[2][2] - a[0][2] * a[2][0] * invdet;
a1[1][2] = -a[0][0] * a[2][1] + a[0][1] * a[2][0] * invdet;
a1[2][0] = a[0][1] * a[1][2] - a[0][1] * a[1][1] * invdet;
a1[2][1] = -a[0][0] * a[1][2] + a[0][2] * a[2][0] * invdet;
a1[2][2] = a[0][0] * a[1][1] - a[0][1] * a[1][0] * invdet;
}
}
//Массив функций
for (j = 0; j < n; j++)
{
f[0] = f1(x, y, z);
f[1] = f2(x, y, z);
f[2] = f3(x, y, z);
}
//Массив переменных
for (j = 0; j < n; j++)
{
K[0] = x;
K[1] = y;
K[2] = z;
}
do {
for (i = 0; i < n; i++)
{
for (j = 0; j < n; j++)
{
b[0] = a[0][0] * f[0] + a[0][1] * f[1] + a[0][2] * f[2];
b[1] = a[1][0] * f[0] + a[1][1] * f[1] + a[1][2] * f[2];
b[2] = a[2][0] * f[0] + a[2][1] * f[1] + a[2][2] * f[2];
}
}
xk = K[0] - b[0];
yk = K[1] - b[1];
zk = K[2] - b[2];
xn = xk - x;
yn = yk - y;
zn = zk - z;
{
if (abs(xn) > abs(yn))
{
if (abs(xn) > abs(zn))
{
max = xn;
}
else
{
max = zn;
}
}
else
{
if (abs(yn) > abs(zn))
{
max = yn;
}
else
{
max = zn;
}
}
}
x = xk;
y = yk;
z = zk;
} while (abs(max) >= esp);
k++;
Edit5->Text = FloatToStr(xk);
Edit6->Text = FloatToStr(yk);
Edit7->Text = FloatToStr(zk);
Edit8->Text = FloatToStr(k);
}
}
она неправильно работает
. На каком наборе данных? Какой результат ожидается/получается? И т.д.