What is the solution of this? 2nd part

  • Thread starter Thread starter Clausius2
  • Start date Start date
Clausius2
Science Advisor
Gold Member
Messages
1,433
Reaction score
7
Hi Max,

Thanks very much, but now I am puzzled solving this problem:

F'''+FF''/2=0

F'(-infinite)=1
F(-infinite)=0
F'(infinite)=0

I am looking for a numerical solution, and I have been trying unsuccesfully to understand the non-linear shooting method, in order to code it with Matlab.

Be a good man, Max, and help an ODE's world wanderer.
 
Physics news on Phys.org
hi

Hi, Clausius2;

> I am looking for a numerical solution, and I have been trying unsuccesfully
> to understand the non-linear shooting method, in order to code it with
> Matlab.

Is it a goal to code it in Mathlab, or you just need to code it somewhere? And why non-linear shooting method? Why not standard RK4 (Runge-Kutta of the 4th order)?
To say the truth, I don't know Mathlab at all. Sorry for the stupid suggestion, but have you tried to find a ready Mathlab program in the Internet for any DE (system of DEs)? Try here for example: http://www.math.umn.edu/~olver/matlab.html.
Then, why not rewrite the original DE as three first-order DEs and solve that system? I did this, but it was 100 years ago and I programmed in QuickBasic.

> Be a good man, Max, and help an ODE's world wanderer.

I would like to, but I'm mostly interested in _analytic_ solutions. And this particular equation is _very_ interesting since they say it cannot be completely integrated analytically (http://mathworld.wolfram.com/BlasiusDifferentialEquation.html).
If I find smth related to your problem, I'll post it here.

Best of luck,
Max.
 
check this

Clausius2,
check this link for using Mathlab for solving of ODEs:
http://math.rice.edu/~polking/odesoft/dfpp.html,
and these ones for numerical solving of ODEs in general:
http://mathlab.cit.cornell.edu/math_software_resources/math_software_resources.html,
http://archives.math.utk.edu/software/msdos/diff.equations/.html.
Best of luck,
Max.

P.S. Well, man. It seems to me that I found the description of integration of the problem similar to yours. Read
Example: Using Continuation to Verify a Solution's Consistent Behavior
right here: http://www.mathworks.com/access/helpdesk/help/techdoc/math_anal/diffeq22.html#709310 .
Input beta=0, 1/2 in front of y*y'' and change the BCs and that seems to be it.
Check this link just in case: http://www.mathworks.com/access/helpdesk/help/techdoc/math_anal/math_anal_example_index.html .
 
Last edited by a moderator:
Last edited by a moderator:
Thanks

Thank you so much, Max.

I'm sure your links will be very useful. Now I am with my final exams, so I will come back to this problem on July.

What a pity you are not a woman, because I would ask you to marry with me. :smile:
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top