MaplePrimes Questions

Hello,

I am working on a problem involving Newton's Law of Cooling in Maple, I want to:

  1. Solve for the cooling constant k,
  2. Find how long it takes for the pizza to cool to 300 F,

Here is the code I am currently using:

with(plots):

T__s := 75;  # Surrounding temperature
ode := diff(T(t), t) = k*(T(t) - T__s);

# General solution
general_solution := dsolve(ode, T(t));

# Initial condition: T(0) = 350
T := t -> 75 + exp(k*t)*c__1;
C := solve(350 = 75 + exp(0)*C__1, C__1);

# Solve for k using T(5) = 340
k := solve(340 = 75 + 275*exp(5*k), k);
k := evalf(k);

# Solve for time to reach T = 300
time_to_reach := solve(300 = 75 + 275*exp(k*t), t);
time_to_reach := evalf(time_to_reach);

# Define T(t) with solved constants
T := t -> 75 + 275*exp(k*t);

# Plot temperature over time
temp := plot(T, t = 0 .. 280, 75 .. 360, title = "Exponential Decay of Pizza Temperature", labels = ["t (minutes)", "T(t) (°F)"], color = blue);

# Plot T__s as a horizontal line
T__sur := plot(T__s, t = 0 .. 280, color = red, linestyle = dash);

# Combine and display plots
display([temp, T__sur]);

The code works, but I am wondering if there is a more efficient or genuine way.

I would greatly appreciate any advice or alternative approaches to solving and visualizing this in Maple.

Thank you in advance for your insights!

loading

Error occurred during PDF generation. Please refresh the page and try again

With the latest version of Maple I get an indenting space when I click "Shift+Enter".
For example 

But I don't like it: suggestion? Thanks

Hint: we can vanish the different power of x and t!


 

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)

pde := diff(w(x, t), `$`(t, 2))-a*(diff(w(x, t), `$`(x, 4)))-b*(diff(w(x, t), `$`(x, 2), `$`(t, 2)))+(diff(w(x, t), t))*(diff(w(x, t), `$`(x, 2)))+2*(diff(w(x, t), x))*(diff(w(x, t), x, t)) = 0

diff(diff(w(x, t), t), t)-a*(diff(diff(diff(diff(w(x, t), x), x), x), x))-b*(diff(diff(diff(diff(w(x, t), t), t), x), x))+(diff(w(x, t), t))*(diff(diff(w(x, t), x), x))+2*(diff(w(x, t), x))*(diff(diff(w(x, t), t), x)) = 0

(2)

Y := w(x, t) = 2*(diff(ln(f(x, t)), `$`(x, 2)))

w(x, t) = 2*(diff(diff(f(x, t), x), x))/f(x, t)-2*(diff(f(x, t), x))^2/f(x, t)^2

(3)

L := eval(pde, Y)

2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))/f(x, t)^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-4*(diff(diff(f(x, t), t), x))^2/f(x, t)^2+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))/f(x, t)^3-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-a*(2*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))/f(x, t)-12*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+60*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-30*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))^3/f(x, t)^4+240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))*(diff(diff(f(x, t), x), x))/f(x, t)^3-20*(diff(diff(diff(f(x, t), x), x), x))^2/f(x, t)^2+720*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4/f(x, t)^5-540*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2/f(x, t)^4+60*(diff(diff(f(x, t), x), x))^3/f(x, t)^3-240*(diff(f(x, t), x))^6/f(x, t)^6)-b*(48*(diff(f(x, t), x))^4*(diff(diff(f(x, t), t), t))/f(x, t)^5-144*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))^2/f(x, t)^4+32*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3+96*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))/f(x, t)^3-48*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))/f(x, t)^4-36*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+32*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2/f(x, t)^3+16*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(diff(f(x, t), t), t), x))*(diff(f(x, t), x))/f(x, t)^3+12*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-48*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^4-240*(diff(f(x, t), x))^4*(diff(f(x, t), t))^2/f(x, t)^6+384*(diff(f(x, t), x))^3*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^5-144*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2/f(x, t)^4+288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))^2/f(x, t)^5-288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^4-72*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2/f(x, t)^4+2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))/f(x, t)+24*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)^3-8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2-12*(diff(diff(diff(diff(f(x, t), t), t), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*(diff(f(x, t), t))/f(x, t)^2-16*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2-12*(diff(diff(diff(f(x, t), t), x), x))^2/f(x, t)^2+4*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2)+(2*(diff(diff(diff(f(x, t), t), x), x))/f(x, t)-2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))/f(x, t)^2-4*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+4*(diff(f(x, t), x))^2*(diff(f(x, t), t))/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), x), x), x), x))/f(x, t)-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+24*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-6*(diff(diff(f(x, t), x), x))^2/f(x, t)^2-12*(diff(f(x, t), x))^4/f(x, t)^4)+2*(2*(diff(diff(diff(f(x, t), x), x), x))/f(x, t)-6*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))/f(x, t)^2+4*(diff(f(x, t), x))^3/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), t), x), x), x))/f(x, t)-6*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))/f(x, t)^2-2*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3-6*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+12*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))/f(x, t)^3-12*(diff(f(x, t), x))^3*(diff(f(x, t), t))/f(x, t)^4) = 0

