## How to plot two functions on different domain in t...

I have the function: f(x)= sqrt(x) for x>=0; f(x)=1/x for x<0 and I would like to plot the function.

How exactly to implement this?

## difficulties with textplot and polygonplot...

restart;
with(plots): with(geometry):
_EnvHorizontalName := 'x':
_EnvVerticalName := 'y':
a := 7:
point(E, 0, a*sqrt(3)/2):
point(B, -a/2, 0):
point(C, a/2, 0):
point(o, 0, a*sqrt(3)/6):
point(A, 0, a/2):
point(H, 0, 0):
R := (3-sqrt(3))*sqrt(2)*a/12:
point(J, 0, a*sqrt(3)/6 - R):
triangle(Tr1, [E, B, C]):
triangle(Tr2, [A, B, C]):
StretchRotation(E1, E, B, Pi/4, clockwise, sqrt(2)/2);
coordinates(E1);
StretchRotation(E2, E, C, Pi/4, counterclockwise, sqrt(2)/2);
coordinates(E2);
triangle(Tr3, [E, B, E1]);
triangle(Tr4, [E, C, E2]);
triangle(Tr5, [B, C, J]);
circle(cir, [point(P1,[0,a*sqrt(3)/6]), R]):
poly := Matrix([[0, a*sqrt(3)/2], [-7/4 + (7*sqrt(3))/4, -7/4 + (7*sqrt(3))/4], [0, a/2], [7/4 - (7*sqrt(3))/4, -7/4 + (7*sqrt(3))/4]], datatype = float);
pol1 := polygonplot(poly, colour = "Magenta", transparency = 0.7, gridlines);
tex := textplot([0.2, a*sqrt(3)/2, "zE"], 'align' = {'above', 'right'});
draw([Tr1(color = cyan),
Tr3(color = green),
tex,
cir(color=blue),
Tr2(color = red),
Tr4(color = grey),
Tr1(color=blue)],'view' = [-5 .. 5, 0 .. 7],
axes = normal, scaling = constrained,size=[800,800]);
Error, (in geometry:-draw) the option must be of type equation or name
line(L1, [B, o]);
Equation(L1);
line(L1, -(7*x*sqrt(3))/6 + (7*y)/2 - (49*sqrt(3))/12 = 0);
reflection(J1, J, L1);
triangle(Tr6, [B, J1, E]);
line(L2, [C, o]);
Equation(L2);
line(L2, -(7*x*sqrt(3))/6 - (7*y)/2 + (49*sqrt(3))/12 = 0);
reflection(J3, J, L2);
triangle(Tr7, [C, J3, E]);
triangle(T1, [E, J1, A]);
triangle(T2, [E, C, E2]);
triangle(T3, [B, H, J]);
triangle(T4, [C, H, J]);

draw([cir(color = orange, filled = true, transparency = 0.1), pol1, Tr6(color = blue, filled = true, transparency = 0.2), Tr5(color = blue, filled = true, transparency = 0.2), Tr7(color = blue, filled = true, transparency = 0.2), T1(color = green, filled = true, transparency = 0.2), T2(color = green, filled = true, transparency = 0.2), T3(color = green, filled = true, transparency = 0.2), T4(color = green, filled = true, transparency = 0.2)], axes = none, scaling = constrained);
Error, (in geometry:-draw) the option must be of type equation or name
NULL;
NULL;

Why these errors messages ? Thank you.

## I want the numerical values of the partial derivat...

 > restart:
 >
 > N:=10
 (1)
 > PDE:=diff(u(y, t), t) = diff(u(y, t), y, y)
 (2)
 >
 > ICBC:= {u(y,0) = 0, u(0,t) = cos(t), u(N, t) = 0};
 (3)
 >
 > sol1:=pdsolve(PDE,ICBC,numeric,spacestep=0.025,timestep=0.00001):
 > u_at_1 := sol1:-value(u(y,t), t=0.1);
 > u_at_1(0);
 (4)

## Degrees package and working with RPMs...

Hello im an amateur using maple for my hidraulics machines course at engineering school. Does anyone know how to use the degrees package on maple 2021, i called it out by using the short form. with(Degrees): and then try to use sind but it doesnt recognice the command.

