#Import statements

NDENSITIES = 100000
densities = np.random.normal(-100, 50 NDENSITIES)
results = np.zeros(shape=n, NDENSITIES))
for density in range(0, NDENSITIES):
    model = Prism(-1000, 1000, -1000, 1000, 2000, 2500, {'density': densities[density]})
    data = np.array(prism.potential(x, y, z, [model]))
    results[:,density] = data