%首先定義非線性函數的M文件mymodel.m:
function yy=model(beta0,x)
a=beta0(1);b=beta0(2);c=beta0(3);d=beta0(4);e=beta0(5);f=beta0(6);
x1=x(:,1);x2=x(:,2);
yy=a.*x1.^2+b.*x1+c.*x2.^2+d.*x2+e*x1.*x2+f;