Текст програми. /*пiдпрограма обчислення визначника*/

#include <stdio.h>

#include <math.h>

/*пiдпрограма обчислення визначника*/

/*методом алгебраiчних доповнень*/

float det(float a[10][10],int n) {

int j,k,l,k1,l1,zn,n1;

float b[10][10],deter,d;

if (n==1) deter=a[0][0];

if (n==2) deter=a[0][0]*a[1][1]-a[0][1]*a[1][0];

if (n>2) {

deter=0;

for (j=0;j<n;j++) {

if (fmod(j,2)==0) zn=1;

else zn=-1;

k=1;l=0;n1=n-1;

for (k1=0;k1<n1;k1++,k++) {

l=0;

for (l1=0;l1<n1;l1++,l++) {

if (l1==j) l++;

b[k1][l1]=a[k][l];

};

};

d=det(b,n1);

deter+=a[0][j]*zn*d;

};

};

return(deter);

}

main () {

float a,mas[10][10];

int i,j,p;

FILE *res;

res=fopen("deter.txt","w");

printf("Введiть порядок визначника ");

scanf("%d",&p);

/*ввiд коефiцiентiв рiвняннь*/

for(i=0;i<p;i++) {

for(j=0;j<p;j++) {

scanf("%f",&a);

mas[i][j]=a;

};

};

fprintf(res,"Матриця\n");

for(i=0;i<p;i++) {

for(j=0;j<p;j++) fprintf(res," %5.2f",mas[i][j]);

fprintf(res,"\n");

};

a=det(mas,p);

fprintf(res,"визначник матрицi det=%5.3f\n ",a);

fclose(res);

}

Результат роботи програми

Матриця

3.00 6.00 4.00 1.00

7.00 0.00 4.00 3.00

2.00 5.00 7.00 9.00

3.00 7.00 4.00 2.00

визначник матрицi det=-308.000

Розв’язок системи рівнянь методом Крамера

Один з методів розв’язку системи лінійних рівнянь - правило Крамера. Використовуючи позначення попереднього розділу, вони формулюються так: якщо система лінійних рівнянь сумісна, то вона має єдиний розв’язок

a1= Текст програми. /*пiдпрограма обчислення визначника*/ - student2.ru

де dj - визначник, одержаний із визначника системи d заміною j-го стовпця стовпцем із вільних членів системи рівнянь.

Представлена програма на мові С реалізує метод Крамера. Визначники обчислюються за тим же правилом, що й у попередньому розділі - з використанням рекурсії. Для порівняння, та ж система рівнянь розв’язана методом Крамера за допомогою пакету MathCad. Результати співпадають.

Наши рекомендации