(4)

numer(lhs(2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))/f(x, t)^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-4*(diff(diff(f(x, t), t), x))^2/f(x, t)^2+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))/f(x, t)^3-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-a*(2*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))/f(x, t)-12*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+60*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-30*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))^3/f(x, t)^4+240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))*(diff(diff(f(x, t), x), x))/f(x, t)^3-20*(diff(diff(diff(f(x, t), x), x), x))^2/f(x, t)^2+720*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4/f(x, t)^5-540*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2/f(x, t)^4+60*(diff(diff(f(x, t), x), x))^3/f(x, t)^3-240*(diff(f(x, t), x))^6/f(x, t)^6)-b*(32*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3+96*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))/f(x, t)^3-48*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))/f(x, t)^4+32*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^3+16*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(diff(f(x, t), t), t), x))*(diff(f(x, t), x))/f(x, t)^3+384*(diff(f(x, t), x))^3*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^5-144*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2/f(x, t)^4+288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))^2/f(x, t)^5-72*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2/f(x, t)^4-288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^4-2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2+48*(diff(f(x, t), x))^4*(diff(diff(f(x, t), t), t))/f(x, t)^5-48*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^4-36*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), t))^2/f(x, t)^3+24*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)^3-16*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2-12*(diff(diff(diff(diff(f(x, t), t), t), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2/f(x, t)^3-240*(diff(f(x, t), x))^4*(diff(f(x, t), t))^2/f(x, t)^6-144*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))^2/f(x, t)^4-8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))/f(x, t)-12*(diff(diff(diff(f(x, t), t), x), x))^2/f(x, t)^2)+(2*(diff(diff(diff(f(x, t), t), x), x))/f(x, t)-2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))/f(x, t)^2-4*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+4*(diff(f(x, t), x))^2*(diff(f(x, t), t))/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), x), x), x), x))/f(x, t)-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+24*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-6*(diff(diff(f(x, t), x), x))^2/f(x, t)^2-12*(diff(f(x, t), x))^4/f(x, t)^4)+2*(2*(diff(diff(diff(f(x, t), x), x), x))/f(x, t)-6*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))/f(x, t)^2+4*(diff(f(x, t), x))^3/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), t), x), x), x))/f(x, t)-6*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))/f(x, t)^2-2*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3-6*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+12*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))/f(x, t)^3-12*(diff(f(x, t), x))^3*(diff(f(x, t), t))/f(x, t)^4) = 0))*denom(rhs(2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))/f(x, t)^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-4*(diff(diff(f(x, t), t), x))^2/f(x, t)^2+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))/f(x, t)^3-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-a*(2*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))/f(x, t)-12*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+60*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-30*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))^3/f(x, t)^4+240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))*(diff(diff(f(x, t), x), x))/f(x, t)^3-20*(diff(diff(diff(f(x, t), x), x), x))^2/f(x, t)^2+720*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4/f(x, t)^5-540*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2/f(x, t)^4+60*(diff(diff(f(x, t), x), x))^3/f(x, t)^3-240*(diff(f(x, t), x))^6/f(x, t)^6)-b*(32*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3+96*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))/f(x, t)^3-48*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))/f(x, t)^4+32*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^3+16*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(diff(f(x, t), t), t), x))*(diff(f(x, t), x))/f(x, t)^3+384*(diff(f(x, t), x))^3*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^5-144*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2/f(x, t)^4+288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))^2/f(x, t)^5-72*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2/f(x, t)^4-288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^4-2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2+48*(diff(f(x, t), x))^4*(diff(diff(f(x, t), t), t))/f(x, t)^5-48*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^4-36*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), t))^2/f(x, t)^3+24*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)^3-16*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2-12*(diff(diff(diff(diff(f(x, t), t), t), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2/f(x, t)^3-240*(diff(f(x, t), x))^4*(diff(f(x, t), t))^2/f(x, t)^6-144*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))^2/f(x, t)^4-8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))/f(x, t)-12*(diff(diff(diff(f(x, t), t), x), x))^2/f(x, t)^2)+(2*(diff(diff(diff(f(x, t), t), x), x))/f(x, t)-2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))/f(x, t)^2-4*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+4*(diff(f(x, t), x))^2*(diff(f(x, t), t))/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), x), x), x), x))/f(x, t)-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+24*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-6*(diff(diff(f(x, t), x), x))^2/f(x, t)^2-12*(diff(f(x, t), x))^4/f(x, t)^4)+2*(2*(diff(diff(diff(f(x, t), x), x), x))/f(x, t)-6*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))/f(x, t)^2+4*(diff(f(x, t), x))^3/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), t), x), x), x))/f(x, t)-6*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))/f(x, t)^2-2*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3-6*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+12*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))/f(x, t)^3-12*(diff(f(x, t), x))^3*(diff(f(x, t), t))/f(x, t)^4) = 0)) = numer(rhs(2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))/f(x, t)^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-4*(diff(diff(f(x, t), t), x))^2/f(x, t)^2+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))/f(x, t)^3-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-a*(2*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))/f(x, t)-12*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+60*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-30*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))^3/f(x, t)^4+240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))*(diff(diff(f(x, t), x), x))/f(x, t)^3-20*(diff(diff(diff(f(x, t), x), x), x))^2/f(x, t)^2+720*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4/f(x, t)^5-540*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2/f(x, t)^4+60*(diff(diff(f(x, t), x), x))^3/f(x, t)^3-240*(diff(f(x, t), x))^6/f(x, t)^6)-b*(32*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3+96*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))/f(x, t)^3-48*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))/f(x, t)^4+32*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^3+16*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(diff(f(x, t), t), t), x))*(diff(f(x, t), x))/f(x, t)^3+384*(diff(f(x, t), x))^3*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^5-144*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2/f(x, t)^4+288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))^2/f(x, t)^5-72*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2/f(x, t)^4-288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^4-2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2+48*(diff(f(x, t), x))^4*(diff(diff(f(x, t), t), t))/f(x, t)^5-48*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^4-36*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), t))^2/f(x, t)^3+24*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)^3-16*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2-12*(diff(diff(diff(diff(f(x, t), t), t), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2/f(x, t)^3-240*(diff(f(x, t), x))^4*(diff(f(x, t), t))^2/f(x, t)^6-144*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))^2/f(x, t)^4-8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))/f(x, t)-12*(diff(diff(diff(f(x, t), t), x), x))^2/f(x, t)^2)+(2*(diff(diff(diff(f(x, t), t), x), x))/f(x, t)-2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))/f(x, t)^2-4*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+4*(diff(f(x, t), x))^2*(diff(f(x, t), t))/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), x), x), x), x))/f(x, t)-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+24*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-6*(diff(diff(f(x, t), x), x))^2/f(x, t)^2-12*(diff(f(x, t), x))^4/f(x, t)^4)+2*(2*(diff(diff(diff(f(x, t), x), x), x))/f(x, t)-6*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))/f(x, t)^2+4*(diff(f(x, t), x))^3/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), t), x), x), x))/f(x, t)-6*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))/f(x, t)^2-2*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3-6*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+12*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))/f(x, t)^3-12*(diff(f(x, t), x))^3*(diff(f(x, t), t))/f(x, t)^4) = 0))*denom(lhs(2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))/f(x, t)^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2/f(x, t)^3-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-4*(diff(diff(f(x, t), t), x))^2/f(x, t)^2+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))/f(x, t)^3-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3-a*(2*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))/f(x, t)-12*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+60*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-30*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))^3/f(x, t)^4+240*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))*(diff(diff(f(x, t), x), x))/f(x, t)^3-20*(diff(diff(diff(f(x, t), x), x), x))^2/f(x, t)^2+720*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4/f(x, t)^5-540*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2/f(x, t)^4+60*(diff(diff(f(x, t), x), x))^3/f(x, t)^3-240*(diff(f(x, t), x))^6/f(x, t)^6)-b*(32*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3+96*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))/f(x, t)^3-48*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))/f(x, t)^4+32*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^3+16*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(diff(f(x, t), t), t), x))*(diff(f(x, t), x))/f(x, t)^3+384*(diff(f(x, t), x))^3*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))/f(x, t)^5-144*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2/f(x, t)^4+288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*(diff(f(x, t), x))^2/f(x, t)^5-72*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2/f(x, t)^4-288*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^4-2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(diff(f(x, t), t), t))/f(x, t)^2-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^2+48*(diff(f(x, t), x))^4*(diff(diff(f(x, t), t), t))/f(x, t)^5-48*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))/f(x, t)^4-36*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), t))^2/f(x, t)^4+4*(diff(diff(diff(diff(f(x, t), x), x), x), x))*(diff(f(x, t), t))^2/f(x, t)^3+24*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))/f(x, t)^3-16*(diff(diff(diff(diff(f(x, t), t), x), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2-12*(diff(diff(diff(diff(f(x, t), t), t), x), x))*(diff(diff(f(x, t), x), x))/f(x, t)^2-4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))/f(x, t)^3+48*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2/f(x, t)^3-240*(diff(f(x, t), x))^4*(diff(f(x, t), t))^2/f(x, t)^6-144*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))^2/f(x, t)^4-8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))/f(x, t)-12*(diff(diff(diff(f(x, t), t), x), x))^2/f(x, t)^2)+(2*(diff(diff(diff(f(x, t), t), x), x))/f(x, t)-2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))/f(x, t)^2-4*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+4*(diff(f(x, t), x))^2*(diff(f(x, t), t))/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), x), x), x), x))/f(x, t)-8*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), x))/f(x, t)^2+24*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2/f(x, t)^3-6*(diff(diff(f(x, t), x), x))^2/f(x, t)^2-12*(diff(f(x, t), x))^4/f(x, t)^4)+2*(2*(diff(diff(diff(f(x, t), x), x), x))/f(x, t)-6*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))/f(x, t)^2+4*(diff(f(x, t), x))^3/f(x, t)^3)*(2*(diff(diff(diff(diff(f(x, t), t), x), x), x))/f(x, t)-6*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))/f(x, t)^2-2*(diff(diff(diff(f(x, t), x), x), x))*(diff(f(x, t), t))/f(x, t)^2+12*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))/f(x, t)^3-6*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))/f(x, t)^2+12*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), x))/f(x, t)^3-12*(diff(f(x, t), x))^3*(diff(f(x, t), t))/f(x, t)^4) = 0))

