MaplePrimes Questions

Hello everyone, 

 

I'm currently taking a Chaos and Fractal Geometry class that did not have a pre-req of computer programming or knowledge of Maple, but I am apparently the only person in the class who has not taken a class/does not have this knowledge so my professor really flew by teaching the material. My sob story aside, I was hoping you great people would be able to help me out with (what are apparently) really basic procedures and plots while I try to read up on how to actaully do it myself (if anyone has a good supplement to recommend, I'm all ears!). Anyway, here are the one's that I'm stumped on:

  • "Write a Maple proc L that takes as input two lists and returns as output a set containing the terms they have in common (not necessarily in the same position). For example, L([1, 3, x, 5],[2, x, 5, 2, 1, 8]) should return [1, x, 5]."
    I'm thinking this is going to use some form of the intersection command, but I cannot for the life of me figure out how to set up the syntax
  • "Write a Maple program that implements the Euclidean algorithm to compute the gcd(x, y) for any natural numbers x, y (except x = y = 0). Your program should not just compute the gcd, but should do so by implementing the Euclidean algorithm directly (i.e. proc(x,y) RETURN(gcd(x,y)) end; is not what I am looking for here! :) )."
    If only he would've let me do the easy one :(
  • "Write a Maple program to plot the Pythagorean spiral. The program should take a single positive integer n as an argument and plot the spiral until the side of length n is reached"
    Not even a clue on how to approach this one. 

 

Thanks again to any and all help/recommendations you can make. I really appreciate it!

Hi there

I use Maple13 for plotting the surface of revolution of y=x^2 around the line y=1 on the interval [1,2]. But the distance between the surface 

and the vertical axis is ignored while there must be 1 unit distance.please specify the correct command.

Thanks for your cooperation

Regards

M. R. Yegan

Hi, 

What I want is quite simple but I have not seen it implemented anywhere yet. 

I have an explicit map F from the plane into itself. (not an ODE but simply a map) 

Then I have a set X, (something like a line segment or a rectangle in the plane) for which I 

really need to keep track of the images $F(X), F^2(X), F^3(X)$, and so on.

Is there a way to plot the image of a set in Maple? 

Thank you very much for your interests in the question. 

 

 

Hello,

 

I am a new fellow of maple software.

 

I am trying to run an example of Poincare package example found in 

http://www.maplesoft.com/support/help/maple/view.aspx?path=examples%2Fpoincare

But when i  wrote the commands found in this site, i got this error:  Error, unable to match delimiters 

The commands  are the following:

> restart;
> with(DEtools, poincare, generate_ic, zoom, hamilton_eqs);
> H := ((p1^2+p2^2)*(1/2))*((e^(2*q2+2*sqrt(3)*q1)+e^(2*q2-2*sqrt(3)*q1)+e^(-4*q2))*(1/24))-1/8;

 

> poincare('' '', stepsize=.05,iterations=5);
Error, unable to match delimiters 

                                        poincare('' '', stepsize=.05,iterations=5)

I would appreciate any helps.

 

Miguel 

if input x^2*y+x*y+x

how to output [x^2*y, x*y, x]

Hello those in mapleprimes,

 

I have a question.

I opened Help and entered diff, for example, to have the explanation of how to use diff.

Whether it is diff or not, is not the problem I want to ask about now. But, what I want to ask is the following.

I cklicked the WS botton at the tool bar of the Help, so that I created the worksheet of the same contents of the help of the diff.

Then, I chose an expression in Calling Sequence. That is,

the first line shows

diff(f,x1, ..., xj)     on the left, and on the right \frac{d^j}{dx_j \cdots dx_I} f (I wrote this in the TeX ways, though

it is written in 2-D math input, actually)

 

I chose an expression on the right side, and from context menu, which appeared though right clicking of my mouse,

I selected convert to and 1-D math input, through which I changed 2-D math of the worksheet of the help file to 1-D math input.

Then, the expression shown is a strange one:

((ⅆ)^j)/(ⅆx[j] `...` ⅆx[`1`]) f;

 

This is still easier one. As for the third line in the Help, it has a more complicated form as

((ⅆ)^r)/(ⅆx[k]^m ⅆx[j] `...` ⅆx[3] ⅆx[2]^n ⅆx[1]^n) f;

 

Do you know how this way to notate is applied in maple worksheet?

Why does maplesoft use this way to notate expressions limited to the Help pages, if it is limited to the help pages.

 

Thanks in advance.

 

taro

how to composite two system of polynomials in maple

for a o b = identity

how to composite two polynomial in maple

i would like to search a o b = 1

 

composite two polynomial equal identity

I am entering into Maple, but it interprets it as

How do I change this partial derivative to a total one. (Both x and y can be assumed invertible functions of parameter u and their derivatives are nonzero)

I need to find all combinations of a binary list of length 3.  For example: [0,0,0] , [1,0,0] , [0,1,0] , [0,0,1] , [1,1,0] , [1,0,1] , [0,1,1], [1,1,1].  Does anyone have a procedure or can write a loop that does this?

I want to extract information  about an ODE.

For example, having the next ODE:

ode := cos(x)*(diff(y[1](x), x, x))+sinh((1/10)*x+1)*y[2](x)+y[1](x)+exp(-x)-1:

extract information as:

the inputs system are :

{y[1](x),y[2](x)}

the independent variable is:

{x}

 

Some suggestion? I am tried with: indets, DEtools[odeadvisor], whattype, evaluating the function, convert it into a table and using the "indices". Probably the answer is related to dsolve/numeric/process_input

Regards

 

restart; with(Physics); with(Tetrads)

0, "%1 is not a command in the %2 package", Tetrads, Physics

(1)

Physics:-Setup(coordinatesystems = {X}, mathematicalnotation = true)

[coordinatesystems = {X}, mathematicalnotation = true]

(2)

Physics:-Define(Ybar[a], Y[a], GammaT[a, b, c], RicciT[b, c] = %d_[c](GammaT[`~a`, b, a])-%d_[a](GammaT[`~a`, b, c])+Physics:-`*`(GammaT[`~m`, b, a], GammaT[`~a`, m, c])-Physics:-`*`(GammaT[`~m`, b, c], GammaT[`~a`, m, a])+Physics:-`*`(GammaT[`~a`, b, m], GammaT[`~m`, c, a]-GammaT[`~m`, a, c]), RiemannT[k, l, m, n] = %d_[k](GammaT[m, n, l])-%d_[l](GammaT[m, n, k])+Physics:-`*`(GammaT[`~a`, m, l], GammaT[a, n, k])-Physics:-`*`(GammaT[`~a`, m, k], GammaT[a, n, l])+Physics:-`*`(GammaT[m, n, a], GammaT[`~a`, k, l]-GammaT[`~a`, l, k]))

RicciT[1, 3]

