restart;
f := 419*x^2 + 116*x*y - 426*x*z + 78*y^2 - 142*y*z + 133*z^2 - 1604*x - 682*y + 1086*z + 2306;
df_dx := diff(f, x);
df_dy := diff(f, y);
df_dz := diff(f, z);
centrum := solve({df_dx = 0, df_dy = 0, df_dz = 0}, {x, y, z});
                      "maple.ini in user"

             df_dx := 838 x + 116 y - 426 z - 1604

              df_dy := 116 x + 156 y - 142 z - 682

             df_dz := -426 x - 142 y + 266 z + 1086

               centrum := {x = 7, y = 11, z = 13}

It is a quadratic surface, with this point : centrum ? : point symmetry around centrum, needs more info.

eqt1 := u(x, y, 0, t) = 4*alpha*(-lambda[1]*t*(lambda[1]^2 - 3*lambda[2]^2)*alpha + (b*lambda[1] + ((r[1] + r[2])*c)/2 + a)*t + beta*(y*lambda[1] + x))*lambda[2]^2*beta/(lambda[2]^2*t^2*(lambda[1]^2 + lambda[2]^2)^3*alpha^3 + 3*t*lambda[2]^2*alpha^2*(2*(b*t + beta*y)*lambda[2]^4/3 + t*c*(r[2] - r[1])*lambda[2]^3*I/3 + 2*lambda[1]*((((r[1] + r[2])*c)/2 + a)*t + x*beta)*lambda[2]^2 - lambda[1]^2*t*c*(r[2] - r[1])*lambda[2]*I - (2*lambda[1]^3*((b*lambda[1] + ((r[1] + r[2])*c)/2 + a)*t + beta*(y*lambda[1] + x)))/3) + lambda[2]^2*((b*t + beta*y)^2*lambda[2]^2 + t*c*(b*t + beta*y)*(r[2] - r[1])*lambda[2]*I + ((b*lambda[1] + c*r[2] + a)*t + beta*(y*lambda[1] + x))*((b*lambda[1] + c*r[1] + a)*t + beta*(y*lambda[1] + x)))*alpha + beta^2)

Can this function be used now ?

@salim-barzani 
interaktietermen_procedure_21-9-2025.mw

@Alfred_F 

There are some corrections  to make for the procedure FastPursuite 
Handling the ode system in Maple for this procedure can be different on two ways.

As you can see the goat is escaping the wolf ..lol,  message (  try  correct this  and add  coordinat axes)



The procedure approach:

The golden ratio ..



cikels-hexacon_gulden_snede_aklfred_mprimes_12-9-2025.mw

How_plot_and_visualization_2D_and_3D_of_imaginary_part-this_idea_mprimes9-9-2025.mw 

@Alfred_F 
A bijective function cannot be periodic ?

@Alfred_F 
alfred_bolpuzzle_mprimes_14-8-2025.mw

linear__selectie_in_pde_parts_mprimes_11-8-2025.mw

@Andiguys 

nusselt_nummer_mprimes_Question_How_Can_I_Generate_and_Verify_the_results_4-8-2025_-versieA.mw

sol := dsolve([de1, de2, bc], numeric, method = bvp[midrich], approxsoln = [y(x) = 1 - exp(-x), z(x) = exp(-x)], abserr = 0.1e-5, maxmesh = 128)

` Solve using BVP with midrich method (collocation-type)`

Angular Motion Over Time

restart;
grid := proc(M::Matrix) local i, j, Ms, m, n, wks; m, n := op(1, M); Ms := map(convert, M, string); wks := XMLTools:-ToString(_XML_Worksheet(DocumentTools:-Layout:-Table(':-alignment' = ':-center', ':-width' = 20, seq(DocumentTools:-Layout:-Column(':-weight' = 3 + max(map(length, Ms[() .. (), j]))), j = 1 .. n), seq(DocumentTools:-Layout:-Row(seq(DocumentTools:-Layout:-Cell(`_XML_Text-field`("alignment" = "centred", "style" = "Text", Ms[i, j])), j = 1 .. n)), i = 1 .. m)))); streamcall(INTERFACE_TASKTEMPLATE(':-insertdirect', ':-content' = wks)); NULL; end proc;
data := [["Case", "Optimal Pc"], ["μ₁ = 0, μ₂ = 0", "Pc₁"], ["μ₁ = 0, μ₂ > 0", "Pc₂"], ["μ₁ > 0, μ₂ = 0", "Pc₃"]];
grid(Matrix(data));

Is this correct ?


 

Download perbiutation_pde_mprimes_21-7-2025.mw