I wish to  plot 2D animate for the soltion of this equation here is the code

restart;
with(PDEtools):
with(ArrayTools):
with(plots):

f:=u->sech(u):
g:=v->sech(v):
h:=1/10:
N:=20:
M:=20:
V:=x->x^2:
psi:=Array(0..N/h+1,0..M/h+1):
for i from 0 to N/h do
psi[i,0]:=evalf(f(i*h)):
od:
for j from 0 to M/h do
psi[0,j]:=evalf(g(j*h)):
od:
for i from 1 to N/h do
for j from 1 to M/h do
psi[i,j]:=-psi[i-1, j-1]+(1-(1/8)*h^2*V((1/2)*h*(j-i-1)))*psi[i, j-1]+(1-(1/8)*h^2*V((1/2)*h*(j- i+1)))*psi[i-1,j]:
od:
od:
ls:=[seq([seq([i*h,j*h,psi[i,j]],i=0..N/h)],j=0..M/h)]:

surfdata((ls),axes=boxed,labels=[`u`,`v`,`psi(u,v)`],shading=zhue,style=patchcontour);

Hello 

 I couldn't make plot3d in my Maple code please see the file I attached it and let me know any idea how can I make plot it in 3D

Thanks

I have here Maple proc, but it doesn't work. I am try to solve this equation 
pde := diff(u(t, x), t, t)-(diff(u(t, x), x, x))+x^2*u(t, x);

with 3 piont known , we have 
stencil:=u(i,j)=-u(i-1,j-1)+(1-h^2/8*v(h/2*(j-i-1))^2)*u(i,j-1)+(1-h^2/8*v(h/2*(j-i+1))^2)*u(i-1,j);

please see the file I attached here

Thank you

Hello

 I wish to solve nonlinear PDE equation numerically  in u & v coordinate . where u=t-x, v=t+x . we note x=(v-u)/2.. 

pde := diff(u(t,x),t,t) - diff(u(t,x),x,x) + x^2*u(t,x);  

 We need only initial condition let f(u)=sech(u) & g(v)=sech(v). So u((i,0)=f(ih), u(0,j)=g(jh) 

we summe step size of time j =  step size of spatial i .

the stencil we use , 

HelloI have Mathematica code(notebook) I want to convert it to Maple code (worksheet)
 Thanks