import numpy as np
import matplotlib
import matplotlib.pyplot as plt

n = 20
malla1d = []
for i in range(n):
    malla1d.append(float(i) / (i + 1))

malla2d = []
for i in range(n):
    for j in range(n):
        malla2d.append([malla1d[i], malla1d[j]])

malla3d = []
for i in range(n):
    for j in range(n):
        for k in range(n):
            malla3d.append([malla1d[i], malla1d[j], malla1d[k]])

# una forma alternativa para 2d, sin usar for
x, y = np.meshgrid(malla1d, malla1d)
plt.scatter(x, y)
plt.show()

# guardar datos para despues graficar, por ejem con gnuplot
np.savetxt("malla1d", malla1d)
np.savetxt("malla2d", malla2d)
np.savetxt("malla3d", malla3d)