I am having problems with using the Newt() subroutine to solve one of the classical problems in fluid dynamics (Flow of Newtonian fluid over a rotating disk, known as Von Karman equations). The problem is, that Newt () uncapable to deal with a will known sensitive boundary value problem. I tried the code on other problems with the same number of equations (five Ode’s, with two unknown values). Despite, the vast amount of publications of different techniques of solving, numerically, these equations, I have a little doubt about that, at least by using some of the widely used techniques. If any one have any information about that or comments about what I am saying, please share with me.Thank you:) :) :)