Code matlab of newton raphson method