Given the second order non linear BVP

(ƒ')^n = 1 + γ θ …..……………………… 1

θ" + (λ +n+1/ 2n +1) ƒ θ' - n (2 λ +1/2n+ 1) * ƒ' θ = 0 ………… 2

Prime in the above eqn’s describe partial differentiation with respect to η

Boundary conditions are

ƒ (0) = 0, θ'(0) = -1

ƒ' (∞) = 0, θ(∞) = 0

where η is a function of x & y given by

Similarity variable, η = x ^ (λ-n/2n+1) * y

Θ is a dimensionless temperature

ƒ is a dimensionless stream function given by

Stream function, ψ = x ^ (λ+n+1/2n+1) * ƒ(η)

how to solve the above system of PD eqn’s

the above system of equations are related to mixed convection of non newtonian fluids