clc;clear;
cd d:\data1;
load d:\data1\shuju.txt;
save shuju;
ziben=shuju(:,1);
laodong=shuju(:,2);
nengyuan=shuju(:,3);
chanzhi=shuju(:,4);
paifang=shuju(:,5);
x=[ziben,laodong,nengyuan];
c=[ones(29,1),x(:,1),x(:,2),x(:,3),0.5*x(:,1).*x(:,1),0.5*x(:,1).*x(:,2)...
,0.5*x(:,1).*x(:,3),0.5*x(:,2).*x(:,1),0.5*x(:,2).*x(:,2),0.5*x(:,2).*x(:,3)...
,0.5*x(:,3).*x(:,1),0.5*x(:,3).*x(:,2),0.5*x(:,3).*x(:,3),chanzhi,0.5*chanzhi.*chanzhi...
,paifang,0.5*paifang.*paifang,chanzhi.*paifang,x(:,1).*chanzhi,x(:,2).*chanzhi...
,x(:,3).*chanzhi,x(:,1).*paifang,x(:,2).*paifang,x(:,3).*paifang];
c=sum(c,1);
A1=-1.*[ones(29,1),x(:,1),x(:,2),x(:,3),0.5*x(:,1).*x(:,1),0.5*x(:,1).*x(:,2)...
,0.5*x(:,1).*x(:,3),0.5*x(:,2).*x(:,1),0.5*x(:,2).*x(:,2),0.5*x(:,2).*x(:,3)...
,0.5*x(:,3)