further work on numerical method