--> //initializations

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

   40000.


--> m=15 //Mg (mass)
 m  = 

   15.


--> l=2 //m (length)
 l  = 

   2.


--> b=2 //m (breadth)
 b  = 

   2.


--> //calculations

--> A=l*b //(area)
 A  = 

   4.


--> wn=sqrt((Cu*(10^3)*A)/(m*10^3)) //(circular frequency)
 wn  = 

   103.27956


--> f=wn/(2*%pi) //(natural frequency)
 f  = 

   16.437452


--> //results

--> printf('natural frequency=%f',f) //cps (hertz)
natural frequency=16.437452
