Dear Maple experts, I have written a not-very-complicated Maple procedure, including eval, if, etc. However,

Download SiS_Quadratic_(goal).mw

it faces an error saying that Maple cannot determine if the procedure is true. Would you please have a look and advise?

I would like to plot graph of 2d function f(x, y) = (-ax/(1+y^2), x+by) where x from -5 to 5 and y from -5 to 5 and parameter bar of a and b. Thank you 

Dear Maple experts

I face an issue in using ‘implicit plot’ or ‘inequal’ commands in the attached file.

I use three Maple procedures to make some simple calculations and draw an output. It is Firm Profit versus ‘alpha’. This graph is correct and makes perfect sense.

Then, I use ‘inequal’ command to compare the same outputs when ‘alpha’ and ‘delta’ change. There is a conflict between these plots and the above diagram. For example, P2 plot says

Download ConflcitInequal.mw

SiS is superior to FC in small ‘alpha’s. But in the initial diagram, it is obvious that it should be the opposite!

Perhaps I am making a mistake in using commands. Would you please help me?

Dear Maple experts, I want to write text in a specific location within a Maple plot like p<c. To do that, I use 'plots:-textplot' as follows:

plots:-textplot([[0.1, 0.80, p < c, color = black]], font = [times, 11], rotation = 'horizontal');

The issue is that in my code, 'p' and 'c' are parameters with specific values. Therefore,  in my plot, their values appear instead of p<c. Even if I use 'p<c' or "p<c" it does not work.

Would you please guide me on how to solve the issue?

I have tried to translate the Mathematica-Code of OEIS (A219954):

with MmaTranslator. In the translator occurs an error message with the IF-Statement ...

I will be happy to have this in Maple. I am interested in the digitCount in the above Code.

Sequence is:

Hello Everyone;

I need help to solve the following system using rk-4 method

restart;
NULL;
NULL;
C := 1.0;
gK := 36.0;
gNa := 120.0;
gL := 0.3;
VK := -77.0;
VNa := 50.0;
VL := -54.0;
III := 20;
alpha_n := 0.01*(v(t) + 55.0)/(1 - exp(-1.0/10.0*v(t) - 11.0/2.0));
beta_n := 0.125*exp((-v(t))/80.0 + (-1)*13.0/16.0);
alpha_m := 0.1*(v(t) + 40.0)/(1 - exp(-1.0/10.0*v(t) - 4.0));
beta_m := 4.0*exp(-1.0/18.0*v(t) - 65.0/18.0);
alpha_h := 0.07*exp((-1)*(v(t) + 65.0)/20.0);
beta_h := 1/(1.0 + exp((-v(t) + 35.0)/10.0));

dsys1 := {diff(h(t), t) = alpha_h*(1 - h(t)) - beta_h*h(t), diff(m(t), t) = alpha_m*(1 - m(t)) - beta_m*m(t), diff(n(t), t) = alpha_n*(1 - n(t)) - beta_n(t), diff(v(t), t) = III - gK*n(t)^4*(v(t) - VK) - gNa*m(t)^3*h(t)*(v(t) - VNa) - gL*(v(t) - VL), h(0) = 0.9996937394, m(0) = 0.02890553447, n(0) = 0.2445865495, v(0) = -70};

Dear Maple experts,

I have a system of several nonlinear equaitons. My code can solve it for a given parameter. But when I want to plot it, it takes too much time with no results. So, I decided to plot it for several given points. I get the answer for the points individually, but I don't know how to apply this to 'plot' command. Would you please help?

Download Reza_Tchno_(v3).mw

Besdies, it seems for some points there is no feasible solution. So, there is no need to have plot for them.

I attach the code. Thank you so much!

Hello Everyone;

Can anybody help me to solve the analytical solution of the following ODE?
AAS := C*diff(y(t), t) + (-B0*y(t)^3 - B1*y(t)^2 - B2*y(t) - B3);y[0] := 1

C,B0,B1,B2 are constants.

Dear Maple experts,

I am running a Maple code with several lines/ curves in a plot.

For two of the curves, I use style=pointline. But in the legend the symbols on the lines are not shown. Can you advise how we can make them appear?

I have attached the file. Thank you so much! 

Download LegendPlot.mw

Hello Everyone;

I need to solve the following nonlinear ODE

C*diff(y(x), x) + (-B0*y(x)^3 - B1*y(x)^2 - B2*y(x) - B3) = 0, y(0)=B4

where B0,B1,B2,B3 and B4 are constants. I am trying in Maple 2021, but receiving solution in the form of integral. Is that any other ways that I will be able exact solution. Maple sheet is atatched. I am waiting for your kind respose.

Thanks

Question1.mw

Download Question1.mw

Download HmaxProc.mw

I am currently unable to type in maple 2021, there is no cursor, nothing shows up when I type etc. I am also unable to save files. I went into document mode, nothing. Some of the dialogs don't show up. Where does maple save documents and how do I change that setting? That might be my problem. Furthermore, there are no tabs visable, like I can't access any documents I open.

HPMthesis.mw

Dear Maple Expeets, In the attached Maple file, I have three procedures that are run smoothly. I am trying to compare outputs from each procedure to see which one if larger than the others. For changes on two parameters (alpha and delta), I used implicitplot and inequal commands to identify each region by a specific color. But the outcome is not reasonable. Would you please help with this? Thanks in advance!

Download Compare_three_regions.mw

the final plot does not look reasonable. Would you please guide me?

aRandStep2D := proc(X0, Y0, dx, dy)
  local X, Y, P, R;
  P := Array(1 .. 2);
  R := rand(1 .. 8)();
  if R = 1 then X := X0 - dx; Y := Y0 + dy; end if;
  if R = 2 then X := X0; Y := Y0 + dy; end if;
  if R = 3 then X := X0 + dx; Y := Y0 + dy; end if;
  if R = 4 then X := X0 - dx; Y := Y0; end if;
  if R = 5 then X := X0 + dx; Y := Y0; end if;
  if R = 6 then X := X0 - dx; Y := Y0 - dy; end if;
  if R = 7 then X := X0; Y := Y0 - dy; end if;
  if R = 8 then X := X0 + dx; Y := Y0 - dy; end if;
  P[1] := X; P[2] := Y;
  return P;
end proc 

SetStart := proc(b)
  local alpha, R, P;
  P := Array(1 .. 2);
  alpha := rand(1 .. b)();
  P[1] := alpha*b;
  P[2] := alpha*b;
  return P;
end proc 

RandomFactTpq := proc(N, pb, dx, dy)
  local alpha, X, Y, f, P, counter, B, n, T;
  P := Array(1 .. 2);
  counter := 0; f := 1;
  B := floor(evalf(sqrt(N))); #Set maximal searching steps
  T := floor(evalf(sqrt(N))); #For SetStart's use
  P := SetStart(T);
  X := P[1]; Y := P[2];
  while f = 1 and counter < B do  #loop
    n := pb - X - Y;
    f := gcd(N, n);
    if f > 1then break; end if;
    P := aRandStep2D(X, Y, dx, dy); #A random move
    X := P[1]; Y := P[2];
    if X < 1 or Y < 1 or N - pb - 1 < X or X <= Y then
      P := SetStart(T);       # Restart when out of borders
      X := P[1]; Y := P[2];
    end if;
    counter := counter + 1;    #Counting the searched steps
  end do;
  if  f>1  then print(Find at point (X, Y), found divisor = f, searching steps = counter);
  else print(This*time*finds*no*result, test*again!); end if;
end proc

wxbRandWalkTpqNew4.pdf