MaplePrimes Questions

I have a set of four equations along with specific parameter values provided in the dataset. I need to create a graph where I can dynamically adjust the values of parameters g2 and y within the specified range given in the sheet.
question_plot_slider.mw

Could you provide a maple code that allows for this flexibility, ideally with interactive controls (like the slider buttons shown in the example figure)

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 ?

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