Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: question about ndsolve
Replies: 2   Last Post: Sep 17, 2013 9:30 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Hagwood, Charles R.

Posts: 12
Registered: 12/31/09
question about ndsolve
Posted: Sep 10, 2013 3:24 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


I know the following differential equation has a solution, but cannot get NDSolve to find it.


dgamma[s_,a_]:=1+9.6*(s-1)^2*s+9.6*s^2*(s-1)
q2[x_]:=Sqrt[2*Pi]*{-Sin[2*Pi*x],Cos[2*Pi*x]}
q1[x_,a_]:=Sqrt[2*Pi*dgamma[x,a]]*{-Sin[2*Pi*gamma[x,a]],Cos[2*Pi*gamma[x,a]]}


h1[x_]:=q2[x].q2'[x]
h2[x_]:=q2[x].q2[x]


sol=NDSolve[{-6*y'[x]^2*(q1[x,.3].q2'[y[x]]) +6*y'[x]^(5/2)*h1[y[x]]+y''[x]*(q1[x,.3].q2[y[x]])-y''[x]*y'[x]^(1/2)*h2[y[x]]==0,y[0]==0,y[1]==1},y,{x,0.01,1},Method->{"Shooting","StartingInitialConditions"->{y[0]==0}}]

I get the error NDSolve::ntcs: Cannot solve constraint equations for initial conditions. <http://reference.wolfram.com/mathematica/ref/message/NDSolve/ntcs.html>


The solution is

y[x_]:=x+16*a*x^2*(x-1)^2
a=.3;

The following expression is zero for all values of x

-6*y'[x]^2*(q1[x,.3].q2'[y[x]]) +6*y'[x]^(5/2)*h1[y[x]]+y''[x]*(q1[x,.3].q2[y[x]])-y''[x]*y'[x]^(1/2)*h2[y[x]]/.x->.9


Does anyone know what is going wrong with NDSolve?

Charles Hagwood








Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.