%d_[3](GammaT[`~1`, 1, 1]+GammaT[`~2`, 1, 2]+GammaT[`~3`, 1, 3]+GammaT[`~4`, 1, 4])+GammaT[`~4`, 1, 2]*GammaT[`~2`, 3, 4]+GammaT[`~1`, 1, 3]*GammaT[`~3`, 3, 1]+GammaT[`~2`, 1, 3]*GammaT[`~3`, 3, 2]+GammaT[`~4`, 1, 3]*GammaT[`~3`, 3, 4]-GammaT[`~2`, 1, 3]*GammaT[`~1`, 2, 1]-GammaT[`~2`, 1, 3]*GammaT[`~2`, 2, 2]-GammaT[`~2`, 1, 3]*GammaT[`~4`, 2, 4]-GammaT[`~3`, 1, 3]*GammaT[`~1`, 3, 1]-GammaT[`~3`, 1, 3]*GammaT[`~2`, 3, 2]-GammaT[`~3`, 1, 3]*GammaT[`~4`, 3, 4]-GammaT[`~4`, 1, 3]*GammaT[`~1`, 4, 1]-GammaT[`~4`, 1, 3]*GammaT[`~2`, 4, 2]-GammaT[`~4`, 1, 3]*GammaT[`~4`, 4, 4]-GammaT[`~1`, 1, 3]*GammaT[`~2`, 1, 2]-GammaT[`~1`, 1, 3]*GammaT[`~4`, 1, 4]+GammaT[`~1`, 1, 4]*GammaT[`~4`, 3, 1]+GammaT[`~2`, 1, 4]*GammaT[`~4`, 3, 2]+GammaT[`~4`, 1, 4]*GammaT[`~4`, 3, 4]+GammaT[`~1`, 1, 1]*GammaT[`~1`, 3, 1]+GammaT[`~2`, 1, 1]*GammaT[`~1`, 3, 2]+GammaT[`~4`, 1, 1]*GammaT[`~1`, 3, 4]+GammaT[`~1`, 1, 2]*GammaT[`~2`, 3, 1]+GammaT[`~2`, 1, 2]*GammaT[`~2`, 3, 2]+GammaT[`~1`, 3, 3]*GammaT[`~3`, 1, 1]+GammaT[`~2`, 3, 3]*GammaT[`~3`, 1, 2]+GammaT[`~3`, 1, 4]*GammaT[`~4`, 3, 3]-GammaT[`~3`, 4, 3]*GammaT[`~4`, 1, 3]-GammaT[`~1`, 1, 1]*GammaT[`~1`, 1, 3]-GammaT[`~1`, 1, 3]*GammaT[`~3`, 1, 3]-GammaT[`~2`, 1, 3]*GammaT[`~3`, 2, 3]-%d_[4](GammaT[`~4`, 1, 3])+%d_[1](GammaT[`~1`, 1, 3])+%d_[3](GammaT[`~3`, 1, 3])+%d_[2](GammaT[`~2`, 1, 3])

(3)

for a to 4 do for b to 4 do RicciT[a, b] end do end do

Error, (in index/PhysicsTensor) expected summation indices of type symbol, received: 1

 

 

Now, if I type the RicciT from (3) it displays the same result. However,......

 

RicciT[1, 3]

%d_[3](GammaT[`~1`, 1, 1]+GammaT[`~2`, 1, 2]+GammaT[`~3`, 1, 3]+GammaT[`~4`, 1, 4])+GammaT[`~4`, 1, 2]*GammaT[`~2`, 3, 4]+GammaT[`~1`, 1, 3]*GammaT[`~3`, 3, 1]+GammaT[`~2`, 1, 3]*GammaT[`~3`, 3, 2]+GammaT[`~4`, 1, 3]*GammaT[`~3`, 3, 4]-GammaT[`~2`, 1, 3]*GammaT[`~1`, 2, 1]-GammaT[`~2`, 1, 3]*GammaT[`~2`, 2, 2]-GammaT[`~2`, 1, 3]*GammaT[`~4`, 2, 4]-GammaT[`~3`, 1, 3]*GammaT[`~1`, 3, 1]-GammaT[`~3`, 1, 3]*GammaT[`~2`, 3, 2]-GammaT[`~3`, 1, 3]*GammaT[`~4`, 3, 4]-GammaT[`~4`, 1, 3]*GammaT[`~1`, 4, 1]-GammaT[`~4`, 1, 3]*GammaT[`~2`, 4, 2]-GammaT[`~4`, 1, 3]*GammaT[`~4`, 4, 4]-GammaT[`~1`, 1, 3]*GammaT[`~2`, 1, 2]-GammaT[`~1`, 1, 3]*GammaT[`~4`, 1, 4]+GammaT[`~1`, 1, 4]*GammaT[`~4`, 3, 1]+GammaT[`~2`, 1, 4]*GammaT[`~4`, 3, 2]+GammaT[`~4`, 1, 4]*GammaT[`~4`, 3, 4]+GammaT[`~1`, 1, 1]*GammaT[`~1`, 3, 1]+GammaT[`~2`, 1, 1]*GammaT[`~1`, 3, 2]+GammaT[`~4`, 1, 1]*GammaT[`~1`, 3, 4]+GammaT[`~1`, 1, 2]*GammaT[`~2`, 3, 1]+GammaT[`~2`, 1, 2]*GammaT[`~2`, 3, 2]+GammaT[`~1`, 3, 3]*GammaT[`~3`, 1, 1]+GammaT[`~2`, 3, 3]*GammaT[`~3`, 1, 2]+GammaT[`~3`, 1, 4]*GammaT[`~4`, 3, 3]-GammaT[`~3`, 4, 3]*GammaT[`~4`, 1, 3]-GammaT[`~1`, 1, 1]*GammaT[`~1`, 1, 3]-GammaT[`~1`, 1, 3]*GammaT[`~3`, 1, 3]-GammaT[`~2`, 1, 3]*GammaT[`~3`, 2, 3]-%d_[4](GammaT[`~4`, 1, 3])+%d_[1](GammaT[`~1`, 1, 3])+%d_[3](GammaT[`~3`, 1, 3])+%d_[2](GammaT[`~2`, 1, 3])

