Program to square the values of a Gaussian Cubefile

The cubefile must be formatted and it must contain only one dataset (one orbital, etc).

The dimension of the cube is limited to 100