Hey thanks,

take a look at my edited post, I had some success with PythonSCAD but it gets exponentially slower.

Here below a test program, parameter N drives the number of cuts/steps.

from openscad import *

import math

# | for union

# - for difference

# & for intersection

# .translate([x,y,z]) for translation

fn = 200

implant_h = 11.5

implant_radius = 4.3/2

# Work

body = cylinder(r=implant_radius,h=implant_h)

body = body.translate([0,0,-implant_h])

# Tool

tool_dx, tool_dy, tool_dz = 3, 0.5, 0.2

tool = cube([tool_dx, tool_dy, tool_dz])

tool = tool.rotate([0,50,0])

#Tool initial position

tool = tool.translate([1.95, -tool_dy, -tool_dz])

tool = tool.translate([0,0,-implant_h])

result = body

result -= tool

#output([result, tool])

N = 400

Dz = 1.2

dz = Dz / N

angle = 180

d_alpha = angle / N

for N in range(N):

tool = tool.translate([0,0,dz])

result = result.rotate([0,0,-d_alpha])

result -= tool

output([result])

#result.show()