Matlab code goes here!