%========================================================================
% Calcula a integral numérica da expressão para cálculo do tempo de 
% retomada.
% Ref. Nicolazzi, L.C.. "Uma introdução à modelagem quase-estática 
%           de automóveis".  Publicação interna do GRANTE, Departamento
%           de Engenharia Mecânica da UFSC, 2012.
%
%                                                criada em 26/09/2017
%                         
%========================================================================

function tempo = func_tempo_retomada(const1,const2,const3,const4,marcha,ic,etam,delta,vr1,vr2,vincr)

soma=0;
icaixa=ic(marcha);
deltacaixa=delta(marcha);
K=const1*icaixa;
A0=-1.278493e+05;
A1=2.195991e+02;
A2=-1.248258e-01;
A3=6.622842e-05;
A4=-2.172773e-08;
A5=3.559221e-12;
A6=-2.289778e-16;
for vel=vr1:vincr:vr2
    PLiq=etam*(A0+A1*K*vel+A2*K^2*vel^2+A3*K^3*vel^3+A4*K^4*vel^4+...
        +A5*K^5*vel^5+A6*K^6*vel^6)-const3*vel-const2*vel^3;
    soma = soma+vel/PLiq*vincr;
end
tempo=const4*(1+deltacaixa)*soma;


