At y=0 the function 'psi' in attached file and its second derivative are vanished. Also by taking G as a constant amount, 'G,y' is zero. Therefore the problem has solution.
We haven't the integration constant. By assuming this constant equal to zero, I try to solve the problem (it is a guess!)
My questions are
1- The MAPLE can not solve the last equations in 1.mw by fsolve. What is the reason? The problem hasn't answer or some codes like the Newton's iterative method must be used instead of fsolve?
2- Is it possible to minimize 'phi' by using two coupled nonlinear equations in 2.mw? If yes then which commands can be used?
3- Please suggest a method to solve the problem, if it is possible.
Thanks a lot