exm05063_2.m

来自「这本书是matlab经典书籍」· M 代码 · 共 7 行

y=dsolve('y=x*Dy-(Dy)^2','x')							%求微分方程解
clf,hold on,ezplot(y(2),[-6,6,-4,8],1)				%画奇解
cc=get(gca,'Children');										%取奇解曲线的图柄	<3>
set(cc,'Color','r','LineWidth',5)						%把奇解画成粗红线	<4>
for k=-2:0.5:2;ezplot(subs(y(1),'C1',k),[-6,6,-4,8],1);end%画通解
hold off,title('\fontname{隶书}\fontsize{16}通解和奇解')