% modulo cisalhamento complexo material
% viscoelastico modelo 4 parametros
% af fator de deslocamento
% pv(1)=g0 modulo cisalhamento inferior
% pv(2)=goo modulo cisalhamento superior
% pv(3)=bt inclinacao
% pv(4)=b1 deslocamento
% ww frequencia de trabalho (rad/s)
function r=g_cmplx(af,pv,ww)
  % inicio
  nr=(pv(1)+pv(2)*pv(4)*(1i*af*ww)^pv(3));
  dn=(1+pv(4)*(1i*af*ww)^pv(3));
  r=nr ./ dn;
  % fim
end