private double PotentialExpectationValue_fm()
{
Complex expectationValue = 0;
for (int j = 1; j < Param.StepNumber; j++)
{
expectationValue += WaveFunctionNorm[j] * Potential_fm.Value(Radius_fm[j]);
}
expectationValue += 0.5 * WaveFunctionNorm[Param.StepNumber]
* Potential_fm.Value(Radius_fm[Param.StepNumber]);
return(StepSize_fm * ComplexMath.Abs(expectationValue));
}
private void UpdatePotential()
{
Potential_fm.UpdateAlpha(AlphaSoft);
}