I am using MAPLE for quantum computation.
It appears to me that the Physics[Expand] and Physics[Symplify] functions do not operate as I expected on
Kets with multiple quantum numbers. The functions no not consider the non commutative nature of tensor
products of Kets. (neither of bras incidently). It would be very useful if thos could be adjusted in some way.

Thank you for your help

LL

Download Expand_Simplify-of-Kets.mw

The initial conditions are U[13.75]=0.7787 and U'[13.75]=0.344037

Hello together,

I´m a new member and I used in Maple 2017 the function "is prime" for the largest known Mersenne prime (2^77,232,917 − 1) with the command "isprime (2^77,232,917 − 1)" to test how much time the programm needs that it returns "true".

1.) Has anybody experiences with the function "is prime" in relation with such a large Mersenne prime ( 2^77,232,917 − 1) or similar Mersenne primes like ( 2^57,885,161 − 1), ( 2^74,207,281 − 1), etc... ? How many days it takes to get the confirmation in my case ?

My processor: Intel i5-4590 CPU @ 3.30 GHz

System type: 64 Bit  

I started the function isprime(2^77,232,917 − 1) before 5 days. The programm is still evaluating, next to "evaluating" the point changes his colors (black and white) continuous, so I think the programm is still working, but the memory stopped at 2440.33 M and the time stopped at 12326.65 s

2.) Has anybody an idea why Maple has stopped to count the time and the memory in this case ? 

I hope that somebody can answer my 2 questions....thx... 

Good morning sirs,

I have been trying to compute this problem to its given me the error message "Error, invalid input: rhs received {f[0](eta) = -(1/50)*eta^2+(1/5)*eta+1/5, theta[0](eta) = -eta+5}, which is not valid for its 1st argument, expr".

The problem is attach below.

Please, anyone with useful informations should help.

Thanks

HPM_4.mw

Download Difficult-to-explain.mw

I am interested and intrigued by optical illusions, such as that entitled 'Drawing A Hole in Line Paper

 - 3D Trick Art' by artist Jonathan Stephen Harris. Check out this YouTube video  to appreciate the following. Not being good at art, I tried to replicate this illusion using Maple..  The program below is merely a start towards this goal.  I have a procedure which draws the outline of the letter A.  Different letters have different peculiarities: I chose the letter A since it is relatively easy in that it does not have curves.  It does have a triangular region, so I created a separate procedure for that.  The whole letter A can be rotated through an angle phi.   Also I have drawn a grid of parallelograms for use as the lines drawn on the paper - but this I think has bugs. .

     The method I have used is rather long and messy, just using coordinate geometry, for the vertex points of the A, then rotating them through an angle phi.   Also in the program is a grid of parallelograms, for use in drawing the lines across the page.  That's as far as I've got.  Major problems I foresee is seeing/calculating where these lines meet the outline of the letter A, where vertical lines would be drawn.  Also, the artist draws in shading - how can that be done in Maple?  My attempt at the task has brought up some issues that there must be a simpler, better method of doing this.  eg first put the coords in a vector, and use matrix multiplication to calculate the new coords.  My method is long and error prone. 

   I'd appreciate some feedback about any attempts similar to this.  In my program I tried to fill the color of the polygon which was the boundary of the A, but it filled in the base trapezoidal region as well.  On top of that I failed to color the small triangular region a different color.  As always, any help or suggestions would be gratefully received. 

.   

restart:

printlevel:=0:

with(plots):

with(plottools):

# nrows= Number of rows ncols is one more than nrows

nrows:=8:ncols:=9:

x0:=0:y0:=0:theta:=Pi/12:psi:=Pi/3:
#Width, w, and length, l, of sides of the parallelograms/polygons c[i,j]

w:=2:l:=3:

for i from 0 to nrows do

  for j from 0 to ncols do

#c[i,j] := rectangle([i+1,j], [i,j+1], color=red):

c[i,j] := polygon([[x0,y0],[x0+i*w*cos(psi),y0+i*w*sin(psi)],[x0+i*w*cos(psi)+j*l*cos(theta),y0+i*w*sin(psi)+j*l*sin(theta)], [x0+j*l*cos(theta),y0+j*l*sin(theta)]], filled=true,color=red):

  end do:

