Data:

My name is Darin Ohashi and I am a senior kernel developer at Maplesoft. For the last few years I have been focused on developing tools to enable parallel programming in Maple. My background is in Mathematics and Computer Science, with a focus on algorithm and data structure design and analysis. Much of my experience with parallel programming has been acquired while working at Maplesoft, and it has been a very interesting ride.

In Maple 13 we added the Task Programming Model, a high level parallel programming system. With the addition of this feature, and a few significant kernel optimizations, useful parallel programs can now be written in Maple. Although there are still limitations and lots more work to be done on our side, adventurous users may want to try writing parallel code for themselves.

To encourage those users, and to help make information about parallel programming more available, I have decided to write a series of blog posts here at Maple Primes. My hope is that I can help explain parallel programming in general terms, with a focus on the tools available in Maple 13. Along the way I may post links to sites, articles and blogs that discuss parallel programming issues, as well as related topics, such as GPU programming (CUDA, OpenCL, etc).

My next post, the first real one, I am going to explain why parallel programming has suddenly become such an important topic.

I was wondering how I can get the equation for Probability Density Function (PDF) for a variable x which follows

a normal distribution with mean 0 and stdev 10 but also has a stop at -10 (ie all x's that are more negative

than -10 becomes -10). I can plot such a function with the below code but I dont know

how to extract the equation for the PDF...humm

I'm trying to plot directly the result of a solve, using the following syntax:

plot(solve(diff(E,N)=0,N),Z=0..20);

( E is dependant of Z and N )
 

it works for an equation that has only 1 solution. But, when using an equation with two solutions, it only plot one of the solutions. Of course, it is always the solution I don't want.

I've tried solving this problem by forcing N to be above 0 ( the unwanted, but always plotted solution is negative ):

plot(solve({Ep-E=0,N>0},N),Z=0..20);

Hello,

I am using the LinearAlgebra package, and I am dealing with very big square matrices (almost a hundred rows and columns).  They have a lot of eigenvalues, with high multiplicity.  

I am performing operations on the eigenspaces.  So what I would need is a matrix with as columns a basis with eigenvectors for a very specific eigenvalue.

I can of course use SubMatrix(Eigenvectors[2],....) but that would take ages and Eigenvectors gives the eigenvalues in an ever changing order.

Hi all,

I just recently started using maple to supplement my math work. When I was just trying some random equations and graphing them, I tried to plot some 3D graphs. When I tried a simple equation, x+y from -2..2, I recieved this image. http://picasaweb.google.com/lh/photo/wnW4UIkdJck46smuTAnb4Q?feat=directlink

Recently I came across a page that was working with ifactor and it seems op now handles the operations a bit different now. 

ifactor(3^43+1);

       (2) (82064241848634269407)

a:=op(2,%):

ifactor(3^41+1);

       (2)(270547105429567)(33703)

b:=op(2,%):

c:=a*b;

     (82064241848634269407) (271547105429567)

I am using fsolve in maple 13 to find the crossing of a pair of implicitly defined functions in a triangular region.  I have been  using the startpoint option, but sometimes fsolve wanders out of bounds.  I have tried specifying the range as {x=0..1, y = 0..x}, but fsolve complains of an invalid range.  Any suggestions?

Thanks!

Michael

Maple experts,

(I am using Maple10)

I have the following system of equations, but I cannot obtain the solution. (I am able to solve by hand, but I want to verify and to solve the same sorte of equation but with more variables).

solve(    {z1s- p*y1s-(1-p)*y1=0, z2s- p*y2s-(1-p)*y2=0, z3s- p*y3s -(1-p)*y3= 0, y1s+y2s+y3s=1,   y1s-thetas1*(1- z1s-z2s) -(1-thetas1)*z1s=0 } , [y1s]);

What is my mistake?

Thanks in advance for the help.

Gabriel

Pretty simple --- in Maple 13, I try plot(x^2, view = 1.995..2.005,3.99..4.01) and I get the below.  Any idea what's going on?

I use Maple 10 to calculate eigenvalues and eigenvectors:

K:=Matrix(n,n,datatype=float,storage=sparse):
M:=Matrix(n,n,datatype=float,shape=diagonal):
# code to assign the non-zero elements of K and M
(lambda,v):=Eigenvectors(K,M):
 

For n>1000 the calculation takes several minutes (on 2.4 GHz processor, 2 GB RAM). Maple calculates all eigenvalues, but I only need the ten smallest eigenvalues (and the corresponding eigenvectors). Is there a way to increase the speed of the calculation?

I need to use a function that transform a 3 coordonates point (x,y,z) into a 2 coordonates point (x,y) on a specified projection plane, which I called fg. I have the function, but I need to apply it to over a thousand points. Is there a more efficient way than what I did below to do so?

Help creating audio file from Array tia sal22

Greetings All
 
I have a text file with two columns the 1st column which goes from 0 to 1 
represents time in seconds, the 2nd column which goes from 0 to 0.9 represents amplitude.

 I would like to create a wave file from the data and export it to a stereo file at 44100 (cd quality).  How can I go about doing this?  I looked up the maple audio tools and they don't have a place to use the x (time) values they only want the y (Amplitude) values.