public double Get_Equilibrium_Chemical_Potential()
{
OneVariableFunction charge_func = new OneVariableFunction(new Func <double, double>(Get_ChargeDensity));
// initialise root finder and search for a root for the chemical potential
RiddersRootFinder finder = new RiddersRootFinder();
double result = finder.Find(charge_func, -1.0 * band_gap, band_gap);
return(result);
}
public double Get_Equilibrium_Chemical_Potential()
{
OneVariableFunction charge_func = new OneVariableFunction(new Func<double, double>(Get_ChargeDensity));
// initialise root finder and search for a root for the chemical potential
RiddersRootFinder finder = new RiddersRootFinder();
double result = finder.Find(charge_func, -1.0 * band_gap, band_gap);
return result;
}
