Question: plotting a system of ode?

eq1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)+beta*(1-diff(f(eta),eta)^2)=0;

eq2:=diff(theta(eta), eta, eta)+pr*f(eta)*(diff(theta(eta), eta))+pr*Ec*(diff(f(eta), eta, eta))^2

+(2-beta)*pr*Q*theta(eta) = 0;

bc:=theta(9)=0,D(theta)(0)=-Nc*(1-theta(0)),f(0)=1,D(f)(0)=lambda,D(f)(9)=1};

Para:={beta=1,Ec=0.5,Q=-0.5,Nc=5,pr=3}:

How to plot -D(theta)(0) vs lambda?

Thanks

Please Wait...