function area = S(x)
r = x(1);
hcyl = x(2);
hcon = x(3);
Scyl = 2*pi*r*hcyl;
Scon = pi*r*sqrt(r^2+hcon^2);
area = Scyl+Scon;
