MaplePrimes Questions

The documentation for the option AllSolutions for int says that the results are always valid for all real parameter values (in the endpoints). That seems like a pretty major claim. Each of these three is already wrong for a=-1/2, b=1/2:

int(1/ln(t), t = a .. b, AllSolutions);
    piecewise(ln(a) < ln(b), piecewise(And(1 < b, a < 1), undefined, piecewise(a = 1, infinity,
    Ei(1, -ln(a)))+piecewise(b = 1, -infinity, -Ei(1, -ln(b)))), ln(b) = ln(a), 0, ln(b) < ln(a),
    -piecewise(And(1 < a, b < 1), undefined, piecewise(b = 1, infinity, Ei(1, -ln(b)))+
    piecewise(a = 1, -infinity, -Ei(1, -ln(a)))))

int(sqrt(t^2-1+I*t), t = a .. b, AllSolutions);
    piecewise(a < b, (1/2)*sqrt(b^2-1+I*b)*b+I*sqrt(b^2-1+I*b)*(1/4)-3*ln(-2*signum(0, -b, 1)^2*
    b^2+2*sqrt(b^4-b^2+1)*signum(0, -b, 1)^2+4*b*sqrt(2*sqrt(b^4-b^2+1)+2*b^2-2)+2*
    signum(0, -b, 1)^2-2*signum(0, -b, 1)*sqrt(2*sqrt(b^4-b^2+1)-2*b^2+2)+6*b^2+2*
    sqrt(b^4-b^2+1)-1)*(1/16)-3*ln((-I*(signum(0, -b, 1)*sqrt(2*sqrt(b^4-b^2+1)-2*b^2+2)-1+I*
    sqrt(2*sqrt(b^4-b^2+1)+2*b^2-2)+(2*I)*b))*(1/sqrt(-2*signum(0, -b, 1)^2*b^2+2*sqrt(b^4-b^2+1)*
    signum(0, -b, 1)^2+4*b*sqrt(2*sqrt(b^4-b^2+1)+2*b^2-2)+2*signum(0, -b, 1)^2-2*signum(0, -b, 1)*
    sqrt(2*sqrt(b^4-b^2+1)-2*b^2+2)+6*b^2+2*sqrt(b^4-b^2+1)-1)))*(1/8)-(1/2)*sqrt(a^2-1+I*a)*a-I*
    sqrt(a^2-1+I*a)*(1/4)+3*ln(-2*signum(0, -a, -1)^2*a^2+2*sqrt(a^4-a^2+1)*signum(0, -a, -1)^2+
    4*a*sqrt(2*sqrt(a^4-a^2+1)+2*a^2-2)+2*signum(0, -a, -1)^2-2*signum(0, -a, -1)*sqrt(2*
    sqrt(a^4-a^2+1)-2*a^2+2)+6*a^2+2*sqrt(a^4-a^2+1)-1)*(1/16)+3*ln((-I*(signum(0, -a, -1)*sqrt(2*
    sqrt(a^4-a^2+1)-2*a^2+2)+I*sqrt(2*sqrt(a^4-a^2+1)+2*a^2-2)-1+(2*I)*a))*(1/sqrt(-2*
    signum(0, -a, -1)^2*a^2+2*sqrt(a^4-a^2+1)*signum(0, -a, -1)^2+4*a*sqrt(2*sqrt(a^4-a^2+1)+2*
    a^2-2)+2*signum(0, -a, -1)^2-2*signum(0, -a, -1)*sqrt(2*sqrt(a^4-a^2+1)-2*a^2+2)+6*a^2+2*
    sqrt(a^4-a^2+1)-1)))*(1/8), b = a, 0, b < a, -(1/2)*sqrt(a^2-1+I*a)*a-I*sqrt(a^2-1+I*a)*(1/4)+
    3*ln(-2*signum(0, -a, 1)^2*a^2+2*signum(0, -a, 1)^2*sqrt(a^4-a^2+1)+4*a*sqrt(2*sqrt(a^4-a^2+1)+
    2*a^2-2)+2*signum(0, -a, 1)^2-2*signum(0, -a, 1)*sqrt(2*sqrt(a^4-a^2+1)-2*a^2+2)+6*a^2+
    2*sqrt(a^4-a^2+1)-1)*(1/16)+3*ln((-I*(signum(0, -a, 1)*sqrt(2*sqrt(a^4-a^2+1)-2*a^2+2)+
    I*sqrt(2*sqrt(a^4-a^2+1)+2*a^2-2)-1+(2*I)*a))*(1/sqrt(-2*signum(0, -a, 1)^2*a^2+
    2*signum(0, -a, 1)^2*sqrt(a^4-a^2+1)+4*a*sqrt(2*sqrt(a^4-a^2+1)+2*a^2-2)+2*signum(0, -a, 1)^2-
    2*signum(0, -a, 1)*sqrt(2*sqrt(a^4-a^2+1)-2*a^2+2)+6*a^2+2*sqrt(a^4-a^2+1)-1)))*(1/8)+(1/2)*
    sqrt(b^2-1+I*b)*b+I*sqrt(b^2-1+I*b)*(1/4)-3*ln(-2*signum(0, -b, -1)^2*b^2+2*sqrt(b^4-b^2+1)*
    signum(0, -b, -1)^2+4*b*sqrt(2*sqrt(b^4-b^2+1)+2*b^2-2)+2*signum(0, -b, -1)^2-
    2*signum(0, -b, -1)*sqrt(2*sqrt(b^4-b^2+1)-2*b^2+2)+6*b^2+2*sqrt(b^4-b^2+1)-1)*(1/16)-
    3*ln(-(I*signum(0, -b, -1)*sqrt(2*sqrt(b^4-b^2+1)-2*b^2+2)-I-sqrt(2*sqrt(b^4-b^2+1)+2*b^2-2)-
    2*b)/sqrt(-2*signum(0, -b, -1)^2*b^2+2*sqrt(b^4-b^2+1)*signum(0, -b, -1)^2+4*b*sqrt(2*
    sqrt(b^4-b^2+1)+2*b^2-2)+2*signum(0, -b, -1)^2-2*signum(0, -b, -1)*sqrt(2*sqrt(b^4-b^2+1)-
    2*b^2+2)+6*b^2+2*sqrt(b^4-b^2+1)-1))*(1/8))

