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
