А я как-то ради прикола написал прогу для расчёта определителя
матрицы по теореме старика Лапласа. Вручную это полный заё...
Матрица 4х4, интересно то, что помогли квадратные уравнения.
Текст Laplas.cpp:
#include <iostream>
#include <stdlib>
#include <math>
#define PR cout<<
#include <John>
short main(void)
{short deti; PR " Hello body!\n\n Matrix:\n"; PR John(deti); return 0;}
Текст John.h:
short John (short deti)
{ short string1=0, string2=0, help1=0,det=0, help2=0, C=0,B=0,
scherta1=0, scherta2=0, hcherta1=0, hcherta2=0, D=0,
helpt2=0, helpt1=0, A[4][4], H[]={5,3,1,5,6,1,3,3,0,2,0,1,4,1,1,2};
for (help1=0; help1<=3; help1++)
{cout<< "\n"; for (help2=0; help2<=3;help2++)
{A[help1][help2]=H[help2+help1*4]; PR A[help1][help2]<<" "; }}
PR"\n\n Number first worker line ";
cin>>string1;
PR"\n Number second worker line ";
cin>>string2;
if(string1<1||string1>4||string1==string2||string2<1||string2>4||string1>string2){PR "Invalid parameters are Incorporated"; exit(1);}
C=24/(string1*string2); B=10-(string1+string2); D=abs(B*B-4*C); D=sqrt(D);
help1=abs((B-D)/2);
help2=abs((B+D)/2);
scherta2=(help1>help2)?help1:help2;
scherta1=(help1+help2)-scherta2;
for(help1=1;help1<5;help1++){
for(help2=2;help2<5;help2++){
if(help2<=help1) continue;
C=24/(help1*help2); B=10-(help1+help2); D=abs(B*B-4*C); D=sqrt(D);
helpt1=abs((B-D)/2);
helpt2=abs((B+D)/2);
hcherta2=(helpt1>helpt2)?helpt1:helpt2;
hcherta1=(helpt1+helpt2)-hcherta2;
det=det+(((A[string1-1][help1-1]*A[string2-1][help2-1]-A[string1-1][help2-1]*
A[string2-1][help1-1])*(A[scherta1-1][hcherta1-1]*
A[scherta2-1][hcherta2-1]-A[scherta1-1][hcherta2-1]*
A[scherta2-1][hcherta1-1]))*pow((-1),(string1+string2+help1+help2)));}}
PR"\n\tDeterminante of your matrix = ";
deti=det;
return deti+0; }
|
