Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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

As an enthusiast and frequent user of unit functionalities, I would like to suggest some improvements.

Example 1: 

Assume a quantity equation w = f(z) where we substitute z = 100*`μm` (and probably other parameters with units) to obtain a numerical result.
After some simplifications the result is

w = 0.3365290139e-4*Unit('m')

w = 0.3365290139e-4*Units:-Unit(m)

(1)

 

A result in this form is the most frequent physical quantity equation with units in use:

A physical quantity denoted by the name w that is equal to a numerical value times a physical unit.

In this form it is also the shortest way to define/name scientific constants, scientific and engineering parameters, and results for variables/unknowns.  


To better interpret the above result a conversion to micron is desirable. This unfortunately does not work with the context panel. Alternatively one could think of mapping a command onto the equation.
(Mapping a command is a good practice to threat the left-hand side and the right-hand side of an equation equally.)

map(convert, w = 0.3365290139e-4*Units:-Unit(m), units, micron)

Error, (in convert/units) unable to convert `1` to `micron`

 

This however does not work, since the left-hand side does not contain the dimension length, but only the name w of the physical quantity.  

A way to convert is:

lhs(w = 0.3365290139e-4*Units:-Unit(m)) = convert(rhs(w = 0.3365290139e-4*Units:-Unit(m)), units, micron)

w = 33.65290139*Units:-Unit(micron)

(2)

For a single result this might be acceptable, but for a list of results this not practical.

In any case not being capable to apply to an equation the same operation to both sides bares the risks of introducing errors.

NULL

Example 2: 

More complicated are results of this kind

NULL

w((1/10000)*Unit('m')) = 0.3365290139e-4*Unit('m')

w((1/10000)*Units:-Unit(m)) = 0.3365290139e-4*Units:-Unit(m)

(3)

A way to convert this to w(100*`μm`) = 33.65290139*`μm`:

w(convert(op(lhs(w((1/10000)*Units:-Unit(m)) = 0.3365290139e-4*Units:-Unit(m))), units, micron)) = convert(rhs(w((1/10000)*Units:-Unit(m)) = 0.3365290139e-4*Units:-Unit(m)), units, micron)

w(100*Units:-Unit(micron)) = 33.65290139*Units:-Unit(micron)

(4)

This is too complicated.

Instead, eval with a manually entered conversion can be used for this

eval(w(100*Units:-Unit(micron)) = 33.65290139*Units:-Unit(micron), Unit('m') = 10^6*Unit('`μm`'))

w(100*Units:-Unit(micron)) = 33.65290139*Units:-Unit(micron)

(5)

However, this is old style: The way we had to work before the introduction of the unit package.

 

In other areas Maples unit conversion functionalities are more advanced.

Plotting with units for example has been improved to the point where unit conversion errors of quantity expressions are almost impossible (in contrast to numerical-value expressions): The user can specify the units to display and the numerical values of the physical quantities are converted accordingly.

 

S1: To catch up with that improvement, conversion of the results as in Example 1 in the form of quantity equations would be desirable as well.

 

Other desirable unit conversions of quantity expressions are

 

• 

S2: Right click (or select and right click) a unit in Maple output and conversion to a new unit with the context panel (try it with one of the above outputs: nothing happens).

• 

S3: Replacement of a unit by a desired unit of the same dimension (Example 2).

• 

S4: Replacement of a unit by conversion to a new equivalent unit (e.g.: "⟦(m kg)/(s^3 A)⟧ to V/(m)") matching the pattern of the unit but not necessarily the pattern of the underlying Unit commands (for this subs can be used already).

 

For all suggestions, it is important that the physical quantity to be converted does not change, i.e. if a unit is changed, the numerical value of the physical quantity is changed accordingly.

 

Download Unit_conversion_in_quantity_expressions.mw

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, 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

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.

How we can calculate solution of this ODE by give assumption to the equation , we have otehr case too, like lambda>0&mu<0

First 70 71 72 73 74 75 76 Last Page 72 of 2233