path: root/569/CH4/EX4.23/4_23.sci

// determine the error
clc;
w=6;
wn=4;
u=w/wn;
eta=0.66;
M=1/{[(1-u^2)^2]+(2*eta*u)^2}^0.5;
Error=(M-1)*100;
disp(Error,'error (%)=')