summaryrefslogtreecommitdiff
path: root/260/CH4/EX4.16/4_16.sce
blob: 7952ea328f67ef114feaca8707db9fc6fc5c2e71 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
//Eg-4.16    
//pg-182

clear
clc


//co-efficients of the polynomial
a=[6.3825 -21.52 21.26 -8 1];
//one root
l=3.7;
n=5;
for i=n-2:-1:1
    b(n-1)=a(n);
    b(i)=b(i+1)*l+a(i+1)
end

fdef=poly(b,'x','c');
disp(fdef)
printf('\nOr\n    The coefficient matrix is \n')
disp(b)