I used the script "DielectricConstantNonEq.pl" to calculate the dielectric constant for three times.The structure and parameters are the same, but the results were different for three times. Is there anything wrong with it?

​​​​​​​