@Carl Love I recalled some of the OP's earlier postings using the same basic code layout. Some of his "parameters" are the global names to which he's assigned values, rather than some procedure's parameters having values passed into a call.

ps. It sure does help when people provide their actual code in an attachment to the Question.

@RezaZanjirani It seems to work, respecting the rhs&lhs of the inequality, as well as preserving the symbol c unevaluated, whether one uses,
   Typesetting:-Typeset(p__s < c)
or,
   Typesetting:-Typeset(c > p__s)

Scenario_ac.mw

Scenario_ac2.mw

So I'm not sure I understand what you mean when you state, "But it does not work on my own code".

ps. Your code has the input as c>p__s , by the way, though you indicate elsewhere that you want it as p__s<c.

@sand15 A fundamental danger of your trick is that some denominator involved in the computation for the generic x[i,j] could be mathematically zero when replaced by the B[i,j] expressions. Subsequent evaluation of the final result can then fall prey to the division-by-zero time-bomb.

@FDS To the right of an Answer's title is a thumbs-up icon and (if you posted the Question) a cup icon.

Anyone (with reputation score of at least 10) can vote-up an Answer by clicking the thumbs-up icon.

The person who posts the Question can also accept it as best by clicking the cup icon.

Naturally, you are free to do this for any of the Questions you post in this forum.

@Rouben Rostamian  The worksheet attached by simplevn1967 was saved by him using Maple 17 (released 2013).

@lemelinm The effect of your first suggestion -- to control the color of the surface "grid" lines differently from the color/shading of the surface -- is already mostly possible to achieve, and without having to recompute the surface's data points (eg. z-values).

The key here is to compute (once) the plot values, either as pure surface or pure wireframe style. Then the other can be constructed by overriding that option, ie. without recomputing the values.

For example, (and you could also adjust other qualities of these wire-frame lines, eg. thickness=2 , etc),

with(plots):
P := plot3d(-x^3+y^2,x=-1..1,y=-1..1,style=surface):
display(P,display(P,overrideoption,
                  color="Orange",style=wireframe));

Naturally, that example applied the override only to a single surface. One might well not want to have it as a blanket effect on all compound parts of a VolumeOfRevolution result. (But even there, it can be done selectively, with a bit of care.)

I'd agree that a simple choice of options for the effect you've described would be more user-friendly. On the other hand, there are several other things I'd rather take priority. Eg, separate x/z, y/z, x/y aspect ratios for the axes.

ps. I think that jtreiman's original suggestions (at top) produce a result that is easier to interpret than the default for the command in question.

@Christopher2222 Calls to time() ought to be calls to time[real]().

You could do that something like this:

caring_phase_pemanenan_predator_dgn_parametersesuai_jurnal_ac3.mw

I indexed table T by t-values, only because I don't know how you plan to access the maxima later. You could also index by the ordinal, `found`, etc.

It's not clear from your latest question whether you want the global maximum (over t=0..100, say) or all the local maxima each time around (when diff(x(t),t) is zero).

Here's a simple way to get that global max from the phaseportrait data itself, and the corresponding time, for one of the initial values.
caring_phase_pemanenan_predator_dgn_parametersesuai_jurnal_ac2.mw

There are several more complicated ways (eg, wrapping x(t) in a proc that does extra tracking), or using dsolve/numeric's events options, etc. That kind of approach can also come in useful if you want to save values "each time around", or each time diff(x(t),t)=0, etc. Let us know if you need.

Are x,y,z all to be considered as real?

Why do you use the symbolic option?

Could you provide some examples of input and desired answer?

@AHSAN If you want to move the horizontal line segment,

   plot([[-5, 0], [5, 0]], color=black)

either up or down then simply change the second values in the lists to something else. Do you not understand that the second values in the lists are the vertical (y) coordinates of the line's end-points?

Similarly, if you want to move the vertical line segment,

   plot([[0, ymin-eps], [0, ymax+eps]], color=black)

either left or right then simply change the first values in the lists. Do you not understand that the first values in the lists are the horizontal (x) coordinates of the line's end-points?

I think you could have easily figured this out for yourself.

@AHSAN The rangesonly option makes getdata return the ranges for all the dimensions of the plot.

A 2D plot has two dimensions. So in your case it's returning two ranges. The indexing, [2], accesses the second of those ranges, which is the vertical range.

I think that you could have figured that out.

Sorry, I don't understand what you mean or want, by, "show these line according to desire position".

@AHSAN That operator/procedure assigned to r is used to construct the columns, as rectangles.

The second parameter of that operator, y, is used to specify the vertical end (ie. height, whether positive or negative) of the column.

The other end of the column, vertically, is at 0, ie. the x-axis.

ps. I think that you could have figured that out.

@AHSAN I have modified the procedure so that it is more flexible, accepting additional plotting options and the numerical format for the column labels.

Help_Bar_Graph_acc.mw

And now it be called to produce these, say,

@Nicole Sharp Regarding your statement about extra parentheses being required in Maple, you could also enter that as,

   evalf(2*Pi^5*1.380649e-23^4/(15*299792458^2*6.62607015e-34^3))