Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi, here is my question.

test:=v^3;

plot3d(test,u=0..1-v, v=0..1);

How can I  transform this plot in the right triangle into an equilateral triangle? 

So that the three edges of the plot are the same length.

I wanna, squish it, or shear?

any tools i can use? Or can I do it using some tricks to change the function?

 

 

Ever wonder how to show progress updates from your executing code without printing new lines each time?

One way to do this is to use a TextArea component and the DocumentTools package. The TextArea could be inserted from the Components Palette in Maple, or programmatically like so:

restart;

with(DocumentTools):

with(DocumentTools:-Components):

with(DocumentTools:-Layout):

s := "0": #initial text value

T := TextArea(s, identity = "TextArea0"):
xml := Worksheet(Group(Input(Textfield(T)))):

insertedname:=InsertContent(xml)[1,1]: #find the inserted component name in case changed

for i to 10 do #start the demonstration procedure
   Threads:-Sleep(1);
   SetProperty(insertedname,value,sprintf("%d",i),refresh=true);
end do:

Maplets:-Examples:-Message("Done");


Download text-area-update-progress.mw

why the output wont appear eventhough I think the working is right?

y'' + 1/x * u' + u = 6 - 9x + x^2 - x^3

with baundary condition :

u(0) = 0 , and u(1) = 0

 I run the attached code, but after waiting for a few minutes, there is no response, does anyone know the reason?

 It only shows Evaluating!!!!

 sy01.mw

eq1:=( d)/(dt)u+(d^(2))/(dy^(2))u + s*( d)/(dy)u + delta * theta = 0;

eq2:=( d)/(dt)theta + (d^(2))/(dy^(2))theta + s*Pr*( d)/(dy)theta +lambda* exp(theta/(1 +(epsilon*theta))) = 0; 

initial and boundary conditons   

t <=0; u = theta = 0, for 0 <= y  <= 1   

t> 0;  u =0, theta = 0   at  y = 0;  

t> 0;  u =1, theta = 0   at   y = 1  ;

where, s, epsilon, Pr, lambda, delta are arbitrary parameters

Assuming I have an expression

you can know that the derivative with respect to y will vanish when y is a constant y0.

How can we use the mathematical software Maple to simplify this expression?

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.

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