Unanswered Questions

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

Is package SumTools supposed to be a superset of the sum command in Maple? Or is it supposed to be a complementary to it? When should one use SumTools vs. sum?

The help on SumTools says

The SumTools package contains commands that help find closed forms of definite and indefinite sums. The package consists of three commands and three subpackages.
 

But does not make it clear when to use it vs. sum or what is its relation to sum command.

I am asking because on one example, sum is able to give an answer, while SumTools can't. I have expected the other way, where if sum failed, then one will try SumTools.

restart;
sum(1/n*sin(n*x),n=1..infinity) assuming 0<x,x<Pi

SumTools:-IndefiniteSummation(1/n*sin(n*x),n=1..infinity) assuming 0<x, x<Pi

 

 

Hello,

I am trying to solve CVRPSD by integer L shape method. I have some problems.

d=(0, 20,20,25,15,10): Demand
p=(1,1,1/3,1/2,1/4): the probability demand

K=50
1) How to write properly Recourse Policy (Q) in maple (I attached my problem) ? Recourse (Penalty cost) is applied when the demand exceed the vehicle capacity. We know that in CVRPSD, we don't know the demand until vehicle arrive at the node(customer).
(Look at my program)

Hopefully anybody could help me to solve this.

Thank you any way.
Exact_CVRPSD.mw


for example, there is an expression exp1=(a+b)*(a+c),  I want to find all the position in DAG of a specific operand, such as 'a' ; also I want to find the position in the DAG of  sub-expression which type is '+' .
 is there any function in maple which can return the  integers indicating position of an operand or type in an expression? i.e. with this function, assume its name is 'find',then
find(exp1, a) will return  [1,1],[2,1];
find(exp1,'+') will return  [1,0],[2,0];

hi

i got a problem to solve 2 nonlinear  pde equation that are not coupled but there is 8 coupled boundary condition. i can use pdsolve in maple.anyone can help me in numerical solving of these pdes?

This is using Maple 2018.2.1 and using Physics version MapleCloud 289 on windows 10, 64 bit with lots of RAM (64 GB). The call to timelimit hangs. Is there a workaround?

restart;

#and in new execution cell

pde :=  diff(w(x,y),x)-((k+1)*x^k*y^2- a*x^(k+1)*tan(x)^m*y + a*tan(x)^m )*diff(w(x,y),y) = 0;
cpu_time := timelimit(180,CodeTools[Usage](assign('sol',pdsolve(pde,w(x,y))),output='realtime'))

 

 

 

 

Hello,

I'm a long-time Maple user and have used the CodeGeneration package in the past, but with a recent download of Maple 2018, I'm now getting an error I have never gotten before. Even when I do something simple (see below, for example), I get an error that reads "Error, (in clear) argument `` is incorrect or out of order". And I get this error no matter which language I choose. Any help would be much appreciated.

> with(CodeGeneration)

> R(exp(x))

Thanks!

 

 

Suppose that A is an nxn matrix over the finite field Z:=GF(2,q) for some q. I wan to get the smitform of A over Z. First I used the package  

with(LinearAlgebra[Generic]) 

and after that I applied the command 

S := SmithForm[Z](A)

but the mentioned command made some errors. In fact, I do not how to define commands igcdex, iquo, irem, sign and abs for SmithForm over finite fields.

Thanks for any suggestions 

how I can find solutions for non linear equations.

I want to find non zero solution.

thank you

ZAH.mw



Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ZAH.mw .
 

Download ZAH.mw

I get this problem in a paper published in 2011. there is a table which compares the solution of a differential equation and gives the range of initial guesses. I got confused when saw maple and approximated maple, please any help me 

I'm attaching the link of paper. Please have a look at table 3. 

https://www.sciencedirect.com/science/article/pii/S0098135411000056

Dear all,
How can I input different spacesteps in numerical solution of PDE (Heat equation) with pdsolve of Maple?

For example, the x range is x=0..L,
and I'd like to solve the PDE with spacestep1=L/100 for x=0..a and spacestep2=L/10 for x=a..L.

Thank you in advance!

Why Maple doesn't calculate this and only rewrites it??

 

with(inttrans);

invlaplace(exp((0.2500000000e-1-2.500000000*sqrt(0.116e-3+.8*p-3.2*10^(-10)/(p+0.2e-4)))*x)/p, p, t);
     
 

iqt + aqxy + ibq (qq*x − q*qx) = 0. write this equation in maple

how i can solve the eqution by ritz methode codes?


I have an object in 6d I'd like to visualise. The region of 6d space I am interested in is described by these equations:

{f[10] = -(.2000000000*(5.*f[21]*f[20]*f[22]-5.*f[20]*f[22]^2+20.*f[20]*f[21]-20.*f[20]*f[22]+135.*f[20]+46.*f[21]))/(f[21]*(f[21]-1.*f[22])),
f[11] = -1.*f[22]-4.,
f[12] = -(1.*(f[22]^2+4.*f[22]-27.))/f[21],
f[20] = f[20],
f[21] = f[21],
f[22] = f[22]}

clearly the first three variables are dependant, and the latter three are independant. I'd like to graph the first three as the latter three vary between bounds and then colour the points on the output based on where they came from in the input, so i can get some intuition about what these equations mean.
 

Good day everyone.

Please can help me with this code on the taylor series expansion involving Fractional Differential Equation (FDE)? Particularlly, the lines highlighted blue and green respectively in relation to FDE.

Thank you and kind regards

#k=2
restart:q:=n*h:
P:=sum((a[k]*x^(k))/GAMMA(k+1-alpha), k=0..3):
assume(alpha>0,alpha < 1):
Q:=fracdiff(P,x,alpha):
e1:=simplify(eval(P, x=q+h))=y[n+1]:
e2:=simplify(eval(Q,x=q))=f[n]:
e3:=simplify(eval(Q,x=q+h))=f[n+1]:
e4:=simplify(eval(Q,x=q+2*h))=f[n+2]:
var:=seq(a[i], i=0..3):
M:=e||(1..4):

Cc:=eval(<var>, solve(eval({M}),{var}) ):
for i from 1 to 4 do
	a[i-1]:=Cc[i]:
end do:
Cf:=P:
E:=collect(Cf, [y[n+1], f[n], f[n+1],f[n+2]], recursive):
print():
s2:=y[n+2]=simplify(eval(Cf, x=q+2*h)):
collect(%, [y[n+1], f[n], f[n+1],f[n+2]], recursive):

s1:=y[n]=simplify(eval(Cf, x=q)):
collect(%, [y[n+1], f[n], f[n+1],f[n+2]], recursive):

Y[n+1]:=convert(taylor(y(x+h),h=0,12),polynom):
F[n]:=convert(taylor((D(y))(x), h = 0,12), polynom):
F[n+1]:=convert(taylor((D(y))(x+h), h = 0,12), polynom):
F[n+2]:=convert(taylor((D(y))(x+2*h), h = 0,12), polynom):


W:=asympt(expand(eval(rhs(s2),[y[n+1]=Y[n+1],f[n]=F[n],f[n+1]=F[n+1],f[n+2]=F[n+2]])),h,6);
X:=convert(taylor(y(x+2*h),h=0,12),polynom)-W;
lte:=convert(asympt(X,h,8),polynom);

 

 

First 9 10 11 12 13 14 15 Last Page 11 of 247