y=dsolve('Dy-3*y=6','x')
figure(1);hold on
for C2=-10:.1:10
x=[0:.1:1]
plot(x,C2*exp(3*x)-2)
end
或者直接用画向量场的命令.