Poiseuille flow in Elliptical channel

Hi, i want to plot 3D velocity field for Poiseuille flow in an elliptical channel Can you help me to plot it.

W(x,y)=a^2*(dp/dx)(1-(x/a)^2-(y/b)^2)/2*mu*(1+(a/b)^2)
Mu=1
a=2
b=1
Dp/dx=1

[Walter Roberson]

What do the underscores indicate here? Subscripts? 1 subscript (x/a)^2 ?

[Garland Garland]

I'm sorry That was aTypographical mistake It's (-)

Answer 1

[Torsten]

r = linspace(0,1,20);
phi= linspace(0,2*pi,36);
[R,PHI] = meshgrid(r,phi);
X=a*R.*cos(PHI);
Y=b*R.*sin(PHI);
V = -1/(2*mu)*dpdz*a^2*b^2/(a^2+b^2)*(1-(X/a).^2-(Y/b).^2);
surf(X,Y,V)

[Garland Garland]

Grateful