# SECTION 5.1de1 := diff(x(t),t) = y(t)/12 - x(t)/3;de2 := diff(y(t),t) = x(t)/3 - y(t)/3;solve(de2,x(t));subs(x(t)=%,de1);expand(%); denew := lhs(%)-rhs(%)=0;# Solve by the methods of Section 4.2:P := 3*r^2+2*r+1/4;rsol := solve(P=0,r);r1 := rsol[1]; r2 := rsol[2];y1 := exp(r1*t); y2 := exp(r2*t);ysol := c1*y1+c2*y2;xsol := 3*diff(ysol,t)+ysol;subs(t=0,xsol);eq1 := % = x0;subs(t=0,ysol);eq2 := % = y0;solve({eq1,eq2},{c1,c2});subs(%,{xsol,ysol});