@Samir Khan Thanks!

(Grist for the efficiency mill.)

Please forgive me if I've missed it in some hidden code block or region in that uploaded worksheet, but I don't see any code to reproduce these images in Maple.

I don't see how anyone could properly assess whether Maple can produce such plots or images using a reasonable amount of time and memory resources, without source code or a worksheet which successfully reproduces the results.

acer

This is very good news for Maple.

The first few areas I would suggest are,

• Advanced plotting: The toolset is quite powerful, but since there will always be many custom needs there will never be a canned solution for every task and problem. Showing how the plotting facilities can do wonders, when combined with even a little programming (not a dirty word), would be a great addition.
  
  
• Applied DEs: Applied problems with symbolic DE solutions, or at least some measure of mathematically driven symbolic analysis. (One topic that comes to mind is Control. Another is Delay DEs.)

acer

@epostma Did you intend to include the option 'output'='residualsumofsquares' in the `Fit` call in procedure `minsumsq`?

Like this. say,

minsumsq := c -> Statistics:-ExponentialFit(ln~([3.05, 3.1, 3.75] -~ c),
                                            [.74e-4, .1806e-3, .584e-4],
                                            output=residualsumofsquares):

plot(minsumsq,-10..3.04);

Just because a fitting problem does not have finite optimal parameter values does not mean that there is no limiting curve for the fit that it continuous between the original independent data points. For example, the constant function x -> 0.00009207055709 may approximately "fit the bill" in such a limiting sense.

All pretty moot, anyway, for 3 data points and only conjecture as to what the poster wanted.

ps. If the error weighting were not by exponential fitting then there may be finite optimal parameter values here.

@epostma Did you intend to include the option 'output'='residualsumofsquares' in the `Fit` call in procedure `minsumsq`?

Like this. say,

minsumsq := c -> Statistics:-ExponentialFit(ln~([3.05, 3.1, 3.75] -~ c),
                                            [.74e-4, .1806e-3, .584e-4],
                                            output=residualsumofsquares):

plot(minsumsq,-10..3.04);

Just because a fitting problem does not have finite optimal parameter values does not mean that there is no limiting curve for the fit that it continuous between the original independent data points. For example, the constant function x -> 0.00009207055709 may approximately "fit the bill" in such a limiting sense.

All pretty moot, anyway, for 3 data points and only conjecture as to what the poster wanted.

ps. If the error weighting were not by exponential fitting then there may be finite optimal parameter values here.

I suspect that what's Axel is driving at is,

expr:=R*D^2 / ( R*D^2 + s*L + s^2*R*L*C );

                                           2         
                                        R D          
                        expr := ---------------------
                                   2          2      
                                R D  + s L + s  R L C


1/expand(1/expr);

                                      1        
                              -----------------
                                          2    
                                  s L    s  L C
                              1 + ---- + ------
                                     2      2  
                                  R D      D   

acer

I suspect that what's Axel is driving at is,

expr:=R*D^2 / ( R*D^2 + s*L + s^2*R*L*C );

                                           2         
                                        R D          
                        expr := ---------------------
                                   2          2      
                                R D  + s L + s  R L C


1/expand(1/expr);

                                      1        
                              -----------------
                                          2    
                                  s L    s  L C
                              1 + ---- + ------
                                     2      2  
                                  R D      D   

acer

Which operating system are you running?

When this problem occurs, is it always the case that there is a worksheet which has been already opened with Maple 15 and minimized to the tray? Is so, then can you maximize that? If you cannot maximize that with the mouse-pointer, and if you are on MS-Windows, then can you maximize if using the <Alt>-hold and <tab> key press?

acer

@PatrickT Please do follow the link in the Comment by pagan. It shares themes with that fast-complex-argument-images Post -- the fastest way to assemble high quality, high density images of mathematical functions in Maple is to act inplace upon hardware datatype rtable with the Compiler.

Other techniques for computing fractal images used in some Application Center uploads (`option hfloat`, or Threads or Grid -- often submitted to use some method that was in vogue) are alone not nearly as fast or efficient.

It may be possible to combine Threads/Task and compiled routines. But since it's not obvious that the Compiler's external-call runtime is thread-safe, it might be necessary to use a Poor Man's Compiler [here & here]. Different techniques are needed for two cases: when the user has already pre-Compiled the mathematical function which must then be used to populate the rtable (eg. for regular plot creation), and when the function and rtable-populating are combined into one single procedure (eg. fractal image creation).

I'd like to finish and submit a posting on the first of those two: fast generation of high density plots/images using a pre-Compiled procedure.

One day this may all be automatic... and people reading these esoteric methods may feel a strange sensation similar to what we can now get re-reading old books on specialized non-digital exposed-film-camera techiques.

None of this relates directly to your projection problem. I asked about curves-vs-surfaces because I was immediately worried about the overlay. I haven't thought of an effective solution to that issue.

acer

@PatrickT Please do follow the link in the Comment by pagan. It shares themes with that fast-complex-argument-images Post -- the fastest way to assemble high quality, high density images of mathematical functions in Maple is to act inplace upon hardware datatype rtable with the Compiler.

