Unanswered Questions

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

I am calling a proc inside a module, where this module had local module variable initialized to false;

This works fine when calling the proc normally.

When calling it from Grid:-Run() it says it is not seeing the module local variable at all.

Why does this happen and is there a workaround this? Here a simple worksheet showing this.

I hope there is a way to make the proc inside the module see the local module variables when using it from Grid node. Otherwise, this whole thing will not work for me, as I have few variables initialized at modules level in number of places.

restart;

A := module()
  local DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",DEBUG_MSG);
  end proc;
end module;

 

_m1690230418240

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

"My flag is ", DEBUG_MSG

 

 

Download grid_with_module_local_variable_feb_7_2025.mw


Update 2/7/2025

I tried to change the access to the module local variables from the module local procs, by adding explicit A:- to each variable name where A here is the module name.

This works for normal calls, but not when using Grid to call the proc.

This workaround fail, it gives error that module does not export `%1`

May be I have to redo all my code not to use local module variables. But in some places I have to do this, in order to detect if something has happened before or not. I use module local variables to store state the presist after call is completed.

It looks like Grid is not meant to be used for calling proc() that uses/lives inside modules. But this makes Grid not very useful then for large application.

restart;

A := module()
  local DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",A:-DEBUG_MSG);
  end proc;
end module;

 

_m1690230418240

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) module does not export `%1`

 

 

Download grid_with_module_local_variable_feb_7_2025_V2.mw

I also tried to make the module variable as export instead of local, and that also did not work. I am starting to run out of ideas what else to try...
 

restart;

A := module()
  export DEBUG_MSG::truefalse:=false;
  export work_proc:=proc(n::integer)
     print("My flag is ",A:-DEBUG_MSG);
  end proc;
end module;

 

_m1690230418208

#THis works as expected
A:-work_proc(0)

"My flag is ", false

#this does not work. THe call works OK, but inside A:-workproc() it does not see module local variables.

Grid:-Set(A:-work_proc):
Grid:-Run(0,A:-work_proc,[ 0 ]);
Grid:-Wait();

Error, (in work_proc) `%1` does not evaluate to a module

 


 

Download grid_with_module_local_variable_feb_7_2025_V3.mw

 

 

scmch.mw

I can't get a graph. Is this code is correct.Please help.

How to integrate eq (4)? Since 'a', 'b', and 'c' are constant. 

restart

with(DEtools)

declare(z(x), y(x))

declare(z(x), y(x))

(1)

eq1 := (1/2)*(-z(x)^3-2*c*z(x))*(diff(diff(y(x), x), x))-((z(x)^2+2*c)*(diff(y(x), x))+b*z(x)^2+c*k+a)*(diff(z(x), x)) = 0

(1/2)*(-z(x)^3-2*c*z(x))*(diff(diff(y(x), x), x))-((z(x)^2+2*c)*(diff(y(x), x))+b*z(x)^2+c*k+a)*(diff(z(x), x)) = 0

(2)

eq2 := simplify(z(x)*eq1)

-z(x)*(z(x)*((1/2)*z(x)^2+c)*(diff(diff(y(x), x), x))+((z(x)^2+2*c)*(diff(y(x), x))+b*z(x)^2+c*k+a)*(diff(z(x), x))) = 0

(3)

eq3 := eval(int(lhs(eq2), x))

int(-z(x)*(z(x)*((1/2)*z(x)^2+c)*(diff(diff(y(x), x), x))+((z(x)^2+2*c)*(diff(y(x), x))+b*z(x)^2+c*k+a)*(diff(z(x), x))), x)

(4)

NULL

Download integration.mw

The modified Liouville equation

How to solve this pde for a general solution ?

The general solution in this form exist.

restart;

with(PDEtools): declare(u(x,t)); U:=diff_table(u(x,t));
PDE1:=U[t,t]=a^2*U[x,x]+b*exp(beta*U[]);
Sol11:=u(x,t)=1/beta*ln(2*(B^2-a^2*A^2)/(b*beta*(A*x+B*t+C)^2));
Sol12:=S->u(x,t)=1/beta*ln(8*a^2*C/(b*beta))
-2/beta*ln(S*(x+A)^2-S*a^2*(t+B)^2+S*C);
Test11:=pdetest(Sol11,PDE1);
Test12:=pdetest(Sol12(1),PDE1);
Test13:=pdetest(Sol12(-1),PDE1);

u(x, t)*`will now be displayed as`*u

 

table( [(  ) = u(x, t) ] )

 

diff(diff(u(x, t), t), t) = a^2*(diff(diff(u(x, t), x), x))+b*exp(beta*u(x, t))

 

u(x, t) = ln(2*(-A^2*a^2+B^2)/(b*beta*(A*x+B*t+C)^2))/beta

 

proc (S) options operator, arrow; u(x, t) = ln(8*a^2*C/(b*beta))/beta-2*ln(S*(x+A)^2-S*a^2*(t+B)^2+S*C)/beta end proc

 

0

 

0

 

0

(1)

The Soll11 can be plotted with a Explore plot in this form of soll11 with th eparameters , but suppose i try to get the general solution in Maple ?

infolevel[pdsolve] := 3

pdsolve(PDE1, generalsolution)

ans := pdsolve(PDE1);

What solvin gstrategy to follow ? : the pde is a non-linear wave eqation  with a exponentiel sourceterm
It seems that the pde can reduced to a ode? :

 

with(PDEtools):
declare(u(x,t));

# Stap 1: Definieer de PDE
PDE := diff(u(x,t), t,t) = a^2 * diff(u(x,t), x,x) + b * exp(beta * u(x,t));

# Stap 2: Definieer de transformatie naar karakteristieke variabelen
# Nieuw: x en t uitgedrukt in ξ en η
tr := {
    x = (xi + eta)/2,
    t = (eta - xi)/(2*a)
};

# Pas de transformatie toe op de PDE
simplified_PDE := dchange(tr, PDE, [xi, eta], params = [a, b, beta], simplify);

# Stap 3: Definieer de algemene oplossing
solution := u(x,t) = (1/beta) * ln(
    (-8*a^2/(b*beta)) *
    diff(_F1(x - a*t), x) * diff(_F2(x + a*t), x) /
    (_F1(x - a*t) + _F2(x + a*t))^2
);

# Stap 4: Controleer de oplossing (optioneel)
pdetest(solution, PDE);  # Moet 0 teruggeven als correct

u(x, t)*`will now be displayed as`*u

 

diff(diff(u(x, t), t), t) = a^2*(diff(diff(u(x, t), x), x))+b*exp(beta*u(x, t))

 

{t = (1/2)*(eta-xi)/a, x = (1/2)*xi+(1/2)*eta}

 

a^2*(diff(diff(u(xi, eta), xi), xi)-2*(diff(diff(u(xi, eta), eta), xi))+diff(diff(u(xi, eta), eta), eta)) = a^2*(diff(diff(u(xi, eta), xi), xi))+2*a^2*(diff(diff(u(xi, eta), eta), xi))+a^2*(diff(diff(u(xi, eta), eta), eta))+b*exp(beta*u(xi, eta))

 

u(x, t) = ln(-8*a^2*(D(_F1))(-a*t+x)*(D(_F2))(a*t+x)/(b*beta*(_F1(-a*t+x)+_F2(a*t+x))^2))/beta

 

0

(2)

missing some steps here : solution u  without  the pde reduced ?
there is a ode ?

# Definieer de ODE # vorige stappen ontbreken van de reduktie
ode := (v^2 - a^2) * diff(f(xi), xi, xi) = b * exp(beta * f(xi));

# Algemene oplossing zoeken
sol := dsolve(ode, f(xi));

(-a^2+v^2)*(diff(diff(f(xi), xi), xi)) = b*exp(beta*f(xi))

 

f(xi) = ln((1/2)*c__1*(tan((1/2)*(-c__1*a^2*beta+c__1*beta*v^2)^(1/2)*(c__2+xi)/(a^2-v^2))^2+1)/b)/beta

(3)

 

, ,

Question : how do i arrive on Soll11   in Maple  ?

 

Download liouville_reduced_2-2-2025_mprimes_vraag.mw

Hi All,

Maple is changing fast. It is not possible to run some older codes. 

Is it possible those who have a valid Maple license to have the old versions free of charge?

I have Maple 7, 2018, 2021 licenses but still have problem running older codes.

restart; with(PDEtools); declare(F(x, t), G(x, t), H(x, t))

F(x, t)*`will now be displayed as`*F

 

G(x, t)*`will now be displayed as`*G

 

H(x, t)*`will now be displayed as`*H

(1)

q := 1-(diff(diff(log(F(x, t)), x), t)); r := G/F; s := H/F

1-(diff(diff(F(x, t), t), x))/F(x, t)+(diff(F(x, t), x))*(diff(F(x, t), t))/F(x, t)^2

 

G/F

 

H/F

(2)

r1s1 := r*s; r1s1der := diff(r1s1(x, t), x)

qt := diff(q(x, t), t)

eq1B := F(x, t)^3*(qt+r1s1der) = 0; eq12B := simplify(expand(eq1B))

-F(x, t)^3*(diff((diff(diff(F(x, t), t), x))(x, t), t))/(F(x, t))(x, t)+F(x, t)^3*(diff(diff(F(x, t), t), x))(x, t)*(diff((F(x, t))(x, t), t))/(F(x, t))(x, t)^2+F(x, t)^3*(diff((diff(F(x, t), x))(x, t), t))*(diff(F(x, t), t))(x, t)/(F(x, t))(x, t)^2-2*F(x, t)^3*(diff(F(x, t), x))(x, t)*(diff(F(x, t), t))(x, t)*(diff((F(x, t))(x, t), t))/(F(x, t))(x, t)^3+F(x, t)^3*(diff(F(x, t), x))(x, t)*(diff((diff(F(x, t), t))(x, t), t))/(F(x, t))(x, t)^2+F(x, t)*(diff(G(x, t), x))*H(x, t)-2*G(x, t)*H(x, t)*(diff(F(x, t), x))+F(x, t)*G(x, t)*(diff(H(x, t), x)) = 0

(3)

D_x_x_G_F := (diff(G(x, t), x, x))*F(x, t)-2*(diff(G(x, t), x))*(diff(F(x, t), x))+G(x, t)*(diff(F(x, t), x, x)); D_t_t_F_F := F(x, t)*(diff(F(x, t), `$`(t, 2)))-2*(diff(F(x, t), t))^2

(diff(diff(G(x, t), x), x))*F(x, t)-2*(diff(G(x, t), x))*(diff(F(x, t), x))+G(x, t)*(diff(diff(F(x, t), x), x))

 

F(x, t)*(diff(diff(F(x, t), t), t))-2*(diff(F(x, t), t))^2

(4)

NULL

rxt := diff(diff(r(x, t), x), t)

eq2B := -2*q*r+rxt = 0

eq22B := simplify(expand(eq2B))

((-F*F(x, t)*G(x, t)+2*G*F(x, t)^2)*(diff(diff(F(x, t), t), x))+(diff(diff(G(x, t), t), x))*F*F(x, t)^2+((2*F*G(x, t)-2*G*F(x, t))*(diff(F(x, t), x))-F*(diff(G(x, t), x))*F(x, t))*(diff(F(x, t), t))-(diff(G(x, t), t))*(diff(F(x, t), x))*F*F(x, t)-2*G*F(x, t)^3)/(F*F(x, t)^3) = 0

(5)

sxt := diff(diff(s(x, t), x), t)

eq3B := -2*q*s+sxt = 0

eq32B := simplify(expand(eq3B))

((-F*F(x, t)*H(x, t)+2*H*F(x, t)^2)*(diff(diff(F(x, t), t), x))+(diff(diff(H(x, t), t), x))*F*F(x, t)^2+((2*F*H(x, t)-2*H*F(x, t))*(diff(F(x, t), x))-F*(diff(H(x, t), x))*F(x, t))*(diff(F(x, t), t))-(diff(H(x, t), t))*(diff(F(x, t), x))*F*F(x, t)-2*H*F(x, t)^3)/(F*F(x, t)^3) = 0

(6)

"#`# How to simplify Eqs. (3), (5) and (6) and write in terms of following bilineat operators` by using (4)"?""

NULL

NULL

Download BE.mw

Hi,

Experiencing the following problem.  One of our servers was cloned, the GUID was replaced and then rejoined to the domain.  All applications are working exept Maple.  The application launches but then closes right away.  No error messages provided so not sure where else to look for possible fixes to this problem.  The application is runing on Server 2022.

Thank you.

restart

with(PDEtools)

with(LinearAlgebra)

NULL

with(SolveTools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

declare(u(x, y, z, t))

u(x, y, z, t)*`will now be displayed as`*u

(2)

pde := diff(u(x, y, z, t), `$`(t, 2))+diff(u(x, y, z, t), `$`(x, 2))-(diff(u(x, y, z, t)^2, `$`(x, 2)))-(diff(u(x, y, z, t), `$`(x, 4)))+diff(diff(u(x, y, z, t), y)+diff(u(x, y, z, t), z)+diff(u(x, y, z, t), t), x)+2*(diff(u(x, y, z, t), y, t))+diff(u(x, y, z, t), `$`(y, 2)) = 0

diff(diff(u(x, y, z, t), t), t)+diff(diff(u(x, y, z, t), x), x)-2*(diff(u(x, y, z, t), x))^2-2*u(x, y, z, t)*(diff(diff(u(x, y, z, t), x), x))-(diff(diff(diff(diff(u(x, y, z, t), x), x), x), x))+diff(diff(u(x, y, z, t), x), y)+diff(diff(u(x, y, z, t), x), z)+diff(diff(u(x, y, z, t), t), x)+2*(diff(diff(u(x, y, z, t), t), y))+diff(diff(u(x, y, z, t), y), y) = 0

(3)

declare(v(t))

v(t)*`will now be displayed as`*v

(4)

declare(f(x, y, z, t))

f(x, y, z, t)*`will now be displayed as`*f

(5)

Q := u(x, y, z, t) = 6*(diff(ln(f(x, y, z, t)), `$`(x, 2)))

LL := diff(diff(diff(diff(diff(f(x, y, z, t), x), x), x), x), x)-(diff(diff(diff(f(x, y, z, t), x), x), x))-(diff(diff(diff(f(x, y, z, t), t), t), x))-(diff(diff(diff(f(x, y, z, t), t), x), x))-2*(diff(diff(diff(f(x, y, z, t), t), x), y))-(diff(diff(diff(f(x, y, z, t), x), x), y))-(diff(diff(diff(f(x, y, z, t), x), x), z))-(diff(diff(diff(f(x, y, z, t), x), y), y)) = 0

diff(diff(diff(diff(diff(f(x, y, z, t), x), x), x), x), x)-(diff(diff(diff(f(x, y, z, t), x), x), x))-(diff(diff(diff(f(x, y, z, t), t), t), x))-(diff(diff(diff(f(x, y, z, t), t), x), x))-2*(diff(diff(diff(f(x, y, z, t), t), x), y))-(diff(diff(diff(f(x, y, z, t), x), x), y))-(diff(diff(diff(f(x, y, z, t), x), x), z))-(diff(diff(diff(f(x, y, z, t), x), y), y)) = 0

(6)

S22 := f(x, y, z, t) = 1+exp((-(1/2)*k[1]-l[1]+(1/2)*sqrt(4*k[1]^4-3*k[1]^2-4*k[1]*s[1]))*t+k[1]*x+l[1]*y+s[1]*z)+exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*sqrt(4*k[2]^4-3*k[2]^2-4*k[2]*s[2]))*t)+B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*sqrt(4*k[1]^4-3*k[1]^2-4*k[1]*s[1]))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*sqrt(4*k[2]^4-3*k[2]^2-4*k[2]*s[2]))*t)

f(x, y, z, t) = 1+exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)+exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)+B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)

(7)

NULL

R11 := eval(LL, S22)

k[1]^5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)+k[2]^5*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)+B[1]*(k[1]+k[2])^5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^3*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^3*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^3*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))^2*k[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))^2*k[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))^2*(k[1]+k[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*k[1]^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*k[2]^2*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*(k[1]+k[2])^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-2*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*k[1]*l[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-2*(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*k[2]*l[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-2*B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*(k[1]+k[2])*(l[1]+l[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^2*l[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^2*l[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^2*(l[1]+l[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^2*s[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^2*s[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^2*(s[1]+s[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]*l[1]^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]*l[2]^2*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])*(l[1]+l[2])^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t) = 0

(8)

L4 := collect(%, [x, y, t], 'distributed')

k[1]^5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)+k[2]^5*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)+B[1]*(k[1]+k[2])^5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^3*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^3*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^3*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))^2*k[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))^2*k[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))^2*(k[1]+k[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*k[1]^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*k[2]^2*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*(k[1]+k[2])^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-2*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*k[1]*l[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-2*(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*k[2]*l[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-2*B[1]*(-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*(k[1]+k[2])*(l[1]+l[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^2*l[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^2*l[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^2*(l[1]+l[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]^2*s[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]^2*s[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])^2*(s[1]+s[2])*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[1]*l[1]^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z)-k[2]*l[2]^2*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-B[1]*(k[1]+k[2])*(l[1]+l[2])^2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t) = 0

(9)

indets(%)

{t, x, y, z, B[1], k[1], k[2], l[1], l[2], s[1], s[2], (4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2), (4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2), exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t), exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z), exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)}

(10)

eq2 := algsubs(exp((-(1/2)*k[1]-l[1]+(1/2)*sqrt(4*k[1]^4-3*k[1]^2-4*k[1]*s[1]))*t+k[1]*x+l[1]*y+s[1]*z) = X, L4)

-(1/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])*k[2]-exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[2]^2*s[1]-exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[2]^2*s[2]+5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^4*k[2]+10*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^3*k[2]^2+10*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^2*k[2]^3+5*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]*k[2]^4-exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^2*s[1]-exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^2*s[2]-(9/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*k[2]^2*k[1]-(9/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*k[1]^2*k[2]-(1/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])*k[1]-(1/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])*k[2]-(1/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])*k[1]+exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]^5+exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[2]^5-(3/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*k[1]^3-(3/4)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*k[2]^3-(1/4)*k[1]*X*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])-(1/4)*k[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])-k[1]^2*s[1]*X-(3/4)*k[2]^3*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)+k[2]^5*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)-k[2]^2*s[2]*exp(k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)+k[1]^5*X-(3/4)*k[1]^3*X-2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]*k[2]*s[1]-2*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*B[1]*k[1]*k[2]*s[2]-(1/2)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2)*k[1]-(1/2)*B[1]*exp((-(1/2)*k[1]-l[1]+(1/2)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2))*t+k[1]*x+l[1]*y+s[1]*z+k[2]*x+l[2]*y+s[2]*z+(-(1/2)*k[2]-l[2]+(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2))*t)*(4*k[1]^4-3*k[1]^2-4*k[1]*s[1])^(1/2)*(4*k[2]^4-3*k[2]^2-4*k[2]*s[2])^(1/2)*k[2] = 0

(11)

eq3 := simplify(eq2)

-(1/2)*(k[1]+k[2])*B[1]*((k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-8*k[1]*k[2]^3-12*k[2]^2*k[1]^2+(-8*k[1]^3+3*k[1]+2*s[1])*k[2]+2*s[2]*k[1])*exp((1/2)*t*(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*(-k[1]-k[2]-2*l[1]-2*l[2])*t+k[1]*x+k[2]*x+l[1]*y+l[2]*y+z*(s[1]+s[2])) = 0

(12)

indets(eq3)

{t, x, y, z, B[1], k[1], k[2], l[1], l[2], s[1], s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2), exp((1/2)*t*(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*(-k[1]-k[2]-2*l[1]-2*l[2])*t+k[1]*x+k[2]*x+l[1]*y+l[2]*y+z*(s[1]+s[2]))}

(13)

eq4 := algsubs(exp((1/2)*t*sqrt(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))+(1/2)*t*sqrt(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))+(1/2)*(-k[1]-k[2]-2*l[1]-2*l[2])*t+k[1]*x+k[2]*x+l[1]*y+l[2]*y+z*(s[1]+s[2])) = V, eq3)

-(1/2)*(k[1]+k[2])*B[1]*(-8*k[2]*k[1]^3-12*k[2]^2*k[1]^2-8*k[1]*k[2]^3+(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+3*k[1]*k[2]+2*s[2]*k[1]+2*s[1]*k[2])*V = 0

(14)

indets(eq4)

{V, B[1], k[1], k[2], s[1], s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}

(15)

eqs := {coeffs(collect(numer(normal(lhs(eq4))), {V}, 'distributed'), {V})}; nops(%); indets(eqs)

{-(k[1]+k[2])*B[1]*(-8*k[2]*k[1]^3-12*k[2]^2*k[1]^2-8*k[1]*k[2]^3+(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+3*k[1]*k[2]+2*s[2]*k[1]+2*s[1]*k[2])}

 

1

 

{B[1], k[1], k[2], s[1], s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}

(16)

vars := indets(eqs); ans := solve(eqs, vars)

{B[1], k[1], k[2], s[1], s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}

 

Warning, solving for expressions other than names or functions is not recommended.

 

{B[1] = B[1], k[1] = -k[2], k[2] = k[2], s[1] = s[1], s[2] = s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2) = (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2) = (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}, {B[1] = 0, k[1] = k[1], k[2] = k[2], s[1] = s[1], s[2] = s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2) = (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2) = (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}, {B[1] = B[1], k[1] = k[1], k[2] = k[2], s[1] = (1/2)*(8*k[2]*k[1]^3+12*k[2]^2*k[1]^2+8*k[1]*k[2]^3-3*k[1]*k[2]-2*s[2]*k[1]-(k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))/k[2], s[2] = s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2) = (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2) = (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}

(17)

case2 := ans[1]

{B[1] = B[1], k[1] = -k[2], k[2] = k[2], s[1] = s[1], s[2] = s[2], (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2) = (k[1]*(4*k[1]^3-3*k[1]-4*s[1]))^(1/2), (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2) = (k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)}

(18)

FF := subs(case2, S22)

NULL

F11 := eval(Q, FF)

pdetest(F11, pde)

-6*k[2]^2*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(B[1]*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+56*k[2]^4*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-6*k[2]^2*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+B[1]*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-24*B[1]*k[2]^2*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+2*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+4*k[2]*s[1]*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-4*k[2]*s[2]*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+2*k[2]*s[2]*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*k[2]*s[1]*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-224*k[2]^4*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-4*k[2]*s[2]*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+4*k[2]*s[1]*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+6*B[1]*k[2]^2*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-8*B[1]*k[2]^4*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+3*B[1]*k[2]^2*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-28*B[1]*k[2]^4*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-28*B[1]*k[2]^4*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+24*k[2]^2*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+16*B[1]*k[2]*s[1]*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-16*B[1]*k[2]*s[2]*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+8*B[1]*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-4*B[1]*k[2]*s[1]*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+4*B[1]*k[2]*s[2]*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*B[1]*exp(-t*k[2]+2*t*l[1]+4*k[2]*x+2*l[2]*y+2*s[2]*z+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-2*B[1]*k[2]*s[1]*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+2*B[1]*k[2]*s[2]*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*B[1]*k[2]*s[1]*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*B[1]*k[2]*s[2]*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*k[2]*s[1]*B[1]^2*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*k[2]*s[2]*B[1]^2*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*k[2]*s[2]*B[1]^2*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*k[2]*s[1]*B[1]^2*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-4*B[1]*k[2]*s[1]*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+4*B[1]*k[2]*s[2]*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-2*B[1]*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-2*B[1]*k[2]*s[2]*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+32*B[1]*k[2]^4*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+2*k[2]*s[2]*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-2*k[2]*s[1]*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*B[1]*k[2]*s[1]*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+3*k[2]^2*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+3*k[2]^2*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-28*k[2]^4*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-28*k[2]^4*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+3*B[1]*k[2]^2*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-3*B[1]*k[2]^2*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+4*B[1]*k[2]^4*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+2*B[1]*k[2]*s[1]*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-2*B[1]*k[2]*s[2]*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+B[1]^2*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+B[1]^2*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-B[1]*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-B[1]*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+4*B[1]*k[2]^4*exp(-(1/2)*t*k[2]+2*t*l[1]+t*l[2]+3*k[2]*x+l[2]*y+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+4*B[1]^2*k[2]^4*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+16*k[2]*s[2]*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-16*k[2]*s[1]*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-8*exp(t*l[1]+t*l[2]+2*k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)-6*k[2]^2*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+56*k[2]^4*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+6*B[1]*k[2]^2*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-8*B[1]*k[2]^4*exp(t*k[2]+2*t*l[2]+2*l[1]*y+2*s[1]*z+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-3*B[1]*k[2]^2*exp((1/2)*t*k[2]+t*l[1]+2*t*l[2]+k[2]*x+l[1]*y+s[1]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))-3*B[1]^2*k[2]^2*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))-3*B[1]^2*k[2]^2*exp(t*l[1]+3*k[2]*x+l[1]*y+2*l[2]*y+s[1]*z+2*s[2]*z-(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+4*B[1]^2*k[2]^4*exp(t*l[2]+k[2]*x+2*l[1]*y+l[2]*y+2*s[1]*z+s[2]*z+(1/2)*t*k[2]+t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2)))/(B[1]*exp(k[2]*x+l[1]*y+l[2]*y+s[1]*z+s[2]*z+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2)+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+exp(t*l[1]+2*k[2]*x+l[2]*y+s[2]*z-(1/2)*t*k[2]+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]-4*s[2]))^(1/2))+exp(t*l[2]+l[1]*y+s[1]*z+(1/2)*t*k[2]+(1/2)*t*(k[2]*(4*k[2]^3-3*k[2]+4*s[1]))^(1/2))+exp(t*l[1]+t*l[2]+x*k[2]))^4

(19)
 

NULL

Download hard_parameters.mw

If I do 

      Grid:-Run(0,foo,[n])

in worksheet, where foo is proc() defined in same worksheet before, which just prints something, but I do not see anything on the screen printed when running. It seems to run, since the Grid:-Wait() seems to take 1-2 second to finished, so I am sure the proc foo() was called. It is just the print output seem not to show.

When changing the above to

      Grid:-Run(0,foo,[n],'assignto'='ans0')

now I see the print on the screen but only when evaluating ans0. Actually the print does not show, had to evaluate ans0 to see the print string and also the return value.

In my code, I do not need to return anything, I just want to print to the screen. 

How to make the first one display on the screen? Is output from foo proc going somewhere else with Grid?

I made sure to use Grid:-Set(foo): to define foo for the grid as the help says to do. 

Here is worksheet which shows this problem. I simply want to use print inside foo and see the result on the screen when running foo on a node. 

Is this a bug?

restart;

interface(version);

`Standard Worksheet Interface, Maple 2024.2, Windows 10, October 29 2024 Build ID 1872373`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1843 and is the same as the version installed in this computer, created 2025, January 25, 22:5 hours Pacific Time.`

Why task running on node do not print anything?

 

restart;

foo:=proc(n)
   print("Entered foo, n=",n);    
end proc;

proc (n) print("Entered foo, n=", n) end proc

Grid:-Set(foo):
Grid:-Setup("local",numnodes=1); #just use one node
n:=1:
Grid:-Run(0,foo,[n]):
Grid:-Wait();

#nothing is printed from foo

 

To make it print, why assignto is needed?

 

restart;

foo:=proc(n)
   print("Entered foo, n=",n);    
   RETURN("done");
end proc;

proc (n) print("Entered foo, n=", n); RETURN("done") end proc

Grid:-Set(foo):
Grid:-Setup("local",numnodes=1);
n:=1:
Grid:-Run(0,foo,[n],'assignto'='ans0'):
Grid:-Wait();
ans0

"Entered foo, n=", 1

"done"

 

 

Download grid_question_jan_27_2025.mw

Update

I found why, but I do not understand why this fixes it. Removing numnodes =n  where n is number of nodes needed, now it works, i.e. proc prints on the screen ok as expected !

restart;

foo:=proc(n)
   print("Entered foo, n=",n);    
end proc;

proc (n) print("Entered foo, n=", n) end proc

Grid:-Set(foo):
Grid:-Setup("local"); # remove numnodes = ... , then it works. why??
n:=1:
Grid:-Run(0,foo,[n]):
Grid:-Wait();

"Entered foo, n=", 1

 

 

Download grid_question_jan_27_2025_V2.mw

But why? why is adding numnodes =1 or any other number makes it not print? Help says

So I do not see why adding this option makes it not print.

This has to be a bug, right? If not, then what is the logical explanation for this?

in a lot of my equation i have such problem and really i don't know how fix this also i try to put : in end and sometime is work and i keep contionues  but sometime not there is any way for solve this problem?

limit.mw

Let us suppose I open Maple, write a worksheet, save it in a folder A and then quit Maple.

Now I run a new session by double clicking on some mw file located in folder B, let us say //B/test_1.mw.
Once opened I do some modifications and decide to save this worksheet into a new file, let us say test_2.mw in the same folder B test_1.mw belongs to (which means I use Save As from the menu bar)

I'm regularly fooled by the fact that the default folder is not B, but the folder A I used in the previous session.

I find this very unpleasant.
Is this a Maple (2015) issue or something related to my operating system (Mac OSX Catalina)?
In case it is a Maple issue which is still present in more recent Maple versions, Would it be possible to set the default backup folder to be the folder to which the active worksheet belongs?

Thanks in advance

It's 2024 and this is still something that doesn't exist? I'd just like to swap the Enter/Shfit+Enter behaviors since I find myself writing a lot of multi-line and custom procs and boy howdy it'd be nice if I could make Maple behave at least the littllest bit like, I dunno, every other product I own and use.

i found solution of PDE but there is some different from my solution and paper solution so there is must be a mistake becuase he solved by maple too he mentioned in the paper i try to figure out but i can't see any mistake from my solution can anyone watch where i did mistake, i change some letter in finding parameter but they are same like p=k&h=A&n=p&w=n

here is paper solution 

parameter-different.mw

I'm calcuating an endomorphism in 2d dimensions. It is contructed out of a tensor contraction, for example, in 6-dimensions, the endomorphism is

K[mu,~nu] = LeviCivita[~alpha,~beta,~gamma,~delta,~upsilon,~nu]*C[alpha,beta,gamma]*C[delta,upsilon,mu]

I appreciate that, in terms of computation, this gets big quickly: it's something like O(exp) in time to sum over repeated indices in each matrix entry. Therefore, I thought, instead of putting the above expression in Define, I could make a matrix with unsummed entries, and then do the sums in parallel using Threads[Map](SumOverRepeatedIndices,...) but looking at my CPU usage and comparing execution times, it doesn't appear that this is working.

Is there any way I can more efficiently calculate these matrix entries?

This isn't the first time that I've seen a question that doesn't seem to have received a comment or answer, but which, when I click on the question title, turns out to have received one.
Here's the opposite phenomenon: a question appears to have two comments or answers, but none of them exist (9:38 GMT+1)

Screen capture from the main page

Screen capture from the question page

1 2 3 4 5 6 7 Last Page 1 of 353