for just the first equation, I have successfully implemented the Crank-Nicolson method (euler method is not accurate enough) but I run into difficulties when I try to do the second one - I don't know how to deal with having both Cr and Cp... is this even possible?
Any help would be greatly appreciated.