developing a Matlab software to get the natural frequency , and dynamic parameters (modal shape for each mode)

i am trying to develop a matlab program to calculate the natural frequency anf the model shapes for beams

using the finite element method , as i will need to change the stiffness in the beam sections to model cracks.

 i am sure there an easy ready method to perform this , i am happy if some one provided me with the method or the basic code