% polfig1.m   96/07/21
%
% Create a figure in polar coordinates with
% a<=r<=b and c<=theta<=d
%
    con=pi/180;
    c=20*con
    d=70*con
    theta=c+(d-c)*(0:20)/20;
    a=1;
    b=2;

    hold off;
    plot([a*cos(c) b*cos(c)],[a*sin(c) b*sin(c)],'w-');
    hold on
    axis([-.5 2.5 -.5 2.5]);
    axis square;
    plot([a*cos(d) b*cos(d)],[a*sin(d) b*sin(d)],'w-');
    plot([-.5 5],[0 0],'w-');
    plot([0 0],[-.5 5],'w-');
    xa=a*cos(theta); ya=a*sin(theta);
    xb=b*cos(theta); yb=b*sin(theta);
    plot(xa,ya,'w-')
    plot(xb,yb,'w-')

% Text

%    xlabel('x');
%    ylabel('y');
    text(1.5,.4,'theta=c');
    text(1.4,1.6,'r=b');
    text(0.5,0.6,'r=a');
    text(.05,1.5,'theta=d');
    text(2,-.15,'x');
    text(-.2,2,'y');
