ea520.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 14 行

% 《工程线性代数（MATLAB版）》第五章例题5.20程序ea520
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 二次型方程插值问题
%
x=[-1;0;2;2;0];y=[0;2;3;-2;-3]; 
A=[x.^2,x.*y,y.^2,x,y],
b=-ones(5,1)
K=A\b
pause
[X,Y]=meshgrid(-4:0.1:4);
ezplot('2*X.^2-X.*Y+Y.^2-4*X+Y-6*ones(size(X))=0'),hold on
plot(x,y,'s'),grid on
axis([-4 4 -4 4])
set(gcf,'color','w')