Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

I'm currently trying to find a numerical solution to a differential equation of the form:

D[W[X], {X, 4}] ==(-(1/(delta + G - (G X)/L)^2) + 1/(delta + (G X)/L)^2)

Subject to the boundary conditions:

W[0] == 0, W[L] == 0, W'[0] == 0, W'[L] == 0

NDSolve appears to be choking on a singularity or a stiffness constraint. Any ideas on how I might go about dealing with this?

Thanks in advance!

share|improve this question
Welcome! Providing actual code is always very useful to understand/reproduce your problem. You can also format your code, have a look at mathematica.stackexchange.com/help/formatting. –  Yves Klett May 13 '14 at 15:30

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.