Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

The purpose of this document is:

a) to correct the physics that was used in the document "Minimal Road Radius for Highway Superelevation" recently submitted to the Maple Applications Center;

b) to confirm the values found in the manual for the American Association of State Highway and Transportation Officials (AASHTO) that engineers use to design and build these banked curves are physically sound. 

c) to highlight the pedagogical value inherent in the Maple language to distinguish between assignment ( := )  and equivalence (  =  );

d) but most importantly, to demonstrate the pedagogical value Maple has in thinking about solving a problem involving a physical process. Given Maple's symbolic mathematics capabilities, one can implement a top-down approach to the physics and the mathematics, working from the general principle to the specific example. This allows one to avoid the types of errors that occur when translating the problem into a bottom up approach, from specific values of the example to the general principle, an approach that is required by most other computational systems.

I hope that others are willing to continue to engage in discussions related to the pedagogical value of Maple beyond mathematics.

I was asked to post this document to both here and the Maple Applications Center

[Document edited for typos.]

Minimum_Road_Radius.mw

Hallo;
the following MAPLE code generates 2-vectors. To collect them into a set K unfortunately does not work. K contains vectors two-fold etc. So K is not a set which I hoped to get. What is wrong?

restart:m:=5;#m Module
with(LinearAlgebra[Modular]):
K:={};
M:=Matrix([[3, 4*168], [4,3]]);
L:=Matrix([[3],[4]]);
#with(LinearAlgebra[Modular]):
for s from 0 to 10 do
a:=Mod(m,(M^s).L,integer);
K:=K union {a};
od:
K:=K;

Gerd

Hi experts

with(plots):
#epsilon:=5*Pi/2:
gammaa:=0.4:
gammao:=0.1:
gammab:=0.002:
omegab:=100:
X:=60:
Omega:=6*Pi:
g:=10:
sigmag:=5*gammaa:
a:=1/(sigmag*sqrt(2*Pi)):

