Discretizzazione eq differenziale alle derivate parziali matlab

Ciao a tutti,
come procedereste per discretizzare e implementare in matlab nel modo più fattibile un'equazione differenziale con tale aspetto $(del)/(delx)(a(x,y)(delF)/(delx)+b(x,y)(delF)/(dely))+(del)/(dely)(c(x,y)(delF)/(dely)+d(x,y)(delF)/(delx))=(dele(x,y))/(delx)$, su un dominio quadrato $[0,1]xx[0,1]$ e con $F$ nulla agli estremi come condizione?
Cerco più che altro idee su come poter fare nel modo migliore, nulla di troppo particolareggiato
Grazie a chiunque voglia aiutarmi

Per un problema di questo tipo, se i coefficienti non sono particolarmente irregolari, io userei il metodo delle differenze finite senza pensarci due volte.