end do:

plots[display](seq(seq(c[i,j], j=0..ncols),i=0..nrows), scaling=constrained);

#t1:=textplot([2,3,`David`]):

#for i from 0 to nrows do

#  for j from 0 to ncols do

   #c[i,j] := rectangle([j+1,i], [j,i+1], color=white):

#  end do;

#end do;

#pl1:=plots[display](seq(seq(c[j,i], i=0..ncols),j=0..nrows), scaling=constrained):

pl1:=plots[display](seq(seq(c[i,j], i=0..ncols),j=0..nrows), scaling=constrained):

plots[display](pl1);

# To draw the letter A

# Using proc poly_out:  also to put in initial coords (x0,y0) of lower left foot of A, then rotate the letter A through an angle of phi

#thet:=Pi/3:phi:=Pi/2:

# thet is angle the left "diagonal" makes with the horizontal.

# l is the length of the diagonals - ie the left and right hand sloping sides of the letter A

# w is width of the feet of the letter A.  (Both equal width)

# topl is the "top length" of the horizontal top part of the letter A.

# (x0, y0) are the coords of the bottom left point of the letter A.

# phi is the anti-clockwise angle of rotation about (0,0)

poly_out:=proc(thet,l,w, topl, x0, y0, phi)

local outA, outAr,trig, trigr,D,E,F,G, corr, eps:

#l:=10:thet:=Pi/3:w:=4:

#topl:=3:

eps:=3*w/4:corr:=.3*w:

trig:=polygon([[x0+(2*l*cos(thet)+topl)/2, y0+l*sin(thet)-eps],[x0+w+(l/3+corr*w)*cos(thet), y0+(l/3+corr*w)*sin(thet)],[x0+2*l*cos(thet)+topl-w-(l/3+corr*w)*cos(thet), y0+(l/3+corr*w)*sin(thet)]]):

