In this previous post, an example is shown that demonstrates the potential problems that can arise following symbolic conversions such as from sqrt(x^2)  to x^(1/2).

This provides a Maple solution to compute the bivariate normal distribution by recursions for numerical inputs. It works even for extreme cases and handles situations, where usual integration with Maple has serious problems (even after reducing to dimension 1), it seems to be reliable and fast and works in 'arbitrary' precision.

To use it call N2_as_sum(1.0, 2.0, 0.8,  200) to compute the BVN for x = 1.0, y = 2.0 and correlation rho = 0.8 with at most 200 recursion steps (it will stop earlier, if no more improvements can be seen).

The calling sequence in ?complex states only the two argument form 'Complex(x, y)', but later in the section "Description", in a bit confusing way, two rules for the single argument form are stated:

I need to do the following. Suppose I have an array, c[ ]. I need to be able to specialize the elements of this array, and then later "clear" the elements, making them indeterminates again. I attempted to do this by creating a temporary holding array

i:='i': for i from 1 to 5 do d[i]:=c[i] end do;
i:='i': for i from 1 to 5 do c[i]:=i end do;
i:='i': for i from 1 to 5 do c[i]:=d[i] end do;
 

I expected and want the output to be

c[0] c[1] c[2] c[3] c[4] c[5] 1 2 3 4 5 c[0] c[1] c[2] c[3] c[4] c[5]

Here's a little test double loop

for i from 1 to 2 do

for j from 1 to 2 do

A[i,j]:=Binomial[i+j,i];

end do;

end do;

I would like to get ANYTHING to output. I get NOTHING, no matter WHAT variations I make. I can output when the loop is just a SINGLE nest. But, as soon as I embed that nest inside a second nest - NOTHING outputs.

Are there any commands which will allow me to create zero-based arrays in Maple?

I have a linear recursion, in 4 variables,

c[i,j,k,l]

To the members of the MaplePrimes community,

I am posting this message in response to some comments that have been made about Maplesoft’s commitment and involvement in the MaplePrimes community.  Let me assure you that we at Maplesoft take these comments seriously, and I would like to take this opportunity to personally address this issue.

Hi, I am facing a problem to copy the Maple output into SWP5.0. When I paste the output, generated by 

print(expression,latex);

in SWP it remains mere ASCII. Kindly guide me. Thanx.

Why does the following give me an imaginary number?

 (-.008)^(1/3);

In a concurrent thread I posted the following simplification procedure

  Tryhard:= proc(expr)
    global E_in_Tryhard;
    subs(pow= `^`,
      codegen[optimize](subs(E_in_Tryhard= expr, ()-> E_in_Tryhard), tryhard))()
  end proc;

and now put into an extra blog post, as it might be helpful for others. Note, that this does not work on all constructs in Maple and certainly the package is a bit dated.