ψ =  ψa - -F(ψa)-
           F ′(ψa )

until |ψ - ψa| < ϵ (accuracy desired)