function henger()
    [r m] = meshgrid(0:0.05:4, 0:0.05:4);
    v = r .* r * pi .* m;
    a = 2 * r * pi .* m + 2 * r .* r * pi;
    v = v .* (a > 70);
    surf(r, m, v);
end