Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015

After manually working out answer for problem 4-4 in Mathews & Walker's Mathematical Methods of Physics , I tried to check my solution with maple2015. Briefly the problem involves inputs periodic with period T, being transformed into outputs, through a kernal G.  The net result is that all input frequencies omega periodic in T are multiplied by (omega_0/omega)^2, except for constant frequency which transforms to zero.  The problem asks to evaluate the kernal G.

Maple2015 correctly evaluated the integral for a constant input, a cosine input, and a sine input, but gave undefined when I tried an exponential(i*x) input which is just a linear combination of the two previous inputs.  I found this interesting because the integral is finite, well defined, and only has an absolute function (in the kernal), which may cause Maple problems, as it correctly evaluated integral when I split it into two regions.  Interestingly if instead of working with a period of T, I used 2*pi, and redfined my G function accordingly, Maple evaluated the exp input integral without any problems.  So the problem appears to be with the T variable, but I correctly used assumptions of T>0, and 0<t<T, so I am not sure why it would work correctly when I use T=2*pi, but failed when using a general period T.  Any help would be welcome.

 

 

restart

assume(T > 0)

assume(0 < t and t < T)

about(T)

Originally T, renamed T~:

  Involved in the following expressions with properties
    T-t assumed RealRange(Open(0),infinity)
  is assumed to be: real
  also used in the following assumed objects
  [T-t] assumed RealRange(Open(0),infinity)

 

about(t)

Originally t, renamed t~:

  Involved in the following expressions with properties
    T-t assumed RealRange(Open(0),infinity)
  is assumed to be: RealRange(Open(0),infinity)
  also used in the following assumed objects
  [T-t] assumed RealRange(Open(0),infinity)

 

assume(n::integer, n > 0)

about(n)

Originally n, renamed n~:

  is assumed to be: AndProp(integer,RealRange(1,infinity))

 

G := proc (x) options operator, arrow; (1/2)*omega0^2*T^2*((1/6)*Pi^2-(1/2)*Pi*abs(2*Pi*x/T)+Pi^2*x^2/T^2)/Pi^2 end proc

proc (x) options operator, arrow; (1/2)*omega0^2*T^2*((1/6)*Pi^2-(1/2)*Pi*abs(2*Pi*x/T)+Pi^2*x^2/T^2)/Pi^2 end proc

(1)

(int(G(t-tp), tp = 0 .. T))/T

0

(2)

(int(G(t-tp)*sin(2*Pi*n*tp/T), tp = 0 .. T))/T

(1/2)*T^2*omega0^2*cos(t*Pi*n/T)*sin(t*Pi*n/T)/(Pi^2*n^2)

(3)

(int(G(t-tp)*cos(2*Pi*n*tp/T), tp = 0 .. T))/T

(1/4)*T^2*omega0^2*(2*cos(t*Pi*n/T)^2-1)/(Pi^2*n^2)

(4)

(int(G(t-tp)*exp((I*2)*Pi*n*tp/T), tp = 0 .. T))/T

undefined/T

(5)

(int(G(t-tp)*(cos(2*Pi*n*tp/T)+I*sin(2*Pi*n*tp/T)), tp = 0 .. T))/T

undefined/T

(6)

simplify((int(G(t-tp)*exp((I*2)*Pi*n*tp/T), tp = 0 .. t))/T+(int(G(t-tp)*exp((I*2)*Pi*n*tp/T), tp = t .. T))/T)

(1/4)*omega0^2*exp((2*I)*t*Pi*n/T)*T^2/(Pi^2*n^2)

(7)

assume(0 < t and t < 2*Pi)

G2 := proc (x) options operator, arrow; 2*omega0^2*((1/6)*Pi^2-(1/2)*Pi*abs(x)+(1/4)*x^2) end proc

proc (x) options operator, arrow; 2*omega0^2*((1/6)*Pi^2-(1/2)*Pi*abs(x)+(1/4)*x^2) end proc

(8)

(int(G2(t-tp)*exp(I*n*tp), tp = 0 .. 2*Pi))/(2*Pi)