int(arctan(t+2*I), t = a .. b, AllSolutions);
   piecewise(a < b, piecewise(a < 0, I*arctan(4*a/(a^2-3))*(1/2)+(1/4)*ln(a^2+1)+(1/4)*ln(a^2+9)-
   (2*I)*arctan(2*I+a)-arctan(2*I+a)*a+I*Pi*(1/2), a = 0, -I*Pi+3*ln(3)*(1/2), 0 < a, I*arctan(4*a/
   (a^2-3))*(1/2)+(1/4)*ln(a^2+1)+(1/4)*ln(a^2+9)-(2*I)*arctan(2*I+a)-arctan(2*I+a)*a-I*Pi*(1/2))+
   piecewise(b < 0, -I*arctan(4*b/(b^2-3))*(1/2)-(1/4)*ln(b^2+1)-(1/4)*ln(b^2+9)+(2*I)*
   arctan(2*I+b)+arctan(2*I+b)*b-I*Pi*(1/2), b = 0, -I*Pi-3*ln(3)*(1/2), 0 < b, -I*arctan(4*b/
   (b^2-3))*(1/2)-(1/4)*ln(b^2+1)-(1/4)*ln(b^2+9)+(2*I)*arctan(2*I+b)+arctan(2*I+b)*b+I*Pi*(1/2))+
   piecewise(And(0 < b, a < 0), -(2*I)*Pi, 0), b = a, 0, b < a, piecewise(b < 0, -I*arctan(4*b/
   (b^2-3))*(1/2)-(1/4)*ln(b^2+1)-(1/4)*ln(b^2+9)+(2*I)*arctan(2*I+b)+arctan(2*I+b)*b-I*Pi*(1/2),
   b = 0, I*Pi-3*ln(3)*(1/2), 0 < b, -I*arctan(4*b/(b^2-3))*(1/2)-(1/4)*ln(b^2+1)-(1/4)*ln(b^2+9)+
   (2*I)*arctan(2*I+b)+arctan(2*I+b)*b+I*Pi*(1/2))+piecewise(a < 0, I*arctan(4*a/(a^2-3))*(1/2)+
   (1/4)*ln(a^2+1)+(1/4)*ln(a^2+9)-(2*I)*arctan(2*I+a)-arctan(2*I+a)*a+I*Pi*(1/2), a = 0,
   I*Pi+3*ln(3)*(1/2), 0 < a, I*arctan(4*a/(a^2-3))*(1/2)+(1/4)*ln(a^2+1)+(1/4)*ln(a^2+9)-(2*I)*
   arctan(2*I+a)-arctan(2*I+a)*a-I*Pi*(1/2))-piecewise(And(0 < a, b < 0), -(2*I)*Pi, 0))

