I do not know why int() on this integrand fails always first time, and works second time it is called. seems like something is not loaded correctly first time?

integrand:=(((-3*x^2-18*x-27)*exp(2)^2+(30*x^3+330*x^2+1170*x+1350)*exp(2)-75*x^4-1200*x^3-7050*x^2-18000*x-16875)*ln(x)+(12*x^2+54*x+81)*exp(2)^2+(-120*x^3-1106*x^2-3510*x-4050)*exp(2)+225*x^4+3560*x^3+20990*x^2+54000*x+50625)/((3*x^4+18*x^3+27*x^2)*exp(2)^2+(-30*x^5-330*x^4-1170*x^3-1350*x^2)*exp(2)+75*x^6+1200*x^5+7050*x^4+18000*x^3+16875*x^2):

print("First time");
int(integrand,x);

print("second time");
int(integrand,x);

Worksheet attached.

Update

Here is a movie. it is few minutes long. This happens by random and not each time. This movie shows the command starting from "restart" are repated 4 times. First two times, no error. Then the error shows up.  So it is random. Maybe it depends if Maple is busy with other things or not. I have each worksheet set to use its own server though. So I have no idea why this happens sometimes and not other times.

issue_int_nov_11_2021.mw

I think there is a problem here

restart;
the_integrand:=(((-2*x^2+2)*exp(exp(exp(3)))^2+(4*x^3-4*x)*exp(exp(exp(3)))-2*x^4-6*x^2+8)*exp(ln(x)-x^2)*ln(exp(exp(exp(3)))^2-2*x*exp(exp(exp(3)))+x^2+4)+((2*x^2-2)*exp(2)*exp(exp(exp(3)))^2+((-4*x^3+4*x)*exp(2)-2*x)*exp(exp(exp(3)))+(2*x^4+6*x^2-8)*exp(2)+2*x^2)*exp(ln(x)-x^2))/(exp(exp(exp(3)))^2-2*x*exp(exp(exp(3)))+x^2+4);

int(the_integrand,x,method=_RETURNVERBOSE)

#try MeijerG
int(the_integrand,x,method=MeijerG);

gives

I do not know if this known or not.

Maple 2021.1 on windows

This is Maple 2021.1 on windows.

===================
restart;
expr:=x^(6+1/3);
res:=series(expr,x=0,6);
==============

gives
        O(x^(19/3))

But type of the above is not series:

===========
type(res,'series');
        false
===========

Yet convert(res,polynom) works

================
convert(res,polynom)
          0
================

Which is correct conversion. But help says that

"convert/polynom
convert a series to polynomial form"

Notice, it says "series" there.

So the input must be type series. But Maple says
O(x^(19/3)) is not type series. I think this is wrong. The
type returned should be series. Now the type returned is 'function'
from the series command.

What Am I overlooking here?

I want to calculate the numerical value of the following integral at a given (r,z) where r and z are real numbers:

g(r,z)=int(exp(-sqrt(k^2 - k0^2)*z)*BesselJ(0, k*r)*k/sqrt(k^2 - k0^2), k = 0 .. infinity) where k0 is a real number.

In other words, it is the hankel transform g(r,z) of the function f(k):

f(k)=exp(-sqrt(k^2 - k0^2)*z)/sqrt(k^2 - k0^2)

The theoretical result is exp(-I*k0*R)/R with R=sqrt(r^2+z^2)

I tried to use the commande evalf(Int(exp(-sqrt(k^2 - k0^2)*z)*BesselJ(0, k*r)*k/sqrt(k^2 - k0^2), k = 0 .. infinity)) but it takes hours to evaluate and I did not get any results yet. I tried to calculate this integral with another software which uses a vectorized adaptive quadrature algorithm and I got the result in a few seconds.

How could I calculate this integral with Maple ?

Thanks.

Hello. Please help me. I need to calculate the integral (3). This integral has many singular points at which there is convergence in the sense of the principal Cauchy value. The Maple integral itself does not count. I don't understand how to find automatically all the singular points on the integration area. Then, perhaps, it would be possible to split the integral into the sum of integrals by regions, as I roughly wrote in the picture. I want to automate this process, because in fact it is necessary to calculate many integrals of the form (4), where f(x,y) are arbitrary functions that can oscillate strongly, so I don't want to write banal quadrature formulas. I would like to use the means of Maple, since the accuracy will be greater and faster, but we need to somehow bypass the special points. I will be glad of any help. Thank you very much

Download Integrate.mw

fyi;

Maple 2021.1 on windows 10.

Screen in worksheet displays this (correct)

But latex generated when compiled using lualatex from TEXLIVE 2021 shows this

code 

