Axel Vogt

5408 Reputation

20 Badges

17 years, 226 days
Munich, Germany

MaplePrimes Activity

These are replies submitted by Axel Vogt

@acer sorry & thank you, I used 2017.3


@Mariusz Iwaniuk you might give a reference

@vv : how you came up to involve the Catalan constant?

@Carl Love 

Oh ... that is worth to be posted as an extra thread.

Yes, multiply numerator and denominator by the conjugate here is the same as 'rationalize'.

Then the "complex division" vanishes and one proceedes as usual (guessing the x is some real), if one knows how to multiply complex numbers:

1/z = (1*conjugate(z)) / (z*conjugate(z)) = complex/real

1/(5*x-I): %= rationalize(%);

                            1       5 x + I
                         ------- = ---------
                         5 x - I       2
                                   25 x  + 1


Maple should stop that cancellations (except spam).

Please consider to introduce moving threads to some trash folder, but do not delete.

If one uses MultiSeries:-asympt(f, q, 4) for f = integrand the constant 8*p prevents convergence in infinity

... the following is meant (I use 1-D Math input [finding it better to handle] and do not use labels)

restart; interface(version);

`Standard Worksheet Interface, Maple 2021.2, Windows 7, November 23 2021 Build ID 1576349`


DepVars := [u(x, y, t), U(xi, eta)];
alias(u = u(x, y, t));

[u(x, y, t), U(xi, eta)]




xi1 := 1/2*(x^2+y^2);
xi2 := t;
u := (h(t)+(x^2+y^2)*(1/2))*arccos(x/sqrt(x^2+y^2))/t+U(xi1, xi2);





(h(t)+(1/2)*x^2+(1/2)*y^2)*arccos(x/(x^2+y^2)^(1/2))/t+U((1/2)*x^2+(1/2)*y^2, t)


dd := '(diff(u, x))*(diff(u, y))';

(diff(u, x))*(diff(u, y))


'dd'= ``;
eval(dd, [t = eta]):
simplify(%) assuming x::real, 0 < y:
``=algsubs(1/2*(x^2+y^2) = xi, %):
simplify(%) assuming 0<xi: # else nothing is known about xi

%, xi = 1/2*(x^2+y^2), eta=t;

dd = ``


`` = (1/2)*(arccos((1/2)*x*2^(1/2)/xi^(1/2))*x*xi+(D[1](U))(xi, eta)*eta*x*xi-(1/2)*y*(xi+h(eta)))*(2*arccos((1/2)*x*2^(1/2)/xi^(1/2))*xi*y+2*(D[1](U))(xi, eta)*eta*xi*y+x*(xi+h(eta)))/(eta^2*xi^2), xi = (1/2)*x^2+(1/2)*y^2, eta = t






@mskalsi try it using "eval" instead of "subs"

For example: ... assuming 0<p  and read the help for that.

Using a search engine gives his publications (where I do not want to register) Downloading a paper should give you his concurrent contact information. Likewise you can try similar for his co-author Fauser. Good luck.

Edit: or contact him using if you are a member of Linkedin

Why does one need "Operators" ?


You seem to expect that the usual rules are correct if working with floating point calculation (this is what tomleslie means).

But this is only true with rounding errors and what you see is just that, that's all.

Here is a simple example where (a+b) + c and a + (b+c) are different:

restart; interface(version); Digits:=10;

`Classic Worksheet Interface, Maple 2017.3, Windows, Sep 27 2017, Build ID 1265877`

Digits := 10


eps:=10^(-Digits); evalf(%);

eps := 1/10000000000




(1+eps) - eps;



1+eps; evalf(%);
% - eps;

1 - %;







Download example_rounding_error.mws

Edit: and here an example for multiplication giving an error in the last decimal place, (a*b) * c and a * (b*c) are different



May be you search for some lecture notes about Complex Analysis. It is not important which one, the contents are similar (though the style varies). Do not only try to see what Maple provides, some Math would help you for that.


That list is certainly not complete (square roots and friends) - but what do you *actually* want to know / achieve?

1 2 3 4 5 6 7 Last Page 2 of 199