Solve Beams using Finite Element Method

These function files solve any engineering beam problem using Finite Element Method.

There are multiple function files. One finds the total stiffness matrix for a beam. Since this is a 2-D beam solver which means each of the nodes in this Euler Bernoulli beam has 2 DOF only (uy and phi), the order of the total stiffness matrix is number of nodes times 2.
Another function file finds the equivalent nodal force vector due to the distributed loads.
solve_beam.m is the main function to be called (like shown in the example file beam_problem.m) that calls the these two functions to solve the beam.
draw_frame and animate functions draw the beam and its displacement at the names suggest.