clear all
close all
clc

C=11;
R=5;

NP = 361;           
qmin = 0*pi/180;   
qmax = 360*pi/180; 

for i = 1:NP
    q(i)=qmin+(i-1)*(qmax-qmin)/(NP-1); 
    A(i)=atan(R*sin(q(i))/(C+R*cos(q(i))));
    MFL(i)=(cos(A(i))^3*(C*R*cos(q(i))+R^2))/((C+R*cos(q(i)))^2);
end
plot(q*180/pi,MFL)
grid on
xlabel('q (^o)')                        % rótulo do eixo x
ylabel('M/FL')                    % rótulo do eixo y          <**>
title('Posição do elo de saida')        % título do gráfico
xlim([qmin*180/pi qmax*180/pi])