blob: fcc9a97430afd19e03bbb18424b3b92a0452ebb0 (
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
//Eg-4.20
//pg-198
clear
clc
a=[1 -10 50 -120 144];
b0=0;
b01=1;
g=-1;
d=1;
for i=1:60
b1=a(2)-g*b01-d*b0;
b2=a(3)-g*b1-d*b01;
r1=a(4)-g*b2-d*b1;
r2=a(5)-d*b2;
dnew=a(5)/b2;
gnew=(a(4)-dnew*b1)/b2;
d=dnew;
g=gnew;
end
disp("Hence gamma delta b1 b2 are respectively")
disp(g);
disp(d);
disp(b1);
disp(b2);
m=[8 -4 1];
f=poly(m,'x','c');
p=roots(f);
disp("two roots are")
disp(p);
m=[18 -6 1];
f=poly(m,'x','c');
p=roots(f);
disp("other two roots are");
disp(p);
disp("second part")
a=[1 -8 27 -50 50];
b0=0;
b01=1;
g=-1;
d=1;
for i=1:60
b1=a(2)-g*b01-d*b0;
b2=a(3)-g*b1-d*b01;
r1=a(4)-g*b2-d*b1;
r2=a(5)-d*b2;
dnew=a(5)/b2;
gnew=(a(4)-dnew*b1)/b2;
d=dnew;
g=gnew;
end
disp("Hence gamma delta b1 b2 are respectively")
disp(g);
disp(d);
disp(b1);
disp(b2);
m=[5 -2 1];
f=poly(m,'x','c');
p=roots(f);
disp("two roots are")
disp(p);
m=[10 -6 1];
f=poly(m,'x','c');
p=roots(f);
disp("other two roots are");
disp(p);
|
