>> syms t x(t) y(t); uu=diff(sin(t))+3*sin(t);
>> y1=diff(y);y2=diff(y,2);y3=diff(y,3);y4=diff(y,4);y5=diff(y,5);
>> z1=dsolve(y5+5*y4+12*y3+16*y2+12*y1+4*y==uu,y(0)==0,y1(0)==0,y2(0)==0,y3(0)==0,y4(0)==0);
>> z1=simplify(z1)
原文：https://www.cnblogs.com/Math-Nav/p/13375208.html