Questions - MaplePrimes

Error, (in …) integer overflow detected in compile...

Asked by: sursumCorda 1174

December 28 2024

1 14

Recently, I read an old article about how to call external subroutines to speed up the evaluation tremendously. The core of that article is:

Evidently, the is equivalent to when b＞0, and the `∸`(a, b) is equivalent to .
The 𝙲 code in that article was composed separately. However, now that it is possible to generate external functions, compile and access them using a single Compiler:-Compile command, it is better to complete the entire process on the fly (in order to reduce oversights and typographical errors while coding). (Besides, the original 𝙼𝚊𝚙𝚕𝚎 implementation appears verbose and less elegant.) So I write:

JonesM1 := Compiler:-Compile((n::posint) -> add(max(1 - max(add(ifelse(j = 0, max(j - 1, 0)!**2, irem(max(j - 1, 0)!**2, j)), j = 0 .. i) - n, 0), 0), i = 0 .. n**2), 'O' = 2): # literal translation of the above formula

But unfortunately, I then get:

JonesM1(1);
                               2

JonesM1(2);
                               3

JonesM1(3);
                               5

JonesM1(4);
Error, (in JonesM1) integer overflow detected in compiled code

How do I get rid of this error message?

Note. Of course there exist a ready-to-use built-in function , but the chief aim is not to find the n-th prime number as quickly as possible; I want to see if 𝙼𝚊𝚙𝚕𝚎 can directly create a compiled version of the one-liner that implements the so-called Jones' algorithm instead of writing separate 𝙲 code like that article. (Don't confound the means with the ends.)

Bug report: Error, (in dsolve) invalid subscript s...

Asked by: nm 10431

December 27 2024

0 2

I do not remember if I reported this before or not. Can't find it. Just in case, I am posting this.

If someone find it is duplicate, feel free to delete this. But this is in latest Maple 2024.2. May be this can be fixed in time by Maple 2025 version.

>

restart;

>	interface(version);

>	Physics:-Version();

>

libname;

$"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"$

>	ode := diff(y(x),x)/y(x)-(3(4x^2+y(x)^2+1))/(2x(4x^2+y(x)^2-2-2x))=0;

(diff(y(x), x))/y(x)-(3/2)*(4*x^2+y(x)^2+1)/(x*(4*x^2+y(x)^2-2-2*x)) = 0

>	DEtools:-odeadvisor(ode);

>	dsolve(ode,y(x));

Error, (in dsolve) invalid subscript selector

>

restart;

>	infolevel[dsolve]:=5;

>	ode := diff(y(x),x)/y(x)-(3(4x^2+y(x)^2+1))/(2x(4x^2+y(x)^2-2-2x))=0:

>	dsolve(ode,y(x));

Methods for first order ODEs:

--- Trying classification methods ---

trying a quadrature

trying 1st order linear

trying Bernoulli

trying separable

trying inverse linear

trying homogeneous types:

trying Chini

differential order: 1; looking for linear symmetries

trying exact

Looking for potential symmetries

trying inverse_Riccati

trying an equivalence to an Abel ODE

equivalence obtained to this Abel ODE: diff(y(x),x) = 3/2*(4*x^2+1)/x/(2*x^2-x-1)*y(x)-(x^2+2*x+3)/x/(2*x^2-x-1)^2*y(x)^2+3/8*(2*x+3)/(2*x^2-x-1)^3/x*y(x)^3

trying to solve the Abel ODE ...

The relative invariant s3 is: -1/432*(8*x^4+40*x^3+45*x^2-270*x+135)/x^3/(x-1)^6/(2*x+1)^4

The first absolute invariant s5^3/s3^5 is: 729/16*(128*x^8+1152*x^7+3696*x^6+1744*x^5+8148*x^4-31500*x^3+6615*x^2-5670*x+8505)^3/(2*x+1)^4/(8*x^4+40*x^3+45*x^2-270*x+135)^5

The second absolute invariant s3*s7/s5^2 is: 1/3*(8*x^4+40*x^3+45*x^2-270*x+135)*(10240*x^12+133120*x^11+697600*x^10+1710080*x^9+3358592*x^8-1701568*x^7+6692592*x^6-18182448*x^5+2088072*x^4-7938000*x^3+2525985*x^2+1786050*x+2679075)/(128*x^8+1152*x^7+3696*x^6+1744*x^5+8148*x^4-31500*x^3+6615*x^2-5670*x+8505)^2

