Dear @C_R , please find the following file

help.mw

@C_R 

Shooting+RK4+Secant_Method.mw

Sending you again @C_R 

Shooting+RK4+Secant_Method.mw

Thank you @C_R for your response. I am sending you my complete code for better understanding.

Shooting+RK4+Secant_Method.mw

Please see if the code not stoped when both conditions fulfilled.

@sand15 many thanks for your answer. I will avoid such phrases next time.

Thank you @sand15! For your reply. I don't want to plot log[10] data. I want to plot surface plot of all original values of MAE which are less than 1.

In the attached file

Help.mw

I constructed sequences like RK11, RK22, RK33, and RK44 to plot surface plot, but don't get the appropriate one, for instant, the obtained surface plot for RK11 is

 Here some lines are directed toward zero which is not needed. The same problem happens with RK22, RK33 and RK44 (see Help.mw). Moreover, I also plot the density plot for the original data.

I am waiting for your positive response. Please take care.

Dear Sir @acer, I hope you are doing well. I have plotted the contours of absolute errors of four approximate solutions. The contours of the first three solutions are reasonable since the same behavior has been obtained in 3d plots. But the fourth contour plot (all blue) is inappropriate as the 3d behavior is different. Therefore, kindly help me to fix the problem

help_cont.mw

I shall be waiting for your positive response. Please take care

@acer Dear Professor, I benefited a lot from the file shared above. But, now I am trying to plot coutourplot of the following continuous function

sin(x)*sin(y)-x*y+(1/6)*x*y^3-(1/120)*x*y^5+(1/5040)*x*y^7-(1/362880)*x*y^9+(1/6)*x^3*y-(1/36)*x^3*y^3+(1/720)*x^3*y^5-(1/30240)*x^3*y^7+(1/2177280)*x^3*y^9-(1/120)*x^5*y+(1/720)*x^5*y^3-(1/14400)*x^5*y^5+(1/604800)*x^5*y^7-(1/43545600)*x^5*y^9+(1/5040)*x^7*y-(1/30240)*x^7*y^3+(1/604800)*x^7*y^5-(1/25401600)*x^7*y^7+(1/1828915200)*x^7*y^9-(1/362880)*x^9*y+(1/2177280)*x^9*y^3-(1/43545600)*x^9*y^5+(1/1828915200)*x^9*y^7-(1/131681894400)*x^9*y^9

and want to include the same options

- the colored density-plot
 - the black contours
 - those same contours, but colored to match the density-plot
 - the colored legend items

Kindly guide me on how I can use these options for contourplot of a continuous function. I am waiting for your response. Please take care

@acer thank you so much, dear Professor!

@acer thank you, Sir! What do I have to change in the above code If I don't want to show the contour lines only, the legend should be there. 

@sand15 Many thanks, the axis was changed as per your suggestions. About "structure Tecplot can manage", I really have no idea about it. 

@sand15 

I applied finite difference method with mesh size `&Delta;x` := L/Mx; `&Delta;t` := T/N; which divides the domain 0<=t<=2 into N+1 number of points and the domain 0<=x<=1 into Mx+1 number of points.

The more accurate solution obtain as taking the higher values for N and Mx thats whats why I took N=Mx=20.

Next, the professional software means in term of plotting like I want to plot it in tecplot as tecplot results have highe resolution than maple plot.

@acer , many thanks!

@acer , Thank you I resolve it, Now just want to ask how I can show the x-axis value and t-axis value as mentioned bellow:

@acer , Now I want to plot the contour plot of the sequence for in x and t as I tried above. 

xData := [0., 0.04995839572*sin(x), 0.09966865249*sin(x), 0.1488899476*sin(x), 0.1973955598*sin(x), 0.2449786631*sin(x), 0.2914567945*sin(x), 0.3366748194*sin(x), 0.3805063771*sin(x), 0.4228539261*sin(x), 0.4636476090*sin(x), 0.5028432109*sin(x), 0.5404195003*sin(x), 0.5763752206*sin(x), 0.6107259644*sin(x), 0.6435011088*sin(x), 0.6747409422*sin(x), 0.7044940642*sin(x), 0.7328151018*sin(x), 0.7597627549*sin(x), 0.7853981634*sin(x)];

where the first output is against when t=0, the second when t=0.05, the third when t=0.1 and so on the last is for t=1.