Title: Numerical calculation method and program implementation