...checking Abel class AIL (45)

...checking Abel class AIL (310)

...checking Abel class AIR (36)

...checking Abel class AIL (301)

...checking Abel class AIL (1000)

...checking Abel class AIL (42)

...checking Abel class AIL (185)

...checking Abel class AIA (by Halphen)

...checking Abel class AIL (205)

...checking Abel class AIA (147)

...checking Abel class AIL (581)

...checking Abel class AIL (200)

...checking Abel class AIL (257)

...checking Abel class AIL (400)

...checking Abel class AIA (515)

...checking Abel class AIR (1001)

...checking Abel class AIA (201)

...checking Abel class AIA (815)

Looking for potential symmetries

... changing x -> 1/x, trying again

Looking for potential symmetries

The third absolute invariant s5*s7/s3^4 is: 243/16*(10240*x^12+133120*x^11+697600*x^10+1710080*x^9+3358592*x^8-1701568*x^7+6692592*x^6-18182448*x^5+2088072*x^4-7938000*x^3+2525985*x^2+1786050*x+2679075)/(2*x+1)^4*(128*x^8+1152*x^7+3696*x^6+1744*x^5+8148*x^4-31500*x^3+6615*x^2-5670*x+8505)/(8*x^4+40*x^3+45*x^2-270*x+135)^4

-> ======================================

-> ...checking Abel class D (by Appell)

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {4/27} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class B (by Liouville)

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {1, 4, 1/4} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class A (by Abel)

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {0, -1/4} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class C (by Abel)

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {2, -11676447873119/75975070592769, 9/5, 15632211369872/75439744512117, 46273613050865/52325357771027, 75312059745574/25138886548531} ***

-> Step 3: looking for a solution F depending on x

_____________________________

C = 9/5 leads to a useless solution (F does not depend on x)

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIL 1.6

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {-4, 16} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIL 1.8

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {0, -116457391291688/45108305127449, -96869842492381/35485755507516, -36964550865207/94238117721032, -32286830321303/11596568583712, 32286830321303/11596568583712, 36964550865207/94238117721032, 96869842492381/35485755507516, 116457391291688/45108305127449} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIL 1.9

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {-2/9, -1/9} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIA 1.51

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {0, -94917840318055/84247876515289, -85939756880989/51399391393709, -82210125508529/36853933366676, -74381886667083/82545981233858, -41168492684238/33804146399567, -15658703496425/19275443365317, -9175348901453/101481647952193, 3/4, 15/4, 5568553686203/113599855351490, 12774469621703/63437040534358, 17836021821409/102823494563886, 39657708622139/74009717243016, 82495450887526/27663991325651, 86656182727564/45157560524183, 90074893410229/54954593917906, 100200889070747/32282555481919, 113612565327585/103754255779069} ***

-> Step 3: looking for a solution F depending on x

_____________________________

C = 15/4 leads to a useless solution (F does not depend on x)

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIA 1.5

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {-1, 1, -113553630998996/78694251194667, -112790344818825/35834119404842, -104905620984375/18860524785743, -95409943222181/78810323073434, -77648002983645/31218435062578, -67259194033608/9576982470445, -46892223838816/86694928762723, -45901561561111/29768419326991, -34674701564566/6522678435631, 26154715634141/21099761863911, 42841215778132/81925179545457, 52638927823233/15127919203723, 54069389554571/5444364811188, 54445812264368/10328928623117, 56815569067370/40738034746481, 75614540760757/62881656939350, 76459718737483/64786816765621, 85896394925571/88677987470966, 90623073438172/24246571690325, 103628692054633/17857341616628, 117754725919014/60191028908095} ***

-> Step 3: looking for a solution F depending on x

_____________________________

C = -1 leads to a useless solution (F does not depend on x)

_____________________________

C = -113553630998996/78694251194667 leads to a useless solution (F does not depend on x)

_____________________________

C = -112790344818825/35834119404842 leads to a useless solution (F does not depend on x)

_____________________________

C = -104905620984375/18860524785743 leads to a useless solution (F does not depend on x)

_____________________________

C = -95409943222181/78810323073434 leads to a useless solution (F does not depend on x)

_____________________________

