function tanxInvxPlotAnon(a,xmin,xmax,n)
% tanxInvxPlotAnon  Plot y = tan(a*x) - 1/x

if nargin<1, a = 1.0;       end
if nargin<2, xmin = -2*pi;  end
if nargin<3, xmax =  2*pi;  end
if nargin<4, n = 500;       end

% -- Define anonymous function
tanfun = @(x,a) tan(a*x) - 1./x;

x = linspace(xmin,xmax,n);
y = tanfun(x,a);
plot(x,y)

end