blob: 276b32e4d6f8650f57e1786155115078b07b5d73 (
plain)
1
2
3
4
5
6
7
|
s=%s
p=s^4+s^3+2*s^2+2*s+3
h=routh_t(p)
disp(h)
disp("The element of 1st column 7th row is zero.......its replaced by epsilon(very small quantity)")
disp("for epsilon<0 seventh element is negative of the 1st column giving 2 sign change implying two positive roots (real) ")
disp("for epsilon>0 tenth element is zero giving same conclusion")
|