257/CH7/EX7.10/example_7_10.sce


1
2
3
4
5
6
7
8
9
10
11
12
13

syms k1 k2 s
T= k1/(s^s+k1*k2*s+k1)

Mp=25               //given
zeta=0.4037
Tp=4

omegaN=%pi/(Tp*(sqrt(1-zeta^2)))
disp(omegaN,"omegaN = ")
k1=omegaN^2
disp(k1," k1 = ")
k2=2*zeta/(sqrt(k1))
disp(k2," k2 = ")