function S = povrch_kelimku(r,R,V)
     if r<0||R<0||V<0
         error('Kuzel nemuze mit zaporny polomer zadne z podstav, ani objem')
     end
  if V==0||(r==0&&R==0)
       S=0;
  else
        h=3*V.*((pi.*(R^2+R*r+r^2))^-1);
        S=pi.*(r^2)+pi.*(r+R)*sqrt(h^2 + (R-r)^2);
  end
end