function [c,ceq] = bincon(x)
r = x(1);
hcyl = x(2);
hcon = x(3);
angdeg = 20.4;
angrad = angdeg/180*pi;
c = atan(r/hcon)-angrad;
volc = 10;
ceq = V(x)-volc;
