Симметричной матрицы
Вычисление собственных значений и собственных векторов
Пример программы (лабораторная работа 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