@Carl Love dear Sir, I have seen anything wrong with the epsilon value as I mentioned above.

@mmcdara thank you, but you used different initial conditions. Can you examine the solution with my ICs?

Num_PDE.mw

@Carl Love thank you! I need the results for spacestep=1/200 and timestep=1/200. But can't find it. Can see my file and resolve my issue to get the result. Thanks

Num_PDE.mw

@acer thank you for your response.

@Rouben Rostamian thank you for your response. Actually, I used MatrixEquation1 for higher M (for M>50) it takes a lot of time on the other hand MatrixEquation2 produces the solution in less time. That's why I want to use MatrixEquation2 (to save time). But the problem is it gives a lesser number of solutions. Kindly guide me on how I can obtain the same number of solutions by using MatrixEquation2 as I obtained in MatrixEquation1. I am waiting for your positive response. Please take care.

@Rouben Rostamian thank you for your response. Now, I want to compute the eigenvalues of matrix B (highlighted in red color in the attached file) taking its entries as integers (mentioned as MatrixEquation1) and taking its entries as decimals (mentioned as MatrixEquation2). The problem is that in both ways I obtained a different number of eigenvalues (obtained lesser values in the case of using MatrixEquation2). I want the same values in both ways since I am using the same matrix (B). Kindly guide me in this matter. I am waiting for your positive response. Please take care and thanks

EIGEN_VALUES.mw

@tomleslie that's great and many thanks. Kindly let me know where we adjust this sequence in the above shared code.

@tomleslie in addition, the actual domain of the problem is 0<=eta<=1 and 0<=X<=1, therefore I need the output of contour plot in this domain instead of [0,4]. 

@tomleslie thanks for your response. I need the dat file because I want to use that files to plot the contour in tecplot for better resolution. 

@tomleslie thank you for your response. I have extend your idea but not attained satisfactory results. Please see the highlighted (red) in the attached file

Buitin.mw

@acer thanks for your response

@Carl Love I want to convert the following relation

[Vector(5, {(1) = 1, (2) = 2, (3) = 3, (4) = 4, (5) = 5}), Vector(5, {(1) = 1, (2) = 1/2, (3) = 1/3, (4) = 1/4, (5) = 1/5}), Vector(4, {(1) = -1/2, (2) = -1/6, (3) = -1/12, (4) = -1/20}), Vector(3, {(1) = 1/3, (2) = 1/12, (3) = 1/30}), Vector(2, {(1) = -1/4, (2) = -1/20}), Vector(1, {(1) = 1/5})]

into the matrix form as,

[[[1,1, , , , ],[ , ,-1/(2), , , ],[2,1/(2), ,1/(3), , ],[ , ,-1/(6), ,-1/(4), ],[3,1/(3), ,1/(12),,1/(5)],[ , ,-1/(12), ,-1/(20), ],[4,1/(4), ,1/(30), , ],[ , ,-1/(20), , , ],[5,1/(5), , , , ]]]

Then how can I convert it into the following table using printf/other command

Thanks in advance.

@Carl Love Thanks. Let me check it!

@tomleslie thank you for your answer. Yes, I need the results as you obtained in the vector form. But, it is more efficient if you use the following logic to make the table

printf("_______________________________________________________________________________\\n");  
for i from 0 by 1 while i<=2 N do:  
if irem(i,2)=0 then
printf("%2.1f%16.7f%16.7f\\n",i,x[i/(2)],y[i/(2)]);
else  
printf("%2.1f%16.7f%16.7f\\n",i,,); fi;od; printf("_______________________________________________________________________________\\n")
I tried to make the table but find some errors. Kindly have a look. Thanks

@Preben Alsholm thanks for highlighting the mistake. I am sending you modify version. But facing same problem. Please have a look and give me some positive response.

Compare_coeff.mw