Matlab code for project work