function [y,yd,H]=fun7(x) x1 = x(1); x2 = x(2); x3 = x(3); y= 2500*(sqrt(4+x1.^2))+1200*(sqrt(1+(x2-x1).^2))+2000*(sqrt(0.5^2+(x3-x2).^2))+1000*(10-x3); if nargout>1 yd(1,1) =(600.*(2.*x1 - 2.*x2))./((x1 - x2).^2 + 1).^(1/2) + (2500.*x1)./(x1.^2 + 4).^(1/2); yd(2,1) = (1000.*(2.*x2 - 2.*x3))./((x2 - x3).^2 + 1/4).^(1/2) - (600.*(2.*x1 - 2.*x2))./((x1 - x2).^2 + 1).^(1/2); yd(3,1) = -(1000.*(2.*x2 - 2.*x3))./((x2 - x3).^2 + 1/4).^(1/2) - 1000; end if nargout>2 H(1,1) = 2500./(x1.^2 + 4).^(1/2) - (300.*(2.*x1 - 2.*x2).^2)./((x1 - x2).^2 + 1).^(3/2) + 1200./((x1 - x2).^2 + 1).^(1/2) - (2500.*x1.^2)./(x1.^2 + 4).^(3/2); H(1,2) = (300.*(2.*x1 - 2.*x2).^2)./((x1 - x2).^2 + 1).^(3/2) - 1200./((x1 - x2).^2 + 1).^(1/2); H(1,3) = 0; H(2,1) = (300.*(2.*x1 - 2.*x2).^2)./((x1 - x2).^2 + 1).^(3/2) - 1200./((x1 - x2).^2 + 1).^(1/2); H(2,2) = 1200./((x1 - x2).^2 + 1).^(1/2) - (500.*(2.*x2 - 2.*x3).^2)./((x2 - x3).^2 + 1/4).^(3/2) - (300.*(2.*x1 - 2.*x2).^2)./((x1 - x2).^2 + 1).^(3/2) + 2000./((x2 - x3).^2 + 1/4).^(1/2); H(2,3) = (500.*(2.*x2 - 2.*x3).^2)./((x2 - x3).^2 + 1/4).^(3/2) - 2000./((x2 - x3).^2 + 1/4).^(1/2); H(3,1) = 0; H(3,2) = (500*(2*x2 - 2*x3)^2)/((x2 - x3)^2 + 1/4)^(3/2) - 2000/((x2 - x3)^2 + 1/4)^(1/2); H(3,3) = 2000/((x2 - x3)^2 + 1/4)^(1/2) - (500*(2*x2 - 2*x3)^2)/((x2 - x3)^2 + 1/4)^(3/2); end