No, it would be very unwise.

acer

@dipamilo One alternative is to keep `x` local to procedure `posi4` and return unapply(objf1,x) as the last line of `posi4`.

posi1_C_modif2.mw

@dipamilo One alternative is to keep `x` local to procedure `posi4` and return unapply(objf1,x) as the last line of `posi4`.

posi1_C_modif2.mw

@kolmogorov You might have a look at this page. The error message may not match yours, but it likely covers the same general problem, namely, "when a procedure that requires numeric arguments is called with symbolic arguments."

When you call plot(p(a,b), a=.., b=..) then Maple will try and evaluate p(a,b) for nonnumeric names `a` and `b`. Your procedure `p` may not have been written to deal with symbolic names for those arguments.

@kolmogorov You might have a look at this page. The error message may not match yours, but it likely covers the same general problem, namely, "when a procedure that requires numeric arguments is called with symbolic arguments."

When you call plot(p(a,b), a=.., b=..) then Maple will try and evaluate p(a,b) for nonnumeric names `a` and `b`. Your procedure `p` may not have been written to deal with symbolic names for those arguments.

@Drahque Yes, each Embedded Plot Component has to be inserted manually, using the Components palette (left sidebar, in the Standard GUI). And then one has to right-click on it to check its name, and then use that in the DocumentTools:-Do call. A bit too much effort.

That's why I mentioned that it could be possible to write a procedure that, when passed a plot command, would insert the fully formed, final Plot Component (already populated with the plot reult). It could act just like regular `plot`, but stick the plot into the component and automagically embed it. But it's time consuming to write a procedure such as that, and it just saves a few mouse-clicks of an only partially convincing workaround to a printing bug that ought to be fixed anyway.

@Drahque Sorry if I wasn't clear enough, but I was thinking about (possible) workarounds for embedded images (as in Erik's posted example) rather than for plots.

I don't even know if using a Label would always help for embedded images -- it seemed to work for the sheet I uploaded, but that by itself is just one example.