Also when using RPMs the units that should be giving me m/s appear as m2/m*s(radious). I know that the software recognice it as angular speed but is there anyway to eliminate the (radious) so i can work with the speed as lineal?

## ellipse in a rectangle...

with(geometry):
_EnvHorizontalName := x:
_EnvVerticalName := y:
a := 7:
b := a*(1/2 + 1/6*sqrt(45 - 24*sqrt(3)))^2:
r := b*sqrt(b)/(sqrt(a + b) + sqrt(a)):
point(A, -a, b): point(B, -a, -b):
point(C, a, -c): point(F, a, b):
Sq := square(Sq, [A, B, C, F]):
circle(C1, [point(P1, [r, 0]), r]):
circle(C2, [point(P2, [(1 + sqrt(3))*r, r]), r]):
circle(C3, [point(P3, [(1 + sqrt(3))*r, -r]), r]):
ellipse(E, x^2/a^2 + y^2/b^2 = 1, [x, y]):
solve({Equation(C1), x^2/a^2 + y^2/b^2 = 1}, {x, y}):
point(T, [5.349255162, 2.829908743]):
IsOnCircle(T, C1);
draw([E(color = cyan), C1(color = yellow, filled = true), T(symbol = solidcircle, symbolsize = 20, color = red), Sq, C2(color = red), C3(color = red),Sq(color=blue)], axes = normal, view = [-a .. a, -b .. b], scaling = constrained);
square: (196+(7*(1/2+(1/6)*(45-24*3^(1/2))^(1/2))^2+c)^2)^(1/2)-(196+196*(1/2+(1/6)*(45-24*3^(1/2))^(1/2))^4)^(1/2) = 0
square: (7/2)*(1/2+(1/6)*(45-24*3^(1/2))^(1/2))^2-(1/2)*c = 0
Error, (in geometry:-square) not enough information to define a square
false

Error, (in geometry:-draw) cannot determine the vertices for drawing .Why all these errors ? Thank you.

## inner bisector equation...

EqBIS := proc(P, U, V)
local a, eq1, M1, t, PU, PV, bissec1;
a := (P - U)/LinearAlgebra:-Norm(P - U, 2) + (P - V)/LinearAlgebra:-Norm(P - V, 2);
M1 := P + a*t;
eq1 := op(eliminate({x = M1[1], y = M1[2]}, t));
RETURN(op(eq1[2])); end proc;
EqBIS*([4, 5], [11, 7/3], [11, 5]);
why such a procedure gives no result Thabk you.

## Error, `MapleSim` does not evaluate to a module...

I have maplesim installed and I get the error "Error, `MapleSim` does not evaluate to a module" when I run A:=MapleSim:-LinkModel(); in a Maple worksheet.

Does anyone know how to solve this error?

## Plot inequality system...

Hi guys,

I know how to plot inequality system through using with(plots) and inequal term. however, I couldn't plot following system of inqulity equations:

alpha <= 0.0002500000000*(-18000.*m^2 + 47271.*m + 39514. + sqrt(3.24000000*10^8*m^4 - 1.701756000*10^9*m^3 - 4.266980559*10^9*m^2 - 3.036299412*10^9*m - 6.95987804*10^8))/(9.*m^2 + 12.*m + 4.), 0.00005000000000*(-90000.*m^2 + 237237.*m + 198158. + sqrt(8.100000000*10^9*m^4 - 4.270266000*10^10*m^3 - 1.069976858*10^11*m^2 - 7.612670111*10^10*m - 1.744924704*10^10))/(9.*m^2 + 12.*m + 4.) <= alpha, -0.6666666667 < m, m < -0.6665522013

please let me know how we can plot it.

with best

## Why can't Maple solve the below code for some par...

Hello everybody.

I wrote the attached code. this code gives a good result when "H" parameter is equal to 0.5. Even in this condition if I change other parameters, I can get a true result as well. But when I change "H" parameter to 0.6,0.7,0.8 or 0.9 the code can't be executed anymore and I can't get a result.

based on the paper, for H=0.6,0.7,0.8 or 0.9. we have a result.

what's the problem?

