exm056_2.m

来自「这是一本介绍MATLAB的教程」· M 代码 · 共 10 行

%exm056_2.m

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