Matlab code for newton raphson method