function y = PrC1c(x)
% Matlab Code by A. Hedar (Nov. 23, 2005).
% Parameters
P = 6000; L = 14; E = 30e+6; G = 12e+6;
t_max = 13600; s_max = 30000; d_max = 0.25;
M = P*(L+x(2)/2); R = sqrt(0.25*(x(2)^2+(x(1)+x(3))^2));
J = 2/sqrt(2)*x(1)*x(2)*(x(2)^2/12+0.25*(x(1)+x(3))^2);
P_c = (4.013*E/(6*L^2))*x(3)*x(4)^3*(1-0.25*x(3)*sqrt(E/G)/L);
t1 = P/(sqrt(2)*x(1)*x(2)); t2 = M*R/J;
t = sqrt(t1^2+t1*t2*x(2)/R+t2^2);
s = 6*P*L/(x(4)*x(3)^2);
d = 4*P*L^3/(E*x(4)*x(3)^3);
% Constraints
y(1) = t-t_max;
y(2) = s-s_max;
y(3) = x(1)-x(4);
y(4) = 0.10471*x(1)^2+0.04811*x(3)*x(4)*(14.0+x(2))-5.0;
y(5) = d-d_max;
y(6) = P-P_c;
% Variable lower bounds
y(7) = -x(1)+0.125;
y(8) = -x(2)+0.1;
y(9) = -x(3)+0.1;
y(10) = -x(4)+0.1;
% Variable upper bounds
y(11) = x(1)-10;
y(12) = x(2)-10;
y(13) = x(3)-10;
y(14) = x(4)-10;
% *************************************************************************
y=y'; 