[SciPy-User] How to use Scipy to solve a Two-Point Boundary Value problem, of a nth system Nonlinear Second-Order Differential Equation

[freeblue]

Hello to everyone here, 

I am trying to find the root of the next equation using Newton method:
http://www.wolframalpha.com/input/?i=d%2Fdx%28du%2Fdx%29+%3D+%28-3%2F%28k1%29[x]%29*%28k4[x]-%28k2[x]%2Bk3[x]%29*u
^
%281%2F3%29%29+*%28%28du%2Fdx%29^%282%2F3%29%29
which has 2 boundary conditions (u(x=0)=0, u(x=n(max) = m (constant)) ...
k1, k2, k3, k4 are arrays, calculated on each x gridpoint

Please, do someone knows if it is possible to solve this nth nonlinear 
second-order differential equation using scipy.integrate.odeint? (nonlinear 
ODE BVP 1-D.
I tried to use Central differencing to "simplify" it but the power number 
(2/3) makes its hard to use later Newton method.

Could be possible then to use scipy?  Do you know or see somewhere any 
relative example ? I am trying to find and write a solution on this the 
last 3 weeks, so any help will be more than welcome!!!!!

thank you in advance,
Kas

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20130330/9ec3c4f2/attachment.html>