#deltao:=1:
D7:=evalf(a*int(g*Omega*(gammao/2-I*Delta)*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2),Delta=-infinity..infinity)):
D8:=evalf(0.5*a*int(g^2*gammao*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
D9:=evalf(a*int(g^2*Delta*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
b1 := (deltao-2*X^2*omegab*Y/(gammab^2+omegab^2)-D9)^2:

 f:=Y*((gammaa+D8)^2+(b1))=Im(D7)^2+(epsilon-Re(D7))^2:
#P1:=implicitplot(f,deltao=-70..50,Y=0..1,numpoints=100,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[4, 3],linestyle=1);

implicitplot3d(f,Y = 0 .. 1,deltao=-40..40,epsilon=0..0.4,numpoints=100,labels=[Y,deltao,E],tickmarks=[3,3,3], style = surface, color = "Niagara Azure");

The above code to plot the implicitplot3d ,

first the figures are not smooth

second thing is how to make contourplot at epsilon= 5*Pi/2.

I appreciated any comments


 

i wanna solve a system of equations using ''fsolve'', I have found this example:

fsolve({2*x-x*y=0,-y+3*x*y=0},{x,y}, {x=0.1..5,y=0..4});

but i wanna understand how does this work,and what is the signification of term 

thank you for your help.

This exact value of this serie must be calculated with help of a integral and another serie (Maclaurin)

Note: there is Taylor serie and Maclarin is a special case of this Taylor serie. 

The serie is also numerical in Maple calculated , to get a idea of the exact value 

foroum_vraag_3_inleiding.mw 

I enclosed some book taks pdf about alternating series /infinite series and products 
to get a idea what is asked for  

 

blz_62.pdf

blz_63.pdf

blz64.pdf

blz66.pdf

blz_67.pdf

blz68.pdf

 

Let say the procedure MyProc returns two variables as

A,B:= MyProc(…);

However, for programming purpose, I do not want to create an unused variable. Assuming one only needs B, not A, is there a way not to create and assign A?

For information, in Matlab, it is a commonly used technique. The syntax is

[~,B]= MyProc(…);

In Maple, the meaning of ~ is different, but I am looking for something similar.

 

 Today I read the  Problem 5.21 from book Maple and Mathematica, A Problem Solving Approach for Mathematics:

The book uses following codes:

with(LinearAlgebra): n:=9; v:=Vector(n); T:=v->add(v[i]*x^(i-1),i=1..n); P:=Vector(n,symbol=p); Q:=Vector(n,symbol=q); Z1:=T(a*P+b*Q); Z2:=a*T(P)+b*T(Q);
expand(T(a*P+b*Q)-(a*T(P)+b*T(Q)));
evalb(expand(Z1)=expand(Z2));

But I think  these codes just verify rather than prove(n is always specific number). This example is not difficult to prove by hand. but I'd like to know Maple how to prove it. It seems not easy.

 

 

 

 

 

 

 

 

 

 

 

Write a Maple command to define the function f that is given by

f(x)=x^5+2*x^3+3 .

In the box below, enter the part of your Maple command that is to the right of the ":=".

(Do not include the semi-colon (";") at the end.)

 

The answer I typed in is x^5 + 2*x^3 + 3 but it is wrong. Can someone help me what am I missing. 

Hi,

How to plot polygon between solutions of Complex equation

Thanks

ComplexQuestion.mw

Is it possible to connect Maple to

altibase 

I am trying to use task programming but getting incorrect results and taking more time when using a single core.  Please help

TASK_PARALLEL3.mw

Hello everyone,

I would have a question about how to colour a region of complicated inequalities.

I'm curious about the region where my 'lambdas' are all negatives regarding the parameters 'pm2' and 'am2', but I always find difficulties even if I try with implicitplot or inequal commands. Also, I would like to put into the same graph with different colours when they are separately smaller than 0.

I've already tried plenty of things, I hopefully attached some pics by snipping tool about them.  And here is an example of my problem:

p61w:= implicitplot([lambda61 <= 0], pm2 = 0 .. 1, am2 = 0 .. 1)

in61 := inequal(lambda61 <= 0, pm2 = 0 .. 1, am2 = 0 .. 1)

 

display(in61,p61w)

 

 

So, the main issue that it is not just coloured where my condition is true, but elsewhere too. In plot 'in61' as you can see, there is a region approx. pm2=0..0.6 where it shouldn't be coloured since my condition is not true there. In the graph; 'display(in61, p61w)',  it is visible that under the red line should be just coloured however my 'inequal' command coloured it outside too. 

Also, I tried with implicitplot, for example:

p6 := implicitplot([lambda61 <= 0], pm2 = 0 .. 1, am2 = 0 .. 1, filled = true, coloring = [blue, white], transparency = .5)

But the result is similarly wrong. And I'm not sure why it happens and how can I fix it.

(It's okay if for example for these particular lambdas (lambda61,lambda62) never negatives together, but I have more systems and I don't think it works properly)

I read similar questions too but I couldn't find an answer so, I would very appreciate any suggestions.

This is my first post in Maple, so please let me know if I have mistaken something.

Suppose I have an expression or a defined function that uses basic operations(typical stuff)... I have to go and add % to every operation to prevent it expanding so I can see the expression in detail. This really dirties up the formula.

 

Is there any way to automate this?

 

E.g., simple example

 

f := x->3%*x %+ x%^4;

 

 

Also %^ seems to have invalid order of operations, this seems like a bug in maple? 2%^x does not behave like 2^x in complex expressions, I have to put parathesis around it to avoid it grouping with other terms. It seems to do this with multiplication too. This then requires using tons of paranethesis to get normal expressions to behave correctly.

 

What would be real cool if one could automate all this and also have some way to control the expansion amount. E.g., some parameter than one can set to control from no expansion to full expansion(it would work by making inner most expressions inert then working to outer).

Also, it seems when I use the inert operators and make sure the precendence is correct maple evaluates the wrong value! I have used inert before and it seemed right but now I get very large values. There is something seriously buggy with %op.

If I type a differential equation in Maple software, I am getting a numerical value as in the below image. Can anybody help me to rectify the error? Thanks in advance.


 

 

 

I am experimenting with using units especially for when there are awkard conversion factors. 

The document is for motor gearbox torques inertia ration and frequency.

I have documented my specific questions in the worksheet.

I am interested is see different ways of setting this up. Like I would rather setup the formulas at the start, then supply figures.


 

restart

``

with(Units[Natural])

UseSystem('SI')

UseSystem(SI)

(1)

Can I set the units for mass, accel, etc. so I don't have to re enter the default units if I change the figures later?

``

Mass := 3000*'kg'

3000*Units:-Unit(kg)

(2)

Accel := 1.5*'m'/'s'^2

1.5*Units:-Unit(m/s^2)

(3)

Rat := 5

5

(4)

``

Radius := (1/2)*(0.8282e-1-0.6e-2)*'m'

0.3841000000e-1*Units:-Unit(m)

(5)

Torque := Mass*Accel*Radius

172.8450000*Units:-Unit(J)

(6)

````

MotorTorque := Torque/Rat

34.56900000*Units:-Unit(J)

(7)

Jload := Mass*Radius^2

4.425984300*Units:-Unit(kg*m^2)

(8)

Jmtr := 1.42*10^(-2)*'kg'*'m'^2

0.1420000000e-1*Units:-Unit(kg*m^2)

(9)

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

12.46756141

(10)

stiff := convert(48, 'units', 'N'*'m(radius)'/'arcmin', 'N'*'m(rarius)'/'rad', 'symbolic')

518400/Pi

(11)

/rad

518400/Pi

(12)

"(->)"

0.16501e6

(13)

``

Freq := sqrt(stiff/Jload)/(2*Pi)/s   i.e. Hertz

30.73072153*Units:-Unit(1/(kg^(1/2)*m))

(14)

stiff := 4567

4567

(15)

Freq

30.73072153*Units:-Unit(1/(kg^(1/2)*m))

(16)

``


 

Download Units_questions.mw

First 12 13 14 15 16 17 18 Last Page 14 of 1696