optM                    optM function
plot_optM               plot_optM function
