A short Matlab implementation for P-1 and Q(1) finite elements (FE) is provided for the numerical solution of 2d and 3d problems in linear elasticity with mixed boundary conditions. Any adaptation from the simple model examples provided to more complex problems can easily be performed with the given documentation. Numerical examples with postprocessing and error estimation via an averaged stress field illustrate the new Matlab tool and its flexibility.