Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Eqn1 = f'''[x] + f[x] f''[x] + ((2 n)/(n + 1))(1 - f'[x] f'[x])- M f'[x]==0
- Eqn2 = T''[x] + Prf[x] T'[x]-Pr ((2 p)/(n + 1)) f'[x] T[x] + (2 /(n + 1))[A
- f'[x] + B T[x]] == 0
- BC1 = f[0] == 0;
- BC2 = f'[0] == λ + β f''[0];
- BC3 = f'[inf1] == 1;
- BC4 = T[0] == 1 + σ T'[0];
- BC5 = T[inf1] == 0;
- param1 = {n -> 0, M -> 0, Pr -> 1, p -> 5.29387, A -> -0.05,
- B -> -0.05, λ -> 0.5, β -> 0.5, σ -> 0.5};
- inf1 = 1.5;
- Sol1 = NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5} /. param1, {f,
- T}, {x, 0, inf1},
- Method -> {"Shooting","StartingInitialConditions" -> {f[0]==0,f'[0] == 0,
- f''[0] == 0, T[0] == 0, T'[0] == 0}}];
- param2 = {n -> 1, M -> 0, Pr -> 1, p -> 5.29387, A -> -0.05,
- B -> -0.05, λ -> 0.5, β -> 0.5, σ -> 0.5};
- Sol2 = NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5} /. param2, {f,
- T}, {x, 0, inf1}, Method -> {"Shooting","StartingInitialConditions" -> {f[0] == 0, f'[0] == 0,
- f''[0] == 0, T[0] == 0, T'[0] == 0}}]
- param3 = {n -> 1.5, M -> 0, Pr -> 1, p -> 5.29387, A -> -0.05,
- B -> -0.05, λ -> 0.5, β -> 0.5, σ -> 0.5};
- Sol3 = NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5} /. param3, {f,
- T}, {x, 0, inf1},
- Method ->{"Shooting","StartingInitialConditions" -> {f[0] == 0,f'[0] == 0,
- f''[0] == 1, T[0] == 0, T'[0] == 0}}]
- Plot[{f'[x] /. Sol1, f'[x] /. Sol2, f'[x] /. Sol3, f'[x]}, {x, 0,
- inf1}, PlotRange -> All, AxesLabel -> {η, f' (η)},
- PlotStyle -> {Black, Red, Green, Blue, Yellow}, Frame -> True,
- FrameStyle -> Directive[Black, Bold, 12], PlotRange -> All,
- Axes -> False, FrameLabel -> {η, f'}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement