Hello again, To be more specific, here is the problem I have been willing to solve with a high accuracy:

tmax=25;
xmax=6*Sqrt[tmax];
sol=NDSolve[{D[c[x,t],t]==D[c[x,t],{x,2}]-c[x,t]^2,c[x,0]==0,Derivative[1,0][c][0,t]==-1,c[xmax,t]==0},c,{x,0,xmax},{t,0,tmax}]
Plot3D[Evaluate[c[x,t]/.sol],{x,0,xmax},{t,0,tmax}]

Unfortunately, this code fails with a message that initial/boundary conditions are inconsistent, and the plot shows just a zero solution, implying that the relative error is 100%. From searching internet I see that I am not the first person who encounters this difficulty, but the workarounds suggested here and there are rather incomprehensible to me. Can anybody suggest how to effectively handle this problem? I have to comment that the apparent limitation of MATHEMATICA to PDEs accompanied by consistent initial/boundary conditions is highly over-restrictive, in my opinion. In the area I have been working for nearly 40 years now, actually all PDE-based models are of this sort, but people obtain useful solutions using finite difference or finite element methods. So, it would be beneficial if MATHEMATICA did not impose such restrictions, or at least allowed one to optionally remove them. Leslaw