2*(diff(diff(diff(diff(f(x, t), t), t), x), x))*f(x, t)^6-4*(diff(diff(f(x, t), t), x))^2*f(x, t)^5-144*(diff(f(x, t), t))*(diff(f(x, t), x))^6+288*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*b-48*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*b-96*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*b-32*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))*(diff(diff(diff(f(x, t), x), x), x))*b-32*f(x, t)^4*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))*b-48*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))*b-240*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*a-16*f(x, t)^4*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*b+144*f(x, t)^3*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2*b+48*f(x, t)^3*(diff(f(x, t), t))^2*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*b+72*f(x, t)^3*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2*b-288*f(x, t)^2*(diff(f(x, t), t))^2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2*b-384*f(x, t)^2*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^3*b+88*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))+8*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*f(x, t)^5*(diff(f(x, t), x))*b+4*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*f(x, t)^5*(diff(f(x, t), t))*b+30*f(x, t)^5*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))*a+12*f(x, t)^5*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), t), t), x), x))*b+2*f(x, t)^5*(diff(diff(f(x, t), t), t))*(diff(diff(diff(diff(f(x, t), x), x), x), x))*b+16*f(x, t)^5*(diff(diff(f(x, t), t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))*b+12*f(x, t)^5*(diff(f(x, t), x))*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*a+8*f(x, t)^5*(diff(diff(diff(f(x, t), t), t), x))*(diff(diff(diff(f(x, t), x), x), x))*b-4*f(x, t)^4*(diff(f(x, t), t))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*b-12*f(x, t)^4*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))*b-48*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2*b-60*f(x, t)^4*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*a-24*f(x, t)^4*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))*b+36*f(x, t)^3*(diff(f(x, t), t))^2*(diff(diff(f(x, t), x), x))^2*b+540*f(x, t)^3*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2*a+144*f(x, t)^3*(diff(diff(f(x, t), t), x))^2*(diff(f(x, t), x))^2*b+48*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))*b+240*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), x), x), x))*a-720*f(x, t)^2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4*a-48*f(x, t)^2*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^4*b+240*f(x, t)*(diff(f(x, t), t))^2*(diff(f(x, t), x))^4*b-40*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))-4*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))+16*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))*f(x, t)^4-24*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(diff(diff(f(x, t), x), x), x))-24*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))-8*f(x, t)^4*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))+120*f(x, t)^3*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2+8*f(x, t)^3*(diff(f(x, t), t))*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))+96*f(x, t)^3*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))+80*f(x, t)^3*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^2*(diff(diff(diff(f(x, t), x), x), x))-216*f(x, t)^2*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2-96*f(x, t)^2*(diff(f(x, t), t))*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), x), x), x))-288*f(x, t)^2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^3+360*f(x, t)*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4-2*(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))*f(x, t)^6*b-2*f(x, t)^6*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))*a+12*f(x, t)^5*(diff(diff(diff(f(x, t), t), x), x))^2*b+20*f(x, t)^5*(diff(diff(diff(f(x, t), x), x), x))^2*a-60*f(x, t)^4*(diff(diff(f(x, t), x), x))^3*a+240*f(x, t)*(diff(f(x, t), x))^6*a-4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*f(x, t)^5+4*f(x, t)^5*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))-2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*f(x, t)^5-4*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))*f(x, t)^5+8*f(x, t)^5*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))-12*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), x), x))^2+4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*f(x, t)^4-8*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(diff(f(x, t), x), x), x))^2+4*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))*f(x, t)^4-12*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2*f(x, t)^3+12*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))^3+16*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(diff(f(x, t), t), x), x), x))-72*f(x, t)^2*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))^4+144*f(x, t)*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^5 = 0

