in here i am stuck and i just can check when my function in sinh cosh and sin and cos  but for  other function like tanh,coth,csch,sech,... i  can't apply and do changing variable need a special technique here which  i can't do it i try a lot to change but changing not happened i don't know and the other function will change if  change function in equation (M) and change for sech and other it will make a big problem also for this function already exist  i dont know how hundle this 

p-f.mw

I have generated the regional plot, but the boundaries between the regions are not clearly marked. How can I highlight the boundaries using a black line or a black dashed line, similar to the example image? What syntax should I use? I have attached my generated image file.

Regional_new.mw

how we solve x^(4)-12*x-12=0, 

mapple sent after solving: RootOf(_Z^4 - 12*_Z - 12, index = 1), RootOf(_Z^4 - 12*_Z - 12, index = 2), RootOf(_Z^4 - 12*_Z - 12, index = 3), RootOf(_Z^4 - 12*_Z - 12, index = 4)

sol:=ln( (y-1)^(1/3)* (y^2+y+1)^(1/3) ) - ln(y) = 2/5* ln(t^2+1)+_C1;
solve( sol,y);

in real domain is fine also. But all my attempts failed. I waited 3-4 minutes each time and stopped it.

Any one can find a trick? Below worksheet showing my attempts and also solution by Mathematica which took 0.3 seconds

Make sure to save all your work first. This problem is known to crash Maple !

Download solve_problem_march_7_2026.mw

I have attached a MapleFlow file that shows a result I don't quite understand.  For some reason the FinalValue result that I get using the DynamicSystems:-StepProperties is different than the steady state value from a DynamicSystems:-ResponsePlot to a unit step using a DifferentialEquation object.

You will notice that I played around a bit, including transforming the DiffEquation object to a TransferFunction object, to see if I could make sense of it.  All of the results are the same as the StepProperties result, it appears that only the ResponsePlot for the DifferentialEquation object is different.

Am I doing something wrong?  Does anyone have any ideas?

Dynamic_Systems_Final_Value.flow

How to solve the following system of equations?

