function y=punish1(x,r)
syms x1 x2 x3;
z=1000-x1^2-2*x2^2-x3^2-x1*x2-x1*x3+r*((x1^2+x2^2+x3^2-25)^2+(8*x1+14*x2+7*x3-56)^2);
y=[diff(z,x1);diff(z,x2);diff(z,x3)];
x1=x(1);
x2=x(2);
x3=x(3);
y=subs(y);