Photometric stereo Matlab code