(5)

L1 := %*(1/2)

(diff(diff(diff(diff(f(x, t), t), t), x), x))*f(x, t)^6-2*(diff(diff(f(x, t), t), x))^2*f(x, t)^5-72*(diff(f(x, t), t))*(diff(f(x, t), x))^6+144*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*b-24*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*b-48*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*b-16*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))*(diff(diff(diff(f(x, t), x), x), x))*b-16*f(x, t)^4*(diff(f(x, t), t))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))*b-24*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))*b-120*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*a-8*f(x, t)^4*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*b+72*f(x, t)^3*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*(diff(f(x, t), x))^2*b+24*f(x, t)^3*(diff(f(x, t), t))^2*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))*b+36*f(x, t)^3*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^2*b-144*f(x, t)^2*(diff(f(x, t), t))^2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2*b-192*f(x, t)^2*(diff(f(x, t), t))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^3*b+44*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))+4*(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x))*f(x, t)^5*(diff(f(x, t), x))*b+2*(diff(diff(diff(diff(diff(f(x, t), t), x), x), x), x))*f(x, t)^5*(diff(f(x, t), t))*b+15*f(x, t)^5*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))*a+6*f(x, t)^5*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), t), t), x), x))*b+f(x, t)^5*(diff(diff(f(x, t), t), t))*(diff(diff(diff(diff(f(x, t), x), x), x), x))*b+8*f(x, t)^5*(diff(diff(f(x, t), t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))*b+6*f(x, t)^5*(diff(f(x, t), x))*(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x))*a+4*f(x, t)^5*(diff(diff(diff(f(x, t), t), t), x))*(diff(diff(diff(f(x, t), x), x), x))*b-2*f(x, t)^4*(diff(f(x, t), t))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*b-6*f(x, t)^4*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), t))*b-24*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))^2*b-30*f(x, t)^4*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))*a-12*f(x, t)^4*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), t), t), x), x))*b+18*f(x, t)^3*(diff(f(x, t), t))^2*(diff(diff(f(x, t), x), x))^2*b+270*f(x, t)^3*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2*a+72*f(x, t)^3*(diff(diff(f(x, t), t), x))^2*(diff(f(x, t), x))^2*b+24*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), t), t), x))*b+120*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), x), x), x))*a-360*f(x, t)^2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4*a-24*f(x, t)^2*(diff(diff(f(x, t), t), t))*(diff(f(x, t), x))^4*b+120*f(x, t)*(diff(f(x, t), t))^2*(diff(f(x, t), x))^4*b-20*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), x), x), x))-2*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))+8*(diff(f(x, t), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), t))*f(x, t)^4-12*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(diff(diff(f(x, t), x), x), x))-12*f(x, t)^4*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))-4*f(x, t)^4*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))+60*f(x, t)^3*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^2+4*f(x, t)^3*(diff(f(x, t), t))*(diff(f(x, t), x))^2*(diff(diff(diff(diff(f(x, t), x), x), x), x))+48*f(x, t)^3*(diff(diff(f(x, t), x), x))^2*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))+40*f(x, t)^3*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^2*(diff(diff(diff(f(x, t), x), x), x))-108*f(x, t)^2*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))^2*(diff(f(x, t), x))^2-48*f(x, t)^2*(diff(f(x, t), t))*(diff(f(x, t), x))^3*(diff(diff(diff(f(x, t), x), x), x))-144*f(x, t)^2*(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^3+180*f(x, t)*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))*(diff(f(x, t), x))^4-(diff(diff(diff(diff(diff(diff(f(x, t), t), t), x), x), x), x))*f(x, t)^6*b-f(x, t)^6*(diff(diff(diff(diff(diff(diff(f(x, t), x), x), x), x), x), x))*a+6*f(x, t)^5*(diff(diff(diff(f(x, t), t), x), x))^2*b+10*f(x, t)^5*(diff(diff(diff(f(x, t), x), x), x))^2*a-30*f(x, t)^4*(diff(diff(f(x, t), x), x))^3*a+120*f(x, t)*(diff(f(x, t), x))^6*a-2*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), t))*f(x, t)^5+2*f(x, t)^5*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(diff(diff(f(x, t), x), x), x), x))-(diff(diff(f(x, t), x), x))*(diff(diff(f(x, t), t), t))*f(x, t)^5-2*(diff(f(x, t), x))*(diff(diff(diff(f(x, t), t), t), x))*f(x, t)^5+4*f(x, t)^5*(diff(diff(diff(f(x, t), x), x), x))*(diff(diff(diff(diff(f(x, t), t), x), x), x))-6*f(x, t)^4*(diff(diff(diff(f(x, t), t), x), x))*(diff(diff(f(x, t), x), x))^2+2*(diff(diff(f(x, t), x), x))*(diff(f(x, t), t))^2*f(x, t)^4-4*f(x, t)^4*(diff(f(x, t), t))*(diff(diff(diff(f(x, t), x), x), x))^2+2*(diff(f(x, t), x))^2*(diff(diff(f(x, t), t), t))*f(x, t)^4-6*(diff(f(x, t), x))^2*(diff(f(x, t), t))^2*f(x, t)^3+6*f(x, t)^3*(diff(f(x, t), t))*(diff(diff(f(x, t), x), x))^3+8*f(x, t)^3*(diff(f(x, t), x))^3*(diff(diff(diff(diff(f(x, t), t), x), x), x))-36*f(x, t)^2*(diff(diff(diff(f(x, t), t), x), x))*(diff(f(x, t), x))^4+72*f(x, t)*(diff(diff(f(x, t), t), x))*(diff(f(x, t), x))^5 = 0

(6)

T := f(x, t) = P[1]^2+P[2]^2+a[7]

T1 := P[1] = t*a[2]+x*a[1]+a[3]

T2 := P[2] = t*a[5]+x*a[4]+a[6]

L2 := subs({T1, T2}, T)

f(x, t) = (t*a[2]+x*a[1]+a[3])^2+(t*a[5]+x*a[4]+a[6])^2+a[7]

(7)

a[1] := 0; a[6] := 0

0

 

0

(8)

L3 := eval(L1, L2)

4*a[4]^2*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])^2*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4+48*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*a[4]^6-240*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4*a[4]^6*a-2*a[4]^2*(2*a[2]^2+2*a[5]^2)*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^5-8*a[4]^2*a[5]^2*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^5-4608*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*(t*a[5]+x*a[4])^6*a[4]^6-384*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*(2*a[2]^2+2*a[5]^2)*(t*a[5]+x*a[4])^4*a[4]^4*b+1920*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])^2*(t*a[5]+x*a[4])^4*a[4]^4*b+32*(t*a[5]+x*a[4])*a[4]^2*a[5]*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4+1152*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*a[4]^4*a[5]^2*(t*a[5]+x*a[4])^2*b-1152*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])^2*a[4]^4*(t*a[5]+x*a[4])^2*b+288*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*a[4]^4*(2*a[2]^2+2*a[5]^2)*(t*a[5]+x*a[4])^2*b+4608*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])*a[4]^6*a[5]*(t*a[5]+x*a[4])^5-24*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4*a[4]^4*(2*a[2]^2+2*a[5]^2)*b-192*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4*a[4]^4*a[5]^2*b+72*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])^2*a[4]^4*b+4320*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*a[4]^6*(t*a[5]+x*a[4])^2*a-11520*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*a[4]^6*(t*a[5]+x*a[4])^4*a+768*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*a[4]^6*a[5]*(t*a[5]+x*a[4])-1728*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*a[4]^6*(t*a[5]+x*a[4])^2-4608*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*a[4]^6*a[5]*(t*a[5]+x*a[4])^3+5760*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*a[4]^6*(t*a[5]+x*a[4])^4+7680*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])*(t*a[5]+x*a[4])^6*a[4]^6*a+8*(t*a[5]+x*a[4])^2*a[4]^2*(2*a[2]^2+2*a[5]^2)*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^4-24*(t*a[5]+x*a[4])^2*a[4]^2*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])^2*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3+1152*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^3*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*a[4]^4*a[5]*(t*a[5]+x*a[4])*b-3072*((t*a[2]+a[3])^2+(t*a[5]+x*a[4])^2+a[7])^2*(2*(t*a[2]+a[3])*a[2]+2*(t*a[5]+x*a[4])*a[5])*a[4]^4*a[5]*(t*a[5]+x*a[4])^3*b = 0

(9)

eqs := {`~`[`=`](PDEtools:-dcoeffs((lhs-rhs)(L3), {t, x}), 0)}

Error, invalid input: PDEtools:-dcoeffs expects its 2nd argument, V, to be of type {`PDEtools/F(x,...)`, And(function,`ODEtools/F(x)`)}, but received {t, x}

 

solve(identity(L3, {t, x}), {a[2], a[3], a[4], a[5], a[7]})

Error, (in unknown) incorrect use of identity(<expr>,<name>)

 

NULL

NULL


 

Download params.mw

params.mw

Hi,

I've just got interested in MapleFlow (and downloaded the trial). I looked at the previous post but I have seen very little activity for MapleFlow. Is it still alive?

Best regards,

Hi, how to rsolve this recursive sequence auto? 

The manual calculation is taking the reciprocal of both sides of the equation, and then do some simple calculations to get the answer:1/(3·2^n-5).

restart; A := 'A'; A[n] := `assuming`([rsolve({A[1] = 1, A[n+1] = A[n]/(5*A[n]+2)}, A[n])], [n::posint]); A[n] := simplify(%); simplify(subs(n = 1, %))

Error, recursive assignment

 

A

 

Error, (in simplify/table) too many levels of recursion

 

Download ask_recrusive_sequence_of_A[n].mw

I want to solve n nonlinear equations in n unknowns named lambda[1], lambda[2], ..., lambda[n]. For simplicity, I assume n=3. Each of the three equations includes matrices and dot products.

Running solve() takes forever...I assume this is because Maple is trying to find the explicit form. I am not sure, but I have reasons to believe that a solution exists.

Questions:

  1. Did I set up the three equations correctly? Is there a simpler or smarter way to set up my problem?
  2. Is it okay for Maple that lambda[i] sometimes appears isolated and sometimes as a component of a vector?  
  3. Is there a way to keep my quantities (except the lambda vector I guess, since it containts the lambda[i] component I need to solve for) undefined/implicit while making sure Maple correctly deals with matrices and dot products?
    e.g., the solution should display terms like p^%T.R.w as such, rather than the corresponding explicit computation...

Worksheet:

restart;

S:=Matrix([[Sigma[1,1],Sigma[1,2],Sigma[1,3]],[Sigma[2,1],Sigma[2,2],Sigma[2,3]],[Sigma[3,1],Sigma[3,2],Sigma[3,3]]]):

w:=Vector(3,symbol=omega):

r:=Vector([sqrt(S[1,1]/(w^%T.S.w)),sqrt(S[2,2]/(w^%T.S.w)),sqrt(S[3,3]/(w^%T.S.w))]):

R:=Matrix([[r[1],0,0],[0,r[2],0],[0,0,r[3]]]):

p:=Vector(3,symbol=rho):

l:=Vector(3,symbol=lambda):

e[1]:=Vector([1,0,0]):
e[2]:=Vector([0,1,0]):
e[3]:=Vector([0,0,1]):

f:=x->x^2:
w__sq:=f~(w):
l__sq:=f~(l):

 

eq1:='(w[1]*((p^%T.R.w)/(2*(l^%T.w__sq)))*(w^%T.S.e[1]))/(w[1]^2*((p^%T.R.w)^2/(4*(l^%T.w__sq)^2))*(w^%T.S.w)+sigma__delta^2*(w[1]^2*(((l__sq)^%T.w__sq)/(4*(l^%T.w__sq)^2)-lambda[1]/(l^%T.w__sq))+1))';

eq2:='(w[2]*((p^%T.R.w)/(2*(l^%T.w__sq)))*(w^%T.S.e[2]))/(w[2]^2*((p^%T.R.w)^2/(4*(l^%T.w__sq)^2))*(w^%T.S.w)+sigma__delta^2*(w[2]^2*(((l__sq)^%T.w__sq)/(4*(l^%T.w__sq)^2)-lambda[2]/(l^%T.w__sq))+1))';

eq3:='(w[3]*((p^%T.R.w)/(2*(l^%T.w__sq)))*(w^%T.S.e[3]))/(w[3]^2*((p^%T.R.w)^2/(4*(l^%T.w__sq)^2))*(w^%T.S.w)+sigma__delta^2*(w[3]^2*(((l__sq)^%T.w__sq)/(4*(l^%T.w__sq)^2)-lambda[3]/(l^%T.w__sq))+1))';

(1/2)*w[1]*(`.`(p^%T, R, w))*(`.`(w^%T, S, e[1]))/((l^%T.w__sq)*((1/4)*w[1]^2*(`.`(p^%T, R, w))^2*(`.`(w^%T, S, w))/(l^%T.w__sq)^2+sigma__delta^2*(w[1]^2*((1/4)*(l__sq^%T.w__sq)/(l^%T.w__sq)^2-lambda[1]/(l^%T.w__sq))+1)))

 

(1/2)*w[2]*(`.`(p^%T, R, w))*(`.`(w^%T, S, e[2]))/((l^%T.w__sq)*((1/4)*w[2]^2*(`.`(p^%T, R, w))^2*(`.`(w^%T, S, w))/(l^%T.w__sq)^2+sigma__delta^2*(w[2]^2*((1/4)*(l__sq^%T.w__sq)/(l^%T.w__sq)^2-lambda[2]/(l^%T.w__sq))+1)))

 

(1/2)*w[3]*(`.`(p^%T, R, w))*(`.`(w^%T, S, e[3]))/((l^%T.w__sq)*((1/4)*w[3]^2*(`.`(p^%T, R, w))^2*(`.`(w^%T, S, w))/(l^%T.w__sq)^2+sigma__delta^2*(w[3]^2*((1/4)*(l__sq^%T.w__sq)/(l^%T.w__sq)^2-lambda[3]/(l^%T.w__sq))+1)))

(1)

 

solve({eq1-lambda[1],eq2-lambda[2],eq3-lambda[3]},{lambda[1],lambda[2],lambda[3]});

Download dot_products.mw

I'm having some trouble finding how to do this: given some field with spacetime indices, say h[mu,nu] and a Lagrangian L of that field, I want to do a n+1 split to find the canonical momentum pi[~mu,~nu]. I can do the splitting fine, but it's calculating the canonical momentum that is the problem. I've tried a couple of different methods using diff and SubstituteTensorIndices but I can't get it to work: I either get the wrong answer or zero (also wrong). Anyone with any ideas or a working example of this, I'd appreciate to hear from!

The XMLTools:-Validate function can be called with a reference to a schema.

I am a bit unsure where the schema location should be, especially if I can have the .xsd file in a library location.

Have done some tests, and as far as I can see, Maple does not give any error messages if the schema is not found at all. So I'm a bit unsure how to check if the schema defintion actually is read at all.

Something has changed in the past weeks (not sure for how long) and I have not figured out what exactly.

Probably it is me doing something differently while editing an input line. Below is an example where the name data is not italic any more

Anybody else experiencing the same?

My interface seetings:

how do i replace equation in maple for istance i want to replace u_t , u_xx, and u_x in the equation "u_t=-au_xx+bu_x"

Hello to whomever this may concern. I am curently trying to estimate parameters and initial conditions for a system of ODEs. I have taken inspiration from the code given by Robert J. Lopez. I am having issues with minimizing the SSE (sum of squared errors) which I named as a procedure using the NLPSolve command from the Optimization package. It returns an error saying there are too many levels of recursion, which I do not understand where it is coming from or how to solve the issue and I was hoping to get an explination on how to solve it.

Below is the code up until the error. Any help or guidance would be greatly appreciated.

Parameter_estimation.mw

 

Hi, I am searching for an exercise in MapleLearn, which I would like to give to my students in the second year of a mathematics degree. The aim is to make them work only with Maple Learn during 15-20 minutes. So the exercise should be rather short. The original idea that I had was to give them the Fibonacci sequence equation to solve : F(n+2)=F(n+1)+F(n). I would like to transform the sequence using a 2x2 matrix, which is classical. Diagonalize the matrix, compute the powers and so on. This is very easy with Maple, but how can I do that with MapleLearn? I could enter the corresponding  2x2 matrix, compute its eigenvalues but that's it. I could not diagonalize it, compute its powers and so on. Maybe someone knows how to do that with MapleLearn? Else I would be happy to discover other exercises with MapleLearn. 

Apart from introducing the software to french universities, the aim is also to make a research experiment (it is an official collaboration with Maplesoft)

when i replace mu=0 in my equation mu[3] and mu[4] also become zero why ?

restart

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)

NULL

B[1] := 0; mu := 0

0

 

0

(2)

U1 := lambda*(1+sqrt(-lambda)*(cosh(xi*sqrt(-lambda))*B[2]+sinh(xi*sqrt(-lambda))*B[1])*`&+-`(sqrt(lambda/(lambda^2*B[1]^2-lambda^2*B[2]^2-mu^2))))*beta[0]*exp(I*(d*x+e*y+f*t^beta/beta))/(B[1]*cosh(xi*sqrt(-lambda))*lambda+B[2]*sinh(xi*sqrt(-lambda))*lambda+mu)

(1+(-lambda)^(1/2)*cosh(xi*(-lambda)^(1/2))*B[2]*`&+-`((-1/(lambda*B[2]^2))^(1/2)))*beta[0]*exp(I*(d*x+e*y+f*t^beta/beta))/(B[2]*sinh(xi*(-lambda)^(1/2)))

(3)

NULL

V1 := -mu[4]*b*lambda^2*(1+sqrt(-lambda)*(cosh(xi*sqrt(-lambda))*B[2]+sinh(xi*sqrt(-lambda))*B[1])*`&+-`(sqrt(lambda/(lambda^2*B[1]^2-lambda^2*B[2]^2-mu^2))))^2*beta[0]^2/((B[1]*cosh(xi*sqrt(-lambda))*lambda+B[2]*sinh(xi*sqrt(-lambda))*lambda+mu)^2*a*mu[3])

Error, numeric exception: division by zero

 

Download tester.mw

The local command is really horribly treated in the manual. Very difficult to figure out heads from tails.

What follows is Maple code to unprotect the metric tensor g_[], and alter the components. that part I figured out by trial and error where the manual fails to clearly describe how.

But, the problem I now have, is that the unprotected metric does not show the tensor "visual matrix" form anymore if I execute  g_[] . So it makes it all unusable.

How do I enter the unprotected and altered g_[] back into maple so it can be used to calculate the Ricci tensor, Einstein tensor and such ? I guess it will have to be protected again etc. ?

Here is the code


restart;
with(Physics);
g_[sc];
g_[];
_local(g_[]);
g_[1, 2] := 1;
g_[3, 3];
Physics:-g_[1, 1];
g_[1, 2];
Physics:-g_[1, 1];
g_[];
for x to 4 do
    for y to 4 do g_[x, y] := Physics:-g_[x, y]; end do;
end do;
g_[];


So the question is. How do you reverse a localalized metric back to being protected after you altered the metric ?

I am sure that this question must be answered. But when I do a search in the "Search" box, I always receive a page where it  is mentioned that there is an error of generating the page.  So this has been the case for more than a month, or maybe more. I don't know why this is happening.  So that is the reason why I ask a question here.

I would like to put an entry in the file maple.ini so that each time I do the command "plot", the size of the plot in the document will be a size=[300,300] in pixels.  That is supposed to give a plot of 8 cm by 8 cm.  So when I print the document, which will be the size of the plot on the printed page.  I want to do that because if I do not modify the size of the plot by default in the document, it takes up too much space in the page.  I want to do it in maple.ini so I don't have to do it manually each time I plot a function.

Thank you in advance for your help.

First 51 52 53 54 55 56 57 Last Page 53 of 2425