@sand15 for 3D, I have not tried yet, but for the line plot, I am using this command: legend = typeset("We = ", WeVals[i], legendstyle = [location = bottom, font = ["TIMES", "italic", 18] and it appears like this 

Even, sir, if there is any issue regarding showing the legend leave, I will do some manual tricks, but I need to verify this: "The file you uploaded is working accordingly, and even on my side, there is no convergence issue in this file. Sir, I am trying to generate legends at the bottom of the generated figures. Actually, my problem has multiple parameters, and I am trying to determine the impact of these parameters on Nu (Nusselt number). As we generated 3D bar graphs, one axis was for Nu, and the other two were for "We" and Br or M and Br. Sir, let me try to explain in more detail. If we want to check variation in Nu with respect to M = [0.5, 1.5, 2.5, 3.5], then the axis of the graph should be Nu, We=[0.05,0.1,0.3,0.5], and Br=[0.1,0.3,0.5..0.9], and M will appear in the legend according to color or values, like one color for M=0.5, another different color for M=1.5, and so on. Similarly, if we want to check the impact of We= [0.05, 0.1, 0.3, 0.5], then the axes of output will be Nu, M=[0.5,1.5,2.5,3.5], and Br=[0.1,0.3,0.5...0.9], and the "We" values appear in the legends. "

@sand15 Sir, thanks a lot for your help. The file you uploaded is working accordingly, and even on my side, there is no convergence issue in this file. Sir, I am trying to generate legends at the bottom of the generated figures. Actually, my problem has multiple parameters, and I am trying to determine the impact of these parameters on Nu (Nusselt number). As we generated 3D bar graphs, one axis was for Nu, and the other two were for "We" and Br or M and Br. Sir, let me try to explain in more detail. If we want to check variation in Nu with respect to M = [0.5, 1.5, 2.5, 3.5], then the axis of the graph should be Nu, We=[0.05,0.1,0.3,0.5], and Br=[0.1,0.3,0.5..0.9], and M will appear in the legend according to color or values, like one color for M=0.5, another different color for M=1.5, and so on. Similarly, if we want to check the impact of We= [0.05, 0.1, 0.3, 0.5], then the axes of output will be Nu, M=[0.5,1.5,2.5,3.5], and Br=[0.1,0.3,0.5...0.9], and the "We" values appear in the legends. And is it possible to write a general method that automatically calculates and puts relevant values of Nu against the parameter for which we are checking the impact on each bar? Sir, I am uploading your worksheet because if I upload my try, it creates confusion. I apologize for any inconvenience caused by my explanation of my problem. Thanks for your help. 

Bar_view_help_sand15.mw

@sand15 thanks sir, I am reading your file. If there is any issue will ask.

@sand15 Sir, the problem occurs due to the large values of "We". Maybe it will be fine for We = 0.05, 0.1, 0.2, 0.3. Can you please share the worksheet? The file you uploaded is working, but I am confused about the values that appear in the Data matrix. How do we obtain these values? Are these values dummy or generated from my problem? Sir, I have been stuck on this problem for the last couple of days. Kindly help.

Here is the file you uploaded in your previou response.

Bars_(2).mw

@sand15 Sir, it's not a disrespectful attitude. We are here to seek help from experts like you. Actually, you spent a lot of time, but the problem is that I do not see my original system of ODEs in the code. Even when I run the file you uploaded, it totally changes the view of the graphs. I wonder why it is happening to me.

@sand15 Thanks, sir, for your suggestion. I really do not want to waste your time. I tried another approach, and the results are reproduced, but the view is not good. I tried a lot, but I cannot enhance orientation and fix the legends issue. If possible, take a look and help in refining graphs. For reference, see the attached picture.

Bar_view_help.mw

@sand15 All information is available in my Maple worksheet. I am not getting an idea of how to get this information. I have only the values of all parameters, which are already given in the worksheet. I am just trying to change my line graph to a bar graph view or a 3D surface plot

@sand15 How to implement it for my problem, not in general? I have uploaded my worksheet for line plot. kindly help

@sand15 sir thanks for your suggestion. I am looking for almost similar style. The picture I uploaded as a reference have retangle style of bar graph, but how to implement in my case?

@janhardo sir thanks alot your support 

@janhardo sir, can you please help me to fix this issue? I think the problem is in this part. 

#-------------------------compute nusselt num ate y=1--------------------#######
                     

      Nu_results[i, j] := eval(-diff(theta(y), y), dsol(1));
end do;

This diff(theta(y), y), dsol(1))  is not working properly. Besides this, if put 0 at place of 1, the same file works good and give different values for Nuesselt numbers.

Same_values_why_help.mw

@janhardo I changed some values and functions in the previous code that you uploaded, but after this it did not provide any values. I followed all of the steps you wrote. Can you please help me correct the file?

@janhardo  Ok sir, thanks a lot for your help. My last query: I am calculating the Nusselt number, but it gives the same values for all parameter sets. Can you please help me why the code is giving the same results for the Nusselt number?

Same_values_why_help.mw

@janhardo Sir, does the file work on your computer? Because it gives error in my side

• @janhardo Sir, thanks for your help. Is it possible to generate single values of skin friction rather at any point of y=0.1 and k=0.1, other parameters remaining the same, then tables or an Excel sheet?