Other techniques for computing fractal images used in some Application Center uploads (`option hfloat`, or Threads or Grid -- often submitted to use some method that was in vogue) are alone not nearly as fast or efficient.

It may be possible to combine Threads/Task and compiled routines. But since it's not obvious that the Compiler's external-call runtime is thread-safe, it might be necessary to use a Poor Man's Compiler [here & here]. Different techniques are needed for two cases: when the user has already pre-Compiled the mathematical function which must then be used to populate the rtable (eg. for regular plot creation), and when the function and rtable-populating are combined into one single procedure (eg. fractal image creation).

I'd like to finish and submit a posting on the first of those two: fast generation of high density plots/images using a pre-Compiled procedure.

One day this may all be automatic... and people reading these esoteric methods may feel a strange sensation similar to what we can now get re-reading old books on specialized non-digital exposed-film-camera techiques.

None of this relates directly to your projection problem. I asked about curves-vs-surfaces because I was immediately worried about the overlay. I haven't thought of an effective solution to that issue.

acer

@PatrickT Are you asking why that `cast` doesn't work the same whether or not you put it into a call to typeset(), in a plot caption? I think it just needs some tweaking and correcting, to satisfy some pickier and less forgiving parser.

Consider, in a Worksheet in Maple 15 and outside of any typeset() environment for plot captions, etc. With the original `cast` above,

seq(cast(i*Pi/4),i=1..4);

Typesetting:-Typeset('%');

So now let's edit it a bit.

cast:=proc(r)
local n,d,t,top;
  n,d:=numer(r),denom(r);
  if type(n,`*`) then
    top:=cat(`mrow(`,cat(`mn("`,convert(op(1,n),string),`")`),
             seq(cat(`,mn(`,convert(op(t,n),string),`)`),t=2..nops(n)));
  else
    top:=cat(`mrow(`,cat(`mn("`,convert(n,string),`")`));
  end if;
  if d<>1 then
    cat(`#mfrac(`,top,`),`,cat(`mn("`,convert(d,string),`")`),`)`);
  else
    cat(`#`,top,`)`);
  end if;
end proc:

seq(cast(i*Pi/4),i=1..4);

Typesetting:-Typeset('%');

Anyway, you could give it a go.

@PatrickT The help-page for topic `value` appears to already have a help-alias for `ditto`.

In Maple 15 with the Standard GUI, if I enter ?% in a worksheet (or if I search for & in an open Help window) then the Help system opens to the `ditto` help-page, sure. But that Help browser's "search" tab on the left shows two matches. One match is for topic `ditto`, and the other match is for topic `value`. So this aspect is already working.

And only one topic can be the default, which actually gets opened automatically. It's understandable that it's the page of the more commonly used ditto functionality.

Maybe the GUI could learn from Wikipedia here: when an ambiguous query is submitted then a default page is shown, but at the top of that page there is a brief note and link to a disambiguation page for that query. The same could be done with Help aliases: if there is more than one match then a brief line and disambiguation link could be inserted on-the-fly into the displayed help content. Not everyone looks at the search/match results, when there is a page opened automatically in the larger right display pane.

@PatrickT The help-page for topic `value` appears to already have a help-alias for `ditto`.

In Maple 15 with the Standard GUI, if I enter ?% in a worksheet (or if I search for & in an open Help window) then the Help system opens to the `ditto` help-page, sure. But that Help browser's "search" tab on the left shows two matches. One match is for topic `ditto`, and the other match is for topic `value`. So this aspect is already working.

And only one topic can be the default, which actually gets opened automatically. It's understandable that it's the page of the more commonly used ditto functionality.

Maybe the GUI could learn from Wikipedia here: when an ambiguous query is submitted then a default page is shown, but at the top of that page there is a brief note and link to a disambiguation page for that query. The same could be done with Help aliases: if there is more than one match then a brief line and disambiguation link could be inserted on-the-fly into the displayed help content. Not everyone looks at the search/match results, when there is a page opened automatically in the larger right display pane.

@PatrickT ...you wouldn't have to call eta(r,u,a,A,d,r,s). You could call it as either,

  eta(r, u, subs(Parameters1,s));

or, after assigning pars:=subs(Parameters1,s) at the top of your worksheet (using `op` if the expression sequence inside some list/set if needed), just calling it as,

  eta(r, u, pars);

Sorry if it seemed like nit-picking, but the idea was that you might not want to call eval of a proc template each and every time you have actually to call the proc, if it must be done many times.

Joe's suggestion, using globals, also seemed relatively straightforward and simple. I don't know that there's a best way.

@PatrickT ...you wouldn't have to call eta(r,u,a,A,d,r,s). You could call it as either,

  eta(r, u, subs(Parameters1,s));

or, after assigning pars:=subs(Parameters1,s) at the top of your worksheet (using `op` if the expression sequence inside some list/set if needed), just calling it as,

  eta(r, u, pars);

Sorry if it seemed like nit-picking, but the idea was that you might not want to call eval of a proc template each and every time you have actually to call the proc, if it must be done many times.

Joe's suggestion, using globals, also seemed relatively straightforward and simple. I don't know that there's a best way.