>  restart; Here are the graphs of a parabola and a straight line: >  plots:-display(         plot(x^2, x=-1..1),         plot((x+1)/2, x=-1..1), color=["Red","Green"]); Suppose I want to plot the part of the parabola that lies below the straight line, and suppose, just to be nasty, I choose to do it with implicitplot: >  plots:-implicitplot(y=x^2, x=-1..1, y=0..(x+1)/2); That is not a parabola at all.  [And where does the "ynew" label come from?]

This behavior was introduced in Maple 2022.

In Maple 2021 we get the expected result:

plots:-implicitplot(y=x^2, x=-1..1, y=0..(x+1)/2);

Download mw.mw

I am looking for a more eligent way to convert a Vector to a Diagonal Matrix.

Download 2024-12-26_Q_Diagonal_Matrix_from_Vector.mw

Hello,

This is a simple dynamics example illustrating the behavior of a disc as a pendulum under ideal conditions. Source attached.

Simple_Disk_Pendulum.mw

This is an ode from textbook. dsolve gives new error I have not seen before. 

Maple 2024.2 on windows 10.

Download dsolve_factor_dec_24_2024.mw

tracelast;  gives long output with this at end

#(IntegrationTools:-Indefinite:-Polynomial,14): return poly/primitivepart*thisproc(primitivepart,var)
 IntegrationTools:-Indefinite:-Polynomial called with arguments: (8*I)*x*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*x1*2^(1/2)-(3*I)*2^(1/2)*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*y+(7*I)*2^(1/2)*KummerM(((1/8)*I)*2^(1/2)+7/4, 3/2, I*2^(1/2)*x^2)*y+12*x^2*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*y+8*x*KummerM(3/4+((1/8)*I)*2^(1/2), 3/2, I*2^(1/2)*x^2)*x1+4*KummerM(((1/8)*I)*2^(1/2)+7/4, 3/2, I*2^(1/2)*x^2)*y, x1, nofactor = false
 #(IntegrationTools:-Indefinite:-Polynomial,8): newpoly := factor(poly)
Error, (in factor) too many levels of recursion
 locals defined as: p = p, primitivepart = primitivepart, base = base, exponent = exponent, subpolys = subpolys, change = change, newpoly = newpoly, u = u

Also, this error can not be cought using try/catch. 

I'm working on solving a system of ODEs in Maple that models an epidemic scenario, tracking the number of susceptible, infected, and recovered individuals over time. Here's what I've done so far:
Download sir_model.mw

Created a table of results at daily intervals with the following code:

results := [seq([t = tval, s = round(evalf(sol(tval)[2][1])), i = round(evalf(sol(tval)[2][2])),r = round(evalf(sol(tval)[2][2]))], tval = 0 .. 50)];

printf("%-10s %-15s %-15s %-15s\n", "Day", "Infected", "Recovered");
printf("---------------------------------------------\n");
for entry in results do
    printf("%-10d %-15d %-15d %-15d\n", entry[t],entry[s], entry[i], entry[r]);
end do;

but I encountered an error (in fprintf) integer expected for integer format

P1.mw

Maple's Student:-ODEs:-ODESteps solves an ode by doing change of variable on the independent variable, but the resulting ode is wrong and final answer is wrong.

Here is one such example

restart;
ode:=diff(diff(y(x),x),x)*sin(x)^2 = 2*y(x);
Student:-ODEs:-ODESteps(ode):

But this result is wrong. First of all, we can not have both x and t  in the same ode. This is what dchange gives

ode:=diff(y(x),x$2)*sin(x)^2-2*y(x)=0;
tr:={PDEtools:-Solve(t=ln(x),x)};
simplify(PDEtools:-dchange(tr,ode,[t]))

It looks like Student:-ODEs:-ODESteps is trying to solve  sin(x)^2*y'' + 2 y=0 as EULER type ode.

But Euler type ode will look like  x^2*y'' +2 y=0  

It seems to have confused sin(x)^2 with x^2. This change of variable it used only works for EULER type ode with polynomial coefficient, not trig coefficients.

Maple 2024.2 on Windows 10

i already use this method for a lot of equation but this time something not normal hapening what is problem?

Download problem.mw

It is not that this a terribly difficult to work out, but I feel I am probably missing something. I need to check if a 3D point lies on a 3D line. What is a good approach here. I started of with the idea all alpha's are equal. but there are exceptions. See P3 and P4

Download 2024-12-21_Q_3D_point_lies_on_3D_line.mw

Download 1.3-Complex_Numbers_bad.mw

I'm trying to get my problems in standard form  a + bi . Questions 19 - 22 are wrong.

As a Maple beginner, I am now interested in symbolic calculations in Maple. As before, I set a problem from a subject area that interests me in order to learn from professional answers.

Determine all regular square (n;n) matrices (determinant not equal to zero) that are commutable with every regular (n;n) matrix with respect to matrix multiplication.

(I know the solution from long ago.)

2024-12-20_Q_simplification_Question.mw
Solve the general cubic. Apply values and simplify. 

Could someone show how Maple simplifies to the value of X=3? I tried doing it manually and I could not figure it out. 

Also is there a Help assistant to see the setps?

Download 2024-12-20_Q_simplification_Question.mw

In a plane, equilateral triangles D(i) with side lengths a(i)= 2*i−1, i = 1; 2; 3; ... are arranged along a straight line g in such a way that the "right" corner point of triangle D(k) coincides with the "left" corner point of triangle D(k+1) and that the third corner points all lie in the same half-plane generated by g. Determine the curve/function on which the third corner points lie!

Given the center x1 of a circle in R^2 with radius d12, and a point p2 on the circle, so that d12=||p2-x1||, denote the points on the line segment from x1 to p2 as x1(t) = x1+t*v12, with t=0..d12, and v12 =( p2-x1)/d12.  I want to animate the points x1(t) moving along the line segment from x1 to p2 and the corresponding circles of decreasing radius, with center x(t) and radius d12-t, so that p2 remains on the circle.

I can animate the points along the line segment from x1 to p2 using ‘style=point, symbol=solidcircle’.

I would like to use plottools-circle, to plot the circles. I have also tried the following type commands for the circles of decreasing radius.

Plot([x1(1)+t*v12(1)+(d12-t)*cos(theta)*v12(1)+ (d12-t)*sin(theta)*u12(1), x1(2)+ +t*v12(2)+(d12-t)*d12*cos(theta)*v12(2)+(d12-t)*sin(theta)*u12(2), theta=0..2*PI]

where u12 is a unit vector orthogonal to v12.

I have not been able to combine the two plots into an animation. Thank you