function funplot(f, lims)
% funplot(f, [left,right]) graph a function f(x)
% in the domain left <= x <= right
% A poor substitute for the built-in FPLOT
xpts = linspace(min(lims), max(lims), 100);
ypts = feval(f,xpts); 
plot(xpts, ypts);
