Finished morozov implementation and newtons iterations. Seems to give very small alpha values -> bugcheck