Programación de ecuaciones diferenciales de segundo orden (como circuitos) con integración trapezoidal implícita en matlab

gt;gt; claro gt;gt; syms a b c d e; gt;gt; y=dsolve('a*D2y b*Dy c*y=0','y(0)=d', 'Dy(0)=e') y = (2*a*e b*d d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a *e b*d - d*(b^2 - 4*a*c)^(1/2))/(...7795