Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

How to solve an initial value or boundary value fractional differential equation problem by Maple or, a Maple code is required for solving fractional differential initial value or boundary value problem.

 

Hi, I am creating an activity to illustrate various rotations around different axes. My goal is to enhance the display by showing the axes (one portion as a dashed line and the other as a vector) and to insert the texts (Ox, Oy, Oz) in an optimal manner. Any ideas to optimize the rendering? Thank you.

uestionRevolution.mw

Why is a change in the text area component triggering an action three times in a row?

Download ThePostmanKnocksThreeTimes.mw

When using 

with(physics)

I can enter a metric tensor. If it is the Schwarzschild metric, it is easy enough to do

Geodesics(output=solutions)

to get a general geodesic solution. But is it possible, once I have entered the metric, to numerically integrate a specific geodesic and then plot it?

In this paper (see the footnote on the page), the authors pointed out:

We use -v16 as it is faster than the -v2019 for real root isolation on our benchmarks.

Since I don't have Maple 16 at present, I cannot confirm such an assertion personally. However, I'd like to ascertain if old Maple 16 (released in 2012) is still faster in that specialized realm now (i.e., in 2023). Is there any benchmark test that supports this claim (or tells an opposite story instead?)?

Note. Here I just try to justify the authors' statement (which has caused confusion). Incidentally, will this library be a part of Maple in some future version? (I find that some recent papers provide new "powerful"—in accordance with their "Experiment" sections—Maple packages to tackle difficult problems efficaciously, but such functionalities are not built into the Maple kernel as yet. (In fact, if Maple is a single integrated all-in-one technical platform, then I will no longer need to change tools and formats at each stage. :) ))

Is it even possible to use a replacement list as attempted below to convert a z-transform to a finite difference expression?

z-transformed impulse response (`sampletime = 1/20`)

(10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10)*U(z)/z+(exp(1/10)*(-3+10*exp(1/10)-7*exp(1/5))-(10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10))*U(z)/z^2-(-3+10*exp(1/10)-7*exp(1/5))*exp(1/10)*U(z)/z^3 = 8*Y(z)*exp(1/10)*exp(3/10)*exp(1/5)+(8*((-exp(1/5)-exp(3/10))*exp(1/10)-exp(3/10)*exp(1/5)))*Y(z)/z+(8*(exp(1/10)+exp(1/5)+exp(3/10)))*Y(z)/z^2-8*Y(z)/z^3

(10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10)*U(z)/z+(exp(1/10)*(-3+10*exp(1/10)-7*exp(1/5))-(10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10))*U(z)/z^2-(-3+10*exp(1/10)-7*exp(1/5))*exp(1/10)*U(z)/z^3 = 8*Y(z)*exp(1/10)*exp(3/10)*exp(1/5)+8*((-exp(1/5)-exp(3/10))*exp(1/10)-exp(3/10)*exp(1/5))*Y(z)/z+8*(exp(1/10)+exp(1/5)+exp(3/10))*Y(z)/z^2-8*Y(z)/z^3

(1)

deg_mx := 3

3

(2)

Replacements to get finite difference equation

[seq(U(z)/z^i = U[q+1-i], i = deg_mx .. 0, -1), seq(Y(z)/z^i = Y[q+1-i], i = deg_mx .. 0, -1)]

[U(z)/z^3 = U[q-2], U(z)/z^2 = U[q-1], U(z)/z = U[q], U(z) = U[q+1], Y(z)/z^3 = Y[q-2], Y(z)/z^2 = Y[q-1], Y(z)/z = Y[q], Y(z) = Y[q+1]]

(3)

Attempt to replace

algsubs([U(z)/z^3 = U[q-2], U(z)/z^2 = U[q-1], U(z)/z = U[q], U(z) = U[q+1], Y(z)/z^3 = Y[q-2], Y(z)/z^2 = Y[q-1], Y(z)/z = Y[q], Y(z) = Y[q+1]][1], (10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10)*U(z)/z+(exp(1/10)*(-3+10*exp(1/10)-7*exp(1/5))-(10*exp(3/10)-7*exp(1/5)-3*exp(2/5))*exp(1/10))*U(z)/z^2-(-3+10*exp(1/10)-7*exp(1/5))*exp(1/10)*U(z)/z^3 = 8*Y(z)*exp(1/10)*exp(3/10)*exp(1/5)+8*((-exp(1/5)-exp(3/10))*exp(1/10)-exp(3/10)*exp(1/5))*Y(z)/z+8*(exp(1/10)+exp(1/5)+exp(3/10))*Y(z)/z^2-8*Y(z)/z^3)

