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

ify"> if (M [i] [j]! = M [j] [i])

{

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 4);

cout <

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 7);

}

}

cout <

}

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 4);

cout <<"матриця не симетрично!" <

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 7);

}

return 0;

}

// перевірка на симетричний кінець

// перевірка на знакоопределенная началоMore_then_0 (int n, double ** M)

{

double d;

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 9);

cout <<"Наводимо матрицю до трикутного вигляду ..." <

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 7);

for (int k = 0; k

{

for (int i = k +1; i

{

if (M [k] [k] == 0)

{

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 4);

printf ("Сталося ділення на 0! nДля виходу натисніть будь-яку клавішу ... n");

_getch ();

exit (0);

}

else

{

d = M [i] [k]/M [k] [k];

for (int j = 0; j

{

M [i] [j] = M [i] [j]-M [k] [j] * d;

}

}

}

}

int z = 0, p = 0;

double DetM = 1;

double * detM;

detM = new double [n];

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 9);

cout <<"Знаходимо головні мінори ..." <

SetConsoleTextAttribute (GetStdHandle (STD_OUTPUT_HANDLE), 7);

for (int i = 0; i

{

DetM * = M [i] [i];

detM [i] = DetM;

printf ("мінор [% d] =% lf n", i +1, detM [i]);

}

for (int i = 0; i

{

if (detM [i]> 0)

{

z + +;

}