Hello !

I want to plot its steamlines and Isotherms for different parameters. Anyone can help me for this issue.

eq1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)-(diff(f(eta),eta)^2)-M*diff(f(eta),eta)=0:
eq2:=(1+R)*1/Pr*diff(theta(eta),eta,eta)+f(eta)*diff(theta(eta),eta)=0:
bc:=D(f)(0)=1,(f)(0)=0,(D)(f)(N)=0,theta(0)=0,theta(N)=1:
Pr:=1:R:=0.5:M:=0.5:N:=10: