clear all
close all
clc
%-------------------------------------------------------------
%% dados numéricos
R=0.25;     % Comprimento da manivela
C=0;        % Distancia do pistao ao centro da manivela
L=0.45;     % comprimento da biela   
NP=50;
qin = 0*pi/180;
qfin = 360*pi/180;
%-------------------------------------------------------------
%% posição
%-------------------------------------------------------------
% solução numérica (precisa da função "funeqcin")
%-------------------------------------------------------------
% valores iniciais das variáveis secundárias
varsec0 = [0*pi/180; R+L];  
for i=1:NP
    q(i)=(i-1)*(qfin-qin)/(NP-1) + qin;
    const=[R; L; C; q(i)];      % valores constantes e variável primária
    options = optimset('Display','off'); 
    % chamada da função que resolve o sistema de equações de posição
    [varsec,fval]=fsolve(@(varsec)funeqcin(varsec,const),varsec0,options);
    % distribuição do resultado nas variáveis secundárias
    a(i)=varsec(1);
    x(i)=varsec(2);
end

plot(q,a);
