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).

Here x is an unknown symbol. The difficulties include the fact that, while `sqrt` can be smart about simplifying numeric values (eg. integers, rationals) the `^` operator has no such opportunity. Once the conversion from `sqrt`...

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).

W := (81/539)*42^(1/2)*hypergeom([-1/2, 3/8, 169/88], [81/88, 19/8], -1/6): evalf(W,50); 0.99999999999999999999999999999999999999999999999996 Show W is 1

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]

Can posts be input using LateX notation?

If so, what would be the TeX tags or delimiters to enclose the TeX input in?

Thanks,

dc

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.

> q := x^(1/2):

> type(q, `^`);
                                     true
> op(0,q);
                                       ^

> subsindets(q, `^`, f->subsop(0=H,f));
Error, (in unknown) improper op or subscript selector

What's wrong with that last one? It's modelled on the first Example in the ?subsindets help-page.

Are there cases where...

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.

Maple 12, Mac OS 10.5 see "invalid character ..." messages! . . How about other platforms?

Dear all, I am working with Maple 11 and 12 on an Intel Mac running Leopard, and I am having trouble with a very simple script that keeps telling "kernel connection lost"; see below. > restart; with(CurveFitting); > poly2 := proc (x, u) return PolynomialInterpolation(x, u, z) end proc; proc(x, u) ... end; > getArrays := proc (left, right) local x, u; x := [seq(h*i, i = left .. right)]; u := [seq(u[i], i = left .. right)]; return x, u end proc; proc(left, right) ... end; > x, u := getArrays(-2, 2); [-2 h, -h, 0, h, 2 h], [u[-2], u[-1], u[0], u[1], u[2]]

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.