vv

14027 Reputation

20 Badges

10 years, 43 days

MaplePrimes Activity


These are replies submitted by vv

To compute a limit could be a very complex task; sometimes almost an "art".
When the expression involves an essential singularity, the result given by Maple has big chances to be wrong.

Unfortunately the help file does not say a word about it, but I am sure that the designers (or at least some of them) are aware of this problem. In these cases, there should appear many "FAIL" answers but this would not be a strong point in marketing.

Here are some real examples (i.e. without a complex direction):

limit((1/2+sin(x))^(1/x), x=infinity);
            1

limit((1+sin(x))^(1/x), x=infinity);
           1

limit( (1+x*cos(x))/(2+x*cos(x)), x=infinity);
           1

@Markiyan Hirnyk 
It seems that you don't understand.
The existence of this f shows that for a CAS it will be very hard to decide whether e.g. the following simple function
g: = sin(exp(-1/z^2+z))*exp(z)-exp(-z)/z^2;
has an essential singularity at z=0 or not.

@Axel Vogt 

The construction of such an f is very technical. It can be done using a deep theorem in complex analysis:
Rosay J.-P., Rudin W. - Arakelian's Approximation Theorem. Amer. Math. Monthly, 96 (1989), 432-434
(the construction itself is not in the article).

@nm 
For a code snippet click the < > icon (near the A icon).

A computation method used by Maple looks very concrete to me.
But the question was not directed exclusively to you.

In this thread I had a concrete question, not an answer. Wasn't this clear enough?

This is just like saying that Riemann Hypothesis is a piece of cake because everyone knows the definition for the zero of a function.

"Borel" was mentioned only by force of habit (being the general context); of course C^1 or even analytic would be enough.

What too general sentences? I have asked about the method used by Maple to compute Probability(...).
You are free to answer if you know something about it.
[I know to compute it by hand, but Maple may use a completely different method].

It seems that Probability has problems and you have found examples in earlier posts.
I would be interested to know how Probability( h(X,Y) < a )  is computed by Maple,  where X,Y are (independent) random variables and h a Borel function. Unfortunately it is very difficult to see this because of the many objects which are generated by the package.

Maple cannot handle essential singularities.
For example, residue does not work in such a point, e.g. residue(exp(1/z), z=0),  even if an undergraduate student knows it.

But I would not expect too much from a CAS in this direction (actually I don't know how Mathematica e.g. works with these problems).

It is not easy to decide whether a singularity is essential or not.
For example there exist analytic functions f for which 0 is an essential singularity and limit(f(t*z), t=0) = 0   for any complex z.

My impression is that MultiSeries:-limit was not designed for complex, and this option should be removed (just use the usual :-limit(...,complex) which seems to work rather well).
All next answers are strange:

with(MultiSeries):
limit(sin(1/z)*z,z=0, complex);
                               0
limit(sin(I/z)*z,z=0, complex); #ok but lucky
                           undefined
limit(GAMMA(z),z=0, complex);
                            infinity
limit(1/z, z=0, complex);
                            infinity
limit(-1/z, z=0, complex);
                           -infinity
limit(I/z, z=0, complex);
                           infinity I
limit(-2*I/z,z=0, complex);
                          -infinity I

 

The Lebesgue-Schoenberg version works easily in higher dimensions:
(For Hilbert/Peano extra work is needed.)

F:=proc(t)
   local u:=abs(2*round(t/2)-t);
   piecewise(u<=1/3,0, u<2/3, 3*u-1,1)
end:
PN:=4:
Pd:=3:  # number of dimensions
for i to Pd do
 P[i] := unapply( add(F(3^(Pd*k+i-1)*t)/2^(k+1),k=0..PN), t):
od:

plots:-spacecurve([P[1](t),P[2](t),P[3](t)],t=0..1, numpoints=100000);

@Carl Love 

But then e.g. a symbolic integration is out of the question.

P.S. I agree that this is not a bug sensu stricto but Maple should do much better in such situations.

@rlopez 

Isn't 'singsol' = all a more natural option?

First 124 125 126 127 128 129 130 Last Page 126 of 177