sol:=(Vector(2, [x(t),y(t)])) = (Vector(2, [8*c[1]*exp((1/2+1/2*89^(1/2))*t)/(-3+89^(1/2))-8*c[2]*exp((1/2-1/2*89^(1/2))*t)/(3+89^(1/2))+2/11*t^2-3/11*exp(t)-2/121*t+23/1331,c[1]*exp((1/2+1/2*89^(1/2))*t)+c[2]*exp((1/2-1/2*89^(1/2))*t)-1/11*t^2-15/22*exp(t)+12/121*t-17/1331])):

sol:=simplify(sol);

latex(sol)

`Standard Worksheet Interface, Maple 2021.1, Windows 10, May 19 2021 Build ID 1539851`

`The "Physics Updates" version in the MapleCloud is 1105 and is the same as the version installed in this computer, created 2021, November 8, 23:55 hours Pacific Time.`

zip file attached include the latex file and worksheet used.

issue_latex_nov_6_2021.zip

What happened to the graph output?

Just to update a bit to be more specific as what I did.

I have been using Maple since 2007 and never had such problems before 2012, when I did my PhD and switched to Mathematica.

1 - Right clicked on the graph and export to PDF - low quality (not seem like vector graphics)

2 - Print using Adobe PDF Printer from the system - low quality (not seem like vector graphics)

3 - EPS is fine, but WHHHHYYYYY????? Is this a bug??

[code] with(plots): z := x + I*y: inequal([ evalc( abs(z+1-I)    ) <= 1, evalc( argument(z-1) ) <= 3*Pi/4 ], x = -2.0..1.5, y = 0..2.2 , scaling = constrained, transparency = 0, numpoints = 500 ) [/code]

I want to calculate the Noether Symmetries from the Lagrangian. Is there any built-in Maple cammad to  where i just put the lagrangian and get the determining PDEs and calculate the Noether Symmetries?

Hello everyone!

I made an integration and it gave me a really long answer. The answer of the book is the same, as you can see, i "smartplotted" the two answers, and they are the same. But how do i get the nice short answer? I tried a lot of different things, but nothing gave me that short answer. 

I would really like to be able to get nice short answers, while it will really clean things up. 

Thank you!

Greetings,

The Function

Download Question_for_maple_primes_6.mw

I want to know how to remove an error message from the last line of using Maple.  I must have made some kind of error, but I want to move on but I can't get past the error message.  How can I remove the error message.  I have tried using "restart" but nothing happened.  I'm stuck until I clear this error message from the last line of my worksheet.  Can you help please?

Thanks

Dave

PS: as you can guess I'm pretty new to Maple.

Download numboccur.mwnumboccur.mw

Hello everbody! 

Ive got a question about my dutch math book, and why it thinks its answer is better than mine. 

It is question c. the integral that is covers the x value from 4 to 6 for f(x). 

I would think that means the integral of 1 to 6 - the integral of 1-4, and thus 4-7=-3. Now i dont know how they thought about getting a smaller area under the function with a larger domain of x, so i am not sure what is going on. Im not even sure if what they ask is possible. The book gave 3 as an answer. My answer was -3. Looking at it, and thinking about it. It may not even be possible. 

Any way, that is the question. 

Picture 2 shows the rules, picture 3 shows the answers of the book. 

Thank you!

Greetings,

The function

Download Question_for_maple_primes_5.mw

Let f(x) = x³+9x²-9x-1/x⁴+1

Find the equations of the 4 lines having slope 0.5 that are tangent to the graph of 𝑓(𝑥) (first find the values of 𝑥 where 𝑓 ′(𝑥) = 0.5).

Hello there,
I'm currently working on a school-related assignment, where I have to put two functions equal each other to find their crossing point. My issue is, that these two functions are waves, and continiue into infinite.
If I put them equal as this, my document crashes:
f(x) = g(x)

My question is then, if I can put some requirement to Maple, so it only calculates from significant X-intervals

elctroWav.mw have the equation 1 :

eq1 := 2*m*(E + e0^2/(rb*rho))*g(rho, t)*(-rb*rho + R)/(h^2*R) + R*(diff(g(rho, t), rho, rho)*(-rb*rho + R)/R - 2*diff(g(rho, t), rho)*rb/R)/(rb^3*rho) - diff(g(rho, t), t, t)*(-rb*rho + R)/(a^2*c^2*R) = 0;

with iv1:

iv1 := D[1](g)(0, 0) - g(0, 0)*rb/R = rb*R, g(1, 0)*(R - rb)/R = 0, D[2](g)(rb, 0) = a*c;

As in maple document electroWav.mw

The issue is that there is an unsolved part with RootOf how can I rezolve it?

-R*h^2*RootOf(AiryBi(-(R*_Z^3*h^2 + 2*e0^2*m*rb^2)/(_Z^2*R*h^2)))^3

Please advice...