@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

@Kitonum Thanks its very useful. 

@tomleslie thanks for your answer. The piecewise function are not shown as you mentioned before Eq. (1) and (2). The appropriate form of piecewise function that I need presented in the file. Kindly help me to fix the problem. Thanks

pwFunc.mw

@acer Sorry its missing. L is the maximum value of abs(diff(f, y)) for all a ≤ x ≤ b, -∞ < y < ∞ for any function f(x,y) or

abs(diff(f, y))<=L for all a ≤ x ≤ b, -∞ < y < ∞

M is the maximum value of abs(diff(y, x, x, x)) for all a ≤ x ≤ b, -∞ < y < ∞ for any function of y(x) or

abs(diff(y, x, x, x))<=M for all a ≤ x ≤ b, -∞ < y < ∞

For example for f:=y-x^(2)+1; 0 ≤ x ≤ 2, -∞ < y < ∞ and y:=(x+1)^2-.5*exp(x); 0 ≤ x ≤ 2, -∞ < y < ∞

the values of L = 1 and M = 0.5exp(2)

@acer thanks for your reply. I have one more question. In the beginning of the following file I declared Digits:=10. But some of values like Y[1,1,1], Y[1,2,1], Y[1,2,2]... have 9 digits after the point decimal (.) like 1.213233616, 1.213873330, 1.214086568 respectively. I need the 10 digits after point decimal (.) which command I have to use instead of Digits?

Refine_Extrapolation.mw

@Kitonum thanks for your answer. I tried to make the above table (Table 5.17 without character y[1,1],y[2,1],y[3,1] etc.) in the following file can you fix it? Because it is more convenient for me. 

I have one more question. In the beginning of the following file I declared Digits:=10. But some of values like Y[1,1,1], Y[1,2,1], Y[1,2,2]... have 9 digits after the point (.) like 1.213233616, 1.213873330, 1.214086568 respectively. I need the 10 digits after point (.) which command I have to use instead of Digits?

Refine_Extrapolation.mw

@tomleslie , In addition, I want to recall a vector B in file Q3 at highlighted portion whereas B matrix is present in file Q4. Which is the command to recall or import a B matrix in file Q3 at highlighted portion from another file Q4.

I shall be very grateful for your support.

Q3.mwQ4.mw