blob: d10efd0d1347e1654ee7de22578c25b1d4904ed2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
//breakaway points
s=%s
sys=syslin('c',(s+4)/(s*(s+2)))
evans(sys)
syms s
d=derivat(sys)
n=numer(d)
a=roots(n) //a=breakaway points
disp(a,"breakaway points=")
for i=1:2
K=-a(i,1)*(a(i,1)+2)/(a(i,1)+4)
disp(a(i,1),"s=")
disp(K,"K=")
end
printf("if K is positive breakaway point lies on RL or else on CRL")
|
