It always makes me happy to see people using Maple for interesting things.  So I was pleased to see this blog post on Technology Review about this paper on arXiv on quantum randomness.  In this case, they are just comparing random numbers generated from lasers (this is why physicists get better press than mathematicians: LASERS!) with pseudo-random numbers generated using Meresenne Twister implemented in Maple, pseudo-random numbers generated using a Celluar Automata method implemented in another computer algebra system, and then binary digits of π treated as a pseudo-random sequence.  (Spoiler: the lasers win)

While not a particularly interesting use of computer algebra systems, it did inspire me to revisit my old blog post on pseudo-random numbers in Maple and now I am working on a follow up that talks about some of the mathematical and statitical tests used to test the quality of pseudo-random number sequences which I hope to post soon.

Hi, Can anyone help me in solving the following system of ode, i'm new in maple and i faced some problems in using it.

the system:

diff(u2(y), y, y)+m*b*rho*h^2*GR*(c3*y+c4) = 0,

diff(u1(y), y, y)+K*(diff(N(y), y))/(1+K)+GR*(c1*y+c2)/(1+K) = 0,

-2*K*(2*N(y)+diff(u1(y), y))/(2+K)+diff(N(y), y, y) = 0

with the following boundary conditions

u1(-1) = 0, u2(1) = 0, u1(0) = u2(0), (D(u1))(0)+K*N(0)/(1+K) = (D(u2))(0)/(m*h*(1+K)), (D(N))(0) = 0, N(-1) = 0

Recently I wanted to create a simple 2D-plot representing this piecewise expression:

T1:=piecewise(x>=0, 3.5*x*13/3+13, x>=20, 3.4*x*13/3+13, x>=30, 3.25*x*13/3+13);

or alternatively: T1:=t->piecewise(t>=0, 3.5*t*13/3+13, t>=20, 3.4*t*13/3+13, t>=30, 3.25*t*13/3+13);

Unfortunately, no matter whether I entered it as a function or as an expression, I was not able to make Maple plotting this function correctly for a domain D

Hi,

Please take a look on this problem I have about converting a Vector to a List. At the end if the program, the list F contains the elements of the Vector Y, but not in the good order.
It gives a result which is correct for a small size of Vector, but not for h = 1000 !

#differents values of h are tried h := 100 and h:=1000
h:=1000;

E := Vector[1..h];
for u from 1 to h do
Y[u] := u;
end do;
for z from 1 to h do
   E[z...

Is it possible to make a plot like this:

This sites needs better rendering of 2D Math.

Here is a sample from the Wikipedia page for Maplesoft, from the pdes example section. The image on Wikipedia is quite nice (it's alt-tags have LaTeX code, which may be a hint).

Hi

I need to insert a legend in a Pie chart, but I want to do it from the command line:

> with(Statistics);
> T := [algodón = 26, lana = 5, tabaco = 7, fruta = 9, arroz = 13, té = 5, carne = 14, otros = 21];
> P := PieChart(T, sector = 0 .. 360, color = ["red", "Orange", "yellow", "green", "ForestGreen", "blue", "MediumPurple", "LightBlue"], captions = relative);
> plots[display](P, title = ["Distribución de las exportaciones...

I have the following differential equation:

de1:=(v(r)^2*(diff(T(r),r)/(1+T(r)^2)+T(r)/r)-T(r))*(v(r)^2*(diff(T(r),r)/(1+T(r)^2)+T(r)/r)+T(r))=(1+v(r)*diff(v(r),r))*(1-v(r)*diff(v(r),r));

I would like to be able to eliminate T' and obtain T as a function of r, v(r), and v' (and possibly v'').

Is this possible?

Thanks, Geoff

I am having a problem working with a slider and a plot component in Maple 14. The slider specifies the domain of an animation I am trying to show. Whenever the slider's value is changed, I have code that should change the plot property "play" to "true" so that the animation begins. However, it fails to play the animation. After some troubleshooting, I found that the exact same code works in a Toggle Button component, and I also discovered that the Slider...

There is a probem in the Optimization package's (nonlinear programming problem) NLPSolve routine which affects multivariate problems given in so-called Operator form. That is, when the objective and constraints are not provided as expressions. Below is a workaround for this.

Usually, an objective (or constraint) gets processed by Optimization using automatic differentiation tools ...

Hi there. I have two expressions of type Matrix(7x1), one I entered directly and the other one is the result of a  procedure calling. They only differ by the term [1,1], so when I do a comparation of the appropiate submatrix using Equal I get a positive answer:

>Equal(SubMatrix(atataaV(1,1,1,1,u1),2..-1,1),SubMatrix(u1,2..-1,1));

true

Now, the elements of the submatrices depend on the symbol 'n',...

I'm brand new to Maple. How can I plot the Gammafunction with real numbers (not in N and C)?

I can only plot the Gammafunction in the complex plane by using the following:

with(plots);

f:=z->GAMMA(z);

complexplot3d(f);

Is it also possible to plot the Betafunction B(a,b):=int(t^(a-1)*(1-t)^(b-1),t=0..1) in such a way that I can assign special values to a and b?

These solutions will help me so much....Thank you for your answer!!!

I often run into this problem. Say I have defined a new variable in terms of some old variables:

     Eq1 := NewVar = OldVar1 + OldVar2

Then suppose I have an expression in terms of the old variables:

     Eq2 := OldVar6 * OldVar5 = OldVar3 * (OldVar1 + OldVar2)

If I want to substitute Eq1 into Eq2 I have to either type

     subs(OldVar1 + OldVar2 = NewVar1, Eq2)

or I can do this

Hi,
I need your help for this problem, please.

$\mathop  {\sum {...\sum }}\limits_{z_r \in {Q_{2}^{*}}}$ such that

$Q_{2}^{*}=((z_1,...,z_r):z_{j}=0,1,\sum_{j=1}^{r}z_{j}=\ell)$, where $\ell = 0, 1, . . . , r$.

For example. If r=3, $\ell=1$; then we are looking for triples with elements 0's or 1's and sum to 1 i.e. $Q_{2}^{*}={(0,0,1),(0,1,0),(1,0,0)}$. So we would like to sum over theses $z_r$'s.

I appreciate your help.
Hasan