Method of lines PDE with Matlab

[Noni]

Hi,

I want to solve a system of PDE but I don't know it's correct because when I put tvec=(0,1500,31) for example in my code , it doesn't work. If someone has an idea, feel free to answer

The file with the pde ( + my reasonning ) and my code is down below.

Thanks