exn527.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 16 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-2-7
% 求空间曲面的极值点
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月

clear,clf,format compact
[x,y]=meshgrid(-4:.1:4);
f = x.^3-y.^3+3*x.^2+3*y.^2-9*x;    % 曲面方程
mesh(x,y,f)                         % 画出三维曲面
xm=[-3,-3,1,1];ym=[0,2,0,2];        % 极点坐标
fm = xm.^3-ym.^3+3*xm.^2+3*ym.^2-9*xm,
fxm = 3*xm.^2+6*xm-9,               % 各极点处的偏导数
fym = -3*ym.^2+6*xm,
fxxm = 6*xm+6, 
fyym = -6*ym+6,fxym = 0,
dis = fxxm.*fyym-fxym.^2		    % 极点特性判别式