function [id, vt] = diodo(tmax, Vmax, f, R)
  % tmax => tiempo máximo
  % Vmax => Voltaje máximo
  % f    => frecuencia de la señal
  % R    => Resistencia

  t = [0:0.01:tmax];         % tiempo
  vt = Vmax*cos(f*t);        % Voltaje variante en el tiempo
  vd = vt;                   % Voltaje inicial en el diodo
  vd(vd>0.7) = 0.7;          % Voltaje en el diodo
  id = (vt - vd)/R;          % Calculo de la corriente
end
