Using NDSolve
while solving system of ODE's. Here is my try
Eqn1 = -C f'''''[x] + f'''[x] + f[x] f''[x] - f'[x] f'[x] + r^2 + λ T[x] == 0
Eqn2 =
T''[x] + p1 f[x] T'[x] +
p2/(3 (1 - tw)) (12 (tw + (1 - tw) T[x] ) + tw)^2 (1 - tw) T'[x] T'[x] +
4 ((1 - tw) T[x] + tw)^3 (1 - tw) T''[x] == 0
BC1 = f[0] == 0
BC2 = f'[0] == 1
BC3 = f'[10] == r
BC4 = f''[10] == 0
BC5 = f'''[10] == 0
BC6 = T'[0] == -B (1 - T[0])
BC7 = T[10] == 0
param = {λ -> 1, p1 -> 10, B -> 10, tw -> 1.5, p2 -> 0.5,
C -> 1, r -> 1}
Sol =
NDSolve[{Eqn1, Eqn2, BC1, BC2, BC3, BC4, BC5, BC6, BC7} /. param,
{f, T}, {x, 0, 10},
Method ->
{"Shooting",
"StartingInitialConditions" ->
{f[0] == 0, f'[0] == 0, f''[0] == 1, f'''[0] == 0, f''''[0] == 0, T[0] = 0}}]
But I'm getting this error
NDSolve::ndode: Input is not an ordinary differential equation.