Hello,

In my the script a determinant of a matrix is set equal to 0 and then solved for. This formula has infinite solutions. So i will find the second root with: Student:-Calculus1:-Roots(Y, N = 0 .. limit, numeric) . This gives me several solutions and then i pick the second one as this is the critical buckling load of the column with the current parameters. 

Everytime i change a parameter i need to rerun a part of the script. I would like to generate a table where i define lets say parameter A and B and then run the script to give me value Ncr. then repeat this a 100 times with diferent parameters to see the influence of the parameters. 

Ive been looking around but i am not able to figure this one out. I think something like this should be possible? Thank you very much in advance!

Has anybody been able to get Maple working on macOS Sonoma 14.4. I know it's not supported yet, but I thought there might be someone who has solved the issue with the app crashing immediately after startup.

The overloaded  procedure here test returns based on 2 lists or 3 lists entered. The two list has a mixed input type with a default value. The default value of "a" can cause a problem if an explicit value in not entered for "a" in foo1. I not sure is the mixed input type is adding to the problem.
By changing the order of the procedures the problem is avoidable here. But this just a simple example. When there are 6 or so procedures it can be very difficult to select a correct ordering.

Is there a way around this, apart from don't have default values?

I could experiment with changing the input order in each proc but that would break up some logical input sequences on me.

Download 2024-02-3_Q_Overload_proc_.mw

I don't why the limit is not providing me the answer; What did I do wrong?

Im trying to replicate a shape in maple but everytime i plot it i get a big red block of line and i don't know why.

g := piecewise(0 <= x and x < 33, -50*x + 33*y, 33 <= x and 38.78 < y, (x - 47.94)^2 + (y - 40.55)^2 - 312.36, 58 <= x and x < 100 and 0 <= y and y <= 38.78, 0.5778*x^2 - 66.08774*x - y + 1887.9165);
implicitplot(g, x = 0 .. 80, y = 0 .. 60, scaling = constrained)

Any help is greatly appreciated

I have been able to sucessfully use the WeibullPlot function of the statistics package.

Is there an output option (similar to summarize = true for other regressions) for the (reference=true) intercept of (X) value (eta) with the 0.632 line and and slope (Beta) of the line, r^2, etc.?

Thanks,
Bill

I am trying to plot the following.

it consist of a line and a circle.

if tried doing it as so:

f := x -> piecewise(0 <= x and x <= 1.588125352, (-1)*0.39*x^2 + 1.459*x - y, 1.588125352 <= x and x < 4, (x - 1.81)^2 + (y - 0.42)^2 + (-1)*0.94^2);

implicitplot(f(x), x = 0 .. 3, y = 0 .. 1.5, scaling = constrained);

every time i get this ugly line in the point where the funktion switches and i can't seem to remove it.

I've tried discont = true, but it doesn't seem to work on implicitplots. 

Hello,

I am using the Weibullplot function and have gotten a good plot out.  The issue is the gradations on the subticks per decade.

If I set it to three (3) I get 0.2, 0.3, 0.4 subgrid correctly, but only 0.6, 0.7, 0.8 between ).5 and 1.0.

If I set it to four (4) - I get the reverse.

Is this a log plot issue or a weibullplot function issue with respct to setting the subticks on the 10ths correcty of a decade span?  I look forward to thoughts.

with(Statistics);
XX := RandomVariable(Weibull(1, 0.6));
AA := Sample(XX, 100);

AA is a sample matrix of 100 Weibull distributed data points.

WeibullPlot(AA, reference = false, style = line, gridlines = true, size = [800, 500], axis = [tickmarks = [default, subticks = 4], color = darkgreen])

WeibullPlot.mw

Thanks,
Bill

Hi. i want to select terms with sum of powers equal to 1 and other terms with sum of powers more than 1 for the following equation.

Download problem.mw

 actually i need an algorithm to select linear terms from nonlinear ones(sum of powers equals to more than one) and vice versa.

It is &ndash; that is not rendered.

Maple 2023

For comparision Maple 2022.

Can someone reproduce this? Anything I can do about it (some settings maybe)?

Both session directly after restart of Maple.

Hi Maple friends,

Id like to create an animation that rotates my 3d plot. Here is the file: https://www.dropbox.com/scl/fi/zop0jpid0ik0a1unxrzzn/Dynamisches_Modell_Final_filled_conts-1.mw?rlkey=eo482ihqw9pw3vn6qolnxlrwa&dl=0

This is the 3d plot:

plot3d([Vth - Vps], E = 0 .. 1000000, T = 0 .. 15, labels = ["E in (EUR)", "T (in Jahren)", "V (in EUR)"], view = 0 .. 180000, labelfont = [times, bold, 12]);

The starting view should be the following (ankle: -137, 74, -2):

After rotation, tId like to see this view (ankle: 180, 0, 180):

Can someone please help me with this issue? That would be so nice!!

Thanks a lot in advance and best regards
Rebekka

 Dear Sir/Madam,

I have faced a problem while using Maple 18. The  warning is the following

Error, (in f) unable to store 'HFloat(1.9648999999999983)-HFloat(500.0)*delta' when datatype=float[8]
Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

