That seems like a very good suggestion. I imposed the BC D(Phi)(0,t) = 0 in the first place to tackle this singularity at the equation. but it seems that when I make these changes it goes away. I have not tried it yet on maple..
However I have a question. Why does maple ask for boundary conditions even for wellposed initial value problems that only needs initial conditions?
Also if you please can look at the original problem I had, of which this is a branch post, I had to make the spacestep really small for it to match the same function at t=0 which I put in as initial condition. When I try to plot the solution at t=0 for large spacestep it differes significantly from the known function at t=0,, I find it bizzare why maple would change it in the process of solving it
Thanks very much for your attention. Looking forward to your suggestions