private static double polarisationFactor(double electricField/*in kV/cm*/)
{
double polFactor = 0.00001386974812686904
+ 0.09020507207607358 * electricField
+ 0.0008134261949342792 * Math.Pow(electricField, 2)
- 0.001365930559363722 * Math.Pow(electricField, 3)
+ 0.00016467328012807663 * Math.Pow(electricField, 4)
- 9.53482 * Math.Pow(10, -6) * Math.Pow(electricField, 5)
+ 2.804041112473679 * Math.Pow(10, -7) * Math.Pow(electricField, 6)
- 3.352604355536456 * Math.Pow(10, -9) * Math.Pow(electricField, 7);
return polFactor;
}