function di=didt(t,i,R,L,v,f,Rsec)

i1=i(1); i2=i(2); A=L; D=L*R; w=2*pi*f;

di(1,1)=A(1,1)*(v*cos(w*t))+A(1,2)*Rsec*i2+A(1,3)*((D(3,1)*i1+D(3,2)*i2...
    -A(3,1)*(v*cos(w*t))-A(3,2)*Rsec*i2)/A(3,3))-D(1,1)*i1-D(1,2)*i2;
di(2,1)=A(2,1)*(v*cos(w*t))+A(2,2)*Rsec*i2+A(2,3)*((D(3,1)*i1+D(3,2)*i2...
    -A(3,1)*(v*cos(w*t))-A(3,2)*Rsec*i2)/A(3,3))-D(2,1)*i1-D(2,2)*i2;

end

% Definição dos parâmetros
R=[2.1581 0 0; 0 0.0114 0; 0 0 0.0556];
L=[8.7492 -8.7828 0.0336; -8.7828 12.3082 -3.5254; 0.0336 -3.5254 3.4918];
f=1e3; v=10e3; Rsec=100; i10=0; i20=0;

% Tempo de simulação
ts=[0 10e-3];

[t,i]=ode45(@(t,y) didt(t,y,R,L,v,f,Rsec),ts ,[i10 i20])
figure(1)
plot(t,i(:,1), t,i(:,2))