exm0663_2.m

来自「matlab程序」· M 代码 · 共 11 行

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