Given a system of 2 first order ode's which are autonomous (i.e. the RHS does not explicitly depend on time), one can make phase plot in Maple using DEtools:-DEplot

I am trying to get my Maple output to look like I get with Mathematica. I am getting close, but there are two issues. First, here is what I get in Maple

restart;

ode1:=diff(x(t),t)=2*x(t)*y(t);
ode2:=diff(y(t),t)=1-x(t)^2-y(t)^2;

DEtools:-DEplot([ode1,ode2],[x(t),y(t)],
                           t=0..1,x = -4..4, y = -4..4,arrows=curve,linecolor=red,
                           numpoints =300,arrowsize='magnitude',scene=[x(t),y(t)]);

Here is the same thing in Mathematica

ode1=x'[t]==2*x[t]*y[t]
ode1=y'[t]==1-x[t]^2-y[t]^2

StreamPlot[{2 x*y,1-x^2-y^2},{x,-4,4},{y,-4,4},
    FrameLabel->{{y[t],None},{x[t],None}},BaseStyle->16,ImagePadding->{{50,10},{50,0}}]

Is it possible to get the arrows to be stream lines similar to how they show up? Also, the axes labels in Maple are internal and hard to see. Is it possible to have on the outside like above so they are easier to see?

On a side question: Should one use DEtools:-DEplot or DEtools:-phaseportrait I could never understand when to use one vs. the other.

Maple 2023.2 on windows

Once the terms have been selected, as below

I only need to keep terms that contain a certain function with generic arguments, for example, tanh(x), in this case, that is, I want from A01

-(m2/4) tanh(2k-2ç)+(m2/16) tanh(4k+2ç)+-(m2/4) tanh(2k+2ç)+(3m2/8) tanh(2ç)

and delete the terms

3/16-(3/16)m2^2

What's the best command to make that? In other words, always exclude terms that do not accompany tanh

I am using Maple to compute the principal part around the point 0, which involves obtaining the series expansion terms with strictly negative powers. This is necessary for my work on evaluating infinite integrals that involve complex functions. While I can calculate these expressions manually, I am exploring whether Maple already offers a convenient tool for this purpose. For instance, Maple's built-in Laurent series expansion can be used to obtain the principal part of functions like BesselK(4, x) as 48/x^4 - 4/x^2. Any assistance on this matter is greatly appreciated. Thank you for your help.

I want to solve the following elliptic integral with symbolic inputs a, b, c, and d, and requires a closed-form solution in Legendre's form in terms of these input variables. I want to know how can I solve this using Maple

  int(1/sqrt((a - t)(t - b)(t - c)(t - d)), t = b .. a)

Condition: a>t>b>c>d

Download maple_query.mw

How do I solve sol1 in terms of DC and T (and tau)?

Download dutycycle.mw

For any given function solely and exclusively given by exponentials, as follows (expression A1), I need to select the argument of each exponential and apply it to any function, for example COS(X)

A1 :=  +5

That is, we would have

A1':=+5

I know this result is zero, but this is a simple example. What is the best way to do this?

I have a linear system of 3 equations in 3 variables and have no issues with solving it using solve().

I am having issues with solving it by imposing assumptions on my parameters. The infolevel[solve] now outputs "Entering solver with 6 equations in 3 variables", so I am likely making mistakes in the syntax solve(Eqs, Vars, UseAssumptions) assuming ... first of all, do I even need to do this given the assume statements that I set up on top of my script? How do I make sure that all the assumptions are preserved throughout all the calculations in my worksheet?

Moreover, I don't know why I get SolutionsLost: setting solutions lost flag.

My script: 061123_solving_with_assumptions.mw

Hi All,

I can enter a definite double integral equation in my own workspace. And I have found the double integral tasks in the task browser. But it's completely unclear how to use or specify the examples in the tasks in my own workspace. How do I add the necessary MultiInt command components shown in the task examples to my own workspace-defined integral equation?

Thanks, mitch

How I put three coordinat node in a graph and capture it. I want to display the coordinat node for the same time with the red grapf for the blue and green graph.

grafik_solusi_AIS_ttk_kritis_8(AIS)tLB5.mw
 

Download grafik_solusi_AIS_ttk_kritis_8(AIS)tLB5.mw

I am trying to print all the equations of sphere in the form (x-a)^2 + (y-b)^2 + (z-c)^2-10^2=0?
my code
restart;
n := 0;
for a to 5 do
    for b to 5 do for c to 5 do mysphere := (x - a)^2 + (y - b)^2 + (z - c)^2 - 10^2 = 0;
if a*b*c <> 0 and igcd(a, b, c) = 1 then n := n + 1; L[n] := [a, b, c];
end if;
end do;
end do;
end do;
L := convert(L, list);
map(mysphere, L);
print(mysphere);

I don't get the correct results. How can get correct results?

How to apply the following relationship in maple commands? In this physics this relation it's utilized in effective field theory

I need some command that simplifies the exponentials (or a series of commands) so that there is no power over power.

What's happening (exp(k))^2

What do I need exp(2k)

What would be the most appropriate method? Below is a photo of my code

Where can I find a tutor on Maples ??? I need a tutor on maple i am on graduate engineering  contat me anyone can help 

So I have a complex-valued function f that depends on two real variables x and y. I want to plot its range, i.e. the set of all values f(x,y), where x and y vary over their entire domain (which is a compact subset of R^2). Since I did not find of a direct way to do this, I decided to just compute the output f(x,y) at a finite set of points and plot these values in the complex plane. This gives me a rough idea of what the range should look like. However, this process is very slow (even though it should not be that slow, in my opinion). Here is the code I have used:

f := x -> exp(2*I*Pi*x);
S1 := (x, y) -> Matrix(3, 3, [[f(x - y), 0, 0], [0, 1, 0], [0, 0, f(-y)]]);
S2 := (x, y) -> Matrix(3, 3, [[f(y), 0, 0], [0, f(y - x), 0], [0, 0, 1]]);
C := Typesetting[delayDotProduct](1/sqrt(3), Matrix(3, 3, [[1, 1, 1], [1, f(1/3), f(2/3)], [1, f(2/3), f(1/3)]]), true);
with(LinearAlgebra);
d := LinearAlgebra[Determinant](C . C);
V := (x, y) -> (((S1(x, y)/d^(1/3)) . C) . (S2(x, y))) . C;
tracevals := [];

for i from 0 to 100 do
    for j from 0 to 100 do tracevals := [op(tracevals), Trace(V(0.01*i, 0.01*j))]; end do;
end do;
with(plots);
P2 := complexplot(tracevals, x = -2 .. 3, y = -3 .. 3);

Note that the specific example (i.e. plotting the trace of a unitary matrix V(x,y)) is not really important to this question. We could instead think of any complex-valued function f(x,y), it should not have any impact on the problem. Maple generates the following output:

Even though this output shows quite well the range that the function will take, it takes almost two minutes to generate. Is there a quicker or more elegant way to compute the range of a function?

We just switched to a new Antivirus program (Bitdefender), and we are getting a bunch of security warnings because of outdated software.

I do have  a couple of those for my PC, e.g. some regarding Apache. I don't have Apache installed on my PC, but I think Maple does use it.

So I wonder which policy Maple has regarding updating its open source libraries.

A quick review on Maple Help - copyright information shows a couple of Open Source licenses code which seems to be delivered in old versions.

Bitdefender classifies them with a high security risk.