@mmcdara Yes, sir you are right; the expression PP is dependent only on x. We have to calculate PP_Max from the expression PP by using the second derivative test, and the procedure is defined in the Maple sheet. Then we need to plot PP_Max against k for different values of epsilon. Further, we need to plot PP/PP_Max against X for different values of epsilon. For reference, I uploaded some pictures for both. Our  results should look like that, not exactly the same.

@mmcdara  I checked, and there is no issue in the plot now. Can you please upload the Maple sheet and how to plot PP\PP_max?

Actually, I am trying to find the maximum relation by using the derivative test. My actual expression is a bit complex and difficult to solve, so I use series from means, but this has again a problem. So could you please have a look at my maple sheet and correct the procedure? I take help only but have problems and do not produce results. I explained the procedure in the MPAL sheet, but you do not know whether it is correct or not.

Maximum_help.mw

The results look like this

@Kitonum Thanks sir

@Kitonum Sir, I have finished the plotting for the ratio of two expressions but am unable to extract expressions for k greater and equal to zero.

Extract_required_solution.mw

@Kitonum sir, Actually, According my problem requiermnets the integration is required of expression "Pre" for x=0..1 and need to find a symbolic expression in the of form of k. 
Besides this, need to plot Expression of Pre_max versus k by taking Pre_max along vertical axis and k along horizontal axis.

Also, need to plot the ratio of Pre / Pre_max versus x for different vakue of k such as 2.14,3,20,4,30.

@Kitonum , Thanks, sir, for your response. Actually, my question has two parts: first, we need to find the maximum expression at the point (critical), where the function is maximum.

and

In the second part, we need to find a symbolic expression. In this part, we obtained 4 answers... [2] When I integrate the given expression, it gives 4 cases. How do we separate the required expression for the k condition?..

I have uploaded the file with more explanation; can you please have a look at it?

Maximum_expression_New_1.mw

@mmcdara Thanks a lot sir

@mmcdara Thanks; when I fix the range for j=1...2 I got the error about Offset=-1/4 being unexpected, so I removed this command, code was run and produced the results, but some curves vanished, and even the values were not aligned on the bars

Error_Matrix_index_out_of_range.mw

@mmcdara thanks sir for time and help. By the way how to fix the error  "Matrix index out of range
"Error_Matrix_index_out_of_range.mw

@mmcdara thanks alot for response and affort. In my work I do not want ignore the negative values as you define new data by using adata := abs~(data)......my opinion is to change only display of colum graphs inside the graphs is look like the view you produced either it is on positive or negative sides. 

@mmcdara Thanks, sir, Is it possible to change the bar chats, Means look like rectangles for the example display in the below picture? Besdie, there is a box with a white line in graphs, it possible to remove them? Also, I tried by adding more option but have error, please correct sir.

another_try.mw

look like the graphs given below

@Carl Love pp is notation and stand for pressure profile.  Below is the picture that I take from a paper of this.

In Mathematica file (pdf) written a simplified form after putting u=diff(psi,y) and its simplified form is gievn in equation 62 and its BCs are also given Mathematica file.

firstly Need to solve for  eq 62 for psi against y and then use the solution of psi to get the solution for u against and similarly use the solution of psi to plot diff(p,x) against x and integrate the solution of diff(p,x) and find solution for p against x .  In mathematica file pp stand for solution of p. I am not understanding mathematics code and how to transform that in Maple?

@acer , sir, thanks for the correction. Is it possible to create one more column for the answer check? Means residual