#` Part(a)`  deSh := diff(Sh(t),t) = (alpha(h)*Nh) +(omega*Rh(t))- (((betahb*Iv(t))/(Nh) +p+ mu(h))*Sh(t));   deEh :=  diff(Eh(t),t) =  (((betahb*Iv(t))/(Nh) )*Sh(t))-((mu(h)+gamm)*Eh(t));   deIh := diff(Ih(t),t) = (gamm*Eh(t))-(theta+mu(h)+delta)*Ih(t);  deRh := diff(Rh(t),t) = (delta*Ih(t))-((mu(h)+omega)*Rh(t));  deSv := diff(Sv(t),t) =  (alpha(m)*Nv*(1-Nw/(Nm)))-(((betamb*Ih(t))/(Nh) +mu(m))*Sv(t));  deEv:=diff(Ev(t),t) =  (((betamb*Ih(t))/(Nh) )*Sv(t))-((mu(m)+delta(m))*Ev(t));  deIv:= diff(Iv(t),t) = (delta(m)*Ev(t))-(mu(m)*Iv(t));   mu(h):= .000045 :   Nh:=5071126:  betahb:=.750:  p:=.09 :  omega:=.001:  alpha(h):=.5:  gamm:=.1667 :  delta :=.328833:  theta:= .0000002:  alpha(m) :=.1167:    Nv :=40570:  Nw:= 10000:  Nm:=50570:   betamb := .375 :  mu(m) := .02941:  delta(m):= .1428:     u:= .5:   F := dsolve( {deSh, deEh,deIh,deRh,deSv,deEv,deIv,  Sh(0)=8500,Eh(0)=100,Ih(0)=500,Rh(0)=0,Sv(0)=3500,Ev(0)=500,Iv(0)=570},     { Sh(t),Eh(t),Ih(t),Rh(t),Sv(t),Ev(t),Iv(t) } , numeric );  with(plots):  odeplot( F, [[t,Sh(t)],[t,Eh(t)] ,[t,Ih(t)],[t,Rh(t)],[t,Sv(t)],[t,Ev(t)],[t,Iv(t)]], t=0..30, numpoints=100, colour=[black,blue,purple, red,green,yellow,orange] , legend = ["suspectible host", "exposed host", "infected host", "recovered host","suspectible  vector","exposed vector"," infected vector"]);
 

How can a variable roller ensure that the nip is within the arc of contact between the web and the roller? Every time I run the simulation, it keeps failing.

Hello everyone, I am working on a Maple code where I am dealing with a matrix X and performing several operations, so I can subsitute a set of lists. However, I feel Im doing some redundant things and I wonder if there are ways to make it more concise and efficient. More specifically,

1. 'vals' and 'ecs' :
   'vals' gives me a set of lists with the values I'll use in matrix X , e.g., vals := {[-1, 1, 0], [0, -1, -1], [0, -1, 1]}.
   'ecs' creates a set of lists of equations that I can use in 'eval' (or 'subs') function, e.g., ecs := {[[X3_4 = -1, X3_3 = 1, X2_3 = 0]], [[X3_4 = 0, X3_3 = -1, X2_3 = -1]], [[X3_4 = 0, X3_3 = -1, X2_3 = 1]]}.
   But clearly, I'm getting an extra pair of brackets for each list so I will need to use double index in the eval command. Is there a way to avoid this?
2. In the subsitution step, I've tried with both 'eval' and 'subs' and many modifications to avoid what I did in step 1 but without success. Also I need this extra 'Nelem' to get a list of indices so I can substitute all possible lists of values in 'ecs'.
3. Finally, I'm wondering if what I'm getting from 'BoolEigs' is actually what I want. My final goal is to check if the eigenvalues of all possible solution matrices I got in 'Xs'  are nonnegative. So I need to avoid numerical computations and perform this step exactly. Is my code correct for this?
    

X := Matrix([[1, -X3_3/2 - 1/2, 0, -X2_3], [-X3_3/2 - 1/2, -2*X3_4 - 1, X2_3, 0], [0, X2_3, X3_3, X3_4], [-X2_3, 0, X3_4, 1]]);
vars := [X3_4, X3_3, X2_3];

w := A^3 - A;
rootz := RootOf(w, A);
Pols := [(-A^2 + 1)/(3*A^2 - 1), (-A^2 - 1)/(3*A^2 - 1), A*(3*A^2 - 1)*1/(3*A^2 - 1)];

vals := {allvalues(eval(Pols, A = rootz))};
ecs := map(x -> convert(vars = x, listofequations), vals);

Nelem := [seq(k, k = 1 .. numelems(vals))];
Xs := map(x -> eval(X, ecs[x][1]), Nelem);#double index for eval
Xz := map(x -> subs(ecs[x][1], X), Nelem);#same for subs

with(LinearAlgebra);
Eigs := map(x -> Eigenvalues(x), Xs);
BoolEigs := map(x -> map(y -> is(Im(y) = 0) and is(0 <= Re(evalf(y))), x), Eigs);
evalf(Eigs);