C = -77648002983645/31218435062578 leads to a useless solution (F does not depend on x)

_____________________________

C = -67259194033608/9576982470445 leads to a useless solution (F does not depend on x)

_____________________________

C = -46892223838816/86694928762723 leads to a useless solution (F does not depend on x)

_____________________________

C = -45901561561111/29768419326991 leads to a useless solution (F does not depend on x)

_____________________________

C = -34674701564566/6522678435631 leads to a useless solution (F does not depend on x)

_____________________________

C = 26154715634141/21099761863911 leads to a useless solution (F does not depend on x)

_____________________________

C = 42841215778132/81925179545457 leads to a useless solution (F does not depend on x)

_____________________________

C = 52638927823233/15127919203723 leads to a useless solution (F does not depend on x)

_____________________________

C = 54069389554571/5444364811188 leads to a useless solution (F does not depend on x)

_____________________________

C = 54445812264368/10328928623117 leads to a useless solution (F does not depend on x)

_____________________________

C = 56815569067370/40738034746481 leads to a useless solution (F does not depend on x)

_____________________________

C = 75614540760757/62881656939350 leads to a useless solution (F does not depend on x)

_____________________________

C = 76459718737483/64786816765621 leads to a useless solution (F does not depend on x)

_____________________________

C = 85896394925571/88677987470966 leads to a useless solution (F does not depend on x)

_____________________________

C = 90623073438172/24246571690325 leads to a useless solution (F does not depend on x)

_____________________________

C = 103628692054633/17857341616628 leads to a useless solution (F does not depend on x)

_____________________________

C = 117754725919014/60191028908095 leads to a useless solution (F does not depend on x)

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIA 1.52

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {-5, -4, -3, 0, 1, 2, -3/2} ***

-> Step 3: looking for a solution F depending on x

*** No solution F of x was found ***

-> ======================================

-> ...checking Abel class AIA 1.53

-> Step 1: checking for a disqualifying factor on F after evaluating x at a number

Trying x = 2

*** No disqualifying factor on F was found ***

-> Step 2: calculating resultants to eliminate F and get candidates for C

*** Candidates for C are {-3, -1, 1, 2, -3/2, -2/3, -1/2} ***

-> Step 3: looking for a solution F depending on x

_____________________________

C = -3 leads to a useless solution (F does not depend on x)

_____________________________

C = -3/2 leads to a useless solution (F does not depend on x)

*** No solution F of x was found ***

trying to map the Abel into a solvable 2nd order ODE

...checking Abel class AIA 2-parameter, reducible to Riccati

Error, (in dsolve) invalid subscript selector

>

restart;

>	ode := diff(y(x),x)/y(x)-(3(4x^2+y(x)^2+1))/(2x(4x^2+y(x)^2-2-2x))=0:

>	dsolve(ode,y(x));

Error, (in dsolve) invalid subscript selector

>	tracelast;

Error, (in dsolve) invalid subscript selector

Download dsolve_invalid_subscript_dec_27_2024.mw

How to Arrange Multiple Plots in a Grid Layout...

Asked by: Paras31 170

December 27 2024

1 4

I am working on visualizing the results of an SIR model in Maple. I’ve created multiple plots for different parameter values using dsolve and odeplot, and I’d like to arrange them in a grid (e.g., 2 rows x 3 columns) for easier comparison. Here’s what I’ve tried:

Download grid.mw

Unfortunately, the insequence=true and gird=[2, 4] options didn’t work as I expected. Instead of arranging the plots in a grid, they appear like an animation

How can I correctly arrange multiple plots into a grid layout in Maple? Are there specific options or steps I’m missing? I would appreciate any guidance or examples.

Trouble Using Export Function in mpl Script...

Asked by: JoyDivisionMan 15

December 27 2024

0 2

I am doing some work in Maple. Specifically, I run the following commands in Maple worksheet mode:

with(Statistics);
L := abs(RandomVariable(UniformDistribution(0, 1)) - RandomVariable(UniformDistribution(0, 1)));
Export("LineDensity.jpg", DensityPlot(L));
PDF(L, x);

Everything works fine. When I export these commands into an MPL script, I get the following:

with(Statistics);
L := abs(RandomVariable(UniformDistribution(0, 1)) - RandomVariable(UniformDistribution(0, 1)));
Export("LineDensity.jpg", DensityPlot(L));
PDF(L, x);
NULL;

