Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

There appears to be a bug in Maple 2023 where constants beginning with a Greek letter cannot be saved.  Assigning a constant to a Greek letter works in the Maple interface but does not work when reading an MPL file encoded as UTF-8.  Example:

read "D:/Nikki/docs/Maple/constants.mpl" :

MPL file contents:

alias(asec=arcsec) : # arcsecant

ζX1 := asec(1) : # first airmass zenith angle
ζX2 := asec(2) : # second airmass zenith angle

The constant "ζX1" can be assigned/saved in Maple but shows as unassigned when imported from an MPL file encoded as UTF-8.

Any assistance in how to use Greek letters in MPL files is much appreciated.

Hello everybody, 

So, i managed to use the examples of https://www.mapleprimes.com/questions/233343-Drawing-Of-Complex-Numbers-On-A-Complex-grid- to get to a figure that works for me. If there is a better way, i would love to hear about it. 

I believe the question does not need translation because it is very clear due to the other examples that have been done before example 3. 

The real question is: how to combine these parameters to form the right figure? The parameters are in the red box in the 3rd picture. 


 

#Toegepaste Wiskunde Deel 2
#Hoofdstuk 6
#Paragraaf 1
NULL

#Voorbeeld 1

plots:-inequal([2 <= x and x <= 7, 1 <= y and y <= 4], x = 0 .. 7, y = 0 .. 7)

 

plots:-inequal([1 <= y and y <= 4, 2 <= x and x <= 7], x = 0 .. 7, y = 0 .. 7)

 

#Voorbeeld 2

plots:-inequal([0 <= x and x <= 4, 0 <= y and y <= sqrt(x)], x = 0 .. 5, y = 0 .. 5)

 

plots:-inequal([0 <= y and y <= 2, y^2 <= x and x <= 4], x = 0 .. 5, y = 0 .. 5)

 

#Voorbeeld 3

plots:-inequal([0 <= x and x <= 2, (1/2)*x <= y and y <= 3*x*(1/2)], x = 0 .. 4, y = 0 .. 4)

 

plots:-inequal([0 <= y and y <= 1, 2*y*(1/3) <= x and x <= 2*y], x = 0 .. 5, y = 0 .. 5)

 

Download Toegepaste_Wiskunde_Deel_2_Hoofdstuk_6.mw

Hi all,

I would like to use the command ConvertIn() from the Galois Field package to convert polynomials to field elements.

This is what I've done:

G:=GF(2,4)
                        G := &Fopf;[16]

b:=G:-random()
                          b := T mod 2

s:=G:-ConvertIn(T)
Error, (in ConvertIn) only integer polynomials in T can be converted

Why is this error showing up?? My argument of the command ConvertIn IS an integer polynomial in T.

How can I avoid this error?

Thanks for your help!

BR David

Hello,

I am generating a code to describe how drug dosages (starting from 1 until dosage 2) are metabolized by the body.

The plot looks weird becasue it does not start from C0 = 4. I don't know how to set t=L so that C0 starts at 4. I suspect that I am misssing something in the piecewise function?  

My code is the follwoing:

C0 := 4;
k := 0.2;
L := 1;
C := t -> C0*(1 + exp(-k*L))*exp(-k*(t - L));
L := piecewise(L <= t and t < 2*L, 0, C(t));

plot(C(t), t = 1 .. 10) 

Thanks,

Andre

I'm trying to do a scatter diagram together with a line graph in one. But the function I putting in is not taking it. Is there anyone that can help me with this problem. And thank you.

I know trig functions are hard to simplify in general but was wondering if there was any simplification command that comes in handy for trig functions of the following form.

sin(theta)^(A - 2)*cos(theta)^2 - sin(theta)^(A - 2) + sin(theta)^A

I am trying to make a very rough animation of a rotating pulsar magnetosphere that would very roughly resemble this gif.
I was able to make a sphere and the dipolar field lines. I was wondering if anyone could provide any possible suggestion regarding.

1. How do i animate the rotation such that the rotation axis is not z but slightly misaligned to z.

2. Is there a better way to make the field likes from the poles (jets).

restart

[1.1, cos(theta), theta]

[1.1, cos(theta), theta]

(1)

sph := [.5, theta, phi]; c1 := [A*sin(phi)^2, 0, phi]; c2 := [A*sin(phi)^2, (1/6)*Pi, phi]; c3 := [A*sin(phi)^2, (1/3)*Pi, phi]; c4 := [A*sin(phi)^2, (1/2)*Pi, phi]; c5 := [A*sin(phi)^2, 2*Pi*(1/3), phi]; c6 := [A*sin(phi)^2, 5*Pi*(1/6), phi]

[A*sin(phi)^2, (5/6)*Pi, phi]

(2)

NULL

NULL

NULL