restart;
el1 := a[0] = 1;
el2 := a[0]+a[1]+a[2]+a[3]+a[4]+a[5]+a[6]+a[7]+a[8]+a[9]+a[10]+a[11]+a[12]+a[13]+a[14]+a[15] = 3/2+sinh(1);
el3 := a[1] = 1;
el4 := a[1]+2*a[2]+3*a[3]+4*a[4]+5*a[5]+6*a[6]+7*a[7]+8*a[8]+9*a[9]+10*a[10]+11*a[11]+12*a[12]+13*a[13]+14*a[14]+15*a[15] = 1+cosh(1);
el5 := 24*a[4]-(2*(1+c))*a[2]+c*a[0] = -1;
el6 := 24*a[4]+120*a[5]+360*a[6]+840*a[7]+1680*a[8]+3024*a[9]+5040*a[10]+7920*a[11]+11880*a[12]+17160*a[13]+24024*a[14]+32760*a[15]-(1+c)*(2*a[2]+6*a[3]+12*a[4]+20*a[5]+30*a[6]+42*a[7]+56*a[8]+72*a[9]+90*a[10]+110*a[11]+132*a[12]+156*a[13]+182*a[14]+210*a[15])+c*(a[0]+a[1]+a[2]+a[3]+a[4]+a[5]+a[6]+a[7]+a[8]+a[9]+a[10]+a[11]+a[12]+a[13]+a[14]+a[15])-(1/2)*c = -1;
el7 := 120*a[5]-(6*(1+c))*a[3]+c*a[1] = 0;
el8 := 120*a[5]+720*a[6]+2520*a[7]+6720*a[8]+15120*a[9]+30240*a[10]+55440*a[11]+95040*a[12]+154440*a[13]+240240*a[14]+360360*a[15]-(1+c)*(6*a[3]+24*a[4]+60*a[5]+120*a[6]+210*a[7]+336*a[8]+504*a[9]+720*a[10]+990*a[11]+1320*a[12]+1716*a[13]+2184*a[14]+2730*a[15])+c*(a[1]+2*a[2]+3*a[3]+4*a[4]+5*a[5]+6*a[6]+7*a[7]+8*a[8]+9*a[9]+10*a[10]+11*a[11]+12*a[12]+13*a[13]+14*a[14]+15*a[15])-c = 0;
el9 := 720*a[6]-(24*(1+c))*a[4]+2*c*a[2]-c = 0;
el10 := 720*a[6]+5040*a[7]+20160*a[8]+60480*a[9]+151200*a[10]+332640*a[11]+665280*a[12]+1235520*a[13]+2162160*a[14]+3603600*a[15]-(1+c)*(24*a[4]+120*a[5]+360*a[6]+840*a[7]+1680*a[8]+3024*a[9]+5040*a[10]+7920*a[11]+11880*a[12]+17160*a[13]+24024*a[14]+32760*a[15])+c*(2*a[2]+6*a[3]+12*a[4]+20*a[5]+30*a[6]+42*a[7]+56*a[8]+72*a[9]+90*a[10]+110*a[11]+132*a[12]+156*a[13]+182*a[14]+210*a[15])-c = 0;
el11 := 5040*a[7]-(120*(1+c))*a[5]+6*c*a[3] = 0;
el12 := 5040*a[7]+40320*a[8]+181440*a[9]+604800*a[10]+1663200*a[11]+3991680*a[12]+8648640*a[13]+17297280*a[14]+32432400*a[15]-(1+c)*(120*a[5]+720*a[6]+2520*a[7]+6720*a[8]+15120*a[9]+30240*a[10]+55440*a[11]+95040*a[12]+154440*a[13]+240240*a[14]+360360*a[15])+c*(6*a[3]+24*a[4]+60*a[5]+120*a[6]+210*a[7]+336*a[8]+504*a[9]+720*a[10]+990*a[11]+1320*a[12]+1716*a[13]+2184*a[14]+2730*a[15]) = 0;
el13 := 40320*a[8]-(720*(1+c))*a[6]+24*c*a[4] = 0;
el14 := 40320*a[8]+362880*a[9]+1814400*a[10]+6652800*a[11]+19958400*a[12]+51891840*a[13]+121080960*a[14]+259459200*a[15]-(1+c)*(720*a[6]+5040*a[7]+20160*a[8]+60480*a[9]+151200*a[10]+332640*a[11]+665280*a[12]+1235520*a[13]+2162160*a[14]+3603600*a[15])+c*(24*a[4]+120*a[5]+360*a[6]+840*a[7]+1680*a[8]+3024*a[9]+5040*a[10]+7920*a[11]+11880*a[12]+17160*a[13]+24024*a[14]+32760*a[15]) = 0;
el15 := 362880*a[9]-(5040*(1+c))*a[7]+120*c*a[5] = 0;
el16 := 362880*a[9]+3628800*a[10]+19958400*a[11]+79833600*a[12]+259459200*a[13]+726485760*a[14]+1816214400*a[15]-(1+c)*(5040*a[7]+40320*a[8]+181440*a[9]+604800*a[10]+1663200*a[11]+3991680*a[12]+8648640*a[13]+17297280*a[14]+32432400*a[15])+c*(120*a[5]+720*a[6]+2520*a[7]+6720*a[8]+15120*a[9]+30240*a[10]+55440*a[11]+95040*a[12]+154440*a[13]+240240*a[14]+360360*a[15]) = 0;

solve({el1, el2, el3, el4, el5, el6, el7, el8, el9, el10, el11, el12, el13, el14, el15, el16}, {a[0], a[1], a[2], a[3], a[4], a[5], a[6], a[7], a[8], a[9], a[10], a[11], a[12], a[13], a[14], a[15]});

i need it in  numerical format

I never used Maplets before.

Been learning Explore, which is OK, but the UI does not look good. Too much white spaces between sliders.  So been looking at Maplets to use instead.

Should one just use Maplets to make interactive demos with sliders, buttons, popup menus and so on or use Explore? Which is better?

I do not waste more time learning Explore more if Maplets is a better choice.

If someone here knows Maplets, here is something I just wrote in Explore.  Could this be coded in Maplets to see how it will look like. It took me 30 minutes to make it using Explore and I am no expert in Explore, so hopefully it should not take someone who knows Maplets much time to produce same thing as Maplet. I wanted to see if the UI will look better or same issues with too wasted spaces between sliders.  I assume with Maplets, there is same concept as Explore, with range of variables, and initial values and so on...

I understand one can run Maplet using Maple viewer without needing to have Maple installed on the PC, so this is an advantage.

Download lotka_volterra.mw

Screen shot of the UI

How would the above look using Maplets instead of Explore?

I am not sure whether this should be in "create a post" or "ask a question." Let me know if this is more appropriate in "ask a question"

The standing wave equation is given by:

PDE:=diff(u(x,t),t,t)=c^2*diff(u(x,t),x,x)

IBC := {u(x, 0) = A0*cos(Pi*x/L), D[1](u)(0, t) = 0, D[1](u)(L, t) = 0, D[2](u)(x, 0) = 0}

pdsolve(PDE,IBC,numeric)

In my problem, the material is a magneto-elastic material where c, the speed of the acoustic wave, is a function of a magnetic field H.  The material is nonlinear and saturable.  I define it by a 3 segment piecewise nonlinear function of H.  The material response is a result of a sinusoidal H field.  I am interested in solving u(x,t).  

With that, c in the PDE has to be rewritten as c(H(t)), pdsolve gives an error as PDE has to be expressed as a function of u,t, or x.  So I redefine c(H(t)) as c(t).  I ran into another error in pdsolve as the piecewise has to be based on t or x and not H.  

The problem is that depending on H I can go through all 3 segments and back in one cycle and I have to find the corresponding t's for the piecewise.  Now I am driving the material with 100 cycles, I have to find and list all those piecewise transition points which is hardly practical.

Is there any other ways to approach and solve this problem?

Hello Friends.

I have used Maple to create a Polynomial Regression model.  The model is called "PRModel."  It works fine.  The independent variable is "X" and the dependent variable is "Y."  Both X and Y are of the vector variety.  The model is as follows:

PRModel:=PolynomialFit(10, X, Y, summarize = true);

I would like an estimate of Y for each value of X.  I would like the estimates of Y to be in a variable called "estY."   I have not been successful with this.  I have tried many different variations of the following, but have not been successful.  

estY := eval(PRModel, X);

Any suggestions as to how I can capture the estimates of Y? 

Thank you.

i recently noticed that i can no longer access maple transactions, usual link https://mapletransactions.org/index.php/maple/index automatically opens https://scholarsportal.info/, and in that web site, some pages do not open, stating

"The page you are trying to access is restricted to IP addresses from our member universities. Please try accessing this page from on-campus or through a VPN."

this seems like a geo ip restriction, i search for an announcement for this change but i couldn't any, could you please clarify situation?

thanks in advance..

On my journey of discovery through the Maple world, I now want to try out Maple's convenient features in the complex plane, something that used to be laboriously worked out and demonstrated on the blackboard with chalk. I couldn't find a suitable introduction in the help text. I'm interested in whether a package needs to be loaded and how to handle polynomials, series, and line integrals (I have a reasonable understanding only of the theory).

once i founded but  i lost the technique and i can't remember how i can reach the point how to find thus parameter and find the solution of pde

t1.mw

I just found a rather incredible error in versions 2025.1 and 2025.2.

The function works fine in version 2022.

The drag function in solve does not work, even for solving a very simple problem.

For example, in the equation 1 + 5x = 3x - 1, if I try to do 3x from right to left, a box opens with a good suggestion (subtract   1+2x = -1)

# subtract 3*x from both sides
((1)) + ~(3 * x)
we can see that the correct expression should be ((1)) - ~(3 * x)

I finally get:

1 + 8x = 6x -1 (it added and did not subtract 3x from both sides)

Please correct this...

I would like to plot multiple curves in a single graph with proper legends for different values of θ, similar to the sample figure. However, I am encountering an error while generating the legend. What modifications should I make to correct this? Additionally, how can I display labels such as θ=0.4, θ=0.5 etc., directly inside the plot as shown in the sample?

Graph_legend_error.mw

These are the original set up to solve for second order partial differetial equations:

E := 1.456*10^11;
rho := 7900;
a := sqrt(E/rho);
L := 0.03639;

f0 := 1/(2*L)*a;
f := f0 - 50;
A0 := 11.75*L/2*10^(-6);

d := 0.8*1000;
T0 := 5/f;
NULL;
PDE := diff(u(x, t), t, t) = a^2*diff(u(x, t), x, x) - d*diff(u(x, t), t);
IBC := {u(x, 0) = A0*cos(Pi*x/L), D[1](u)(0, t) = 0, D[1](u)(L, t) = 0, D[2](u)(x, 0) = 0};
pds := pdsolve(PDE, IBC, numeric, timestep = 2*10^(-7));

 I know it solves properly because I can plot it using pds:-plot without problem.  The question is how to extract the value at specific x and t.  I started to try some of the more intuitive way such as u(0,0) without getting an answer.  I tried a number of eval expression such as eval(pds,[x=0,t=0]), eval(u(x,t),[x=0,t=0]) and so on without success.  

Then I tried using :-value as mentioned in an example from the internet, also resulting in an error.

val:=pds:-value(x=0,t=0)

It seems to be simple.  I guess I am too new with maple so be gentle.  Thank you.