The first one probably has the correct answer inside, but it has conditions like ln(a)<ln(b), so that case never gets selected when the values are complex.

I was told that the following workout was done in Maple.  I have tried to read material about how to do it but I am completely lost.  Can someone indicate me where I can read in oder to do what the image says or give me some tips please?  

where all functions are dependent on the variables (u,v).

Observation: subscripts means partial derivatives of the function while superscripts are just for naming different functions,i.e Gamma^1 and Gamma^2 are two functions.

 

Sergio

 

Sergio

 

f := (z, t) -> ln(t)^2/((t^2+1)*(t-z));

int(f(z, t), t = 0 .. infinity) assuming Im(z) > 0;
       int(ln(t)^2/((t^2+1)*(t-z)), t = 0 .. infinity)

int(f(a + I*b, t), t = 0 .. infinity) assuming a::real, b > 0; # 0*infinity
       -sqrt(a^2+b^2-2*b+1)*signum(I*arctan(b, a)-I*arctan(-b, -a)-I*Pi)*
       infinity/((I*b-I+a)*(I*b+I+a))

int(f(z, t), t = 0 .. infinity) assuming Re(z) > 0, Im(z) > 0;
       int(ln(t)^2/((t^2+1)*(t-z)), t = 0 .. infinity)

int(f(a + I*b, t), t = 0 .. infinity) assuming a > 0, b > 0;
      -((3*I)*Pi^3*b+3*Pi^3*a-(16*I)*Pi^2*arctan(b/a)-(6*I)*Pi*ln(a^2+b^2)^2+(24*I)*
      Pi*arctan(b/a)^2+(6*I)*ln(a^2+b^2)^2*arctan(b/a)-(8*I)*arctan(b/a)^3-8*Pi^2*
      ln(a^2+b^2)+24*Pi*ln(a^2+b^2)*arctan(b/a)+ln(a^2+b^2)^3-12*ln(a^2+b^2)*
      arctan(b/a)^2)/(24*(-b^2+(2*I)*a*b+a^2+1))

So it looks like the first three can be made to work as well (and the result in terms of z will be much neater).

 

the Eigenvalues are showing with I and am not expecting a complex eigenvalues so what is that I stand for? Can you please help? Thank You

Hi

I've got this list:

L:=[[TC,DB], [], [TD,JK], [IW,CM], [], [KJ,DJ]]

What command to remove the 'null sets', leaving :=[[TC,DB],[TD,JK], [IW,CM], [KJ,DJ]]

this doesn't work:

remove(has, L, 0)

a::real and b::real;
Error, type `real` does not exist

`and`(a::real, b::real);
                       a::real and b::real

Since a::real by itself is fine, why does the conjunction give an error?

These two are evaluated differently as well:

