Advanced Search

ablmf

25 Reputation

4 Badges

12 years, 345 days
Contact ablmf

MaplePrimes Activity


These are questions asked by ablmf

How to get the name part of an indexed n...

ablmf 25
evalf
January 29 2018
1 1

For example, if we have

 

 

How can get the 0.312 part of a? I have tried evalf and op without success.

A bug in MultiSeries?...

ablmf 25 Maple
series multiseries
January 16 2018
2 6

I am trying to use MultiSeris package to expand a function (nu3 in the code).

The function has two variables, w and u. I want Maple to treat w as a constant and u as the variable.

As you can see from the code, somehow if I replace w by any rational number and do the expansion, I get an expansion with a constant * u term. However, if I replace w with zeta(3) (or pi), then some how the u term vanishes. I think this cannot be correct. Is this a bug?
 

nu3 := -sqrt(-2*sqrt(w)*sqrt(w+2)*u^4+2*u^4*w+4*sqrt(w)*sqrt(w+2)*u^2+u^4-4*u^2*w-2*sqrt(w)*sqrt(w+2)-2*u^2+w^2+2*w+1)+sqrt(-2*sqrt(w)*sqrt(w+2)*u^4+2*u^4*w+4*sqrt(w)*sqrt(w+2)*u^2+u^4-4*u^2*w-2*sqrt(w)*sqrt(w+2)-2*u^2+2*sqrt(-2*sqrt(w)*sqrt(w+2)*u^4+2*u^4*w+4*sqrt(w)*sqrt(w+2)*u^2+u^4-4*u^2*w-2*sqrt(w)*sqrt(w+2)-2*u^2+w^2+2*w+1)+2)+w-1

-(-2*w^(1/2)*(w+2)^(1/2)*u^4+2*u^4*w+4*w^(1/2)*(w+2)^(1/2)*u^2+u^4-4*u^2*w-2*w^(1/2)*(w+2)^(1/2)-2*u^2+w^2+2*w+1)^(1/2)+(-2*w^(1/2)*(w+2)^(1/2)*u^4+2*u^4*w+4*w^(1/2)*(w+2)^(1/2)*u^2+u^4-4*u^2*w-2*w^(1/2)*(w+2)^(1/2)-2*u^2+2*(-2*w^(1/2)*(w+2)^(1/2)*u^4+2*u^4*w+4*w^(1/2)*(w+2)^(1/2)*u^2+u^4-4*u^2*w-2*w^(1/2)*(w+2)^(1/2)-2*u^2+w^2+2*w+1)^(1/2)+2)^(1/2)+w-1

 (1)

with(MultiSeries)

[AddFunction, FunctionSupported, GetFunction, LeadingTerm, RemoveFunction, SeriesInfo, asympt, limit, multiseries, series, taylor]

 (2)

series(subs(w = 7/13, nu3), u = 0, 3)

series((-(-(2/13)*7^(1/2)*33^(1/2)+400/169)^(1/2)-6/13)+((1/13)*26^(1/2)*((2*7^(1/2)*33^(1/2)*(-(2/13)*7^(1/2)*33^(1/2)+400/169)^(1/2)+2*7^(1/2)*33^(1/2)-27*(-(2/13)*7^(1/2)*33^(1/2)+400/169)^(1/2)-27)/(-(2/13)*7^(1/2)*33^(1/2)+400/169)^(1/2))^(1/2))*u-((1/13)*(2*7^(1/2)*33^(1/2)-27)/(-(2/13)*7^(1/2)*33^(1/2)+400/169)^(1/2))*u^2+O(u^3),u,3)

 (3)

series(subs(w = 1, nu3), u = 0, 3)

series(-(-2*3^(1/2)+4)^(1/2)+(2^(1/2)*((2*3^(1/2)*(-2*3^(1/2)+4)^(1/2)+2*3^(1/2)-3*(-2*3^(1/2)+4)^(1/2)-3)/(-2*3^(1/2)+4)^(1/2))^(1/2))*u-((-3+2*3^(1/2))/(-2*3^(1/2)+4)^(1/2))*u^2+O(u^3),u,3)

 (4)

series(subs(w = zeta(3), nu3), u = 0, 2)

series((-(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2)^(1/2)+zeta(3)-1)-((2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2)^(1/2)*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)-2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2)^(1/2)*zeta(3)+2*zeta(3)*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)-(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2)^(1/2)+(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2*zeta(3)+1)/((-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+2*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)+2)^(1/2)*(-2*zeta(3)^(1/2)*(zeta(3)+2)^(1/2)+zeta(3)^2+2*zeta(3)+1)^(1/2)))*u^2+O(u^4),u,4)

 (5)

``

``


 

Download bug-in-MultiSeries.mw

How to simplify formula like e^(ln(2*Pi)...

ablmf 25 Maple
simplification exp convert
December 08 2017
3 0

I have tried

assume(n >= 0, n::integer);

But still,

simplify(exp(-n*ln(2*Pi)))

does nothing.

I want to convert it to (2*Pi)^(-n).

How to stop maple from simplify an integ...

ablmf 25 Maple
integral int inert
June 28 2017
2 1

If I just want to define an itegral and do not want maple to simplify it to a closed form, what should I do?

For example, I want define

s := int(exp(-x^2)*cos(2*x*y), x = 0 .. infinity).

Maple automatically simplify s to

(1/2)*sqrt(Pi)/exp(y^2)

But I want to keep s in integral form.

How to define a function with result of ...

ablmf 25
simplify
March 07 2017
1 0


 

I need a function like this:

 

    f := proc (n) options operator, arrow; sum(i, i = 1 .. n) end proc

proc (n) options operator, arrow; sum(i, i = 1 .. n) end proc

 (1)

But I know that the sum actually simplifes.

simplify(sum(i, i = 1 .. n))

(1/2)*n^2+(1/2)*n

 (2)

So I actually want define my function as

f := proc (n) options operator, arrow; (1/2)*n^2+(1/2)*n end proc

proc (n) options operator, arrow; (1/2)*n^2+(1/2)*n end proc

 (3)

Is there any way to do this with programming instead of copying the result of simplify manually?


 

Download define-functions.mw

Page 1 of 1