Question: numerical solution of PDE with MOL - help with assignment of RHS of ODEs in dsolve procedure

Hello, I am trying to solve this PDE using the Method of Lines and so I turn it into N ODEs. The following is how I plan on solving it; p1:=proc(N,x,Y,YP) FD:=proc(i,N,Y) if (i=1) then "use this stencil (like (Y[1] - Y[2]) / 60*dx); elif (i=2) then "another stencil"; ... ... else "last, general stencil"; end if; end proc; SD:=proc(i,N,Y)
Please Wait...