Download trip_int_ex2.mw

You cannot link just the character length or length of the 1D code representation of a general expression with the width the takes to pretty-print it as 2D Output. It generally needs a bit more finesse, according to type of expression.

However, I've seen at least one quite decent attempt at estimating 2D rendering width for algebraic expressions. I'd have to go hunting to find it.

But the elision fade-out was done using Typesetting constructs, in the instance I'm recollecting. I'm not 100% sure that the further/repeat 2D type-setting by Tabulate would work ok. If not then the DocumentTools:-Layout constructors might be used to get the same kind of Table result.

Alas it is the end of the day for me.

@nmacsai These two ways of injecting the COLBAR substructure into the cylinder's PLOT3D structure end up with the same result.

Download colbar_M2023_2.mw

I used contourplot to get its accompanying color-bar. You could alternatively use densityplot and then pass that your choice of colorscheme option. Eg,

colbar_M2023_3.mw

[edit] The constructed color-bar is itself just a 2D densityplot with aligned values shown. There are doubtless more efficient ways to construct it, and I just chose one that was reasonably straightforward. The trick is to wrap it in a COLBAR call, and then make it the new last operand of the PLOT3D structure. You could likely use other kinds of 2D plots. The GUI recognize the COLBAR within the PLOT3D and renders it (only) on the right side.

@Dkunb evalf(Int(...)) uses option remember,system so doing duplicates in rapid succession will utilize "remembered" results.

So your idea to do the numeric integration once (each, followed by substitution) is a well-meaning idea but not beneficial if all's done properly.

What can make a big effect here is using method=_d01akc (periodic), or method=_d01ajc for you evalf(Int(...)) calls. Do that after substituting numeric parameter values, otherwise an inadvertant evalf will strip off that option.

You latest code is too convoluted and perhaps not robust. The combination of Expand and then multiplication by I looks nonrobust. Best would be if you provided the slow code in full, including the float results. Otherwise fixing it up is problematic.

Did you intend exp(x^2) instead of e^(z^2) ?

Download rho-analysis_acc.mw

@C_R

IVPOrder is a local of the Student:-ODEs module, not an export.

A few ways get what you're after,

restart;
kernelopts(opaquemodules=false):
stopat(Student:-ODEs:-IVPOrder1,14);

        [ODEs:-IVPOrder1]

restart;
stopat(Student:-ODEs::IVPOrder1);

        [ODEs:-IVPOrder1]

restart;
stopat((Student::ODEs)::IVPOrder1);

        [ODEs:-IVPOrder1]

I am surprised to see such a mistake in a Library routines. It's pretty common to do instead something like,

   indets(expr, And(name,Not(constant)))

or,

   indets(expr, name) minus {constants};

@MaPal93 Don't spawn more wholly sepatate new threads on this.

In your second to last example you first asked, "What proportion of the population needs to be vaccinated in order to prevent an epidemic...", and then you set p=0.5 where initial condition is R(0)=p, and then ended with the statement,
   "The required vaccination coverage is around 50% ."

Is that supposed to show some deduced computation? Clearly the "required vaccination coverage" was not computed there. We might accept that 50% is adequate (the R(t) plot is flat), but why wouldn't slightly less also work? The word "required" makes it seem as if no less would do.

The same thing happens in the last example. You start by setting p=0.83, and then end with the statement that 83% is required. Why is that much required? What's the computational criterion to accept that 0.83 was the smallest value that was adequate?

@dharr Is this at all related to procedure f in the 3rd attachment in this Answer? (tracking largest positive RootOf, to help with "switch" when, say, roots coincide at a point, or number of real roots changes, etc).

There I used fnormal and simplify(..,zero) and doubled-Digits to help with imaginary artefacts.

Plot problem reported, regression in Maple 2024.0,

plot(Pi/2-arccos(h/sqrt(h^2+2)),h=0..3),
plot(Pi/2-arccos(h/sqrt(h^2+2)),h=0..3,adaptive=true);

plot(arccos(1/sqrt(x))),
plot(arccos(1/sqrt(x)),adaptive=true);

@minhthien2016 The result is correct.

But you can get rid of the warning with,

   solve(temp = Pi/3, h, useassumptions);

What happens if you replace sum by add?

note: It'd be more useful if you would provide the complete code-to-reproduce.

@nmacsai If there were a way to change the color of the edge of a POLYGONS structure then we'd likely give it in an answer.

I believe that a POLYGONS plotting structure can have either a style where the edge is shown and black, or a style where the edge is not shown.

A POLYGONS plotting structure only has one COLOUR substructure, and that relates only to the color on the interior face.

@AHSAN 

You should,
1) define properly (clearly and in full detail) what you mean by, "check the convergence or validity of the numeric method".

2) define properly (dlearly and in full detail) what you mean by, "streamline pattern" in this context. Simply restating it as, say, "stream functions" is not adequate.

If I don't respond further it's likely because I didn't see adequate details.