Write a Matlab code that computes the shaded relief of a DEM

I am in need for a Matlab code that calculates the shaded relief of a DEM (Digital Elevation Model). The task should be done in Matlab and should use a set of equations that calculates shaded relief. I can provide more information about the task, the formulas and the DEM file upon request. More information on shaded relief can be found here:

http://wiki.gis.com/wiki/index.php/Shaded_Relief
https://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/how-hillshade-works.htm