Advertisement
Guest User

Untitled

a guest
Sep 22nd, 2017
53
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.75 KB | None | 0 0
  1. s = NDSolve[{f'''[eta] + 0.5*f[eta]*f''[eta] == 0.0, f[0] == 0.0,
  2. f'[0] == 0.0, f'[Infinity] = 1.0}, f, {eta, 0, 1}];
  3. Plot[Evaluate[f[eta] /. s], {eta, 0, 1}, PlotRange -> All]
  4.  
  5. NDSolve::deqn: Equation or list of equations expected instead of 1.` in the first argument
  6.  
  7. s = NDSolve[{Derivative[3][f][x] + 1/2 f[x] Derivative[2][f][x] == 0,
  8. f[0] == 0, f'[0] == 0, f'[#] == 1}, f, {x, 0, 1}] & /@ Range[1, 50, 5];
  9.  
  10. Plot[Evaluate[f[eta] /. s], {eta, 0, 1}, PlotRange -> All,
  11. PlotLegends -> (ToString[#] & /@ Range[1, 50, 5])]
  12.  
  13. sol = NDSolve[{f'''[η] + 0.5 f[η] f''[η] == 0,
  14. f[0] == f'[0] == 0, f'[10] == 1}, f, η]
  15. Plot[f[η] /. First[sol], {η, 0, 10}]
  16. Plot[f'[η] /. First[sol], {η, 0, 10}]
  17. Plot[f''[η] /. First[sol], {η, 0, 10}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement