Maple Questions and Posts

modified Bessel function of the third kind with in...

Asked by: sand15 1067

15 hours ago

1 6

I'm not sure of the Maple function I have to use to represent the modified Bessel function of the third kind with index 1.

The Bessel function of the third kind is also named the Hankel function.

I mistakenly started to use HankelH1 or HankelH2 but quickly realized they are both complex-valued functions while the paper I'm working on uses a real-valued modified Bessel function of the third kind: finally, does this latter exist in Maple 2015?

Thanks in advance

Why dsolve do not resolve constant of integration ...

Asked by: nm 11483

October 04 2025

0 3

Given an ode with IC. When solution is explicit, Maple resolves the constant of integration as expected and returns solution with no c__1 in it.

But when asked for implicit solution, also with same IC, it now returns solution with c__1 still there.

Is this by design or a bug? Should not constant of integration be resolved using IC in both cases? If unable to solve for c__1 because solution is implicit, should it then not return solution all?

Does this happen in earlier versions of Maple?

>

restart;

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

>	IC:=D(y)(0)=0; ode:=diff(y(x),x)^2+x*y(x)=0: sol:=dsolve([ode,IC]);

>	sol:=dsolve(ode,'implicit');

y(x) = 0, -x^2/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, x^2/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))-c__1 = 0

>	#WHY did not resolve constant of integration here?? sol:=dsolve([ode,IC],'implicit');

x^2/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(-x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, -x^2/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))+3*(-x*y(x))^(1/2)/((x^3+9*y(x))*(x^2+3*(-x*y(x))^(1/2)))-c__1 = 0, y(x) = 0

>

Download why_C_still_in_solution_maple_2025_1_oct_4_2025.mw

why PDEtools:-Solve crashes Maple server?...

Asked by: nm 11483

October 03 2025

0 0

For some reason, solve() works Ok, but PDEtools:-Solve causes kernel crash.

Attached worksheet.

Anyone knows why and can someone try to reproduce this? I am using Linux CaucyOS (arch linux based) and Maple 2025.1

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

>

eq:=[-3/2*_C6^2+4*(60588*_C6^3+6659802*_C6^2-7113718800*_C6+168667884119)/(40392*_C6-971364)-1/8*(-6*_C6+1)^2-4*(1406664*_C6^2-1186536946*_C6+27728821045)/(26928*_C6-647576)+5*_C6-125/6+9*(-1/24*(-6*_C6+16)^2-4*(1406664*_C6^2-1186536946*_C6+27728821045)/(26928*_C6-647576)+(60588*_C6^3+6659802*_C6^2-7113718800*_C6+168667884119)/(40392*_C6-971364))*(-49/2+_C6) = 4, -1/24*(-6*_C6+49)^2+191*(1406664*_C6^2-1186536946*_C6+27728821045)/(26928*_C6-647576)-(60588*_C6^3+6659802*_C6^2-7113718800*_C6+168667884119)/(40392*_C6-971364)+1/12*(-6*_C6+9801)^2 = 0];
Cs:={_C6};

${_C6}$

>	solve(eq,Cs);

${_C6 = _C6}$

>	PDEtools:-Solve(eq,Cs);

Download maple_crash_Solve_oct_3_2025.mw

Also reported to Maplesoft support.

Physics:-Setup(assumingusesAssume = true): causes ...

Asked by: nm 11483

October 02 2025

0 2

I found very strange behaviour of Maple 2025.1 on Linux.

Same exact code. Calling timelimit(sol,ode) twice on two different solutions. If I do not add

Physics:-Setup(assumingusesAssume = true)

At the start, then both timelimits finish OK. But once Physics:-Setup(assumingusesAssume = true) is added at the start, the second timelimit hangs.

I waited 2 hrs and nothing happens. Maple just freezes. Can't even stop the server from worksheet by clicking on the red button at lower left corner.

This is using latest Physics. Does anyone know why this happens? It seems due to some memory cache issue?

Make sure to save all your work before trying this just in case you have to kill Maple application.

The strange thing, unable to stop the server by clicking on red button or clicking on RESTART KERNEL icon at top, or even clicking on the debuger icon at lower left corner. Only way was to kill Maple itself from Linux command line.

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

Example 1. Not using Physics:-Setup(assumingusesAssume = true): gives NO hang

>

sol_1:=y(x) = 1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2
+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-\
8*3^(1/2))^(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^3*3^(1/2)+
9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*
3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)+36*I*_C2+36*I*x+12*_C2*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(I*3^(1/2)-1):
ode:=diff(y(x),x)-y(x)^3 = 8:
timelimit(30,odetest(sol_1,ode));

Error, (in simplify/ln/relations) time expired

>

sol_2:=y(x) = -1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)
*sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^
2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2
-8*3^(1/2))^(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^3*3^(1/2)
+9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8
*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)-36*I*_C3-36*I*x+12*_C3*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(1+I*3^(1/2)):
timelimit(30,odetest(sol_2,ode));

Error, (in collect) time expired

>

Example 2. using Physics:-Setup(assumingusesAssume = true): second timelimit always hangs

>

restart;

>	Physics:-Setup(assumingusesAssume = true):

>

sol_1:=y(x) = 1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2
+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-\
8*3^(1/2))^(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^3*3^(1/2)+
9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*
3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C2*3^(1/2)+12*3^(1/2)*x+36*I*_C2+36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)+36*I*_C2+36*I*x+12*_C2*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(I*3^(1/2)-1):
ode:=diff(y(x),x)-y(x)^3 = 8:
timelimit(30,odetest(sol_1,ode));

Error, (in expand) time expired

>	#this below will now hang

>

sol_2:=y(x) = -1/2*(3*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)
*sec(_Z)^2-8*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^
2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2
-8*3^(1/2))^(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^3*3^(1/2)
+9*tan(RootOf(2*3^(1/2)*ln(2*3^(1/6)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8
*3^(1/2))^(1/3))-3^(1/2)*ln(4*3^(1/3)-2*3^(1/6)*(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*
sec(_Z)^2-8*3^(1/2))^(1/3)+(9*tan(_Z)*sec(_Z)^2+9*3^(1/2)*sec(_Z)^2-8*3^(1/2))^
(2/3))+12*_C3*3^(1/2)+12*3^(1/2)*x-36*I*_C3-36*I*x+6*_Z))^2+3*3^(1/2)*tan(
RootOf(2*3^(1/2)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(1/3)+2)-3^(1/2
)*ln((3*3^(1/2)*tan(_Z)*sec(_Z)^2-8+9*sec(_Z)^2)^(2/3)-2*(3*3^(1/2)*tan(_Z)*sec
(_Z)^2-8+9*sec(_Z)^2)^(1/3)+4)-36*I*_C3-36*I*x+12*_C3*3^(1/2)+12*3^(1/2)*x+6*_Z
))+1)^(1/3)*(1+I*3^(1/2)):
timelimit(30,odetest(sol_2,ode));

Download hangs_timelimit_with_physics_maple_2025_1_oct_2_2025.mw

Could someone try to reproduce this on their Maple? If so, I will send bug report.

why dsolve gives Error, (in dsolve) this is a bug ...

Asked by: nm 11483

October 02 2025

1 0

Never seen such a message before

Here is MWE, Trace shows it comes from (SolveTools:-PolynomialSystemSolvers:-PseudoResultant:-AttemptFactorization,2)

>

restart;

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>	ode:=diff(y(x),x$2)+3y(x)+D(y)(3)=0; IC:=y(1)+4y(2)+y(4)=0;

diff(diff(y(x), x), x)+3*y(x)+(D(y))(3) = 0

>	dsolve(ode)

x+(1/3)*3^(1/2)*(tan(3^(1/2)*x)-arctan(tan(3^(1/2)*x)))-((1/3)*(D(y(x)))(3)*(1+tan(3^(1/2)*x)^2)^(1/2)+(1+tan(3^(1/2)*x)^2)^(1/2)*y(x))*c__1-c__2 = 0

