function dy = rigid_diodo_tunel(t,y)
%Definição das Variáveis
c=2; r=1.5;l=5;
u=1.2;
dy = zeros(2,1);    % a column vector
hy1=17.76*y(1)-103.79*y(1).^2+229.62*y(1)^3-226.31*y(1)^4+83.72*y(1)^5;

dy(1)= (1/c)*(-hy1+y(2));
dy(2)=(1/l)*(-y(1)-r*y(2)+u);