257/CH7/EX7.38/example_7_38.sce


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

Mp=0.5/2 *100;
zeta=0.4036;

for(n=0:3)
    T=n*%pi/omegaD
    if(n==2)
        T=0.2
        omegaN=2*%pi/(T*(sqrt(1-zeta^2)))
        disp(omegaN," omegaN = ")
    end
end