A=[y x1 x2]
y=A(:,1);x1=A(:,2);x2=A(:,3);x =[x 1 x2];
fun=inline('c(1)+c(2)。*x(:,1)+c(3)。*x(:,2)',' c ',' x ');
c=lsqcurvefit(fun,rand(1,3),x,y);
c =
25.3689655159132 1.68965517250985 -0.965517241460331
y = 25.3689655159132+1.68965517250985 * x 1-0.965517241460331 * X2