AOA... There are three question

1. I want to convert exp(Iota*theta) into ternometric function i.e., 

exp(Iota*theta) = cos(theta)+Iota*sin(theta)

Is there any comand pl help...

2. Also i want to rationalize the complex number...

3. I want to seprate real and imaginary parts of a comaplex numbers

I try to sort a polynomial using the graded reverse lexicograpic order. According to the Documentation this is achieved via tdeg.
So here is what i tried:

with(Groebner):
sort( x+y+z, order = tdeg(x,y,z));

or

sort(x+y+z, [x,y,z], tdeg);

In both cases maple returns "x+y+z" instead of the expected "z+y+x". What am i doing wrong?

Hello,

I have to simplify a piecewise function and Maple gets a more complicated solution than needed.




I don't know how to handle this kind of problems with Maple?
I don't understand why Maple doesn't see this?
Am I doing something wrong?

Thanks in advance for your help / advice.


# the code of my example
restart:
Mf(x):=piecewise(x<=L/2,1/2*x*F,x>1/2*L,1/2*x*F-F*(x-1/2*L));
# Make a dimensionless function:
# -    Mf(x):= Mf(xi)*F*L
# -    variable ξ  ( xi:=x/L )
eq[1]:=Mf(xi)*F*L=Mf(x);
Mf(xi):=solve(eq[1],Mf(xi));
Mf(xi):=subs(x=xi*L,Mf(xi));
# F is the Force and L is the Length of the beam:
Mf(xi):=simplify(Mf(xi)) assuming F>0,L>0;
print("When I simplify this function by hand it will be");
Mf(xi):=piecewise(xi<=1/2,1/2*xi,xi>1/2,-1/2*xi+1/2);

Hello, everyone. I faced some promblems on maple. Hope you can help me.

I want to set up the commutation rules like Pauli sigma matrix in Physics Packages,

but failed to get the correct one.

Here is my maple code:

restart;
with(Physics);
Setup(mathematicalnotation = true);

Setup(quantumop = {Q}, algebrarule = {%AntiCommutator(Q[j], Q[k]) = 2*KroneckerDelta[j, k], %Commutator(Q[j], Q[k]) = 2*I*(Sum(LeviCivita[i, k, l]*Q[l], l = 1 .. 3))});

Commutator(Q[1], Q[2]);
                              
AntiCommutator(Q[1], Q[2]);
                              
AntiCommutator(Q[1], Q[1]);
                              
Commutator(Q[1], Q[1]);

In fact, Commutator(Q[1], Q[2]) give a incorrect result 0 while the correct answer is  2*I*Q[3].

Do anyone know how to resolve this?

Hello,

can I remove tick marks from axis in plots ?

Hi,

thank for any  help.

I need you help to make a movies ( animation)  according to the index k=1..10.  k is used in Vect_T[k], vector of time to make an animation of the true solution, and it's also used in Matrix(M[..,..,k]), h

True solution:ow can I make a seq of plots for different value of k. 

plot3d(True_solution(x,y,Vect_T[k]),x=0..a,y=0..b,axes=boxed,grid= [11,11]);

Analytic solution

plots:-matrixplot(Matrix(M[..,..,k]),axes=boxed,tickmarks= [[1="0", 2=" ", 3=" ", 4=" ", 5=".4", 6=" ", 7=" 0.6", 8=" ", 9=".8", 10=" ", 11="1"],[1="0", 2=" ", 3=" ", 4=" ", 5=".4", 6=" ", 7=" 0.6", 8=" ", 9=".8", 10=" ", 11="1"],default],labels=[x,y,u],axes = boxed);

I tried, this:

S:=seq(plots:-matrixplot(Matrix(M[..,..,i]),labels=[x,y,u],axes = boxed),i=1..10):
plots:-display(S);
 But I get only one figure.

 also,

I tried

BoxTrue:=Null;
for k from 1 to 10 do
BoxTrue:=BoxTrue, plot3d(True_solution(x,y,Vector_time[k]),x=0..a,y=0..b,axes=boxed,grid= [11,11]): od:
display(BoxTrue, insequence=true);
also, no display

