Hi everyone! I have to solve THIS system of coupled nonlinear PDEs with THESE boundary conditions in one dimensional case (domain: [-L,L]). This is a boundary value problem (BVP). There, t is the time-variable whereas u is the space-variable. The functions gamma, beta and tau are nonlinear in lambda (which is function of u) through power law (prop to lambda^-2). How can I solve it using Matlab?