MaplePrimes Questions

Hello everyone,

I need to expand this:


so I have as a final expressioin


This is the property that has to be true so the probability is one.

>assume(u::complex, v::complex);  does not work and the simple commande

>expand() neither.  Any idea?

Let us assume we have a square matrix $A$ with known entries. We know the system is asymptotically stable by inspection of the eigenvalues, hence we know that the Lyapunov equation


will have a symmetric matrix P as a solution with Id being the identity matrix. How can I obtain the entries of the P matrix? I have tried the following code but was not able to solve for the entires of the P matrix.


I would be happy if someone could share a solution to this problem.



A := <<a__11|a__12>,<a__21|a__22>>;
P := <<p__11|p__12>,<p__12|p__22>>;
Id := <<1|0>,<0|1>>;
eqn := Transpose(A).P+P.A = -Id;


I am having a problem with image tools read a file from workbook.

In one document i get

Error, (in ImageTools:-Read) unable to open database file

In another document I cant save the file. 

Please find my error in the doc which produces outputlines always. Is there any otherway (other than view -  show contents and uncheck output which hides all other outputs as well in a programm) by commands, which can suppress these output lines.

Also is there any other way of giving continuously varying inputs (and visible in RotaryGauge ) like sin wave form (0 to 8 to 0 in a period) ?

with(DocumentTools); for i to 12 do SetProperty("RotaryGauge0", value, 6); SetProperty("RotaryGauge1", value, .66*GetProperty("RotaryGauge0", value))*time*[]; 3^100; GetProperty("RotaryGauge0", value); SetProperty("RotaryGauge0", value, 4); SetProperty("RotaryGauge1", value, .7*GetProperty("RotaryGauge0", value))*time[real](); 3^100; SetProperty("RotaryGauge0", value, 2); SetProperty("RotaryGauge1", value, .5*GetProperty("RotaryGauge0", value))*time[real](); 3^100 end do








I am trying to do parallel computing by using task programming. In the attached file, I am having the following error:

Error, (in generate) type 'sizeType' does not exist.

The error comes out only when the size of the computation is large enough and so the Threads piece of my "generate" function is run.

Could you please help me to understand what is wrong?

Here is the link to my .mv file:



I would like to remove asterisk * from a string....

such as...



Any ideas would be appreciated???

Dear All,

I am trying to formulate and solve a linear system of equations to determine the 15 unknowns (link to the script below.) The funcion solve does not seem to return a vlaue due to the lengthy and complex nature of the equations neither do linsolve or LinearSolve. I am wondering whether there exists an argument of another trick that can be employed in order to tackle with such complicated systems. Any help or hint is highly appreciated and desirable.

Thank you,


Federiko script:


I have a functional equation in 2 variables, and I want to write its solution(s) as (Puiseux) power series. When the order is relatively big, the output appears not in the increasing order of powers. For example, I execute 

map(allvalues, algeqtoseries(P, t, u, 60, true))

and get as one of the solutions

O(t^(65/3)) + 13*t^9*(1/3) + 2*t^5*(1/3) + 115*t^15 + 37*t^13 + 5*t^7*(1/3) + 37*t^11*(1/3) + (1/3)*t^3 + 11836*t^21*(1/3) + 1193*t^19 + 1100*t^17*(1/3) [...]

What should I do in order to get the series with terms appearing in the increasing order of powers of t?



Maple will multiply two matrices for me perfectly.  

However, when I change the value in one of the cells and try to re-evaluate the product, Maple starts giving me the sum of something in all of my cells.   


Are you simply not allowed to ever change the initial values you set for a matrix?  This doesn't seem to make sense to me.  

I have attached a screen shot to show an example with a random matrix.  

Dear Maple users,

I recently had to deal with the following sum :


for an arbitrary positive value of k. This sum obviously seems to diverge but Maple gives me 0 as a result for it. It looks very strange to me, but it might be the same sort of result as the sum of all the positive integers is -1/12. I am wondering whether the result Maple gives me is correct, and if yes, where does it come from ?

Any idea ?

Thanks a lot!

for example i have this equation: y=x+c , y and x are variable, but c is constant. i want to plot implicitly with respect to y and x, but the result has no line or something and just two axes. what can i do?

If I have 4 latitude and longitude coordinates in degree.  How can I convert all of them to radian and assign each to a specific variable inside a procedure?


> DBP(36.9,102.04,37.345,101.6104)
Error, invalid loop statement termination



is(k::OrProp(even, odd)) assuming k::integer;

is(k::Or(posint, nonposint)) assuming k::integer; # a synonym for OrProp

is(`or`(k::posint, k::negint, k = 0)) assuming k::integer;

Not including the infix or because that would fail for k=0 for different reasons.


Hey, friends, I have a differential equation to confuse me too much!  The question is shown in the below:


y''(x)-0.00003019*y(x)^0.337=0  ,       0<x<2945


y(0)=0,    y'(2945)=0.0116 


I want to have the numerical solution, I used to plot y(x) and y'(x) , but errors shows me "jacobian is singular at left endpoint, use midpoint method instead."  Thank you for you any help in advance!!!  

Can you give me the code to figure it out in Maple?  Thank you so so so so much!


eliminate({-a+x+y+z, -b+x*y+y*z+x*z, -c+x*y*z, -P+x^4+y^4+z^4}, {x, y, z});
    [{x = (sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2+a*b*c-c^2)/(2*b*c),
      y = -(sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2-a*b*c+c^2)/(2*b*c),
      z = c/b}, {-a^4+4*a^2*b-4*a*c-2*b^2+P}],
    [{x = -(sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2-a*b*c+c^2)/(2*b*c),
      y = (sqrt(c^2*(a^2*b^2-2*a*b*c-4*b^3+c^2)/b^4)*b^2+a*b*c-c^2)/(2*b*c),
      z = c/b}, {-a^4+4*a^2*b-4*a*c-2*b^2+P}]

The expression for P is fine, but how are x and y square roots and z=c/b??


eliminate({y*z, x+y-z}, z);
                       [{z = 0}, {x + y}]

Is this supposed to happen? I was expecting y*(x+y):

Groebner:-Basis({y*z, x+y-z}, lexdeg([z], [x, y]));
                     [x*y + y^2, z - x - y]


First 824 825 826 827 828 829 830 Last Page 826 of 2325