This is a follow-up to a rather old posting about maple.vim.  There is now a maplevim project on google code to helpfully get the community to work on upgrading maple.vim a bit more frequently. 

Right now there is a shipping version of maple.vim checked-in, as well as a set of 'pieces' (in the pieces sub-directory...

While reviewing code the other day, I came across the following snippet (here converted to a procedure).

Ds := proc(V::set, n::posint, t)
local i,v;
    {seq(seq((D@@i)(v)(t), i=1..n), v in V)};
end proc:

The purpose of this is to generate a set of derivatives at a point of a set of unassigned names. For example

 Ds({x,y},2,0);
               ...

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!!!