Hi, 

I'm trying to set up the dirac algebra using the Physics package in maple 18. There are dirac gamma matricies (Dgamma) already specified, but I can't seem to manipulate their commutation relations. 

So I've tried building my own: 

restart;
with(Physics);

Setup(noncommutativeprefix = {gamma});

g[1] := gamma[1]; g[2] := gamma[2]; g[3] := gamma[3]; g[4] := gamma[0];

InverseMetric := rhs(g_[`~mu`, `~nu`, matrix])

Algebra :=  (a, b) -> %AntiCommutator(g[a], g[b]) = 2*InverseMetric[a, b];

Rules := Matrix(4, 4, Algebra);

Setup(algebrarules = Rules);

Error, (in Physics:-Setup) unable to set AntiCommutator(gamma[0], gamma[1]) = 0 because, taking into account {AntiCommutator(gamma[0], gamma[0]) = 2}, we would have gamma[0] and gamma[0] anticommutative and AntiCommutator(gamma[0], gamma[0]) <> 0

It seems like Maple can't handle the Dirac Algebra? Or have I done something obviously wrong?

Any help is appreciated. 

Thanks

Find the smallest number N such that

π^2/6 − Σ (from n=1 to N) 1/n² < 0.001

using Maple commands.

Successively transform the expression x + y + z into x ∗ y ∗ z and [x, y, z], using Maple commands.

Hello,

I am going to model a statically balanced mechanism by using zero free length spring in MapleSim. In path Multibody> Force and Moments>, there is a Translational Spring and Damper which I can use, and I just enter zero length for spring unstreched length. In physical model, zero free length spring never get zero length because of the existing coils, but in MapleSim it reaches zero during simulation. Does anyone here know how can I model zero free length springs?

Thanks

How to find the coefficients of the terms , ,...

how to save the session of maple worksheet so that no need to calculate again after restart computer

after calculated 70 million result for a long time, how to save the result and session so that no need to recalculate again?

Hi,

I use Maple on my laptop computer (windows 8.1) with a small screen. My problem is that the text in the menus, palettes, plots etc. are so small so I can´t read it. I have set Large toolbar icons under the menu tools-options-interface and I use the zoom buttom on the toolbar but this only magnify the text in a document. I also use the magnifying glass in windows, but I find it very inflexible. Is there any better solution to this problem?

Kind Regards

Leif Jonson

Hello

I have done this question again but i didnt get an answer.So does anyone know how i can substitute the following variables to an equation?I tried to create a list but i dont know how the variables can be substituted in the right place.data.docxdata.docxalfaeq.mw

Thanks

Contour lines must be ordinary circles. In fact, we get:

plots[contourplot](1/(x^2+y^2), x=-1..1, y=-1..1);

If we use the additional options, the result is even worse:

plots[contourplot](1/(x^2+y^2), x=-1..1,y=-1..1, numpoints=10000);

Hello i want to sort according to u derivatives (k) system.  And finding determining equations system and solving this system. Thank you very much.  

Download det.eq..mw

I am not sure I know what is going on. I expected to obtain the controllable canonical form in this example, but I am not. May be I am not using it correctly. 

Given A,B, I wanted to transfer the state space to controllable canonical form

http://www.maplesoft.com/support/help/Maple/view.aspx?path=DynamicSystems%2FSSTransformation

In this form, the A matrix will have 1 on the super diagonal, and the last row will have the coefficients of the charaterestic polynomial in reverse order with a minus sign. The B matrix will have all zeros, except for the last entry.  This is what the example on the above page actually shows. 

But when I tried it on my A,B, I do not get this form for the new B matrix. Here is a MWE

restart;
with(DynamicSystems):
A:=Matrix([[0,0,1,0],[0,0,0,1],[-2,-1,0,0],[1,-1,0,0]]);
B:=Matrix([[0],[0],[1],[0]]);
C:=Matrix([[0,0,0,0]]):
D0:=Matrix([[0]]):
sys:=StateSpace(A,B,C,D0): #just to see the polynomial
CharacteristicPolynomial(sys, s);

SSTransformation(A,B,C,D0,form=ControlCanon,output=['A','B']);

The above should be

I am sure I am doing something wrong, but what?

Maple 18.01, windows 7

Below is the function that I have.

f := (t-1)^(1/3)

p:=2;

b[n] := 2/p*(Int(f*sin(2*Pi*n*t/p), t = 0 .. p))

I also included a picture below to show what it is doing. Some help would be greatly appreciated. All I need to know is why maple doesn't want to evaluate bn?

If I have a function such as fk := exp(-x^2)*sin(Pi*x)  and I want to:

Create an animation to illustrate the differences between line thicknesses and then after the thicknesses, illustrate line styles. I also need to add the line thickness and style as text to the animations.

How would I go about doing this?

Hello, sorry for may asking the most simple question ever, but we cant find a solution.

If we do have a simple eqation like this:

U=I*R

how to solve vor R=? using maple?

Okay, we do feel really stupid for asking that and okay we can do that by hand but how to use maple for this, especially when it becomes more complex? Are there any examples of this kind?

Tahnks a lot Eve and Andrew

So if I have a procedure like 

with(GraphTheory)

tneighbors := proc (G::Graph)

local numvertices::integer, i::integer, currentvertex;

numvertices := nops(Vertices(G));

for i to numvertices do

currentvertex := Vertices(G)[i];

if nops(Neighbors(G, currentvertex)) = 2 then print(currentvertex)

end if;

end do;

end proc;

How do I make it so the output gets returned as a set?

For example, if I do  twoneighbors(G); and get
1

7

4

How do I make it so the output is listed as a set like {1,4,7}? Thanks.