ported some cylindrical mathematica code to julia