Items tagged with bug bug Tagged Items Feed

I think I have found a subtle bug in my Maple 9.

If I execute the following code:

###
 B:=[ 3.5, 0.0001, 73.45]; A:=[[seq(1/op(j,B), j=1..3)]];
 writedata("a.txt",A);
###

the output file a.txt contains the line
0.2857142857    1       0.01361470388
(beware: maple has written "1" instead than "10000"!)

Questions:
1) is this bug already known? Is there a patch?
2) does it affect newer versions of Maple?

I have come across this previously but I never had the energy to report it.
The problem is that there exist a bug in Vector[row] when n>10.

For example run the below code:

restart:
with(LinearAlgebra):
a1 := Vector[row]([seq(w[i], i = 1 .. 5)]);
a2 := Vector[column]([seq(w[i], i = 1 .. 5)]);

a3 := Vector[row]([seq(w[i], i = 1 .. 20)]);
a4 := Vector[column]([seq(w[i], i = 1 .. 20)]);

a1.a2;
a3.a4;

   ...

Dear,

======================== This code will give error ============

restart;

with(Statistics);

i:=1:

RandomVariable(ProbabilityTable([1/2,1/2]));

=========================

I  used Maple12 and had a hard time debugging my code to find the error . Finally, I found that when we assign the name i a value, then the RandomVariable(ProbabilityTable([...]) will give error message below:

Error, (in Statistics:-Distribution...

with(Statistics):
Phiinverse := (x) -> Quantile(Normal(0, 1), x) :

evalf(Phiinverse(1-1/2^24));

and the result is 5.691523310-2.147372008*I

Now, I can get a real-valued solution to this by using the equivalent
functionality from stats instead of Statistics.

with(stats):
PhiinvClassic:= (x) -> statevalf[icdf, normald[0, 1]](x):

evalf(PhiinvClassic(1-1/2^24));

test.zip

1) Compile under vc express 2010 debug config. Don't forget point out your location of maplec.lib in configuration properties/Linker/General/Additional library directories. Assure that configuration is at least Debug or all.

2) Run sheet test_bombard.mw from Debug directory.

 

Have permanent crash on BombardByRandoms call. If i would transmit creating of rtable from 'InitInternal'...

In regards to executing a command to output text in color I have stumbled across a minor bug?

Tools -> options -> display -> typsetting level : Extended

Typsetting[Settings](autoatomic) #checking default autoatomic setting (don't know if this really matters though)
                 true

Type in    testing    - select the...

I'm having the problem reported here back in 2008 and more recently in 2010:

http://www.mapleprimes.com/questions/39313-Postscript-Plot-Sizing

http://www.mapleprimes.com/questions/96487-Graphics-Output#comment96501

I'm not sure when the axiswidth, axisheight options of plotsetup got broken for me. I'm not certain, but it looks like I was able to control the width and height of my plots not so long ago (version 14 or 13) or am I mistaken and this bug,...

When I type

cos(x), sin(x)

and start the Plot Builder (using the contextual menu) the first option it gives me is for a 2-D plot.  However, when I type

cos(x), x

and start the Plot Builder the first option is a 2-D polar plot and the "standard" 2-D plot is nowhere on the list.  Is this a known bug?  (I'm using Maple 15.01 on Windows 7.)

 

Thanks,

I have the following Maple worksheet:

restart;
randomize();

select1 := combinat[randperm]([1, 2, 3, 4, 5, 6, 7, 8, 9, 10]);

Because of the "randomize" function, each time I execute this worksheet I should expect a different output. However, the real output is always the same every time. Is this a bug in combinat[randperm]() or am I missing something?

Attached file: randperm1.mw

I'm running simulations in which I produce many plots at once, the purpose being to compare differences as parameters are changed. I'd also like to use some of these plots for publication. I've had problems getting good-looking 3d plots. Here's a strange thing I noticed.

Export the graph by mouse-clicking, export to bitmap:

EDIT : IMAGE POSTED IN MY COMMENT BELOW

Export the garph by executing a plotsetup('bmp') command.

 

This is a continuation branch of the thread "cut/paste" won't work..." See point #1 below.

1. Why won't #$%^ Mapleprimes editor let me reply to the original post? I press "add comment" and I get pointed back to the branch two steps back. Anything I write comes out blank when I submit, so now there are two blank replies credited to me! Does that give me points I don't deserve?

2. Classic worksheet isn't "still there". It says so in the readme that classic worksheet...

The order of 1 in any finite field (that I tried) created by GF is NULL. For example,

F:=GF(3,2):
use F in order(one) end;

It should be 1.

Alec

It looks like for me using shift+tab does not work in linux when one selects 1D as the default input notation and is typing in a worksheet. 

Does anyone know where we should submit bug reports for things like this?


thanks.

Patch request - Patch work on bugs for Maple ceases when a new version of Maple nears release, possibly months before the actual release.  I think it's safe to say no work will ever be done to patch up older versions.  But how can we, the user, introduce our own fix code into maple?  Because we can't fix the built in code how can we temporarily or permanantly make our own patches?

For example prior to Maple 14 a bug exists in a ' for - in ' loop. ...

Earlier this week a colleagues approached me with a problem he was trying to solve with Maple.

The basic problem is to compute the asymptotic expansion of an expression. It's not difficult to see the leading term of this expansion, but Maple's result does not match this expectation. Let me demonstrate:

B := (n,h) -> (n+h)!/r^(n+h) * (1+(P0+h*P1)*r/n);
                           /    (P0 + h P1) r\
          factorial(n + h) |1 + -------------|
                           \          n      /
(n, h) -> ------------------------------------
                         (n + h)             
                        r                    

simplify( asympt( B(n,1)/B(n,0)-1, n, 3 ) );
                            n + O(1)
asympt( simplify( B(n,1)/B(n,0)-1 ), n, 3 );
                                   2                         
 n   1 + r P1 - r   P0 r + r P1 - r  P0 + (-1 - r P1 + r) P0 r
 - + ------------ + ------------------------------------------
 r        r                            r n                   

      /            2   \          
      \-r P1 + P0 r  P1/ P0    /1 \
    + --------------------- + O|--|
                2              | 3|
               n               \n /

 

The real issue is the drastically different results from the previous two commands.  The only difference in the commands is the explicit simplify before or after the call to asympt.

It would be best if this last expression could be seen in the form:

(n+1)/r + (P1-1) + (1-r*P0)*P1/n - r*(1-r*P0)*P0*P1/n^2 + O(1/n^3);
  n + 1            (1 - P0 r) P1   r (1 - P0 r) P0 P1    /1 \
  ----- + P1 - 1 + ------------- - ------------------ + O|--|
    r                    n                  2            | 3|
                                           n             \n /

What's the difference in the two uses of asympt? This simpler example appears to show the essence of what is going on. On one hand,

B := (n,h) -> 1/r^(n+h);
             1   
(n, h) -> --------
           (n + h)
          r      

simplify( asympt( B(n,1)/B(n,0)-1, n, 3 ) );
Error, (in asympt) unable to compute series

while on the other hand

asympt( simplify( B(n,1)/B(n,0)-1 ), n, 3 );
                              r - 1
                            - -----
                                r 
2 3 4 5 6 7 8 Page 4 of 8