public double u01, u02, u03;
public U(double x,double y,double z)
static void Main(string[] args)
double[,] u = new double[3, 1];
double[,] J2u = new double[3, 3];
double[,] J1u = new double[3, 1];
double[,] Jt = new double[3, 3];
double[,] J_1 = new double[3, 3];
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
Console.Write("J''u[{0},{1}]=", i + 1, j + 1);
J2u[i, j] = Convert.ToDouble(Console.ReadLine());
J1u[0, 0] = 10 * u[0,0] + 2 * u[1,0] - 2 * u[2,0] + 2;
J1u[1, 0] = 2 * u[0,0]+ 4 * u[2,0] + 2;
J1u[2, 0] = -2 * u[0,0] + 4 * u[1,0] + 2 * u[2,0] + 1;
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
J = J2u[0, 0] * J2u[1, 1] * J2u[2, 2] + J2u[0, 1] * J2u[1, 2] * J2u[2, 0] + J2u[1, 0] * J2u[2, 1] * J2u[0, 2] - J2u[2, 0] * J2u[1, 1] * J2u[0, 2] - J2u[1, 0] * J2u[0, 1] * J2u[2, 2] - J2u[0, 0] * J2u[1, 2] * J2u[2, 1];
Console.WriteLine("J=" + J);
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
Jt[i, j] = Math.Pow(-1, Math.Abs(0)) * (J2u[1, 1] * J2u[2, 2] - J2u[1, 2] * J2u[2, 1]);
Jt[i, j] = Math.Pow(-1, Math.Abs(1)) * (J2u[1, 2] * J2u[2, 0] - J2u[1, 0] * J2u[2, 2]);
Jt[i, j] = Math.Pow(-1, Math.Abs(2)) * (J2u[1, 0] * J2u[2, 1] - J2u[1, 1] * J2u[2, 0]);
for (int j = 0; j < 3; j++)
Jt[i, j] = Math.Pow(-1, Math.Abs(1)) * (J2u[2, 1] * J2u[0, 2] - J2u[2, 2] * J2u[0, 1]);
Jt[i, j] = Math.Pow(-1, Math.Abs(2)) * (J2u[2, 2] * J2u[0, 0] - J2u[2, 0] * J2u[0, 2]);
Jt[i, j] = Math.Pow(-1, Math.Abs(3)) * (- J2u[2, 0] * J2u[0, 1] + J2u[2, 1] * J2u[0, 0]);
for (int j = 0; j < 3; j++)
Jt[i, j] = Math.Pow(-1, Math.Abs(2)) * (J2u[0, 1] * J2u[1, 2] - J2u[0, 2] * J2u[1, 1]);
Jt[i, j] = Math.Pow(-1, Math.Abs(3)) * (J2u[0, 2] * J2u[1, 0] - J2u[0, 0] * J2u[1, 2]);
Jt[i, j] = Math.Pow(-1, Math.Abs(4)) * (J2u[0, 0] * J2u[1, 1] - J2u[0, 1] * J2u[1, 0]);
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
if (i == j + 1 || i == j - 1)
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
Console.Write("Jt[{0},{1}]={2}\t", i + 1, j + 1, Jt[i, j]);
for (int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
J_1[i, j] = 1.0 / J * Jt[i, j];
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
Console.Write("J^(-1)[{0},{1}]={2}\t", i + 1, j + 1, J_1[i, j]);
double[,] x = new double[3, 1];
for (int i = 0; i < 3; i++)
for (int j = 0; j < 3; j++)
x[i, 0] = J_1[i, 0] * J1u[0, 0] + J_1[i, 1] * J1u[1, 0] + J_1[i, 2] * J1u[2, 0];
Console.WriteLine("{0} po schetu iteracia",o);
for (int i = 0; i < 3; i++)
J1u[0, 0] = 10 * u[0, 0] + 2 * u[1, 0] - 2 * u[2, 0] + 2;
J1u[1, 0] = 2 * u[0, 0] + 4 * u[2, 0] + 2;
J1u[2, 0] = -2 * u[0, 0] + 4 * u[1, 0] + 2 * u[2, 0] + 1;
Console.WriteLine("Norma pervoi proizvodnoi={0}\nu =({1};{2};{3})", k, u[0, 0], u[1, 0], u[2, 0]);
Console.WriteLine("\n\tKolichestvo itogovyh iteracii={0}\nu =({1};{2};{3})", o, u[0, 0], u[1, 0], u[2, 0]);