Im(a) > 0 and a <> 0;
                           0 < Im(a)
`and`(Im(a) > 0, a <> 0);
                      0 < Im(a) and a <> 0

The consequence is that these two will work differently:

is(a <> 0) assuming Im(a) > 0 and a <> 0;
                             FAIL

is(a <> 0) assuming `and`(Im(a) > 0, a <> 0);
                             true

is(a <> 0) assuming Im(a) > 0; # sadly, just Im(a) > 0 is not enough
                             FAIL

It looks like (a and b) and `and`(a, b) just do completely different things:

x := proc() local r; r := rand(); print(r); r end proc;

a > x() and a > x() and a = a;
                          395718860534
                        395718860534 < a

`and`(a > x(), a > x(), a = a);
                          193139816415
                          22424170465
         193139816415 < a and 22424170465 < a and a = a

That is, (a and b) first simplifed the expression and then evaluated a>x() once, but `and`(a, b) evaluated the arguments without doing any simplifications.

Also, should this work (it works with `real`, which supposedly doesn't exist):

Re(a+b) assuming a::realcons, b::realcons;
                           Re(a + b)

 

I've compiled a program with Maple 18. unfortunately it gives me the message "Large sortie de plus de 1000000 noeuds". If any one have any idea to resolve this problem. Thanks you very mutch for your support. 

With Startup Code Editor opened, save and close the parent worksheet. Any unsaved changes in the Code Editor are lost.

 

I tried to solve 4 simultaneous equations but the result is an empty matrix.

Is that because there is no solution or did I make something wrong?

You can download the file by using the link below.

SimEquations.mw

I have noticed that I don't receive e-mails anymore when contributions are submitted to my subscriptions.
I used to.
Has this happened to anyone else?

It is embarrasing to have asked somebody a question and gotten a reply you are not made aware of.

What to do about it?
 

hi 

I have a matrix (for example a 6*6 matrix) and I want to add a row and a column between row and column of number 3 and 4. it means that finally, we have a 7*7 matrix.

tnx

 

The documentation says that Subgroup(c) returns the subgroup of the coset c.

Also, the display of Elements(lc) is garbled in Maple as well, not just in the image.

with(GroupTheory)

sgr := Subgroup(Elements(SymmetricGroup(4)), SymmetricGroup(5))

_m917312091392

(1)

lc := LeftCoset(Perm([[1, 2]]), sgr)

_m917308956640

(2)

Subgroup(lc)

Error, invalid input: GroupTheory:-Subgroup expects its 1st argument, generators, to be of type {list, set, identical(undefined)}, but received _m917308956640

 

Elements(lc)

{_m917312836128, _m917312839712, _m917312840544, _m917312841376, _m917312842272, _m917312843296, _m917312844224, _m917312845152, _m917312846240, _m917312847072, _m917312848768, _m917312849696, _m917312850624, _m917312851552, _m917312852288, _m917312853248, _m917312854144, _m917312855072, _m917312856000, _m917312857472, _m917312858592, _m917312859552, _m917315257248, _m917315263744}

(3)

 

I want to highlight the intersection between 2 graphs; pp(m,a) and the plane m=-0.2

PP := .8707945038*exp(-50.00000000*(m-0.842e-1)^2+(2.745342070*(m-0.842e-1))*(a-2.3722)-.1046792095*(a-2.3722)^2)

How can I do that?

Thank you

Hello every body

I have a plot that contains two curve. I need to chracterize the curves with `min(D_{T})` and `max(D_{E})`. Is it possible to write legend of plot such that when include it in latex file,  T and E be presented in indices?

i used this  commend to plot three equation and it actually work .my question her how to show any point cordinates just i pointed on it ? could i do it by maple or what ?
help please 
implicitplot3d({f[1], f[2], f[3]}, Q[h] = 0 .. 100, S[h] = 0 .. 100, R[h] = 0 .. 100);

1 2 3 4 5 6 7 Last Page 1 of 1499