--> //initializations

--> Cu=4*10^4 //(coefficient of elastic uniform compression)
 Cu  = 

   40000.


--> D=0 // (damping factor)
 D  = 

   0.


--> A=2.5 //m^2 (base area)
 A  = 

   2.5


--> mf=2.5 //Mg (mass of vertical compressor)
 mf  = 

   2.5


--> fo=40 //hertz (operating frequency)
 fo  = 

   40.


--> //calculations

--> ms=0.2*mf 
 ms  = 

   0.5


--> m=mf+ms //(total mass)
 m  = 

   3.


--> fn=(sqrt((Cu*(10^3)*A)/(m*10^3)))/(2*%pi) //(natural frequency)
 fn  = 

   29.057584


--> r=fo/fn //(ratio of frequency)
 r  = 

   1.3765769


--> M=1/sqrt((1-r^2)^2+(4*D^2*r^2)) //(magnification factor)
 M  = 

   1.1173633


--> //results

--> printf('natural frequency=%f',fn) //hertz
natural frequency=29.057584
--> printf('magnification factor=%f',M)
magnification factor=1.117363
