Plot the calculated density
Read the potential from file
Use spline interpolation for nen