(4)

 

If I type a different RicciT, then ...

 

RicciT[2, 3]

Error, (in index/PhysicsTensor) expected summation indices of type symbol, received: 1

 

 

The for loop is changing the tetrad definition of RicciT.

 

NULL

 

Download Question_about_tetrads_in_for_loop.mw

Dear friends,

 

I am trying to do an integration which consists of a list of data; let me explain it more:

I have to do

however, I do not have f(x) in the form of a function. f(x) is a list which is:

also, I have x as

How can I do the integration between 1<x<6 ?

 

 

 

 

Hi,

 

I have a problem with creating a function contains n-dimensional independent variables with random function .. The function looks like:

f(x)=sum(U[i] * abs(x[i])^i,i=1..n)

where U is a random between 0 and 1

 

For 1-dimensional problem, it will be like:

f(x)=U*abs(x)

The plot by MATLAB will be:

 

While, for 2-dimensional problem, it will be like:

f(x[1],x[2])=U[1]*abs(x[1]) - (-U[2]*abs(x[2])^2)

 

and the plot with MATLAB will be:

 

The MATLAB code is:

x1min=-5;
x1max=5;
x2min=-5;
x2max=5;
R=1000; % steps resolution
x1=x1min:(x1max-x1min)/R:x1max;
x2=x2min:(x2max-x2min)/R:x2max;

for j=1:length(x1)

    % For 1-dimensional plotting
    f1(j)=rand*abs(x1(j));

    % For 2-dimensional plotting
    for i=1:length(x2)
        fn(i)=f1(j)+rand*abs(x2(i))^2;
    end

    fn_tot(j,:)=fn;

end

figure(1)
plot(x1,f1);


figure(2)
meshc(x1,x2,fn_tot);

 

 

I have used Maple because of its great graphics and animation, and for the mathematical analysis capability. However, to animate the above problem for 1- and 2-dimensional problems, I always fail!

My attempt has been just successed to plot 1-dimensional problem (without animation) .. And it tooks long-time with many codes!

restart;

X := Statistics:-RandomVariable(('Uniform')(0, 1));
R := Statistics:-Sample(X, 1000);

Vect := abs(Vector[column]([seq(0 .. 10, 10/999)]));

V := `~`[`-`](Vect, 5);

with(LinearAlgebra);
Mat := Multiply(V, R);

with(ArrayTools);
f := ArrayTools:-Diagonal(Mat);

F := abs(f);

Axis := Transpose(V);

f5 := plot(Axis, F, style = line);

 

And finally I got the graph:

 

But the problem is:

How can I animate it? Maybe I need to express it with new Maple code?

How can I plot and animate it when n=2?

 

Thanks

 

I would like to define all the geometry parameters in a multibody model in MapleSim in mm.

Consequently, i define the mm unit for a parameter in 3 places :
1)  in the Parameters area of the subsystem 

2) in the Tab Inspector of the subsystem

 

3) in the Tab Inspector of the rigid body frame of the subsystem

 

The problem is that I receive the warning "Possible double conversion unit on a parameter" and the dimension seems to be as 10^-3 *mm soit too much reduced.

I see that in some examples of the library of Maplesim the parameters of the parameters area of a subsystem are defined as real and consequently, it doesn't cause this kind of problem.

However, how can I do if I want to define my parameters in the parameters area of a subsystem as a position in mm ?

Thanks a lot for your help

 

 

First 1224 1225 1226 1227 1228 1229 1230 Last Page 1226 of 2429