Error while using optimizaton algorithm on PDE function with four parameters

Question

I wrote a program to solve a concentration based PDE in 1D in python (see: opti(eta_m,eta_e,M_e,M_c). The function itself solves a system of equations with a se