A := 2; plot3d([c1, c2, c3, c4, c5, c6, sph], phi = 0 .. 2*Pi, theta = 0 .. 2*Pi, coords = spherical, plotlist = true, scaling = constrained)

2

 

 

NULL

s1 := [sqrt(theta), 1, 1/10]; s2 := [sqrt(theta), 3, -1/10]; s3 := [sqrt(theta), 2, 1/10]; s4 := [sqrt(theta), 1.7, 1/7]

p := plot3d([sph, s1, s2, s3, s4, s5, s6, c1, c2, c3, c4, c5, c6], phi = 0 .. 2*Pi, theta = 0 .. Pi, coords = spherical, orientation = [-75, 58, 10], scaling = constrained, color = ["SkyBlue", "yellow", "yellow", "yellow", "yellow", "yellow", "yellow", "red", "red", "red", "red", "red", "red"])

``

NULL

NULL

plots[animate](plot3d, [[sph, s1, s2, s3, s4, s5, s6, c1, c2, c3, c4, c5, c6], theta = 0 .. Pi, phi = 0 .. 2*Pi, orientation = [B, 45, 52], coords = spherical], B = -180 .. 180)

NULL

 

Download Pulsar_temp.mw

3. Any other suggestions that could help with the aesthetics.

Maple 2023 still can't verify its own series solution. This happens with Frobenius method when the indicial root are repated. I think I mentioned or reported this long time ago but unable to find a link now.

My question is: Since this looks like it will not be fixed by Maplesoft anytime soon, any one could suggest an alternative method to verify this solution?

restart;

ode:=2*x^2*diff(y(x), x, x) + 2*x*diff(y(x), x) - x*y(x) =0;
sol:=dsolve(ode,y(x),'series');
odetest(sol,ode,'series')

Maple's solution is correct, I solved it by hand and get same solution.

Maple 2023 on windows 10.

I could not find a way to check whether or not a given value belongs to an interval expressed as RealRange (or a union of them e.g. RealRange(Open(-10),-5), RealRange(Open(-6),infinity).

Ideally, I would like a Boolean function that takes two arguments (i) a value v and (ii) a (union of) RealRange r, and returns true/false depending on whether or not value v is in/out of the range r.

hallo  every body 
please i have a simple question 
how do i calculate this sum in maple 18 
thank you 

[6600.0*theta(q)-17000.0*theta(q-1)+14400.0*theta(q-2)-4000.0*theta(q-3) = v(q)-.20*theta(q)+.20*theta(q-1)]

[6600.0*theta(q)-17000.0*theta(q-1)+14400.0*theta(q-2)-4000.0*theta(q-3) = v(q)-.20*theta(q)+.20*theta(q-1)]

(1)

isolate([6600.0*theta(q)-17000.0*theta(q-1)+14400.0*theta(q-2)-4000.0*theta(q-3) = v(q)-.20*theta(q)+.20*theta(q-1)], theta(q))

Error, (in isolate) invalid arguments for isolate

 

solve([6600.0*theta(q)-17000.0*theta(q-1)+14400.0*theta(q-2)-4000.0*theta(q-3) = v(q)-.20*theta(q)+.20*theta(q-1)], theta(q))[]

theta(q) = 0.1515105603e-3*v(q)+2.575709827*theta(q-1)-2.181752068*theta(q-2)+.6060422411*theta(q-3)

(2)

It works for h(x)  in help(isolate)


Download no_isolation.mw

I have access to a powerful GPU, a NVIDIA Tesla A100.

How to leverage it to solve my system of 3 nonlinear equations in 3 variables? Script: 160523_stylized_problem_GPU.mw

I don't think that it is just by adding the following two lines right before my solve() block, right?
CUDA:-Enable(true)
CUDA:-IsEnabled()

(https://www.maplesoft.com/products/maple/features/cuda.aspx)

I also read https://mcsr.olemiss.edu/docs/gpus-on-maple/ and https://www.mapleprimes.com/posts/36412-GPU-Programming-CUDA-OpenCL-And-Maple but it's still not clear to me. Thanks!

When trying to load/import/read an MPL file with Maple 2023, I get an error message of "Error, recursive assignment" which only links to a page which states "There is no help page available for this error."

https://www.maplesoft.com/support/help/errors/view.aspx?path=Error,%20recursive%20assignment

Presumably there is an error in the Maple syntax that prevents execution, but I cannot find a way to debug.  Maple Code Editor Diagnostics says "No Errors" but Maple Code Editor Console says "Error: recursive assignment."  The MPL file is 485 lines but I am assuming the error could be anywhere in the file.  Is there any way to fix this?

Update: I found the error by copying just a handful of lines at a time into Maple until the error was reproduced.  It was a typo of just one character out of 16,960 characters.  But this is very frustrating.  Why can't Maple 2023 tell me what line the error is on?  Is there a better way to debug MPL files?

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

First 199 200 201 202 203 204 205 Last Page 201 of 2218