When I execute this script, I get the following error message:

Warning: persistent store makes readlib obsolete
Error: (in Export) exported file LineDensity.jpg could not be created

I need to do the work in MPL script form so that I can have larger problems processed. Does anyone know why I would be having this problem with the Export function?

why taking double derivative with two different in...

Asked by: salim-barzani 210

December 27 2024

0 6

I did a lot of time but this time i don't know why not run any one have idea?

>

(1)

>

(2)

>

(3)

>

(4)

>

(5)

>

Download why.mw

angle triangle calculator...

Asked by: yangtheary 265

December 27 2024

1 5

How factoring two term after substituting solutio...

Asked by: salim-barzani 210

December 27 2024

0 8

I already get the same results, but there's something about the factoring process that I encountered for the first time in this ODE. In the paper, it says that G′ satisfies a certain condition, but I’m not sure exactly what that means. Did the author use it for substitution, or did they change (m+F) into another variable and then solve? I’m not exactly sure what approach was taken. Does anyone have any idea or insight into this?

>

(1)

>

(2)

>

{t = tau, x = (-Zeta*c[3]-tau*c[4]-`Υ`*c[2]+xi)/c[1], y = `Υ`, z = Zeta, u(x, y, z, t) = U(xi)}

(3)

>

(4)

>

(5)

>

(6)

>

(7)

>

F := sum(a[i]*(m+1/(diff(G(xi), xi)))^i, i = -1 .. 1)

a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

(8)

>

a[-1]*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^2

(9)

>

(10)

>

a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

(11)

>

(12)

>

(13)

>

(14)

>

(15)

>

U(xi) = a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

(16)

>

diff(U(xi), xi) = a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

(17)

>

(18)

>

(19)

>

(20)

>

"collect(L,(m+1/(diff(G(xi),xi))))^( )"

Error, (in collect) cannot collect m+1/diff(G(xi),xi)

Download factoring.mw

Error, (in dsolve) type `System` does not exist wh...

Asked by: nm 10431

December 26 2024

1 4

I am trying to see if I can get speed up by using dsolve inside thread.

I made very simple example of global list of two differential equations to start with.

Next, created two threads where each picks one ode from the global list to process. So they should in theory run in parallel. The list of ode's is a global list in the worksheet for now.

But I keep getting error when calling dsolve

Error, (in dsolve) type `System` does not exist

I tried also passing the actual ode to the thread, still, same error.

Next, I did not pass anything, but called dsolve directly from inside thread proc on same ode. The ode is local variable inside the proc. I still get same error.

Does this mean dsolve is not supported by threads?

But when I searched this subject, AI says it works in threads:

Everything works OK when I run dsolve in worksheet outside thread (i.e. normally).

I will show below worksheet showing these cases. I must not be doing something right. But what? Can one not pass data from global worksheet to the thread this way? Or does one needs to load something in each thread to make this work?

>	interface(version);

>	Physics:-Version();

>

libname;

$"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"$

Example 1. Passing index of list to thread

>

restart;

>	g_list:=[sin(t)diff(x(t),t$2)+cos(t)diff(x(t),t)+2x(t)=0, diff(y(x),x)=lambdasin(lambdax)y(x)^2+acos(lambdax)^ny(x)-acos(lambda*x)^(n-1)]:

>	work_func:=proc(i::posint) :-dsolve(g_list[i]): end proc:

>	Threads:-Wait( seq( Threads:-Create( work_func(i)), i=1..2) );

Error, (in dsolve) type `System` does not exist

Example 2. Passing actual ode itself to thread

>

restart;

>	g_list:=[sin(t)diff(x(t),t$2)+cos(t)diff(x(t),t)+2x(t)=0, diff(y(x),x)=lambdasin(lambdax)y(x)^2+acos(lambdax)^ny(x)-acos(lambda*x)^(n-1)]:

>	work_func:=proc(ode::`=`) :-dsolve(ode): end proc:

>	Threads:-Wait( seq( Threads:-Create( work_func(g_list[i])), i=1..2) );

Error, (in dsolve) type `System` does not exist

>

Example 3. Normal processing. No threads

>

restart;

