Minecraft Planer

mc-planer.py

import math
import itertools
from collections import namedtuple

point = namedtuple("point", "x y z")

a,b,c = 10,10,5

dx,dy,dz = 10,10,5

points = set()

for x,y in itertools.product(range(0,dx+1), range(0, dy+1)):
    try:
        z = c*math.sqrt(1-x*x/a/a-y*y/b/b)
        points.add(point(x,y,int(z)))
    except ValueError:
        pass

for z in range(0, dz+1):
    print(f"layer z={z}")
    s = '\n'.join(''.join(('#' if point(x,y,z) in points else '.') for x in range(0, dx+1)) for y in range(0,dy+1))
    print(s)