clear all
close all
clc

Rb=50;
Rmax=60;
NP=100;
for i=1:NP
    R(i)=(i-1)/(NP-1)*(Rmax-Rb)+Rb;
    Af(i)=acos(Rb/R(i));
    B(i)=tan(Af(i))-Af(i);
end
%%
figure('name','Evolvente')
polar(B,R)
grid on