Error, (in algsubs) cannot compute degree of pattern in z

 

NULL

Download replace_Y_of_z_by_Y_of_q.mw

Dear all

I have an equation, I would like to introduce a variable R[0], how can I use subs in my equaiton 

subs_equation.mw

Thank you

This code fails in Maple 2023.  It used to work in earlier versions.  Did I miss something in the release notes?

restart;

kernelopts(version);

`Maple 2023.0, X86 64 LINUX, Mar 06 2023, Build ID 1689885`

with(plots):

setoptions3d(font=[Times,roman,12]);

arrow(<1,1,1>);

Error, (in plottools:-rotate) invalid arguments for 3-D transformation

 

Hi

How can I obtain the following equation

from

under  the condition  and  at  by Maple and what is the form of ? (All symbols are constant with respect to  except -exp(-  .
1.docx

Greetings,

I am trying to display an array of plots. I did a search and found this: "https://www.maplesoft.com/support/help/view.aspx?path=plot%2Farrayplot".

However, when I try it, I get an error message: "Error, (in plots:-display) element 1 of the rtable is not a valid plot structure".

Here's the code identical to the sample page. If I display(A__1) or display(A__2) individually, it works.

with(plots);
A := Array(1 .. 2);
A__1 := plot(sin(x), x = -Pi .. Pi, axes = framed, labels = [x, sin(x)]);
A__2 := plot(cos(x), x = -Pi .. Pi, axes = boxed, style = point);
display(A);
 

Running Maple 2023 on MacBook Pro Monterey V12.6.5 (latest as far as I know).

Thanks for any assistance.

Hello 

How to remove this error 

3d_plots.mw

The following program:

with(RealDomain):
sol1 := solve(x-1 > 0, x);
sol2 := solve(1/x^2 > 0, x);
sol3 := solve(1/(x^2-1) > 0, x);
sol4 := solve(x^2+1 > 0, x);
sol5 := solve(x^2+1 < 0, x);

Gives the following output:

sol1 := RealRange(Open(1),infinity)
sol2 := {x <> 0}
sol3 := RealRange(-infinity,Open(-1)), RealRange(Open(1),infinity)
sol4 := x
sol5 := NULL

I would like to force some consistency in the output of solve because I need to feed the output to a third-party program. For example, what's the deal with sol2 and sol3? Why are they formatted differently?

Ideally, I would like the previous output to be:

sol1 := RealRange(Open(1),infinity)
sol2 := RealRange(-infinity,Open(0)), RealRange(Open(0),infinity)
sol3 := RealRange(-infinity,Open(-1)), RealRange(Open(1),infinity)
sol4 := RealRange(-infinity,infinity)
sol5 := NULL

How do I achieve this?

I find that no matter how I modify the style in the Maple command, the LaTeX output seems to be always fixed. This is a bit frustrating. 

G:=Graph(6,{{1,2},{1,4},{4,5},{2,5},{2,3},{3,6},{5,6}});
DrawGraph(G);
vp:=[[0,0],[0.5,0],[1,0],[0,0.5],[0.5,0.5],[1,0.5]]:
SetVertexPositions(G,vp):
DrawGraph(G);

Latex(G, terminal, 100, 100)

\documentclass{amsart}
\begin{document}
\begin{picture}(100,100)
\qbezier(12.40,87.60)(12.40,50.00)(12.40,12.40)
\qbezier(12.40,87.60)(31.20,50.00)(50.00,12.40)
\qbezier(12.40,12.40)(31.20,50.00)(50.00,87.60)
\qbezier(12.40,12.40)(50.00,50.00)(87.60,87.60)
\qbezier(50.00,87.60)(68.80,50.00)(87.60,12.40)
\qbezier(50.00,12.40)(68.80,50.00)(87.60,87.60)
\qbezier(87.60,87.60)(87.60,50.00)(87.60,12.40)
\put(12.400000,87.600000){\circle*{6}}
\put(10.194,96.943){\makebox(0,0){1}}
\put(12.400000,12.400000){\circle*{6}}
\put(8.726,3.531){\makebox(0,0){2}}
\put(50.000000,87.600000){\circle*{6}}
\put(50.000,97.200){\makebox(0,0){3}}
\put(50.000000,12.400000){\circle*{6}}
\put(50.000,2.800){\makebox(0,0){4}}
\put(87.600000,87.600000){\circle*{6}}
\put(91.274,96.469){\makebox(0,0){5}}
\put(87.600000,12.400000){\circle*{6}}
\put(89.806,3.057){\makebox(0,0){6}}
\end{picture}
\end{document}

 

DrawGraph(G,
stylesheet=[vertexborder=false,vertexpadding=20,edgecolor = "Blue",
vertexcolor="Gold",edgethickness=3], font=["Courier",10],size=[180,180])

Even when I change their colors (any other thing), they always go their own way. I don't know if there is a setting that can make the LaTeX output more flexible and in line with our expectations after adjustment.

Hello, 

Ive followed the examples, and yet not the result the book gave. The question before this one, same type, and i did get the right answer there. 

Would you have a clue why it is not going as planned? 

Translation: "Through a resistancewire with the measured resistance R= (1.5+-0.5)Ohm goes an electric current with the measured strength I= (2.5+-0.05)A, During a timespan of t=(10+-0.5)s develops an amount of heat Q=RI^2*t. Give the estimate of the maximum relative deviation in Q."

Here i got the right answer of .42, and yet when doing the exact same thing in assignment 4, i got the wrong answers. 

Translation: "Two bodies with mass M and m have a massmidpointdistance r. For the gravitational force between these bodies counts: F=y*M*m/r^2, at which y is the gravitational constant. Give an estimation of the maximum realative deviation in F on ground of the following measurementdata: M=(7.0+-0.05)kg, m=(2.0+-0.02)kg and r=(5.0)+-0.05)m."

I had the answer -0.0029. The books says its wrong, but i cant see what i did wrong. 

Translation of the example: "According to the law of Ohm the currentstrength in the stream(/current)circle(circuit) with tension(voltage) U and resistance R is equal to I=U/R. For a given currentcircle((part of a)circuit) is given: U=(30.0+-0.05)Volt and R=(2.5+-0.03)Ohm. Determine the maximum absolute and relative deviation in I. 

Solution

The maximum absolute deviation in I is

This gives the currentstrength:

The maximum realtive deviation in I is: .... or about 1.4%"

Could somebody tell me what i am doing wrong?

Thank you!

Greetings,

The Function

#opdracht 3

R := 1.5; A := 2.5; t := 10; Q := R*A^2*t

1.5

 

2.5

 

10

 

93.750

(1)

NULL

restart

Q := R*A^2*t

R*A^2*t

(2)

diff(R*A^2*t, R)

A^2*t

(3)

diff(R*A^2*t, A)

2*R*A*t

(4)

diff(R*A^2*t, t)

R*A^2

(5)

NULL

R := 1.5; A := 2.5; t := 10

1.5

 

2.5

 

10

(6)

.5*A^2*t+2*R*A*t*0.5e-1+.5*R*A^2

39.6875

(7)

Q := R*A^2*t

93.750

(8)

39.6875/(93.750)

.4233333333

(9)

#Opdracht 4

restart

M := 7; m := 2; r := 5; F := y*M*m/r^2

7

 

2

 

5

 

(14/25)*y

(10)

restart

diff(y*M*m/r^2, M)

y*m/r^2

(11)

diff(y*M*m/r^2, m)

y*M/r^2

(12)

diff(y*M*m/r^2, r)

-2*y*M*m/r^3

(13)

M := 7; m := 2; r := 5; y := 9.81

7

 

2

 

5

 

9.81

(14)

0.5e-1*y*m/r^2+0.2e-1*y*M/r^2+(-2*y*M*m/r^3)*0.5e-1

-0.1569600000e-1

(15)

F := y*M*m/r^2

5.493600000

(16)

-0.1569600000e-1/(5.493600000)

-0.2857142857e-2

(17)

NULL

Download Mapleprimes_Question_Book_2_Paragraph_5.13_Question_4.mw

First 79 80 81 82 83 84 85 Last Page 81 of 2097