Note added: Issue resolved, see my comment below.

I have a sum of several thousands addends each of which is the product of a c-number times a product of 6 Grassmann-odd degrees of freedom, the latter of which does each belong to a set of 24 Grassmannians. The specific numbers are not that important, though.

This sum should equal zero. So I would like to add all c-numbers multiplying the same product of six Grassmannians, taking proper care of anticommutativity, of course. The sum would then be zero if all these sums of c-numbers are zero. Unfortunately, using Physics:-Coefficients is far too slow; actually, it has never succeeded in even completing the calculation.

Therefore, I have tried to loop through all the addends, splitting each one of them using selectremove(), and then adding the c-numbers in an Array (properly indexed), or in a table (associatively indexed, of course). Consistently, by converting the Array and table to two sets, the two methods result in the same set of equations. But solving these equations yields a result that is not stable: it varies from session to session.

I am baffled. Can anyone give me a hint to a safe and reasonably quick method for extracting these c-number-valued equations?, for there has to be something wrong with what I do.

i need to plot tow functions but the legend not in top, bottom, left or right 

Given the sequence defined by the recursive relation a[n+1] = r*a[n](1-a[n])
You need to use the procedure iterate.
Throughout this problem you should choose initial values in the interval 0<a0<1.
(a) Let r=3/2. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit depend upon your choice of initial value? Plot the terms you have calculated
(b) Let r=2.8. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit
depend upon your choice of initial value? Plot the terms you have calculated How does this sequence differ from that in part (a).
(c) Let r=3.2. Calculate a moderate number of terms in the sequence (between 10 and 20). Show that the sequence does not appear to converging. Plot the terms you have calculated and describe how the sequence behaves in this case.
(d) Consider intermediate values between 2.8 and 3.2 to determine more precisely where the transition in behaviour takes place. Provide a few plots (no more than 4) showing the values you have investigated.
(e) Consider the values of r in the range 3.43<r<3.46. Determine as accurately as you can the value of r for which the period of oscillation doubles.
(f) As r increase further period doubling occurs. Try to find the when the sequence appears to oscillate between 8 values.
(g) Let r =3.65 and calculate a considerable number of terms (at least a few hundred) and plot your values.
(h) For r=3.65 choose a0=0.3 and then a0=0.301. Find and plot some terms in the sequence for each initial value. Determine how long the terms in the two sequences remain close together and when they begin to depart significantly from each other.

Hi. I have a commutation table, that I've got as output of Greg Reid algorithm.

Can I use to define Lie algebra for components DifferentialGeometry & Lie Algebras?

Worksheet is attached...

FCT.mw

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?

I have attached the whole worksheet, I hope someone can help! Thanks Detectors_Tutorial.mwDetectors_Tutorial.mw

Even if there was some method of exporting the graph to another medium that would also be very helpful

This is the code  hw2_final.mw

Let me explain it.

I am sure that the mistakes must be about the expresstion of the I1(t) and I2(t). Actually if you delete I1(t) and I2(t) , the whole code works and get the picture at the bottom. 

What I want is to put the expresstion of the I1(t) and I2(t) into 'sol:=dsolve...' and 'plots...' to get the picture of I1(t) and I2(t) with respect to t. Before the t* which subject to Phi(t*)=0 (The blue line in the picture at the bottom is Phi) I want I1(t) and after t* I want I2(t).

I1(t) = (int(sqrt(2*(H(t)+omega*cos(q(t)))), q(t) = q(t)-2*Pi .. q(t), numeric))/Pi.    what I want of this experesstion is to get  'int(sqrt(2*(H(t)+omega*cos(q(t)))' from  'q(t)-2*Pi' to 'q(t)' by numeric method.This q(t) is the solution of the ODE sys.

For example(the number I used is not true,just for example) , at the point t=20, q(t)=30-2*Pi.

so I1(t)= (int(sqrt(2*(H(t)+omega*cos(x))), x = 30-2*Pi .. 30, numeric))/Pi.The I2(t) I want is similar to I1(t).

How can I solve it?

Hello, 

I am pretty new to Maple and Im trying to find the parameters of an equation using some kind of fit routine but I can only find such a routine to fit an expression to data and not vice versa.

my equation is as follows

epsilon(E):= a0 + a1*ln(E) + a2*ln(E)^2 + a3*ln(E)^3

I do have data for E which I imagine I need. 

E:= {121, 244, 344, 411, 444, 778, 867, 964, 1085, 1112, 1212, 1299, 1408} all in keV :)

Any suggestions/help would be much appreciated, although I am new to the program so go easy on me :) 

Cheers

Integral.mw

Hi all

Please help me to solve this Integral...

There have come unwanted lines and marks . I donot know how to remove them. Using doc.block, remove block seems to be little tough to incorporate! Please enlighten me. Modified doc. is most welcome. Thanks. Ramakrishnan V 

Download GausianFinal15Nov2015.mwGausianFinal15Nov2015.mw

There are the complexes (C), quaternions (H or Q), octionions (O), sedenions (S) and the pathions (P). I have found the multiplication tables of them, although according to signs (+ or -) there differents at pathions. The important question is that How can I multiply two bases, i_n and i_m of higher dimensions, like in the routions or in the voudions?

Should I xor the indexes of the bases? Like this way: i_1 * i_2 = i_(1^2) = i_3

What is about the signs?

hello dear
how can i solve this three couplde pde?please help me

I want to solve an ODE from Game Theory, the Cournot competition.

It says

p(q1+r2(q1))+p'(q1+r2(q1))*r2(q1)-c2'(r2(q1))=0

 where, I think,

' means diff(,q1),

c2(q2)=c*q2 for a fixed c in [0,1]

and

p(q)=max(0,1-q).

So c2,p and r2 are functions.r2 goes from [0,inf) to [0,inf).

I look for r2, which should be r2(q1)=(1-q1-c)/2 when correctly solved.

However, the command dsolve says Error in dsolve (divison by 0).

 What is wrong? How do I obtain the solution for r2 in Maple?

How many numbers exist which are: less than 1000, multiple of 18 and also multiple of 10? Use Maple to find the answer 

Hello,

I tried to solve below equation, but it gives me zeros result. Please help me to find their inverse laplace. 

 It will be clearer if was pasted on Maple:

restart

Ps := [P[0], P[1], P[2], P[3], P[4]]:

eqs := [P[0](s) =~ (P[1](s)*mu[1]+P[2](s)*mu[2]+P[3](s)*mu[3]+P[4](s)*mu[4])/(s+lambda[1]+lambda[2]+lambda[3]+lambda[4]), P[1](s) = lambda[1]*P[0](s)/(s+mu[1]), P[2](s) = lambda[2]*P[0](s)/(s+mu[2]), P[3](s) = lambda[3]*P[0](s)/(s+mu[3]), P[4](s) = lambda[4]*P[0](s)/(s+mu[4])];

Ls := solve(eqs, Ps(s))[];

P(t)=~inttrans[invlaplace]~(rhs~(Ls), s, t);

Thank you

I have several *.mw files which I have built seperately. Now, for readability reasons, I would like to create a Master document, which will ues the first files, and then generate a PDF document.
I tried the read command, but it does not recognize the files.

Any help is appreciated.