One wonders what Maple commands you tried.

It's a fascinating question that I'll have to think about. To make sure that the terms of the problem are absolutely clear

1. Is it a right cylinder?
2. Is it a circular cylinder?
3. Is it in 3-space? or higher dimension?
4. Is it placed arbitrarily? That is, can the axis be any line in the space from question 3?
5. How many points are given on the surface?

What is your desired output for n=26, p=2?

@Alejandro Jakubi There's plenty of room to fit all options on the top level.

@emma hassan The matrix/vector multiplication operator is a simple dot (.).

Are the values that you got for vectors U[1] to U[4] in the posted worksheet correct?

@emma hassan The matrix/vector multiplication operator is a simple dot (.).

Are the values that you got for vectors U[1] to U[4] in the posted worksheet correct?

@syn_burz I got no output and no message about the kernel, although it was obvious that the kernel had crashed. I'm using Maple 17, 64-bit, Windows 8.

@syn_burz I got no output and no message about the kernel, although it was obvious that the kernel had crashed. I'm using Maple 17, 64-bit, Windows 8.

What should the output be if the number is 26 and the prime is 3?

@cskoog Thanks. I posted my Reply primarily to find out what you just said.

Check the results of command kernelopts(version).

@Markiyan Hirnyk It would be useful, as a first step to the 3D goal, to have a canonical ordering of the vertices, edges, or faces of polyhedra.

Please re-run the command without the 'evaluate' option and post the results. If the code contains large matrices, please reduce them to 2x2 examples.

 Your initial conditions contain

ur(e,0)=1/ch(e)

What is ch? Do you mean cosh? 

If I change it to cosh, I still get an error, but it's a different error with a more-informative message. But I am not pursuing that error any further unless you say that cosh(e) is what you meant.

@emma hassan Sorry, I don't know enough about finite difference solutions to PDEs to help you anymore with this.