% polfig3.m   96/07/21
%
% Create a figure in polar coordinates with
% h_1(theta)=<r<=h_2(theta) and c<=theta<=d
%
    con=pi/180;
    c=20*con;
    d=70*con;
    theta=c+(d-c)*(0:20)/20;
    h1=1+.3*cos(3*theta);
    h2=2+.3*sin(3*theta);
    lh=length(h1);
    r1=h1(1)+(h2(1)-h1(1))*(0:20)/20;
    r2=h1(lh)+(h2(lh)-h1(lh))*(0:20)/20;

    hold off;
    plot(h1.*cos(theta),h1.*sin(theta),'w-');
    axis([-.5 2.5 -.5 2.5]);
    hold on
    axis square;
    plot(h2.*cos(theta),h2.*sin(theta),'w-');
    plot(r1*cos(c),r1*sin(c),'w-');
    plot(r2*cos(d),r2*sin(d),'w-');

    plot([-.5 5],[0 0],'w-');
    plot([0 0],[-.5 5],'w-');

% Text

%    xlabel('x');
%    ylabel('y');
    text(2,-.15,'x');
    text(-.2,2,'y');
    text(.1,.4,'r=h_1(theta)');
    text(1.5,.4,'theta=c');
    text(1.7,1.6,'r=h_2(theta)');
    text(-.1,1.3,'theta=d');
