y=dsolve('D3y=D2y-Dy-y+t^2','y(0)=0','Dy(0)=1','D2y(0)=-1');
tt=linspace(0,10,1001); %表示t从0到10
for i=1:1001
t=tt(i);
yy(i)=eval(y);
end
plot(tt,yy)亲……求解释…………看不懂…………而且，答案跟用你程序运行出来的图片不一样你的答案t的范围是什么我就取了0到10的1001个等距点然后每个点求y的值然后画出来，有什么问题么？为嘛不能用ezplot(y)呢…………因为y的表达式太长。。。