function f=fun1(Pos,x)
%ÆäÖÐPos(1)=k1,Pos(2)=k2.....;

f=Pos(1).*(2.52125.*10.^(-6).*x(:,1).^(2)-0.00815.*x(:,1))+Pos(2).*(5.11549.*10.^(-6).*x(:,2).^(4)-4.43569.*10.^(-4).*x(:,2).^(3)+0.01157.*x(:,2).^(2)-0.05903.*x(:,2))+Pos(3).*(-66538.82632.*x(:,3).^(6)+68596.56918.*x(:,3).^(5)-25259.909.*x(:,3).^(4)+3778.81796.*x(:,3).^(3)-138.15946.*x(:,3).^(2)-13.21417.*x(:,3))+Pos(4);
end