finish ComputeInitialPotentials.m