half-hexagon in OpenSCAD

height = 5;
radius = 20;

difference() {
  cylinder(height, radius, radius, $fn=6);
  translate([-radius, 0, -height/4]) {            
    cube([radius*2, radius, radius]);
  }
}

It works for these parameters but it needs more work to be really parametric. But hey, it’s my first OpenSCAD thing. I actually exported as .STL and printed it.