omega0^2*exp(I*n*t)/n^2

(9)

 

Download MathewsWalkerProblem4-4.mwMathewsWalkerProblem4-4.mw

 

 

I hoped that Maple would return the value of 1 in all commands (see below). However, introducing a scaling parameter, sigma, yields the unevaluated expression. Why? I still think it should evaluate to the value of 1.

 

kind regards,

Harry (not a mathematician, but a psychologist)

 

 

 

integral.mw

How to increase the size of toolbar icons in Maple 2015 32-bit Classic on Windows?

Hello everyone! I got some trouble in process a list. Hope you can help:

Assume i got a list like this:

 

{{k = k, l = RootOf(_Z^2+_Z*k+k^2-1), o = -k-RootOf(_Z^2+_Z*k+k^2-1)}, {k = k, l = RootOf(_Z^2+_Z*k+k^2+1), o = -k-RootOf(_Z^2+_Z*k+k^2+1)}, {k = 0, l = 1, o = -1}, {k = 0, l = -1, o = 1}, {k = 1, l = 0, o = -1}, {k = 1, l = -1, o = 0}, {k = -1, l = 0, o = 1}, {k = -1, l = 1, o = 0}, {k = RootOf(_Z^2+1), l = 0, o = -RootOf(_Z^2+1)}, {k = RootOf(_Z^2+1), l = -RootOf(_Z^2+1), o = 0}}

 

Now all i want is remove Complex and RootOf from this list, how can i do that?

Thank for your reading adn your help!

I have the following expression

f=u/(sqrt(u*(u-1)))

and I want to simplify it. Eventhough that I tell Maple that u is real and greater than 1 but it does not simplify the expression. What is wrong? Please see the attached file.

Radical.mw

How do you put labels on individual columns in ColumnGraph?

Hello! Hope every is fine. I want to expand the following expression

exp(2*c*t+2*d*n-d)*alpha*c*a[0]*b[1]^2-exp(2*c*t+2*d*n-d)*alpha*c*a[1]*b[0]*b[1]-exp(2*c*t+2*d*n-d)*alpha*a[0]*a[1]*b[1]+exp(2*c*t+2*d*n-d)*alpha*a[1]^2*b[0]+exp(2*c*t+2*d*n)*alpha*a[0]*a[1]*b[1]-exp(2*c*t+2*d*n)*alpha*a[1]^2*b[0]+exp(c*t+d*n)*alpha*c*a[0]*b[0]*b[1]-exp(c*t+d*n)*alpha*c*a[1]*b[0]^2-exp(2*c*t+2*d*n-d)*a[0]*b[1]^2+exp(2*c*t+2*d*n-d)*a[1]*b[0]*b[1]-exp(c*t+d*n-d)*alpha*a[0]^2*b[1]+exp(c*t+d*n-d)*alpha*a[0]*a[1]*b[0]+exp(2*c*t+2*d*n)*a[0]*b[1]^2-exp(2*c*t+2*d*n)*a[1]*b[0]*b[1]+exp(c*t+d*n)*alpha*a[0]^2*b[1]-exp(c*t+d*n)*alpha*a[0]*a[1]*b[0]-exp(c*t+d*n-d)*a[0]*b[0]*b[1]+exp(c*t+d*n-d)*a[1]*b[0]^2+exp(c*t+d*n)*a[0]*b[0]*b[1]-exp(c*t+d*n)*a[1]*b[0]^2

 

like this 

