function y= myfun(beta,x)
b1 = beta(1);
b2 = beta(2);
b3 = beta(3);
b4 = beta(4);
b5 = beta(5);

x1 = x(:,1);
x2 = x(:,2);
x3 = x(:,3);

y = (b1*x2 - x3/b5)./(1+b2*x1+b3*x2+b4*x3); %Ä£ÐÍ·½³Ì

%´ýÄâºÏµÄ×Ô±äÁ¿
x0=[
   470   300    10
   285    80    10
   470   300   120
   470    80   120
   470    80    10
   100   190    10
   100    80    65
   470   190    65
   100   300    54
   100   300   120
   100    80   120
   285   300    10
   285   190   120];
%´ýÄâºÏµÄº¯ÊýÖµ
y0=[
8.55
3.79
4.82
0.02
2.75
14.39
2.54
4.35
13
8.5
0.05
11.32
3.13];
%²ÎÊý³õÖµ
beta0=[1 0.05 0.02 0.1 2]';
%µ÷ÓÃ
beta = nlinfit(x0,y0,@myfun,beta0);
%Ïà¹ØϵÊý
x1=x0(:,1);
x2=x0(:,2);
x3=x0(:,3);
y=(beta(1)*x2 - x3/beta(5))./(1+beta(2)*x1+beta(3)*x2+beta(4)*x3);%¹À¼ÆÖµ
R=corrcoef(y0,y)