No BSD License
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