1.mw

## Equation of an ellpise after rotation...

restart:local D:
with(plots):
with(plottools):

a := 7;
b := a/2;
NULL;

r := b/2;
c := 2*sqrt(10)*r;
A := [c/2, c/2];
B := [c/2, -1/2*c];
C := [-c/2, (-c)/2];
D := [-c/2, c/2];
f := (x, y) -> x^2/a^2 + y^2/b^2 = 1;
Ell1 := rotate(implicitplot(f(x, y), x = -9 .. 9, y = -6 .. 6, color = blue), Pi/4);

g := (x, y) -> x^2/b^2 + y^2/a^2 = 1;
Ell2 := rotate(implicitplot(g(x, y), x = -10 .. 10, y = -10 .. 10, color = blue), Pi/4);

quadri := plot([A, B, C, D, A], x = -6 .. 6, y = -6 .. 6, filled = true, color = red, transparency = 0.7);
Cir1 := disk([c/2 - r, -c/2 + r], r, color = blue);
Cir2 := disk([-c/2 + r, -c/2 + r], r, color = blue);
Cir3 := disk([-c/2 + r, c/2 - r], r, color = blue);
Cir4 := disk([c/2 - r, c/2 - r], r, color = blue);
display([quadri, Ell1, Ell2, Cir1, Cir2, Cir3, Cir4], scaling = constrained, axes = normal);
how to find the equations of ellipses in order to fit them to discs and whiten the interior of ellipses ? Thank you.

## Difficultu with pointplot...

restart;
with(LinearAlgebra);
with(plots);
NULL;

EqBIS := proc(P, U, V) local a, eq1, M1, t, PU, PV, bissec1; a := (P - U)/LinearAlgebra:-Norm(P - U, 2) + (P - V)/LinearAlgebra:-Norm(P - V, 2); M1 := P + a*t; eq1 := op(eliminate({x = M1[1], y = M1[2]}, t)); RETURN(op(eq1[2])); end proc;
NULL;
A := <4, 8>;
B := <4, 2>;
C := <1, 4>;
EqBIS(A, B, C);
NULL;

Cen := proc(M, N, R) local eq1, eq2, sol; eq1 := EqBIS(M, N, R) = 0; eq2 := EqBIS(N, M, R) = 0; sol := simplify(solve({eq1, eq2}, {x, y})); RETURN([subs(sol, x), subs(sol, y)]); end proc;
Cen(A,B,C):

CircleParm := t -> [(-t^2 + 1)/(t^2 + 1), 2*t/(t^2 + 1)];
P1 := Transpose(convert(CircleParm(1/4), Vector));

P2 := Transpose(convert(CircleParm(5), Vector));
P3 := Transpose(convert(CircleParm(-1/10), Vector));
P4 := Transpose(convert(CircleParm(-3/2), Vector));

C1 := Transpose(convert(Cen(P1, P2, P3), Vector));
Pts := [P1, P2, P3, P4, C1];

display(implicitplot([x^2 + y^2 - 1], x = -2 .. 2, y = -4 .. 2, colour = [blue], scaling = constrained), pointplot(Pts, symbolsize = 16));
Error, (in plots:-pointplot) incorrect number of coordinates in points data
NULL;Why this error. Thank you.

## What is generating the R[0] term, and how do I rem...

I'm using Maple 2021, and my Euler-Lagrange equation is coming out with a R[0] term that should not be there. I want to know what is causing this error and how I can fix it.

Lagrangian

 (1.1)

Time Equation

 (2.1)

## Font-size of caption, labels and axes on plots sma...

Maple 2021.
I expect after exporting to pdf view of page with plot like on fig. 1.

It means with font=[TIMES, ROMAN, 12] as in maine text.
But I obtaned this view of page on pdf (fig. 2) with smaller font-size.

In previouse versions there are no the problem like this.
Dear community, how to fix this annoying inconvenience?

fig. 1.  fig. 2.

## Why does Maple give sometimes erroneous results wh...

Consider the following integral, shown below in this image.

>> Link to the Maple sheet: example.mw <<

Why does Maple provide erroneous results? Is there a bug in the software? I use Maple 2021.