Need help in MATLAB programming for coding partial differential equation

I need help for coding coupled partial differential equation in MATLAB. There are 4 main equation and 4 nos of unknown variables. For the other known variables I have already prepared code and file is attached. These equations are varying in single space and time variable. Also there is initial and boundary conditions given for apply.