%===================================
% THE NUMBER a
%===================================
a = input('Enter the number a: ')
if(a==0||a==1||a==2||a==3||a==4)
n = 2;
else
n = 3;
end
%===================================
% THE MATRIX A
%===================================
A = zeros(n,n);
for p = 1 : n
for q = 1 : n
A(p,q) = (a+pi)/(p+q+a);
end
end
%===================================
% THE VECTOR b
%===================================
b = zeros(n,1);
for p = 1 : n
b(p) = 1/p;
end
%===================================
% THE INVERSE MATRIX
%===================================
B = inv(A);
%===================================
% THE SOLUTION OF THE LINEAR SYSTEM
%===================================
x = B * b;
display('The solution: ')
x
display('The norm of x: ')
norm(x)
%===================================
% THE MATRIX C
%===================================
C = A*A+B;
e = eig(C);
display('Eigenvalues: ')
e
display('Their product: ')
prod(e)
Назад към LaTeX, Matlab, Maple, Mathematika...
Регистрирани потребители: Google [Bot]