I want to solve numerically the equation that you can see on the “code” picture. The result that you can observe in “result” picture are not correct as you can see on the plot joined to this question.