image thumbnail
from Transient Pipe Flow by Housam Binous
Compute the transient velocity profile of a flow in a pipe.

f=un(x,n)
function f=un(x,n)

global zer n

G=-8;
nu=0;

f=-G/4*(1-x.^2).*besselj(nu,zer(n)*x).*x;

end

Contact us at files@mathworks.com