Hi, i want to plot stream tubes for the following problem, can anybody help me? Thanks

restart;
with(linalg);
Digits := 15;

velocity=[u(r,z),v(r,z),w(r,z)]

u(r,z)=cr, v(r,z)=r*omega, w(r,z)=0 at z=0

eq1:= 2*f(eta)+diff(h(eta),eta)
eq2 := (1/4)*(diff(h(eta),eta))^2-g(eta)*g(eta)-(1/2)*h(eta)*diff(h(eta),eta,eta)+(1/2)*diff(h(eta),eta,eta,eta);

eq3 := -g(eta)*diff(h(eta),eta)+h(eta)*diff(g(eta),eta)-diff(g(eta),eta,eta);

where f(eta), g(eta) and h(eta) are the solution of above system 

bc:=f(0)=1,f (N)=0,h(0)=0,D(h)(N)=0,D(h)(0)=-2,g(0)=R,g(N)=0; where R is constant