Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Mathematica Not Solving

  1. Jun 16, 2015 #1

    joshmccraney

    User Avatar
    Gold Member

    Hi PF!

    I am wondering if any of you can explain why Mathematica won't solve this ODE? I know an analytic solution exists (it's quadratic). Any help is greatly appreciated! $f0[x]$ is the name of the function.

    Code (Text):
     
    a = -1/5;
    DSolve[{-((2 a f0[x])/(1 + a)) + x f0'[x] +
        2 (f0'[x])^2 + f0[x] f0''[x] == 0,
      f0[1] == 0, f0'[1] == -1/2}, f0[x], x]
     
  2. jcsd
  3. Jun 19, 2015 #2
    Some hints to improve your post.
    Give us the exact differential equation in LaTeX. having to "translate" hardly readable mathematica code doesn't help getting answers.
    You say you know an analytic solution exists, do you know the form as well?

    More on topic, have you tried NDSolve to see if mathematica can work with the DE and gives a reasonable result?
     
  4. Jun 19, 2015 #3

    joshmccraney

    User Avatar
    Gold Member

    $$f f''+2f'^2+x f'+\frac{1}{2}f = 0$$ where ##f## is a function of ##x## subject to ##f(1)=0## and ##f'(1)=-1/2##. The analytic solution is ##f=(1-x^2)/2##.

    I have tried NDSolve and it works great, provided I perturb the boundary condition as ##f(1)=0.00001##.

    Any ideas why it can't figure out the solution?
     
  5. Jun 19, 2015 #4
    So I've been looking at this for a bit.

    First I would get rid of the factor of ##f## from the second derivative term.
    From that I would be led to believe you could consider ##g = f^\prime/\sqrt{f}## as a new variable.

    ##g^\prime = \frac{f^{\prime\prime}}{\sqrt{f}} - \frac{f^\prime}{f}##
    Then ##f^{\prime\prime} = \left(g^\prime+\frac{f^\prime}{f}\right)\sqrt{f}##
    So we can rewrite your DE as

    ##\sqrt{f}\left(g^\prime + g\right) + 2 \left(g^\prime\right)^2+xg+\frac{1}{2} = 0##

    That's where I got stuck, maybe this gives you some ideas.
    I think the key is rewriting the equation such that mathematica recognizes it.
    Also wolfram alpha gave me a solution (for different initial values) although no analytical expression.
    Since wolfram alpha uses mathematica in the back this is promising.
     
  6. Oct 6, 2015 #5
    Are u sure the solution of that differential equation is really [itex] f=(1-x^2)/2[/itex]? I didn't get it because the initial condition of this solution is not same as you mentioned and actually it is [itex] f'(1)=-1[/itex].

    I think the factor [itex]1/2[/itex] is missing. However, still I can't get the solution same as you. That's interesting.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Mathematica Not Solving
Loading...