>	g_list:=[sin(t)diff(x(t),t$2)+cos(t)diff(x(t),t)+2x(t)=0, diff(y(x),x)=lambdasin(lambdax)y(x)^2+acos(lambdax)^ny(x)-acos(lambda*x)^(n-1)]:

>	work_func:=proc(ode::`=`) :-dsolve(ode): end proc:

>	for item in g_list do work_func(item); od:

>

#no error

>

Example 4. do not pass anything. Just call dsolve

>

restart;

>	work_func:=proc(i::posint) local x,t; local ode:=sin(t)diff(x(t),t$2)+cos(t)diff(x(t),t)+2*x(t)=0; :-dsolve(ode): end proc:

>	Threads:-Wait( seq( Threads:-Create( work_func(i)), i=1..2) );

Error, (in dsolve) type `System` does not exist

>

Download error_dsolve_using_threads_dec_26_2024.mw

implicitplot misbehaves in Maple 2024...

Asked by: Rouben Rostamian 8159

December 26 2024

3 5

>

restart;

Here are the graphs of a parabola and a straight line:

>	plots:-display( plot(x^2, x=-1..1), plot((x+1)/2, x=-1..1), color=["Red","Green"]);

Suppose I want to plot the part of the parabola that lies below

the straight line, and suppose, just to be nasty, I choose to do it

with implicitplot:

>	plots:-implicitplot(y=x^2, x=-1..1, y=0..(x+1)/2);

That is not a parabola at all. [And where does the "ynew" label come from?]

This behavior was introduced in Maple 2022.

In Maple 2021 we get the expected result:

plots:-implicitplot(y=x^2, x=-1..1, y=0..(x+1)/2);

Download mw.mw

Parallel processing routines...

Asked by: segfault 130

December 26 2024

1 5

It is not clear to me how to set up parallel processing for let's say 10 completely independent tasks/processes process1,process2 .....process10

Further, no communication between the 10 tasks/processes, process1,process2 .....process10, takes place and they only use previously calculated functions and variables in the worksheet.

So how do you set up the 10 processes to run in parallel using proc or grid. ??

I cannot make out from the examples in the manual how and where the processes are run.

Convert Vector to diagonal Matrix...

Asked by: Ronan 1311

December 26 2024

1 1

I am looking for a more eligent way to convert a Vector to a Diagonal Matrix.

>

restart

>

>	with(LinearAlgebra):

>

>	V:=Vector[column](3, [0.5863730366, 0.1171249270, 0.2965020364])

(1)

>	Vm:=Matrix(3,[[V[1],0,0],[0,V[2],0],[0,0,V[3]]])

(2)

>	Vm1:=Matrix(3,3):

>	for i to 3 do Vm1[i,i]:=V[i]; end do:

>

Vm1

(3)

Download 2024-12-26_Q_Diagonal_Matrix_from_Vector.mw

Bug report: dsolve generates internal error Error,...

Asked by: nm 10431

December 25 2024

1 4

This is an ode from textbook. dsolve gives new error I have not seen before.

Maple 2024.2 on windows 10.

>	interface(version);

>	Physics:-Version();

>

restart;

>

libname;

$"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"$

>	ode:=y(x)^2diff(y(x),x$3)-(3y(x)diff(y(x),x)+2xy(x)^2 )diff(y(x),x$2)+(2diff(y(x),x)^2+2xy(x)diff(y(x),x)+3x^2y(x)^2)diff(y(x),x)+x^3y(x)^3=0;

y(x)^2*(diff(diff(diff(y(x), x), x), x))-(3*y(x)*(diff(y(x), x))+2*x*y(x)^2)*(diff(diff(y(x), x), x))+(2*(diff(y(x), x))^2+2*x*y(x)*(diff(y(x), x))+3*x^2*y(x)^2)*(diff(y(x), x))+x^3*y(x)^3 = 0

>	dsolve(ode);

Error, (in factor) too many levels of recursion

>	infolevel[dsolve]:=5;

>	dsolve(ode);

Methods for third order ODEs:

--- Trying classification methods ---

trying 3rd order ODE linearizable_by_differentiation

differential order: 3; trying a linearization to 4th order

trying differential order: 3; missing variables

trying differential order: 3; exact nonlinear

trying 3rd order, integrating factor of the form mu(y'') for some mu

