Dear,i meet a problem in Physics package.i can't find why i got two different answers.

restart; with(Physics); Setup(mathematicalnotation = true);Setup(op = {C});

Cp := Dagger(C);

alg := %Commutator(Cp, C) = -1;

Setup(alg)

S := sqrt(N-(Dagger(C)+sqrt(N)*conjugate(eta))*(C+sqrt(N)*eta)):

eq1 := collect(expand(op(1, S)), N):

xi = op(1, eq1);

f := map(simplify, map(proc (x) options operator, arrow; x*sqrt(xi) end proc, convert...

Hi

I have problem with solve polynomial...maple can't manage with it.

My polynomial :

1-r+4*r^12*x^12+16*r^6*x^4-8*r^7*x^6+2*r^8*x^8+24*r^6*x^2-36*r^7*x^4+24*r^8*x^6-6*r^9*x^8-24*r^7*x^2+60*r^8*x^4-80*r^9*x^6+60*r^10*x^8-24*r^11*x^10-x+8*r^4*x^2-4*r^5*x^4-16*r^5*x^2+6*r^4+2*r^2-6*r^5-r^7+4*r^6+8*r^8*x^2-28*r^9*x^4+56*r^10*x^6-70*r^11*x^8+56*r^12*x^10-28*r^13*x^12+8*r^14*x^14-r^15*x^16-5*r^3

I'm getting for *good roots* and something "1/4(RootOf(.......

I've tried a few different ways to change the axis label size but am a bit lost. I'm new to the maple environment.

How to transform between permutation group and linear transformation matrix in maple

is there a command like gap

gap> i:=IsomorphismPermGroup(G); 

to find permutation group which has isomorphism

Hello,

I have a small Finite Element program in Maple. I want to accelerate the Matrix assembly. For Matlab this works like described here (http://milamin.sourceforge.net/technical-notes/sparse-matrix-assembly) under "Triplet Sparse Storage". There, two indexing Vectors are used to place the values into a sparse Matrix with the Matlab function sparse. If index pairs are duplicate, they are added.

Now I want to ask if this is also possible in Maple?...

hey friends i know how to write vector column matrix for nth order but what if i want it somehow like the one  attached in file matrix.mw

I have following lists

a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c

I want to write the bare rewrite rules for these so that if any list is given as input these rules produce the next one. It is obvious every third is same as first.

I want to get rid of the unnecessary output lines after the evaluation of the following ocde

restart; with(linalg);
a := 2; b := .29; d := 1.85;
for h from .5 by .1 to 1 do
eq1 := x*(-b*x^2-x+1);
eq2 := y*((a*x*x)/(b*y^2)-d-h*y);
S := solve({eq1, eq2}, {x, y});
SS := solve(subs(S[3], {omega^4+(h*y+x)*omega^2+h^2*x-y}), {omega});
tau := simplify(subs(S[3], subs(SS[3], (b^2*h*y+a*x)/omega)))
end do

Hi all,

Say i have two Vectors:

V1:=Vector(2,[a,b]);

V2:=Vector(4,[d,b,e,a]);

and how do I get

V3:=Vector([d,e]);

?? I currectly use a crude way using 'set' operations:

V3:=convert(V2,set) minus convert(V1,set);

V3:=convert(V3,Vector); # Does not work straightford, so

V3:=convert(V3,list):
V3:=convert(V3,Vector); # Which works,

Also, is there a strightforward way to which elements are the same and their indices?

Let us consider the group of elements {(a,b,c,d): a::integer,b::integer,c::integer,d::integer,a+b+c+d=0, a=c mod 12} with the component-wise addition. How to represent it in Maple? Is it possible? What I know is

restart;

with(GroupTheory):

CyclicGroup(infinity);     

GroupTheory:-CyclicGroup(infinity, form = "fpgroup")

G := DirectProduct(CyclicGroup(infinity), CyclicGroup(infinity), CyclicGroup(infinity), CyclicGroup(infinity));

when solving using output=listprocedure i face this message :

Sys:= {seq}(eq[k], k= 1..18):

fsolve(Sys, {seq}(a[k], k= 1..18),output=listprocedure):

Error, (in fsolve) output = listprocedure is an invalid option

what should i do? how can i list my output to use in latter commands? for example i want to use the output of the problem to plot an expression that the outputs are the coefficients of that expression?

so i got this
> low := 170; high := 195; j := 0; byy := 5; plott := vector((high-low)/byy+1);
with(plots);
for i from low by byy to high do
proc (X) options operator, arrow; X(t) end proc;
proc (Y) options operator, arrow; Y(t) end proc;
ode := diff(X(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvx = VXX[i-1], vvx))), diff(Y(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvy = VYY[i-1], vvy)));
ics := X(0) = Cordx[i-1], Y(0) = Cordy[i-1];
ds := dsolve([ode, ics...

I need three orthogonal unit vectors, k, i, and j. I give the first two vectors, k=[sqrt(2)-1,-1,1] and i=[1,sqrt(2)-1,0], then use crossproduct to produce the third vector j=crossprod(k,i) . But it fails to produce the correct vector since dotprod(j,k) is not 0. 

I have a two-variable matrixfunction A(x,y) with entries A_i,j(x,y).
I want to use mtaylor for each entries and the result to store in a matrix B.
Is it possible to do using the map command?
If not, what methods do exist?
(Instead of typing mtaylor for each entires.)

Sandor

restart:

dsys := {diff(x(t),t)=y(t),diff(y(t),t)=-x(t),x(0)=1,y(0)=0}:

dsn := dsolve(dsys,numeric):

dsn(1);          [t = 1., x(t) = HFloat(0.5403023317785669),  y(t) = HFloat(-0.8414711011553073)]

rhs(dsn(1)[2]);                   HFloat(0.5403023317785669)

After solving  this system I wish to obtain a numeric value to the integral:

int(x(t...