blob: 1a3e2fe4d8b4ece5fe32e31823c0af74ea9a637a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
//Eg-5.14
//pg-252
clear
clc
A=[1 3 6;3 5 7;6 7 4];
k=1;
summ=0;
for i=2:3
summ=summ+A(1,i)^2;
end
sig=summ^.5;
w1=0;
w2=(.5*(1+(abs(A(1,2)/sig))))^.5;
w3=A(1,3)/(2*sig*w2);
w=[w1;w2;w3];
B=eye(3,3)-2*w*w';
A1=inv(B)*A*B;
disp("required tridiagonal matrix is")
disp(A1)
|
