Maple 18 Questions and Posts

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

Hello all,

So far I have been unable to find this question anywhere, but I apologize if it is a duplicate. I'm trying to evaluate the integral of sechq(x), where q is a positive integer. Mathematica is able to tell me the result (a hypergeometric function), but for some reason, Maple seems not to be able to compute this integral, it just gives me back the integral. A higher info-level on the 'int' function reveals a line that says 'Risch d.e. has no solution', but I'm not sure if that has anything to do with my problem. Any suggestions or tips on how to get an answer out of Maple would be greatly appreciated!


I am working with the Maple 18.02 version. I just want want to perform a basic polynomial expansion using the command "expand" and it does not respond as it should according to what Maple Programming Help says it would. For example:

Maple Programming Help says:

I get:


Also, one sees that this isn't even true, as x(x+2) + 1 = x^2 +2x +1, which is not equal to x^2 + 3x +2.

Moreover, maple tells me it is equal..:

What is going on here? I woul like to get the full expanded form (without factors). Also, this is obviously not true, or maybe Maple means something else by x(x+2) +1...

Thank you!

I have a Maple code which generates a matrix, saves it to a .txt file and this is then read in to a C++ program. I have hit a snag with these matrices, in that they are generating absolutely enourmous .txt files. I need to get to a 5000*5000 matrix yet a 200*200 is generating a 100MB file.

The matrix elements contain a lot of algebraic terms which I would like to keep general as these are defined in the C++ code. They also contain hypergeometric functions which in the example below I have left unsimplified (although they are simplified before reaching C++). I have tried various operations/combinations to simplify, but the file sizes still come out very large.

From previous experience expanding the expressions and then simplifying allows Maple to "do more" with it, but it does not seem to work in this instance. I have tried map(options,expr), simplify(expr,options), combine(expr,options), convert (expr,options) etc... The script attached only contains a small example 10*10 matrix on its own without the code which generates it due to the size of the code. What is the best way to simplify these matrices to generate the smallest .txt file?

Any help is appreciated.



To check my arguments in a procedure I need something like

