Solving boundary layer problems using matlab numerically

Actually i am solving the set of nonlinear ode in matlab with bvp4c solver. For example, the governing eqns are $$ s''(\eta)-\theta'(\eta)+Nr*f'(\eta)=0,$$ $$\theta''(\eta)+s(\eta)*\theta(\eta)+Nt*f'(\eta)*\theta'(\eta)+Nb*(\theta'(\eta))^2 =0 $$ $$f''(\eta)-(Nb/Nt)*s'(\eta)=0 $$ with the asymptotic boundary conditions $$s(0)=0, \theta(0)=1, f(0)=1$$ $$s'(\infty)=0,\theta(\infty)=0, f(\infty)=0$$ After solving these equations with bvp solver i plotted the graphs of $s,s',\theta,f $ versus $\eta$. But i need to plot the graphs of $s,s',\theta,f$ versus each parameter $Nr,Nt,Nb$. How will i put the command in bvp4c command for this. I tried in some other way like loop for $Nr$ in that bvp block. The loop takes only first value of $Nr$. If you know to plot some other way help me. thanks.

[Walter Roberson]

To be honest, those equations are unreadable.

[Prakash Dhanapal]

sir, such a type of equations like commands in latex version. The equations are example of some system of nonlinear equations.

$$ s''(\eta)-\theta'(\eta)+Nr*f'(\eta)=0,$$

$$\theta''(\eta)+s(\eta)*\theta(\eta)+Nt*f'(\eta)*\theta'(\eta)+Nb*(\theta'(\eta))^2 =0 $$

$$f''(\eta)-(Nb/Nt)*s'(\eta)=0 $$.

i want to plot $\theta$ versus $Nr$. if you know help me.

[Walter Roberson]

I would have to study LaTeX first before I could understand those equations. How about you convert them to graphical form and post the mathematical typeset version of the equations to a site such as tinypic ?