Axel Vogt

5821 Reputation

20 Badges

20 years, 358 days
Munich, Bavaria, Germany

MaplePrimes Activity


These are answers submitted by Axel Vogt

TAB never worked for me (on WIN, various Maple versions, classic) and use blanks, the editor is not very user friendly

Searching for Groebner will almost immediately show you the keyword (at the 3rd place ...):

?Monomial orders for multivariate polynomials
?tdeg

i am indifferent for that, still ... at least until a reasonable company statement comes up (well, no reason to believe all what may be stated in it then ...)

coefficients:={j = -139.55659155625, k = -3.8511665004543};
select(has,%,j);
op(%);
rhs(%);
                           -139.55659155625

or

eval(j, coefficients);

                           -139.55659155625


Compile as C with __declspec(dllexport) void __stdcall for your function and in case of problems use dependency walker to look into the DLL, whether the function is present.

In case you want C++ separate you usual code and the interfacing into different sources

Here is an old project in MSVC6, just translate to your concurrent MS version:

http://Download 102_MatMult_MSVC6

you can use the Watcom compiler, which comes with Maple, here is an old thread

www.mapleprimes.com/forum/external_callling_and_new_watcom_compiler

 

Maple has a function Li, descrobed in the help: it is afunction of 1 argument.

But Maple has no function Li[a], a = 2/3

Question 1: this is the same as a minimum for  - f , see

www.maplesoft.com/support/help/view.aspx or www.maplesoft.com/support/help/view.aspx

Question 2: no, you should not read it from graphics

you can use indets(%, symbol) or similar to find it out ... but if you want to plot you will have to specify the variable -  what is the actual question?

  restart;
  'Int((3*polylog(5/2, -exp(b*(u-i)))/b^(5/2)+i*polylog(3/2, -exp(b*(u-i)))/b^(3/2))*r^2, r = 0. .. 10.)';
  convert(%, rational); 
  g:=value(%);

       1000 polylog(5/2, -exp(b (u - i)))
  g := ----------------------------------
                       5/2
                      b

                  i polylog(3/2, -exp(b (u - i)))
         + 1000/3 -------------------------------
                                3/2
                               b


Then use your data:

  i := ... # why using Typesetting[delayDotProduct] ?
  i:=convert(i, rational);

  b := 1/43; u := 43; 
  G := 4.3*10^(-3); G:=convert(G,rational); # better 
  M := 10^4; 
  r[m] := 1; # nonsense, just write rm:=1, otherwise you may
             # have problems with the r used for i 
  N := 1; q := M/N; v := sqrt(G*M/r[m]); n := 9/(16*Pi^2*r[m]^3*v^3);
  f1 := (3/10)*q*G*M-(3/5)*q*G*M;  

Finally you get your f, depending on r (and quite lengthy)

  f := -8*Pi^2*n*sqrt((1/2)*Pi)*g:
  indets(f, symbol);
                               {Pi, r}

Even working with Microsoft Word page breaks become dependent of the printer driver, which generally is not part of the application but located at the operating system, thus a preview is only approximative.

If you use Postscript with GSview you can select which pages are omitted on your oddly formatted document and thus can produce a pdf without blank pages - for free.

Interesting how you both crawled into Maple and tortured yourself ...

  b:=6.0*10^13;
  L:=[(.9999999996*cos(1.166666667*10^(-14)*Pi*x)*cos(5.000000000*10^(-15)*Pi*x)-
  1.380952380*sin(1.166666667*10^(-14)*Pi*x)*sin(5.000000000*10^(-15)*Pi*x)), 1];
  plot(L, x = 0 .. 6*b);

The problem is your trigonometric expression, which is L[1] in the above notation:

Numerical it equals (almost) 1:

  L[1]; fnormal(%);
                             0.9999999996

And it gives 1.0 if working with 14 Digits precision (which is where it may come
from, as your exponents 10^(-14) suggest: probably you switched precision through
your code ... it comes from mixing different precisions, yes?).


So what is your actual problem? Do you really want to know the (numerical) values
where the curve intersects the constant line y=1?

There is (possibly) no closed form, so you have to do it numerical.

... look up your manual or help, see plot or www.maplesoft.com/support/help/view.aspx

if you want a sequence, than just use that statement: test2 := j->seq(x[ modp(j-i ,24)], i = 0 .. 8);
First 63 64 65 66 67 68 69 Last Page 65 of 92