myproc := proc(M :: Matrix(square, rational)

                       , N :: Matrix(shape=triangular[lower, unit], datatype = rational
                       , O :: Matrix(shape = square, dimension = 5


end proc;

How does that work in Maple? What is the correct Syntax? I tried many different things that doesn't fit.


I started using Maple recently. The output sometime is not simplified like that in the attached picture. How can I simplify such expressions?


I am trying to write a procedure to see, which (di)-graphs are isomorphic (here represented by there 3*3 adjecency-matrices). When I try the procedure for all 3*3-matrices with entries in {0,1} (there are 512 of them), I get the following error:
"Error, (in GraphTheory:-IsIsomorphic) invalid subscript selector"

Can you possibly say, what I am doing wrong? My code is the following:

getIso3 := proc(liste)
  local i,k,M1,c,d:
  c := 0:
  M1 := [[liste[1]]]:
  for i from 2 to numelems(liste) do
    for k from 1 to numelems(M1) do
        if IsIsomorphic(Digraph([a,b,c], liste[i]),Digraph([a,b,c],M1[k][1])) then
          M1[k] := [op(M1[k]), liste[i]]:
          c := 1:
        end if:
    end do:
    if c=0 then
      M1 := [op(M1), [liste[i]]]:
      c := 0:
    end if:   end do: 
end proc:


My input is a list (JJ) of 512 3*3 matrices constructed the following way :

all9Perm := proc(list)
  local P,i,m,n,A:
  P := list:
  for i from 0 to 9 do
    m:= i:  n:= 9-i:
    A := combinat:-permute([1$n, 0$m]):
    P := [op(P), op(1..numelems(A),A)]:
end proc:
K := []:
L := all9Perm(K):
listoflistsToListofmatrices := proc(liste)
  local M,i:
  M := []:
  for i from 1 to numelems(liste) do
      M := [op(M), Matrix([
          [ liste[i][1] , liste[i][2] , liste[i][3] ],
          [ liste[i][4] , liste[i][5] , liste[i][6] ],
          [ liste[i][7] , liste[i][8] , liste[i][9] ]])
  end do:
end proc:
JJ := listoflistsToListofmatrices(L):


When I run this procedure on some of the 512 matrices it does work, but it crashes somewhere around matrix 350. I have try so split the list of the 512 matrices, and I am able to run the procedure on these splits, but this is very inconvenient :-)

I hope you can help me. Also if this can be done in an easier way - I am new to programming and recieve help with a smile.

Yours, Tomas.


I want to fit a curve to the function that you can see in the picture "maple4.png" (in the appendix) to determine two parameters. This works in Maple 2.
Unfortunately Maple 18 shows me an error. Can you help me with this problem?
Best regards



Hi I am getting this message while soliving my first order non-linear initial value problem.

please share me your thoughts

ode := dsolve({diff(y(x), x)-2*(diff(y(x), x))^3 = -2*x+2, y(0) = 1}, type = numeric, range = 0 .. .24);
Error, (in DEtools/convertsys) unable to convert to an explicit first-order system



Hi everybody

In the attached file, when I run the code an error appears while calculating "omegaL1" parameter. The dimensions of the matrices are correct but the source of error is unknown to me. Can anyone help me?

Thanks in advance 

This must have a simple answer but I have been unable to figure it out after many attempts. 

I am trying to create a Clifford algebra, and then use the results in the multiplication table "MT". The multiplication table elements are correct as displayed, but I don't know how to access the results in the table (i.e. the products of the basis elements).  For example, trying to access the table results as matrix elements like MT[2,3] doesn't work, presumably because it is not a matrix. In other words, I need a matrix that contains the same information as the multiplication table.

DGsetup([x, y, z], M);
I12 := Matrix([[-1, 0, 0], [0, -1, 0], [0, 0, -1]]);
AD3b := AlgebraLibraryData("Clifford(3)", Cl3Q, quadraticform = I12);
DGsetup(AD3b, '[e0, e1, e2, e3, e12, e13, e23, e123]', '[omega]');
MT := MultiplicationTable(Cl3Q, "AlgebraTable");


Hi everyone, i am using Maple 18 and i have a problem in converting a equation to a nice polynomial form (a cubic equation with a form of A*x^3+B*x^2+C*x+D), can anyone please help me on the command? Thanks in advance.

My equation is "  d := s*x*(E*K*q-K*r+K*sigma[1]+r*x)*(1-x*(E*K*q-K*r+K*sigma[1]+r*x)/(K*sigma[2]*L))/(K*sigma[2])+sigma[1]*x-x*(E*K*q-K*r+K*sigma[1]+r*x)/K  " or for simplicity is


Can someone please teach me on the command? Really appreciate the help!

Hello, thanks for read me

I don't know why my code don't work, I'm trying to calculate the magnitude of a complex vector but I get a error in the next image you can view it, 


Can someone help me? thank you

Greetings People!

I am modelling a dynamical system using variational formulation. The final step is to plug in, the energy values, in the Euler-Lagrange Equation. The functional has both explicit dependence on time, and implicit dependence (in the form of the generalised coordinate). Since motive of the exercise is to obtain the governing equations for the generalized coordinates, this implicit dependence is unknown. How can Maple be used to derive the differential equations in this situation?


As an example, consider,

Here, ld and l are variables, which are function of time. In order to calculate the term,

How do I proceed? Thanks in advance.

Good day, I need to 2D plot from points which I get by solving numerically 17-degree equations. Firstly I take an only first solution of the equation as below, and I have a Matrix with 3 column which represents X, Y, and VALUE respectively. On the left side should be the first column ( X ), and right axis Y (2. column). As seen from Matrix X and Y getting the value between 0 and 10. Is there any option that I can plot my data in 2D?  Thanks in advance.

points := seq(seq(seq(Fun[n, i, j], n = 0 .. step), i = 0 .. step), j = 1);
Mat := Matrix((step+1)^2, 3, [points]);
         .                             X                  Y                 VALUE
pointplot3d(Mat, style = point, color = black);

i want to solve the system of equation ( 1 )  , (2)  ,  (3)   under the assumation that x , y have the CDF in (4)  ,  (5)

diff(L(lambda[1], lambda[2], alpha), lambda[1]) = n/lambda[1]+sum(x[i], i = 1 .. n)-(sum(2*x[i]*exp(lambda[1])/(exp(x__i*`λ__1`)-1+alpha), i = 1 .. n))

diff(L(lambda[1], lambda[2], alpha), lambda[1]) = n/lambda[1]+sum(x[i], i = 1 .. n)-(sum(2*x[i]*exp(lambda[1])/(exp(x__i*`λ__1`)-1+alpha), i = 1 .. n))


diff(L(lambda[1], lambda[2], alpha), lambda[2]) = m/lambda[2]+sum(y[j], j = 1 .. m)-(sum(2*y[j]*exp(lambda[2])/(exp(y__j*`λ__2`)-1+alpha), j = 1 .. m))

diff(L(lambda[1], lambda[2], alpha), lambda[2]) = m/lambda[2]+sum(y[j], j = 1 .. m)-(sum(2*y[j]*exp(lambda[2])/(exp(y__j*`λ__2`)-1+alpha), j = 1 .. m))


diff(L(lambda[1], lambda[2], alpha), alpha) = (n+m)/alpha-(sum(2/(exp(x[i]*`λ__1`)-1+alpha), i = 1 .. n))-(sum(2/(exp(y[j]*`λ__2`)-1+alpha), j = 1 .. m))

diff(L(lambda[1], lambda[2], alpha), alpha) = (n+m)/alpha-(sum(2/(exp(x[i]*`λ__1`)-1+alpha), i = 1 .. n))-(sum(2/(exp(y[j]*`λ__2`)-1+alpha), j = 1 .. m))


G(x, lambda[1], alpha) = 1-alpha/(exp(lambda[1]*x)-1+alpha)

G(x, lambda[1], alpha) = 1-alpha/(exp(lambda[1]*x)-1+alpha)


G(y, lambda[2], alpha) = 1-alpha/(exp(lambda[2]*x)-1+alpha)

G(y, lambda[2], alpha) = 1-alpha/(exp(lambda[2]*x)-1+alpha)






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