Group Abstract Group Abstract

Message Boards Message Boards

0
|
4K Views
|
1 Reply
|
0 Total Likes
View groups...
Share
Share this post:

Solve DE by fourth order Runge Kutta code coupled with the shooting method

seema kumari
Posted 10 years ago

how to solve this problem in mathematica, using a forth order Runge Kutta code coupled with the shooting method. A ordinary differential equation is

f''''[y] - 2*a^2*f''[y] - a^2*(a^2 - Ra)*f[y] == 0,

with boundary condition is,

f[0] == 0, 

f'''[0] - a^2*f'[0] == 0, f[1] == 0, f''[1] == 0

?
POSTED BY: seema kumari
1 Reply
Sort By:

You can solve the boundary problem numerically with NDSolve:

With[{a = 2, R = 1}, 
 sol = NDSolveValue[{f''''[y] - 2 a^2 f''[y] - a^2 (a^2 - R a) f[y] == 0, 
   f[0] == 0, f'''[0] - a^2*f'[0] == 0, f[1] == 0, f''[1] == 0}, f, {y,
    0, 1}, Method -> {"Shooting"}];
 Plot[sol[y], {y, 0, 1}]]
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard