Hi all

I have used pointplot3d to plot a set of points using three vecotrs v1, v2 and v3 which all contain many elements and the same number of elements.

I would like to know how to plot a new set of elements for new vectors v4,v5, and v6 on to the same plot as the previous set of values.

Could I then give each of the two sets a different color?

Hello

I'm working on implementing a plot of the Lorenz equations, and I'm trying

  proc(something, something, more something) local something more;

  uses DEtools;

  something;

  end proc

And I get this response:

  Error, `DEtools` is not a module or member

What can I do? "uses" works for plots, LinearAlgebra etc. Why doesn't it work for DEtools?

I have this non-linear ODE in q(p):

ode := diff(q(p), p) = (q(p)-a*p)/(p-a*q(p))*(p^(-b)-1)/q(p)^(-b);

for some real parameters a and b, it may be solved in exact form for the special value a=1.

Is there any hope of getting a closed-form solution or of expressing the solution in a neat way for the case 0<a<1? even some specific value of a would be of interest, e.g. a=1/2.

dsolve(eval(ode,{a=1}));

let us consider the series

y = a_0 - sum(a_i,i=1..n)

here the terms a_i are defined recurrsively as follows

a_i = (a_(i-1) -4*ln(x)/ln(10) + 4/10 + 4*ln(a_(i-1)) )/ ( 1+ 4/(a_(i-1) * ln(10)))

how can i program this in maple. so i can...

Why doesn't 'ln(e)' simplify to '1' in Maple 12.02?

The problem arose from trying to simplify the solution of an ODE that contains a combination of exponentials, logarithms and powers. Maple refuses to simplify, whereas I can easily read the simplified answer. Maple seems to ignore the simplest of algebra concerning power, logarithms and exponentials.

what is the general sum formula for the following seires (we may notice from these few terms it has a pattern)

I so thankful for your help. Please help me expressing this expression compactly is there any pattern here?

I have a mapleTA question that blocks students from starting a mastery assignment until I tell them how to do the problem. But some students are gaming it but finding out the problem from an earlier section. I want to game back by making the answer change depending on the wall clock, that is the section they are in. To do so, I need a clock function to tell the time. Since I can call maple inside mapleTA, a maple solution will also work.

ElasticBeam.mw

Does anyone have any idea why I'm getting this plotting error? This should be working!!! Ahhh.

Hello,

I need some help with plotting the function with varying variables lambda and m ( varying from 0.1 to 1) which subject to the mutual constraint.

My function is:

> d:=(lambda, m) -> (0.1281008386e-3*m*lambda*(20.004936*lambda+21.004936*m+1)^2/(lambda+m)^2+0.1281008386e-3*(21.004936*lambda+21.004936*m-m*(1-lambda))^2/(lambda+m)^2)-(0.2690752038e-2+0.2690752038e-2*m*lambda);

The constraint is:

Is there a way to have interval arithmetics for complex inputs
(maple seems to name that range arithmetics in the help)?

x:=Pi/5;
shake(x);

              INTERVAL(0.628318529462 .. 0.628318531974)

shake(x*I);

    Error, (in evalr/shake) not a real number

Hi all,

I have gone through the help pages in my Maple 12, but cudn't really solve for roots of more than 1 pollynomial equations in parametric form. (I saw roots,parametric pollynomial function etc. & many more but cudnt do it with success).

my equations are coupled like, for example,

axy + bx - c*x*(x+y) = 0,

dxy + ey - f*y*(x+y) = 0,

etc. and I need the roots (x,y...

Dear guides

I have an equation: "eq := Pi-arcsin(sqrt(2E)/sqrt(1+E))-arcsin(sqrt(2E)/sqrt(2+E))-2*sqrt(2E)=0;

when I use "solve" maple say: "Warning, solutions may have been lost" and get a bad result. How can I ontain a better answer?

thanks.

The goal of computing only a select number of eigenvectors of a real symmetric floating-point Matrix comes up now and then. For very large Matrices the memory requirements can be more restrictive than the timing.

The attached worksheet and code computes this, more quickly and with significantly less memory allocation than does the usual task of computing all eigenvectors. By using the supplied Matrix itself as a partial "workspace" the amount of additional workspace and memory allocation for the task is negligible.

For example, having created the very large Matrix in the first place,  essentially no significant further memory allocation is required to compute the largest eigenvalue and its associated eigenvector.

A little about this routine `SelectedEigenvectors` follows.

It only works in hardware double precision. It expects a float[8] datatype Matrix (because you are serious about using minimal memory!). It uses the CLAPACK function dsyevx, using the "wrapperless" version of Maple's external-calling mechanism. It seems to work fine in the systems I've tried so far: Maple 13 and 14 on both 32bit and 64bit Linux and Windows.

Whether it computes and returns the selected eigenvectors (alongside the selected eigenvalues, which are always returned) is controlled by the 'vectors=truefalse' optional argument. By default it uses the Matrix argument as partial workspace and so destroys the original data; but this can be overridden with the 'preserve=true' optional argument. The requested accuracy can be relaxed with the 'epsilon=float' optional argument, which might sometimes speed it up.

The input Matrix is presumed to be symmetric. By default it uses the data in the lower triangle, but this can be changed to be the upper triangle with the 'uplo' optional argument.

The choice of eigenvalues is controlled by the two integer arguments `il` and `iu`. If il=iu=n then only the nth largest eigenvalue is computed. If il=1 and iu=4 then the four smallest eigenvalues are computed.

It returns three things: a Vector of dimension n whose first m entries are the selected eigenvalues, a nxm Matrix whose columns are the m associated eigenvectors, and a Vector of dimension n whose entries indicate whether corresponding eigenvectors failed to converge.

I didn't enable float arguments such as `vl` and `vu` which in principle could allow one to supply a floating-point range in which to find eigenvalues.

I didn't make an optimization of having it do an initial "dummy" external call in which no calculation would be done, but which would instead query for and subsequently utilize the optimal-performance additonal float workspace size.

For reasons mysterious to me, on Windows the 64bit version runs almost exactly half as fast as the 32bit version.

Usually, the workspace for eigen-solving is implemented to be at least O(n^2) for an nxn Matrix. But this routine does only O(m+n) extra workspace allocation to compute the m eigenvectors. And that is linear. Which is the Big Deal.

A 5000x5000 datatype=float[8] (ie. hardware double precision) Matrix takes 200MB of memory. With the preserve=true option, this routine can compute just the largest eigenvector with only about 200MB of additional allocation. And if the original Matrix is no longer required then with the preserve=false option this routine can do that task with less than 1MB further allocation. In comparison, the regular LinearAlgebra:-Eigenvectors command would require about 600MB of additional memory allocation while computing all eigenvectors.

At size 5000x5000 this routine is only about four times faster than LinearAlgebra:-Eigenvectors. I suspect that is because it still has to compute in full the reduction to tridiagonal form.

Download dsyevx.mw