I wasn't trying to convey that using a Label might help for plots. One might try PlotComponents rather than Labels and programatically place plots onto such. That would likely give better visual quality than with Labels and plot-export to images. It's still more awkward than just issuing a plain `plot` call (but at least it doesn't involve having to do file i/o or export to an image). They can be resized, but not just through an easy mouse drag. And there is a "border" to such components, as they appear recessed. But it may get avoid the trailing whitespace, under print/preview at least.

imageprint3.mw

These workarounds are not nearly as convenient as just calling `plot`. It might be possible to construct a nice procedure that got used as conveniently as `plot` but which (behind the scenes) did the above jobs and simply inserted the final, populated Plot component instead of a regular plot. Nicer would be if printing of regular plots were more generally fixed.

Would it be some kind of workaround, to place the image on an Label embedded component? The following seems to Print preview without the trailing whitespace in Maple 16.02 on my Windows 7,

imageprint.mw

That process takes more steps with the mouse-pointer, and I'm not sure how acceptable that would be for you Erik (or for your students). Basically, instead of main menubar->Insert->Image->FileManager it would instead be Component-palette->Label->Properties->Choose-image->FileManager.

On the bright side, I've had some progess in the whole image base64/dotm thing, and I might be able to construct a procedure that can seamlessly embed an image-overlaid Label with just a simple function call (or maybe even palette entry). A bit like this, but with fewer manual steps. I keep trying to find the spare time...

acer

It's tricky but possible to get at the Theta/Phi/Psi in the GUI's internally stored state in a programmatic way, following a rotation of a 3D plot with the mouse-cursor. But I'm not sure that it can be done continuously (repeatedly, for the same set of movements). And it sounds as if you want it to update repeatedly, in any single set of mouse movements over the plot. Is that right? Is that a strict requirement?

If you'd be ok with having three Slider components for each quantity Theta/Phi/Psi, and using only those to do rotations, then sure you can have the current values displayed and updated continuously while any one slider is adjusted. The rotations could be accomplished, in this scenario, by doing on-the-fly updates of the 3D plot in an embedded PlotComponent. This may not be satisfactory for you.

acer

@Carl Love I don't quite see how this guards against the situation that all three roots are close. If they are, then it might generate the 2nd and 3rd "highest-in-t" roots.

I don't see how to compute safely which narrow `t` range to use, either, or how high Digtit has to be at the nth step. Unless there is some more knowledge, such as (making this up now) that the close root pairs get closer from iteration to iteration, the Digits from the nth is always enough to resolve the (n+1)th close pair, the 3rd root is not as close, etc, then I don't see much that is as safe as trying to find (potentially) all three roots.

I haven't studied the modeling though. There may be justifications, based on the math.

@Carl Love I don't quite see how this guards against the situation that all three roots are close. If they are, then it might generate the 2nd and 3rd "highest-in-t" roots.

I don't see how to compute safely which narrow `t` range to use, either, or how high Digtit has to be at the nth step. Unless there is some more knowledge, such as (making this up now) that the close root pairs get closer from iteration to iteration, the Digits from the nth is always enough to resolve the (n+1)th close pair, the 3rd root is not as close, etc, then I don't see much that is as safe as trying to find (potentially) all three roots.

I haven't studied the modeling though. There may be justifications, based on the math.

ps. If the suggested process makes sense, then you could tweak it a little so that it stopped looking when 3 roots were found (at a given iteration). That would save a bit of time. But it would still involve the costly fail-out for any  iteration at which there were only eactly 1 or 2 roots. This is all under the assumption that, as you wrote, there are 3 roots at most for each iteration.

acer

ps. If the suggested process makes sense, then you could tweak it a little so that it stopped looking when 3 roots were found (at a given iteration). That would save a bit of time. But it would still involve the costly fail-out for any  iteration at which there were only eactly 1 or 2 roots. This is all under the assumption that, as you wrote, there are 3 roots at most for each iteration.

acer

Can you use the `avoid` option of fsolve?

acer

As Robert Israel (and at some earlier date, Alec Mihailovs?) has shown, it can be done without recourse to integration.

Is there an easy way to change the color of the shaded region, when using VolumeOfRevolution? I'm not sure.

This could be put into a procedure, which might be able to accept general ?plot/opions more readily than does VolumeOfRevolution. Here's a crude shot at that:

restart:

RegionBetweenCurves:=proc(f1,f2,rng::name=range(numeric))
    local v1, v2;
    v1:=lhs(rng);
    plottools:-transform(unapply([v1,v2+f2],v1,v2))
       (plot(f1-f2,rng,'filled'=true,_rest));
end proc:

RegionBetweenCurves(x^2-1, -x-1, x=-1.5 .. 1.5, color=gold);
RegionBetweenCurves(x^2-1, -x-1, x=-1.5 .. 1.5, color=cyan);
RegionBetweenCurves(sqrt(x-2), x-4, x=4..6);

I do notice a change in behavior between Maple 15 and 16 for this. In Maple 16 one of the boundary curves may be displayed more darkly than is the filled region. This seems to depend on the choice of color. In Maple 15 the shading appears to be uniform for both region and boundary. Compare with color=gold, for example.

It might be nice to get the region shaded exactly as is the boundary. Less nice is to force both curves to display with the same (possibly darker) shading as each other (even if darker than the region).

restart:
RegionBetweenCurves:=proc(f1,f2,rng::name=range(numeric))
    local v1, v2;
    v1:=lhs(rng);
    plots:-display(
      plottools:-transform(unapply([v1,v2+f2],v1,v2))
         (plot(f1-f2, rng, 'filled'=true, _rest)),
      plot([f1,f2], rng, _rest) );
end proc:
RegionBetweenCurves(x^2-1, -x-1, x=-1.5 .. 1.5, color=gold);
RegionBetweenCurves(x^2-1, -x-1, x=-1.5 .. 1.5, color=cyan);
RegionBetweenCurves(sqrt(x-2), x-4, x=4..6);

acer