Trying the formal computation of integrating factors depending on any 2 of [x, y, y', y'']

*** Sublevel 2 ***

Methods for first order ODEs:

--- Trying classification methods ---

trying a quadrature

trying 1st order linear

<- 1st order linear successful

Successful computation of 3 integrating factors: [x*exp(-1/2*x^2*(I*2^(1/2)+1))*KummerM(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2)/y(x), x*exp(-1/2*x^2*(I*2^(1/2)+1))*KummerU(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2)/y(x), x*exp(-1/2*x^2*(I*2^(1/2)+1))*(Int(x*KummerM(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2)*exp(-1/2*x^2*(I*2^(1/2)-1)),x)*KummerU(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2)-Int(x*KummerU(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2)*exp(-1/2*x^2*(I*2^(1/2)-1)),x)*KummerM(3/4+1/8*I*2^(1/2),3/2,I*2^(1/2)*x^2))/y(x)]

Attempting computing related first integrals...

Error, (in factor) too many levels of recursion

Download dsolve_factor_dec_24_2024.mw

tracelast; gives long output with this at end

#(IntegrationTools:-Indefinite:-Polynomial,14): return poly/primitivepart*thisproc(primitivepart,var)
 IntegrationTools:-Indefinite:-Polynomial called with arguments: (8*I)*x*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*x1*2^(1/2)-(3*I)*2^(1/2)*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*y+(7*I)*2^(1/2)*KummerM(((1/8)*I)*2^(1/2)+7/4, 3/2, I*2^(1/2)*x^2)*y+12*x^2*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*y+8*x*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*x1+4*KummerM(((1/8)*I)*2^(1/2)+7/4, 3/2, I*2^(1/2)*x^2)*y, x1, nofactor = false
 #(IntegrationTools:-Indefinite:-Polynomial,8): newpoly := factor(poly)
Error, (in factor) too many levels of recursion
 locals defined as: p = p, primitivepart = primitivepart, base = base, exponent = exponent, subpolys = subpolys, change = change, newpoly = newpoly, u = u

Also, this error can not be cought using try/catch.

No examples for applying kernelopts -T memory al...

Asked by: segfault 130

December 24 2024

1 0

I cannot find any examples of how to raise the software memory limits to maximum with the mentioned -T switch in the manual to avoid memory allocation error popup.

Neither are there any examples posted on the web.

It cannot be a hardware allocation error, as the system memory is large.

So how is the -T switch implimented

How to Properly Display a Table of Numerical Resul...

Asked by: Paras31 170

December 24 2024

0 3

I'm working on solving a system of ODEs in Maple that models an epidemic scenario, tracking the number of susceptible, infected, and recovered individuals over time. Here's what I've done so far:
Download sir_model.mw

Created a table of results at daily intervals with the following code:

results := [seq([t = tval, s = round(evalf(sol(tval)[2][1])), i = round(evalf(sol(tval)[2][2])),r = round(evalf(sol(tval)[2][2]))], tval = 0 .. 50)];

printf("%-10s %-15s %-15s %-15s\n", "Day", "Infected", "Recovered");
printf("---------------------------------------------\n");
for entry in results do
    printf("%-10d %-15d %-15d %-15d\n", entry[t],entry[s], entry[i], entry[r]);
end do;

but I encountered an error (in fprintf) integer expected for integer format

Why in two platform for ODE equation i have two di...

Asked by: salim-barzani 210

December 24 2024

0 17

P1.mw

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post

MaplePrimes Questions

Error, (in …) integer overflow detected in compile...

Bug report: Error, (in dsolve) invalid subscript s...

How to Arrange Multiple Plots in a Grid Layout...

Trouble Using Export Function in mpl Script...

why taking double derivative with two different in...

angle triangle calculator...

How factoring two term after substituting solutio...

Error, (in dsolve) type `System` does not exist wh...

implicitplot misbehaves in Maple 2024...

Parallel processing routines...

Convert Vector to diagonal Matrix...

Bug report: dsolve generates internal error Error,...

No examples for applying kernelopts -T memory al...

How to Properly Display a Table of Numerical Resul...

Why in two platform for ODE equation i have two di...

Save this setting as your default sorting preference?

Ask a Question

Create a Post

Generating PDF…

Save this setting as your default sorting preference?
Note: You can change your preference any time in your account settings
Don't show this again