
function settore_ellittico

theta=linspace(pi/4,3*pi/4,1000);
thetaF=linspace(0,2*pi,1000);

A=[0.5 0]; B=[0 0.2];

Q=[A(1)*cos(theta') A(2)*cos(theta')]+[B(1)*sin(theta') B(2)*sin(theta')];
QF=[A(1)*cos(thetaF') A(2)*cos(thetaF')]+[B(1)*sin(thetaF') B(2)*sin(thetaF')];

grey=[240 240 240]/256;
plot(QF(:,1),QF(:,2),'-','color',grey,'LineWidth',2);
hold on;
Qpoly=[0 0; Q; 0 0];
plot(Qpoly(:,1),Qpoly(:,2),'k-','LineWidth',2);
fill(Qpoly(:,1),Qpoly(:,2),grey);
% axis tight
% axis square
axis off
hold off
