summaryrefslogtreecommitdiff
path: root/1523/CH12/EX12.8/ex12_8.sce
blob: c7b06d159d793a178c9121e88edd5220e8124e60 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
// Network Synthesis : example 12.8 : (pg 12.5)
s=poly(0,'s');
p1=((s^4)+(5*(s)^2)+10);
p2=(5*(s^3)+(4*s));
[r,q]=pdiv(p1,p2);
[r1,q1]=pdiv(p2,r);
[r2,q2]=pdiv(r,r1);
[r3,q3]=pdiv(r1,r2);
printf("\nEven part of P(s) = ((s^4)+(5*(s)^2)+10)");
printf("\nOdd part of P(s) = (5*(s^3)+(4*s))");
printf("\nQ(s)= m(s)/n(s)");
// values of quotients in continued fraction expansion
disp(q);
disp(q1);
disp(q2);
disp(q3);
printf("\nSince two terms are negative, P(s) is not hurwitz");