grossk.m

clear;
%einlesen
A = xlsread('data.xls');
[m,n] = size(A);
%kreis
phi=-pi:pi/100:pi;
x=sqrt(2)*cos(phi);
y=sqrt(2)*sin(phi);
plo=plot(x,y,'k');
hold on;
%grosskreise
for i=1:m;
    n=A(i,1);
    az=-deg2rad(A(i,2))+pi/2;
    a=deg2rad(A(i,3));
    d=-pi/2:pi/100:pi/2;
    x=cos(az-atan(sin(d)./(cos(a).*cos(d)))).*(2-2.*sin(a).*cos(d)).^(1/2);
    y=sin(az-atan(sin(d)./(cos(a).*cos(d)))).*(2-2.*sin(a).*cos(d)).^(1/2);
    plo=plot(x,y,'k');
    hold on;
    %for text set
    x=cos(az)*(2-2.*sin(a)).^(1/2);
    y=sin(az)*(2-2.*sin(a)).^(1/2);
    plo=plot(x,y,'k','Marker','+');
    text(x+0.02,y,num2str(n),'FontSize',13);
    hold on;
end
%format
    %nord
    x=[0;0];
    y=[sqrt(2);1.5];
    plo=plot(x,y,'k');
    text(-.15,1.6,' N','FontSize',20);
axis square;
axis off;
saveas(plo,'grossk','eps');