exp(2*c*t+2*d*n)*exp(-d)*alpha*c*a[0]*b[1]^2-exp(2*c*t+2*d*n)*exp(-d)*alpha*c*a[1]*b[0]*b[1]-exp(2*c*t+2*d*n)*exp(-d)*alpha*a[0]*a[1]*b[1]+exp(2*c*t+2*d*n)*exp(-d)*alpha*a[1]^2*b[0]+exp(2*c*t+2*d*n)*alpha*a[0]*a[1]*b[1]-exp(2*c*t+2*d*n)*alpha*a[1]^2*b[0]+exp(c*t+d*n)*alpha*c*a[0]*b[0]*b[1]-exp(c*t+d*n)*alpha*c*a[1]*b[0]^2-exp(2*c*t+2*d*n)*exp(-d)*a[0]*b[1]^2+exp(2*c*t+2*d*n)*exp(-d)*a[1]*b[0]*b[1]-exp(c*t+d*n)*exp(-d)*alpha*a[0]^2*b[1]+exp(c*t+d*n)*exp(-d)*alpha*a[0]*a[1]*b[0]+exp(2*c*t+2*d*n)*a[0]*b[1]^2-exp(2*c*t+2*d*n)*a[1]*b[0]*b[1]+exp(c*t+d*n)*alpha*a[0]^2*b[1]-exp(c*t+d*n)*alpha*a[0]*a[1]*b[0]-exp(c*t+d*n)*exp(-d)*a[0]*b[0]*b[1]+exp(c*t+d*n)*exp(-d)*a[1]*b[0]^2+exp(c*t+d*n)*a[0]*b[0]*b[1]-exp(c*t+d*n)*a[1]*b[0]^2

i.e., expand exp(2*c*t+2*d*n-d) into exp(2*c*t+2*d*n)*exp(-d) 

waiting your kind response 

hi.please help me for solve this equations

thanks...

Tur.mw

Non dimensionalisation is a vary common task, and I was suprised that I couldn't find a maple tool to automate it . Has anyone developed their own package for it?

I want to automatically do it to the system equations for some Dynamical systems to make some of the other processing I do with them easier.

I was hoping to start with somehting in the form of 

Diff(x[1],t)=f[1](p[1]....p[n],x[1]...x[m])

...

Diff(x[m],t)=f[m](p[1]....p[n],x[1]...x[m])

where each f[i] is some kind of quotient of multivariate polynomials in the variables and parameters:
and end up with something like

Diff(y[1],s)=f[1](q[1]....q[p],y[1]...y[m])

...

Diff(y[m],s)=f[m](q[1]....q[p],y[1]...y[m])

where p<n

I want to cancel some expressions in numerator and denominator of a quotient. But Maple deos not cancel it!

Please see the attached file.


Simplifying_Radicals.mw

in LinearAlgebra Eigenvectors calculation.

Maple 2015 Error

 

 

So the above output startled me.  I have used the Maple Linear Algebra Eigenvalues, Eigenvectors commands many times with no problem.  Can any one explain to me what is going on.  The program correctly calculates the eigenvalues for the matrix which are all distinct for a real symmetric matrix, and thus should have three distinct non-zero eigenvectors, yet the eigenvectore command returns only zeros for the eigenvectors.  I calculated an eigenvector by hand corresponding to the eigenvalue of 1 and obtained (1, -sqrt(2)/sqrt(3), -1/sqrt(3).

 

So this is either a serious bug or I am going completely insane. 

es posible obtener la función del contorno de un dominio convexhull ?

 

 

 

Hi, friends!

I'm not a math =) but it is interesting

How can i solve this equation like the gambler's ruin with Maple's function rsolve

f(n)=0.5*f(n-1)+0.5*f(n+1), f(0)=1, f(6)=0

rsolve({f(n)=0.5*f(n-1)+0.5*f(n+1), f(0)=1, f(6)=0}, {f});

it returns this 

{f(n) = 7 f(5) - f(5) (n + 1)}

I don't understand :( 

for example Wolfram Alpha return the true result

Found a strange behaviour in Mapke 2015 of the sqrt-function after loading the GRTensor package:

the square-root of a non-square integer, e.g. sqrt(5), does not terminate. 5^(1/2) instead works fine.

Can be reproduced with Maple 18, but not with Maple 11.

I consider this a serious bug, as it makes any expressions containing such roots useless.

As it worked with Maple 11 I am inclined to see it as your fault.

 

hi.i trust that attached equation has more answer but fsolve only gain some of them!!! how i can gain another that i know value of them?

another root  that i known, are : 0.165237712988657e-1    and     .103583272213766    and    .290071279318035

thanks 

root.mw

First 55 56 57 58 59 60 61 Last Page 57 of 73