@awass 

Using dsolve/numeric with parameters is not discouraged. Indeed, there is a dedicated and documented mechanism for doing it.

What is now deprecated is calling dsolve/numeric for IVPs that have unassigned global names in them (ie, parameters), and then invoking dsolve's returned solution procedures only after assigning numeric values to those names.

See my Reply below.

@awass Did you see that,

   F(7.0), F(7.1), G(7.0), G(7.1)

return 1D Arrays, whose entries are the numerically solved values y(x) corresponding to the x-values in 1D Array A. (But, for the different parameter values a=7.0 and a=7.1.)

You can always wrap those 1D Arrays in a call to Vector, if you'd much prefer that.

The whole point of the code in my Answer (just as in dharr's) is to compute (for an arbitrary parameter value) a 1D Array of y(x) values that correspond to a forced 1D Array of x-values. We both uses an even spread of x-values.

Otherwise, I don't see what you're requesting.

@sand15 The OP stated that TM1,TM2,TM3 are all positive.

But when Pn=0.2, say, they are all negative for all w=0..1.

Shouldn't the conditions that they are all positive further restrict the regions in your plot?

(...using an extra quadratic, etc, if you are doing it algebraically.).

@AHSAN I cannot say anything about your formula, where you have,

     # Evaluate skin friction at a specific point

etc. How would anyone else know what you mean by your formula, or whether your formula is correct?

You can adjust the 3D plot by zooming in and panning it, in the inlined plotting window. You can do that in the usual way, ie. right-click on the 3D plot, and in the popup menu go down to the "Manipulator" item and toggle it to "Zoom In", or "Pan", end then use the mouse to zoom/pan. Only the actual plot will zoom&pan, so its size&position relative to the color-bar and title can be thusly controlled.

The GUI has its own notion of how large the plot should be, at first, and it depends on the size of the title as well as the presence of the colorbar, and it also depends on the GUI's decision to allow you to rotate the plot without the axes going out-of-view. You can adjust it to your taste, by zooming in/out and panning.

note: zoom effects don't always get retained when the sheet is closed/reopened.
3D_shadow_at_the_base_Help_alt_ac.mw

See the plot3d option axes and axis for choices of the axes-style and location (high, low, etc). Those gray gridlines aren't available with a dedicated option for 3D plots, but you can mimic them easily by adding lines at the edges of your plot's view. You can do that. See also this old posting (or others like it).

You should carefully read the documentation pages relating to the commands and functionality about which you ask.

@Jean-Michel Yes, the subscripts are displayed slightly differently (upright Ronan, versus italic) to provide a visual cue.

It's meant to help perceive what was utilized. Sometimes people use both together, accidentally, while thinking that they're all the same one.

@Jean-Michel After you assign the same expression to both names then your check return true because it's comparing the same assigned values rather than the names to which they're assigned.

@AHSAN Did you try my followup, with a code adjustment for the color-bar's dimensions? I used Maple 2024.2 there, but I'm hoping the effect is the same in your Maple 2025.

Or is your PNG immediately above from my original Answer's code?

@sand15 In my example above I showed colored vertical axes alongside black axis-labels, using a typesetting construction (of which you're already familiar). A similar effect could also be done for the tickmarks, as a different color from the axis-bar.

I don't know of a way to get a linestyle or transparence for the axis-bars, except by faking it with a dedicated line (and supressing the axis proper).

@AHSAN Of course, you could simply get rid of the color-bar altogether, by specifiying colorbar=false instead, in the call to contourplot3d.

But I believe that you could also adjust the original by shortening the color-bar in the substructure, after construction. Eg,

Download 3D_shadow_at_the_base_Help_acc.mw

@sand15 That too is nice and visually easy to understand.

[edit: it's more complicated to do, though, and requires adjustment for other data.]

@Carl Love I don't know why it was done that way. Perhaps the author really preferred consistency, or half-expected that some other module locals might be needed later.

There are other appliable module exports of Student:-Calculus1 which are documented/user-level commands of the package. Some have other locals, and some do not. Eg,

Download SC1_am.mw

The OP has also asked me how I knew that it was such an appliable module, ie. a module with member ModuleApply, which (like a procedure) could be called with arguments.

One answer is simply that having such a module member ModuleApply is the mechanism that allows for that behaviour of a module's being callable like a procedure (ie. which could be applied to arguments).

But I also already knew it because I've worked with that command before.

ps. The printing mechanisms for modules have had a few quirks, over time, as the following shows:

Download module_hmm.mw

@aroche Thanks, Austin!

I guess I was so hung up on avoiding the OP's earlier use of x>=0 that I was willing to accept R::real, unjustified, like a mirage. Thanks for the correction.

@dharr It seems to me that the OP's sentence,

  "Then,  `-\\infty` to `- + \\infty`."

pertains to the results that attain -- undesirably -- when he executes his original attempt.

And then he mentioned what he wanted,

   "How can I replace \infty with +\infty but keep -\infty unchanged?"

@dharr The temporary shunting aside is indeed a classic trick.

But perhaps you intended just,

    s2 := SubstituteAll(%, "#@#","-\\infty");