Given the lagrangian from the pendulum physics example on the help page,

or text format.

L := -sin(omega*t-phi(t))*(diff(phi(t), t))*a*l*m*omega+(1/2)*(diff(phi(t), t))^2*l^2*m+cos(phi(t))*g*l*m

Hi Everybody.

Doing some calculation in quantum mechanics, I stuble on that integral:

I see immediately that the answer is 1/2.  But Maple 18 cannot give an answer other than a limit that he cannot evaluate.  I even try assumption that p and hbar are realcons.  I get infinity.

Any idea?

Thank you in advance for your help.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

Hello,

Can anyone explain to me why there is a difference between:

f := int((A*sin(omega*t+phi)-B*sin(omega*t))^2, t);

and the expanded version:

f := int((A*sin(omega*t)*cos(phi)+A*cos(omega*t)*sin(phi)-B*sin(omega*t))^2, t);

Thanks

Hi,

Is there any difference between

Matrix(4,5,(K,C)->K+C);

and

Array(1..4,1..5,(K,C)->K+C);

Say if I have a very 'complicated' procedure myfunc(K,C) that takes two options, but it runs all dependently.

Matrix(4,5,(K,C)->myfunc(K,C));

and

Array(1..4,1..5,(K,C)->myfunc(K,C));

Which one is more efficient? The final ouput of each run from myfunc is just a integer value.

The reason I am asking it that i think both runs on a 'single' thread (core) as CPU usage is always very low, around 15-20%.

If I look at the task manager, some cores (threads) arent doing anything.

Is there anyway to speed things up?

Thanks,

I am trying to get a Fourier transform of a Gaussian:

so I say

and get

The Fouriertransform of a Gaussian is well known and the result I expect is something like

exp(1/2*sigma^2*omega^2)

ignoring normalizations & other factors. I know that I can add functions to inttrans, but I kind-of expected inttrans[fourier] to know how to transform a Gaussian, it is a commonly used transformation. Even if I set phi0 to 0 it does not produce anything useful.

???

Mac Dude

hello, i'm trying to rewrite the RiemannSum function, it's for my homework, 

i wrote the entire fonction, but looks like my function is not evaluated, i don't understand why!!! if anyone could see what's wrong please let me know, thanks. here is  my piece of code, and i also attached the file

> with(Student[Calculus1]);

> with(plots);

> with(plottools);

> f(x):= (1)/(5.6 sqrt(2 Pi))*exp(-((x-169)^( 2 ))/(2(5.6)^(2)));

SR:=proc(a, b, f, n) description "Les sommes de Riemann";  local Delta, x,SR,xetoile, i;  Delta:= (b-a)/(n);  for i from 0 to n by 1 do  x(i):= a + i*Delta;  end do;  for i from 1 to n by 1 do  xetoile(i):= (1)/(4)(x(i) + 3 x(i-1));  end do;  SR := 0;  for i from 1 to n by 1 do  SR := SR + f(xetoile(i))*Delta;  end do;  return SR  end proc;

if i compute the procedure SR(0, 150, f(x), 10), it doesn't give a result

Hi.

Please, I need your help, to form and solve a linear system.

1) I want to solve the system (13) using "GenerateMatrix" , the system is   evaluated  at  point x[m]=-a+m*h, m=0..2N , In my code I get  equation (15), ( I think it's okay), I juste use subs, x=x[m].
My function, ph(x) approximate by phi[m]=phi(x[m]).  and I want to compute the vecteur [phi[1], phi[2],...,phi[2*N+1]] is my solution.

2) The nodes, are x[k]=-a+k*h, k=0..2*N;  I want to display these nodes in the x-axis, i write the code with the names of each point, but there is an error.

There is some a mistake in my code maple.

Here, my code.

Fred.mw

Hi

coeff dosn't work correctly in maple 17

Here is the code

restart;

E := m1*(diff(x1(t), t, t))+3*k*x1(t)-k*x3(t)-k*x2(t);
coeff(E, x1(t));
Error, unable

the coeff of diff(x1(t),t,t) and x2(t) and x3(t) can be determined and are true, but the coeff of x1(t) can not be determined??!!!!!!!!!!!!