Hello wizards,

I'm given to understand that using add() or if possible Task:-Add() is more efficient than a FOR...DO loop. Today I'm asking about the limits of this generalization. My illustration is probably missing some evalf's, but hopefully the concept is clear:

A computation I'm working on involves between 10E4 and 10E6 computations like

Maple's Statistics package contains many predefined probability distributions; well-known ones such as the normal distribution and lesser-known ones such as the Gumbel distribution. For these distributions, we ship efficient algorithms that can quickly generate a large number of sample points. To generate a sample of size 10⁶ of both of these distributions, and print the time it took to do this (in seconds), you can run the following:

with(Statistics):

Consider the following C code:

Hi, 

How can I create a matrix M by M with maple, where the elements of diagonal comes from b[i]=i+1 and the other elements comes from c[i,k]=i+k.

Amr

Using the Physics package I am trying to define commutation rules between operators. However it looks like when the commutator involves a Dagger, it is not "remembered" properly. Here is a (contrived) example (you might have to run it on your machine in case not clear from copy/paste what's going on):


restart:
with(Physics):

Setup(quantumoperators={a, b, c});
                ...

I've received a suggestion offline saying that I should animate the Pacman. There are two items to animate: Pacman's mouth and pac-dots. To review how pacman's mouth moves, please have a quick look at Google Pacman doodle.

Pacman's mouth only has three stages: open, half open and close. Let's take the pacman that I created in the previous post. Each frame of the animation is done by adjusting the radian values of the "pie" (aka Pacman's body and mouth).

This describes in more detail (with permission and citation) the method posted here of solving problem 2 in the XKCD comic Substitute.

Download RaptorMath.mw

Today is my birthday, and in fact it is also the birthday of at least one other Maplesoft employee (not surprising since more than 23 people work here - considering the generalized birthday problem, I even know of 3 people here who share the same birthday).  Of course, it turns out that birthdays are not evenly distributed through out the year and so I wanted to know if someone with an August birthday is more likely to share than someone with an April birthday. 

I would like to include the current date in the filename when I use writedata.  How can I do that?

I know this doesn't work but I was trying to do something like this

  test:=[1,2,3,4,5,6]
  writedata("c:/data-FomatTime(%B%d%H%M).txt",test)

Look at this article (Its loading can take a few minutes: approximately 5 MB.). All those pictures could be created with Maple.

PS. Working link to  the article.

Motion_of_a_Charged_.mw

Hi,

The procedure I uploaded, was originally written in C++. I made a text file and plot the file with C++ without any problem.

But in the maple version of the same program, when I use the command:

                          ...

I am using Maple to generate PDEs. I have expressions of the form p(dx,dy,dz, etc.) where p is a nonlinear polynomial. I would like to substitute d/dx for dx etc and then later apply the resulting differential equation to a function. I am having multiple problems with this.

1. Maple seems to be unable to handle a derivative unless there is a function for the derivative to act on. (This is not a huge issue, but it would be helpful.)

2. Even if I do specify a function Maple maps dx^3 to diff(f,x...

I would like to know how to plot different series of data on the same graph.

Thanks in Advance

Hello,

I have a grafic, a power curve, and a few values of the curve.

The values are in 0.5 steps. I would like to calculate the values for all the 0.1 steps.

For example:

I know the values for 4, 4.5, 5 etc. and would like to find out the values for 4.1, 4.2, 4.3 etc.

Is this possible using Maple?

Best regards,

Calin

When I plot csc(k*x) for various values of k, the timing varies considerably depending on the value of k.

Here are the timings for k = 1..5 (note the strange case of k = 3):

  seq(time(plot(csc(k*x), x = 0 .. 2*Pi, -3 .. 3, discont = true)), k = 1 .. 5);
               0.031, 0.031, 14.133, 0.046, 0.062

For k = 6, the time is apalling: