blob: f532a2bb5e060a4ceee59dbf735c14e9e9bd9c1f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
//page 34
clear;
close;
clc;
a=rand(1);
b=rand(1);
c=rand(1);
L=[1 0 0;a 1 0;b c 1];
disp(L,'L=');
U=eye(3,3);
disp(U,'U=');
E=[1 0 0;-a 1 0;0 0 1];
disp(E,'E=');
F=[1 0 0;0 1 0;-b 0 1];
disp(F,'F=');
G=[1 0 0;0 1 0;0 -c 1];
disp(G,'G=');
disp('A=inv(E)*inv(F)*inv(G)*U')
A=inv(E)*inv(F)*inv(G)*U;
disp(A,'A=');
disp('When U is identity matrix then L is same as A');
//end
|
