function [aa,bb]=cerchi
% Outputs
% aa,bb: coefficienti angolari delle linee
%        delle caratteristiche

c=[0 0]; % centro

%------- plot dei cerchi
n=10;
r=1:n;    % raggi
x=linspace(0,2*pi,200);
plot(c(1)+r(1)*cos(x),c(1)+r(1)*sin(x),'-r','LineWidth',1)
axis equal
hold on
for i=2:n
  plot(c(1)+r(i)*cos(x),c(1)+r(i)*sin(x),'-b','LineWidth',1)
end

delta=pi/12;
% angoli della suddivisione del vino
a=[pi+delta,pi+2*delta,-2*delta,-delta, pi/2-delta/2,pi/2+delta/2];
% angoli della suddivisione del cibo
b=a+pi;

%--------------------------------------------------
% plot dei segmenti per le caratteristiche del vino
%---------------------------------------------------
% punti iniziali
xin=c(1)+r(1)*cos(a);
yin=c(2)+r(1)*sin(a);

% punti finali
strW=strcat({'alcohol level ';'tannin content ';'sapidity/effer.';'acidity ';...
       'aroma intensity/I.A.P.'; 'sweetness/softness'});

for i=1:length(a)    
 xfin(i)=xin(i)+r(n)*cos(a(i));
 yfin(i)=yin(i)+r(n)*sin(a(i));
 
 %text(xfin(i),yfin(i),strW(i),'FontSize',11);
 
end

xx=[xin; xfin];
yy=[yin; yfin];

ww=plot(xx,yy, '-r','LineWidth',2.5);
axis auto
%--------------------------------------------------
% plot dei segmenti per le caratteristiche del cibo
%---------------------------------------------------
% punti iniziali
xinc=c(1)*ones(1,6);
yinc=c(2)*ones(1,6);

% punti finali
xfinc=xinc+(r(n)+1)*cos(b);
yfinc=yinc+(r(n)+1)*sin(b);

xxc=[xinc; xfinc];
yyc=[yinc; yfinc];

ff=plot(xxc,yyc,'-b','LineWidth',2.5);

%------------------------------------------
legend('Wine Character.','Food Character.',2);

% angols of the lines of the characteristics 
for i=1: length(xx)
 aa(i)=2*atan(yy(1,i)./(sqrt(xx(1,i)^2+yy(1,i)^2)+xx(1,i)));
 bb(i)=2*atan(yyc(2,i)./(sqrt(xxc(2,i)^2+yyc(2,i)^2)+xxc(2,i)) );
end
hold off
end