trigr:=polygon([[(x0+(2*l*cos(thet)+topl)/2)*cos(phi)-(y0+l*sin(thet)-eps)*sin(phi),(x0+(2*l*cos(#thet)+topl)/2)*sin(phi)+(y0+l*sin(thet)-eps)*cos(phi)],[(x0+w+(l/3+corr*w)*cos(thet))*cos(phi)-(y#0+(l/3+corr*w)*sin(thet))*sin(phi),(x0+w+(l/3+corr*w)*cos(thet))*sin(phi)+(y0+(l/3+corr*w)*sin(th#et))*cos(phi)],[(x0+2*l*cos(thet)+topl-w-(l/3+corr*w)*cos(thet))*cos(phi)-(y0+(l/3+corr*w)*sin(th#et))*sin(phi),(x0+2*l*cos(thet)+topl-w-(l/3+corr*w)*cos(thet))*sin(phi)+( #y0+(l/3+corr*w)*sin(thet))*cos(phi)]], color=white,filled=true):

outA:=polygon([[x0,y0],[x0+l*cos(thet),y0+l*sin(thet)],[x0+l*cos(thet),y0+l*sin(thet)],[x0+l*cos(thet)+topl,y0+l*sin(thet)],[x0+2*l*cos(thet)+topl,y0], [x0+2*l*cos(thet)+topl-w,y0],[x0+5*l*cos(thet)/3+topl-w,y0+l*sin(thet)/3],[x0+l*cos(thet)/3+w,y0+l*sin(thet)/3], [x0+w,y0]]):

outAr:=polygon([[x0*cos(phi)-y0*sin(phi),x0*sin(phi)+y0*cos(phi)],[(x0+l*cos(thet))*cos(phi)-(y0+l*sin(thet))*sin(phi),(x0+l*cos(thet))*sin(phi)+(y0+l*sin(thet))*cos(phi)],[(x0+l*cos(thet))*cos(phi)-(y0+l*sin(thet))*sin(phi),(x0+l*cos(thet))*sin(phi)+(y0+l*sin(thet))*cos(phi)],[(x0+l*cos(thet)+topl)*cos(phi)-(y0+l*sin(thet))*sin(phi),(x0+l*cos(thet)+topl)*sin(phi)+(y0+l*sin(thet))*cos(phi)],[(x0+2*l*cos(thet)+topl)*cos(phi)-y0*sin(phi),(x0+2*l*cos(thet)+topl)*sin(phi)+y0*cos(phi)], [(x0+2*l*cos(thet)+topl-w)*cos(phi)-y0*sin(phi),(x0+2*l*cos(thet)+topl-w)*sin(phi)+y0*cos(phi)],

[(x0+5*l*cos(thet)/3+topl-w)*cos(phi)-(y0+l*sin(thet)/3)*sin(phi),(x0+5*l*cos(thet)/3+topl-w)*sin(phi)+(y0+l*sin(thet)/3)*cos(phi)],

[(x0+l*cos(thet)/3+w)*cos(phi)-(y0+l*sin(thet)/3)*sin(phi),(x0+l*cos(thet)/3+w)*sin(phi)+(y0+l*sin(thet)/3)*cos(phi)], [(x0+w)*cos(phi)-y0*sin(phi),(x0+w)*sin(phi)+y0*cos(phi)]], color=grey,filled=true):

plots[display](outAr,trigr, axes=none, scaling=constrained);  # view=[-l..l,-l..l]);

#outA & trig removed from display – these give the `upright` #letter A

end proc:

plot1:=poly_out(Pi/3,15,3,3, 0,0, Pi/15):

plots[display](plot1, pl1);

I'm just unable to see how this could be correct, if someone can take a look it would be much appreciated.

Download mapleHELP.mw

Dear users! I need the help in attached file. Please see and fix it. I am waiting your positive answer.

File_to_help.mw

Dear All

I need to minimize a quadratic function with constraints, but i get

Error, (in Optimization:-QPSolve) unexpected parameters

restart;
with(Optimization);
with(plots);

obj := (c0-x)^2+(c1-(2/3)*x-(1/3)*y)^2+(c2-(1/3)*x-(2/3)*y)^2+(c3-y)^2;
cnsts := 32*x+19*y = 15*c0+18*c1+15*c2+3*c3;
          Minimize(obj, cnsts);

I would like to obtain the optimal parameter ( x,Y) that satisfies the minimization problem with constraints , and the results must depend on c0,c1,c2 and c3. 

Many thanks for your help

hi . i want to solve couple differential  equation with ode analyzer  but i get error

i want to plot y[1] and y[2]. 
for example for this two equation 
diff(y[1](x), x$2)+3*(diff(y[1](x), x$1))+2*y[1](x)+2*y[2](x) = 3, diff(y[2](x), x$2)+diff(y[2](x), x$1)+2*y[2](x)+2*y[1](x) = 8;

How to solve matrix in maple.

Example: Set A := Matrix(3, 3, [[2, -3, 1], [-3, 5, 0], [1, 0, 5]]). Find matrix L such that A = L^T*L.

I can solve by hand but how to use Maple?

in plot line style should to like  and .so on....

line & diamond   ,   line  and circle   ,  line and cross,    line and cross,

 same should be displayed in legends

I have a list L:=[[0,0,0], [1,0,0], [1,1,0], [0,1,0], [0,0,1], [1,0,1], [1,1,1], [0,1,1]] (8 vertices of a cube). How can I select four vertices of the list to make a regular tetrahedron?

Hello,

I am experiencing some trouble with the function Minimize of the Optimization package. I have some trouble understanding what Maple is doing with nice functions already.

f := (x-2)^2+1;
                                 2    
                          (x - 2)  + 1
Minimize(f, {x >= -2, x <= 2});
             [1., [x = HFloat(1.9999999999999998)]]
Minimize(f(x), {x >= -2, x <= 2});
                    [1., [x = HFloat(2.0)]]
g := proc (x) options operator, arrow; (x-2)^2+1 end proc;
            2    
x -> (x - 2)  + 1
Minimize(g, {x >= -2, x <= 2});
Error, (in Optimization:-NLPSolve) constraints must be specified as a set or list of  procedures
Minimize(g(x), {x >= -2, x <= 2});
             [1., [x = HFloat(1.9999999999999998)]]

Hello there!

I have just an elemental question. Given an equation as input, how could I write a procedure to know the variables involved in that equation? 

For example, my input is 2x+3y=0, and I would like to have the output  {x,y}.

Thanks in avance for the help!