function y = sinplot(xl,xh)
% sinplot(xl,xh) plots sin(x) for xl < x < xh
x = (xh-xl).*(0:.01:1);   % 101 data point, including end points
y = sin(xl + x);
plot (x,y);
grid on
xlabel ('Angle (radians)')
ylabel ('Sine of angle')
title  ('Sine plot')
return