>	dsolve([ode,IC])

Error, (in dsolve) this is a bug

>

tracelast

dsolve called with arguments: [diff(diff(y(x), x), x)+3*y(x)+(D(y))(3) = 0, y(1)+4*y(2)+y(4) = 0], arbitraryconstants = subscripted, atomizenames = true, build = false, numeric = false, type = none
#(dsolve,80): error

\`dsolve/IC\` called with arguments: [diff(diff(y(x), x), x)+3*y(x)+(D(y))(3) = 0, y(1)+4*y(2)+y(4) = 0], {y(x)}, skipimplicit = false, skippparticularsolforlinearODEs = true, solution = {}, usesolutions = particular and general
#(\`dsolve/IC\`,53): return procname(_passed,':-usesolutions' = "general and particular")

\`dsolve/IC\` called with arguments: [diff(diff(y(x), x), x)+3*y(x)+(D(y))(3) = 0, y(1)+4*y(2)+y(4) = 0], {y(x)}, skipimplicit = false, skippparticularsolforlinearODEs = true, solution = {}, usesolutions = general and particular
#(\`dsolve/IC\`,57): ans := procname(_passed,':-usesolutions = "general"');

\`dsolve/IC\` called with arguments: [diff(diff(y(x), x), x)+3*y(x)+(D(y))(3) = 0, y(1)+4*y(2)+y(4) = 0], {y(x)}, skipimplicit = false, skippparticularsolforlinearODEs = true, solution = {}, usesolutions = general
#(\`dsolve/IC\`,277): zz := map(op,{\`dsolve/IC/_C\`({ANS[i]},funcs,x,ics)});

\`dsolve/IC/_C\` called with arguments: {y(x) = -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(3^(1/2)*x)^2)^(1/2)-3^(1/2)*tan(3^(1/2)*x)+3^(1/2)*arctan(tan(3^(1/2)*x))+3*_C[2]-3*x)/((1+tan(3^(1/2)*x)^2)^(1/2)*_C[1])}, {y(x)}, x, [y(1) = _C1, y(2) = _C2, y(4) = -_C1-4*_C2]
#(\`dsolve/IC/_C\`,1): ans := \`dsolve/IC/_C/do\`(solns,depvars,t,inits,'evaluated_ans', "default",':-giveup = giveup');

\`dsolve/IC/_C/do\` called with arguments: {y(x) = -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(3^(1/2)*x)^2)^(1/2)-3^(1/2)*tan(3^(1/2)*x)+3^(1/2)*arctan(tan(3^(1/2)*x))+3*_C[2]-3*x)/((1+tan(3^(1/2)*x)^2)^(1/2)*_C[1])}, {y(x)}, x, [y(1) = _C1, y(2) = _C2, y(4) = -_C1-4*_C2], evaluated_ans, default, giveup = giveup, usecansolve = false
#(\`dsolve/IC/_C/do\`,103): csol := [\`ODEtools/Solve/EnvDropMultiplicity\`(eqns,consts)];

\`ODEtools/Solve/EnvDropMultiplicity\` called with arguments: {-(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(3^(1/2))^2)^(1/2)-3^(1/2)*tan(3^(1/2))+3^(1/2)*(3^(1/2)-Pi)+3*_C[2]-3)/((1+tan(3^(1/2))^2)^(1/2)*_C[1]) = _C1, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(2*3^(1/2))^2)^(1/2)-3^(1/2)*tan(2*3^(1/2))+3^(1/2)*(2*3^(1/2)-Pi)+3*_C[2]-6)/((1+tan(2*3^(1/2))^2)^(1/2)*_C[1]) = _C2, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(4*3^(1/2))^2)^(1/2)-3^(1/2)*tan(4*3^(1/2))+3^(1/2)*(4*3^(1/2)-2*Pi)+3*_C[2]-12)/((1+tan(4*3^(1/2))^2)^(1/2)*_C[1]) = -_C1-4*_C2}, {_C1, _C2, _C[1], _C[2]}, keepalreadysolveduntouched = false, removelabel = false
#(\`ODEtools/Solve/EnvDropMultiplicity\`,29): sol := :-solve(ee,X,_rest)

solve called with arguments: {-(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(3^(1/2))^2)^(1/2)-3^(1/2)*tan(3^(1/2))+3^(1/2)*(3^(1/2)-Pi)+3*_C[2]-3)/((1+tan(3^(1/2))^2)^(1/2)*_C[1]) = _C1, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(2*3^(1/2))^2)^(1/2)-3^(1/2)*tan(2*3^(1/2))+3^(1/2)*(2*3^(1/2)-Pi)+3*_C[2]-6)/((1+tan(2*3^(1/2))^2)^(1/2)*_C[1]) = _C2, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(4*3^(1/2))^2)^(1/2)-3^(1/2)*tan(4*3^(1/2))+3^(1/2)*(4*3^(1/2)-2*Pi)+3*_C[2]-12)/((1+tan(4*3^(1/2))^2)^(1/2)*_C[1]) = -_C1-4*_C2}, {_C1, _C2, _C[1], _C[2]}, AllSolutions = FAIL, ConditionalSolutions = FAIL, DropMultiplicity = FAIL, Explicit = FAIL, MaxSols = FAIL, SolveOverReals = FAIL, SymbolicSolutions = true, TryHard = FAIL, UseAssumptions = false, domain = default, parameters = {}, parametric = false, split = false, useunits = FAIL
#(solve,133): _MaxSols := oldmaxsols

Engine:-Main called with arguments: {-(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(3^(1/2))^2)^(1/2)-3^(1/2)*tan(3^(1/2))+3^(1/2)*(3^(1/2)-Pi)+3*_C[2]-3)/((1+tan(3^(1/2))^2)^(1/2)*_C[1]) = _C1, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(2*3^(1/2))^2)^(1/2)-3^(1/2)*tan(2*3^(1/2))+3^(1/2)*(2*3^(1/2)-Pi)+3*_C[2]-6)/((1+tan(2*3^(1/2))^2)^(1/2)*_C[1]) = _C2, -(1/3)*(_C[1]*(D(_Z))(3)*(1+tan(4*3^(1/2))^2)^(1/2)-3^(1/2)*tan(4*3^(1/2))+3^(1/2)*(4*3^(1/2)-2*Pi)+3*_C[2]-12)/((1+tan(4*3^(1/2))^2)^(1/2)*_C[1]) = -_C1-4*_C2}, {}, {_C1, _C2, _C[1], _C[2]}
#(SolveTools:-Engine:-Main,37): sol := [SolveTools:-PolynomialSystem(eqns,vars,ineqs,':- maxsols' = \`if\`(type([_MaxSols],[{integer, infinity}]) and 0 < _MaxSols,_MaxSols,100))]

PolynomialSystem:-ModuleApply called with arguments: {-(1/3)*(_X000003*(D(_Z))(3)*(1+tan(3^(1/2))^2)^(1/2)-3^(1/2)*tan(3^(1/2))+3^(1/2)*(3^(1/2)-Pi)+3*_X000004-3)/((1+tan(3^(1/2))^2)^(1/2)*_X000003) = _X000001, -(1/3)*(_X000003*(D(_Z))(3)*(1+tan(2*3^(1/2))^2)^(1/2)-3^(1/2)*tan(2*3^(1/2))+3^(1/2)*(2*3^(1/2)-Pi)+3*_X000004-6)/((1+tan(2*3^(1/2))^2)^(1/2)*_X000003) = _X000002, -(1/3)*(_X000003*(D(_Z))(3)*(1+tan(4*3^(1/2))^2)^(1/2)-3^(1/2)*tan(4*3^(1/2))+3^(1/2)*(4*3^(1/2)-2*Pi)+3*_X000004-12)/((1+tan(4*3^(1/2))^2)^(1/2)*_X000003) = -_X000001-4*_X000002}, {_X000001, _X000002, _X000003, _X000004}, {}, FAIL, backsubstitute = true, domain = absolute, engine = default, explicit = false, maxsols = 100, preprocess = true, preservelabels = false, tryhard = true
#(SolveTools:-PolynomialSystem:-ModuleApply,15): sol := SolveTools:-PolynomialSystem:-Main(SolveTools:-Utilities:-New( equations,notz,SolveTools:-Utilities:-Intersect(vars,indets(equations ))),':-backsub' = reallybacksub,':-domain' = domain,':-maxsols' = realmaxsols,':-engine' = engine,':-preprocess' = preprocess,':- tryhard' = tryhard);

PolynomialSystem:-Main called with arguments: [{-(1/3)*(_X000003*(D(_Z))(3)*(1+tan(3^(1/2))^2)^(1/2)-3^(1/2)*tan(3^(1/2))+3^(1/2)*(3^(1/2)-Pi)+3*_X000004-3)/((1+tan(3^(1/2))^2)^(1/2)*_X000003)-_X000001, -(1/3)*(_X000003*(D(_Z))(3)*(1+tan(2*3^(1/2))^2)^(1/2)-3^(1/2)*tan(2*3^(1/2))+3^(1/2)*(2*3^(1/2)-Pi)+3*_X000004-6)/((1+tan(2*3^(1/2))^2)^(1/2)*_X000003)-_X000002, -(1/3)*(_X000003*(D(_Z))(3)*(1+tan(4*3^(1/2))^2)^(1/2)-3^(1/2)*tan(4*3^(1/2))+3^(1/2)*(4*3^(1/2)-2*Pi)+3*_X000004-12)/((1+tan(4*3^(1/2))^2)^(1/2)*_X000003)+_X000001+4*_X000002}, {3*(1+tan(3^(1/2))^2)^(1/2)*_X000003 <> 0, 3*(1+tan(2*3^(1/2))^2)^(1/2)*_X000003 <> 0, 3*(1+tan(4*3^(1/2))^2)^(1/2)*_X000003 <> 0}, {_X000001, _X000002, _X000003, _X000004}, {}, true, false, 1, {_X000001, _X000002, _X000003, _X000004}], backsub = true, domain = absolute, engine = default, maxsols = 100, preprocess = true, tryhard = true
#(SolveTools:-PolynomialSystem:-Main,86): sol := [SolveTools:-PolynomialSystemSolvers:- PseudoResultant(SolveTools:-Utilities:- GetEquations(sys),SolveTools:-Utilities:- GetVariables(sys),SolveTools:-Utilities:- GetInequations(sys),':-maxsols' = maxsols,':- tryhard' = tryhard)];

PseudoResultant:-ModuleApply called with arguments: {(_X000003*(D(_Z))(3)*(-2*cos(3^(1/2))+1+sec(3^(1/2)))+3*_X000001*_X000003*(-2*cos(3^(1/2))+sec(3^(1/2)))+3*_X000002*_X000003+3^(1/2)*tan(3^(1/2)))/cos(2*3^(1/2)), -(sec(4*3^(1/2))+sec(3^(1/2)))*_X000003*(D(_Z))(3)+3*(sec(4*3^(1/2))-sec(3^(1/2)))*_X000001*_X000003+12*_X000002*sec(4*3^(1/2))*_X000003+3^(1/2)*Pi+3^(1/2)*tan(4*3^(1/2))-3^(1/2)*tan(3^(1/2))}, {_X000001, _X000002, _X000003}, {_X000003 <> 0}, backsub = true, maxsols = 100, tryhard = true
#(SolveTools:-PolynomialSystemSolvers:-PseudoResultant:-ModuleApply,25): SolveTools:-PolynomialSystemSolvers:-PseudoResultant:-ApplySubstitution ({},[],numer(tm1),notzero union map(xx -> denom(xx) <> 0,tm1), unknowns);

ApplySubstitution called with arguments: {}, [], {_X000003*(D(_Z))(3)*sec(3^(1/2))+3*_X000001*sec(3^(1/2))*_X000003-2*_X000003*(D(_Z))(3)*cos(3^(1/2))-6*_X000001*_X000003*cos(3^(1/2))+3^(1/2)*tan(3^(1/2))+_X000003*(D(_Z))(3)+3*_X000002*_X000003, -_X000003*(D(_Z))(3)*sec(4*3^(1/2))+3*_X000001*sec(4*3^(1/2))*_X000003+12*_X000002*sec(4*3^(1/2))*_X000003+3^(1/2)*Pi+3^(1/2)*tan(4*3^(1/2))-_X000003*(D(_Z))(3)*sec(3^(1/2))-3*_X000001*sec(3^(1/2))*_X000003-3^(1/2)*tan(3^(1/2))}, {1 <> 0, _X000003 <> 0, cos(2*3^(1/2)) <> 0}, {_X000001, _X000002, _X000003}
#(SolveTools:-PolynomialSystemSolvers:-PseudoResultant:-ApplySubstitution,48): eqns := map(SolveTools:-PolynomialSystemSolvers:-PseudoResultant:- AttemptFactorization,eqns,not0,unknowns);

AttemptFactorization called with arguments: _X000003*(D(_Z))(3)*sec(3^(1/2))+3*_X000001*sec(3^(1/2))*_X000003-2*_X000003*(D(_Z))(3)*cos(3^(1/2))-6*_X000001*_X000003*cos(3^(1/2))+3^(1/2)*tan(3^(1/2))+_X000003*(D(_Z))(3)+3*_X000002*_X000003, {_X000003 <> 0, cos(2*3^(1/2)) <> 0}, {_X000001, _X000002, _X000003}
#(SolveTools:-PolynomialSystemSolvers:-PseudoResultant:-AttemptFactorization,2): error \`this is a bug\`

Error, (in dsolve) this is a bug

locals defined as: e = e

>

Download dsolve_error_this_is_bug_2025_1_oct_2_2025.mw

Note that changing the IC to IC:=y(1)+4*y(2)=0; makes it work with no error.

It is only when adding y(4) does the error shows up.

update Bug report emailed to Maplesoft.

Installation Maple and Mathematica...

Asked by: Hajra Abid 5

October 01 2025

0 0

Hi! how I install Maple 2025 and Mathematica in my laptop?
Provide me Complete Set-up.

Incorrect calculation of limits...

Asked by: mthkvv 180

September 27 2025

1 3

Hi there.

There is an issue with following limits in Maple 2025:

Maple returned an incorrect zero answers. Plotting functions we see that non-zero limits do exist:

Thank you.

limits.mw

Enter key to display, without new line?...

Asked by: GFY 65

September 25 2025

0 1

I would like to keep writing commands on the same line without automatically creating a new line break. Is there a way to achieve this in Maple?

Changing column datatype: Excel Import...

Asked by: Motocycle Mayor 15

September 23 2025

2 9

I want to change the datatype of a column of a Dataframe imported from Excel.

The import succeeds, and I can change the data type of numeric columns (from the default 'anything' to, say, float).

But conversion fails for all string columns. For example:

DataNew := SubsDatatype(Data, "UN Region", string);

throws the following error:

Error, (in DataSeries:-SubsDatatype) in datatype conversion; try using the 'conversion' option for data type string

The column in question has values such as "Southern Asia", which I would have thought would convert.

What am I doing wrong?

How to obtain this solution using dsolve?...

Asked by: nm 11483

September 20 2025

0 1

Maple 2025.1 unable to solve this ode. Sympy gives the following two solutions which Maples verifies are correct.

Any trick or option that can help dsolve find these solutions?

>	interface(version);

>

restart;

>	ode:=diff(y(x),x) = (1+cos(x)sin(y(x)))tan(y(x));

>	sol:=dsolve(ode);

>	sol_1:=y(x)=arcsin( 2exp(x) / ( c__1 + sqrt(2)exp(x) * sin(x+Pi/4) ) ) + Pi

y(x) = arcsin(2*exp(x)/(c__1+2^(1/2)*exp(x)*sin(x+(1/4)*Pi)))+Pi

>	odetest(sol_1,ode)

>	sol_2:=y(x)=arcsin( 2exp(x) / ( c__1 - sqrt(2)exp(x) * sin(x+Pi/4) ) ) ;

y(x) = arcsin(2*exp(x)/(c__1-2^(1/2)*exp(x)*sin(x+(1/4)*Pi)))

>	odetest(sol_2,ode)

>

Download How_to_find_solution_sept_20_2025.mw

update:

OK, found out how. Needed transformation u(x)=sin(y(x)). Maple probably did not have this in one of the things to try.

Download How_to_find_solution_sept_20_2025_V2.mw

On format of this solution from dsolve. Why [{...}...

Asked by: nm 11483

September 19 2025

0 0

Why when given IC for this ode, where the IC do not really makes much sense, so was not used. But the question is on the format of the output of the Maple dsolve. It gives solution as [{y(t) = c__1}] instead of y(t) = c__1

Related question. Since Maple did not use the IC, should there have been warning message generated that IC was ignored?

Download strange_format_of_solution_sept_19_2025.mw

Can I use AI to Assist in Creating Maple Documents...

Asked by: ianmccr 170

September 17 2025

0 1

Maple documents appear to be an excellent way of documenting problems and procedures, and notetaking in general, but I have always found them laborious and frequently encountered problems with formatting. One may as well use Tex. This is the reason that I was struck by the ability of the GenerateDocument() procedure in the NaturalLanguage package to create Maple documents.

Is there any practical way to use NaturalLanguage for notetaking and documentation of procedures. By this I mean, the ability to input math and or explanatory text and request(direct) the AI to produce a maple document.

How in combined plot make the movement in both li...

Asked by: salim-barzani 1650

September 14 2025

1 4

In this kind of contour plot i have two line but when i change time variable t just contour of one line wil move the other is not do any movement and is stop how i can make the second plot one second line move too? also there is any way for ploting this kind any other option?

line-2-done.mw

How perform quantum field theory (QFT) integrals?...

Asked by: nmacsai 280

September 09 2025

1 0

It is possible to perform the simplest QFT calculations with second quantization, in Maple? Bosons in a box. See attached example. bosons_in_a_box.mw

Sure any general purpose programming language is capable of performing this task with enough effort. What I am interested in is if the physics tools has a standard way of dealing with these calculations. The general impedement when attempting the calculation is that integrations are perfomed by replacements with delta functions or kronecker delta functions, and its not clear how to force the Maple Physics package to recognize this or if that's possible. Part of the problem is that integrations in maple are defined in one dimension at a time where as in QFT the integration element is almost always atleast three dimensional, d^3x or dxdydzy, the later of which can get extremely cumbersome with even a small number of fields under consideration. I don't find much of what I am refering to mentioned in the help pages and I doubt these types of QFT calculations are possible to perform in Maple without addressing these issues.

bosons_in_a_box.mw

My maple document gives wrong units...

Asked by: Jacobegerup 5

September 09 2025

1 5

Hello everyone!
I have had an issue for the past weeks, where it seems like Maple has a problem identifying the correct units and also sometimes having an issue with defining a variable.

In the picture below you can see I have defined rho, m and tried to Solve V. I get m^3, which is fine but I cannot change the unit in the right bar (see picture).

I even tried to just take square root of my V to see if I could then change units. It seems like it thinks I'm playing with weight.
I even had 2 teachers trying to help me find a solution, without luck.

I have tried executing the whole document and also only bits of the document without luck.

Does anyone know this problem and has a solution?

Thanks in advance!

- Jacob

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