Симметричной матрицы

Вычисление собственных значений и собственных векторов

Пример программы (лабораторная работа 4).

Варианты задания.

,

где ; – номер факультета; – номер группы; – номер студента по журналу.

Текст программы на F90

use dfimsl use linear_operators integer,parameter::n=3 real,parameter::s=3,g=6,s2=2*s,p=2*(s+g),m=s-p,a22=2*p+s real,dimension(n,n)::a=(/p,m,s2,m,a22,m,s2,m,p/),y,lambda(n) a=a/6; print 100,'matrix A',a lambda=eig(a,v=y); print 100,'eigenvalues',lambda; print 100,'eigenvectors',((y(i,j),j=1,n),i=1,n) 100 format(//10x,a/(3f8.3)) end

Результаты счета

matrix A

3.000 -2.500 1.000

-2.500 6.500 -2.500

1.000 -2.500 3.000

eigenvalues

9.000 2.000 1.500

eigenvectors

-0.408 0.707 0.577

0.816 0.000 0.577

-0.408 -0.707 0.577

Press any key to continue

Дата добавления: 2014-01-11