I am trying to use finite difference method to numerically solve a pair of partial differential equation by using Matlab.

So far, this is my working .


Did my working calculation on finite differences correct? For matlab coding (don't known anything) i still learn the basic thing on how to create the coding from youtube.