Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

I know that is easy for you but i am new in this topics.I want to write a code to generate a matrix which have entries consist of married.

L := sum( 1/ln(k), k=2..n ) * ln(n)/n;
limit(L, n=infinity);
# Should be 1

Just curious: in Maple 2017, is it OK?


Hi, i don't know what happened but every time that I run the simulation appears this: 

invalid input: Multibody:-GetMultibodyData expects its 4th argument, lsProbes, to be of type list(list), but received registerMultibodyVariable

If anyone can solve this problem, I am going to be very grateful. Thanks



Dear sir,

in the program boundary conditions D(f)(0)=0 doesn't showing result but when use d(f)(0)=1 it will execute, why is this can you explain this ?program.mw

Pls do we have any Maple code for Adomian Decomposition Method?

So, I'm trying to delelop an algorithm for the method of multiple scales. Starting with a simple ODE:

diff(x(t), `$`(t, 2))+x(t) = 0

After scaling, it should be written in the form:

(d/dT[0]+epsilon*d/dT[1]+epsilon^2*d/dT[2])^2*(epsilon^3*X[3]+epsilon^2*X[2]+epsilon*X[1])+epsilon*X[1]+epsilon^2*X[2]+epsilon^3*X[3] = 0

A proto-algorithm would be:

ode := diff(x(t), `$`(t, 2))+x(t) = 0;
i_ini := 1; i_fin := 3; j_ini := 0; j_fin := 2;
PDEtools:-dchange({t = sum(epsilon^j*T[j], j = j_ini .. j_fin), x(t) = sum(epsilon^i*X[i](T[1]), i = i_ini .. i_fin)}, ode, [{T[0], T[1], T[2]}, {X[1], X[2], X[3]}])

It is not working, though. Could anyone help me out?

Thanks in advance.

hello everyone. i wanna ask. how to do a coding of stability region on diagonally implicit Runge-Kutta for order 3 and order 4? Thank you in advanced :) 


I am trying to solve analytically a simple system of partial differential equations with boundary conditions and I am not able to do it. Even in the very simple case of

pdsolve([diff(u(x, y, t), y, y) = 0, diff(p(x, y, t), y) = 0, u(x, 0, t) = 1, (D[2](u))(x, 1, t) = 0, p(x, 1, t) = 2], [p(x, y, t), u(x, y, t)]);

I don't get any answer.  However if I remove the boundary conditions I get the right answer

pdsolve([diff(u(x, y, t), y, y) = 0, diff(p(x, y, t), y) = 0], [{p(x, y, t), u(x, y, t)}]);
 {p(x, y, t) = _F3(x, t), u(x, y, t) = _F1(x, t) y + _F2(x, t)}

Can maple 2015 solve analytically systems of partial differential equations with boundary conditions? I have not been able to find any example anywhere.

Thanks a lot for your help.



HI Maple primes.  We try to make sense of 'rsolve' in the Maple world.  What is command for source code?





Hi Mapleprimes,

We know that '' rsolve '' is a recurrence equation solver.  It is more than an expression simplifier.

Congratulations to the Maple computer algebra team for creating such a great computer tool.  simply want to know more.


Surely there are many steps to determine the values to place.




I have an important question regarding the symbolic toolbox from maple for matlab

I have declared three symbolic variables:

syms t;

x = sym('x(t)', 'real');

y = sym('y(t)', 'real');

z = sym('z(t)', 'real'); 

f = x + y + z;

Now i want to  calculate the derivation of f regarding the variable t:

df = diff (f, t);

The result should be 

df = dx(t) / dt + dy(t) / dt + dz(t) / dt

but  instead i get the result df = 0. If i use the command diff (f) i get the result

df = dx(t) / dt 

Does anyone know how to fix this problem?


Hi Guys,


got my last Problem solved :) Now i do have a new one...

restart; with(RealDomain); with(CodeGeneration); with(ExcelTools); with(plots);

n := 1; 
for i to 256 do 
for j to 256 do 
r := evalf(sqrt((i-dx1)^2+(j-dy1)^2)); 
Ints := R0[i, j]; 
IntsR[n] := [r, Ints]; 
n := n+1 
end do; 
end do; 
IntsR := [seq(IntsR[i], i = 1 .. n-1)]

The list IntsR consists of unsorted values doublets. Now it would be nice to get some sort of mean value of my Ints over r.  Didnt really find a solution for it until now... do you have any hint?

I saw a question earlier today regarding 3 particles in a gravitational field.  It has since disappeared. 

I see no rhyme or reason why it was deleted by someone, whoever did.  It was Math related and it was Maple related hence no reason for it's deletion.  This happened just before Mapleprimes and Maplesoft went offline for a couple of hours. 

Can someone explain?  Did anyone else see it?

i count the number among group
but when the list a large such as over 1000 records, the count will be over 30,000
use which denominator to find probability?
is there any functions in maple for this case?
correlationlist1 := [[1,2,3],[1,3,5]....]:
PAB := Matrix(50):
for ii from 1 to nops(correlationlist) do
 for jj from 1 to nops(correlationlist[ii]) do
  for kk from 1 to nops(correlationlist) do
   for qq from 1 to nops(correlationlist[kk]) do
    if ii <> kk then
     PAB[correlationlist[ii][jj],correlationlist[kk][qq]] := PAB[correlationlist[ii][jj],correlationlist[kk][qq]] + 1: # group to group relations
    end if:

If there is a list


how to use correlate function?

assume [1,2] and [2,1] count as 2 

when find correlation between 1and 2

and how about 

if [1,2] and [2,1] are different

first position as key

means 1 relate with 2 count as 1 case for 1 only

[2,1] count as 1 case for 2 only


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