I am trying to make a plot that shows the amplitude of a second order linear system relative to an external driving frequency ω, but I am having trouble defining my function. Here's what I have so far.

When I pass the amp function ω=0 and F=4, I don't get back what I expect. I would think that those symbols should evaluate to the values I passed in. I have tried different parameter names etc..

All I am trying to do is plot the function shown above for F=4 and 0≤ω≤5. I know I could manually plug the values in, but I thought I'd ask the proper way to do it. Any help is greatly appreciated. Thanks folks!

Voting is now open for the next individual prize to be awarded as part of the Möbius App Challenge.  The winner will receive an Xbox One Prize Pack! 

Here are the finalist Apps:

• Fixed Point Equation
• Testing Convergence of a Series
• Mod-N-Calulator
• Epsilon-delta Limit Definition

Note that, if you ever have any problems viewing Apps in your browser, or simply want to work offline, you can always download a Möbius App and view it in Maple or the free Maple Player. To download a Möbius App, follow the link to the App and then click on the Download button near the top left of the page.

You can vote for your favorite through our Facebook page or, if you’re not on Facebook, send an email with your vote to Mobius-Project@maplesoft.com.

And remember, we are now accepting entries for the next quarterly prize. You could win a Music Prize Pack, including the 64GB 5th Generation Apple iPod Touch, Sennheiser In-Ear Noise Cancelling Headphones and the Bose SoundLink Bluetooth Speaker III!  See the Möbuis App Challenge for details.

Voting closes April 25^th, 2014.

Voting is now open for the next individual prize to be awarded as part of the Möbius App Challenge.  The winner will receive an Xbox One Prize Pack! 

Here are the finalist Apps:

• Fixed Point Equation
• Testing Convergence of a Series
• Mod-N-Calulator
• Epsilon-delta Limit Definition

Note that, if you ever have any problems viewing Apps in your browser, or simply want to work offline, you can always download a Möbius App and view it in Maple or the free Maple Player. To download a Möbius App, follow the link to the App and then click on the Download button near the top left of the page.

You can vote for your favorite through our Facebook page or, if you’re not on Facebook, send an email with your vote to Mobius-Project@maplesoft.com.

And remember, we are now accepting entries for the next quarterly prize. You could win a Music Prize Pack, including the 64GB 5th Generation Apple iPod Touch, Sennheiser In-Ear Noise Cancelling Headphones and the Bose SoundLink Bluetooth Speaker III!  See the Möbuis App Challenge for details.

Voting closes April 25^th, 2014.

Hi guys!

I want to do a fractalzoom into the mandelbrot fractal. The basic idea was to produce all single plot pictures, and then display/animate the zoom by iterate through the list. But I cant find a solution that really works. Any help would be very appreciated :)

So where am I? A single picture is made like this: (Sorry math tool didnt take it)

plot3d(0, -2 .. 1, -1.5 .. 1.5, orientation = [-90, 0], grid = [250, 250], style = patchnogrid, color = mandelbrot, lightmodel = None):

As u can see, I use the colorfunc property to color the plane. The mandelbrot procedure evalueates how many iterations are needed to go beyond 2 (> 2).

I cant post pictures, so heres the code if u want to reproduce some:

mandelbrot := proc (x, y)

local c, z, iterationCount, position, iterations;

z := 0; iterations := 0; iterationCount := 35;

position := Complex(x, y);

c := evalf(position);

while iterations < iterationCount do

    if 2 <= abs(z) then return iterations;

    else z := z^2+c; iterations := iterations+1;

    end if

od;

return iterations

end:

As you can see i get a list of PLOT3D structures, each with different axis sizes or views. The problem now is when i use the straightforward way:

display(plots, insequence=true);

the animated frames get smaller and smaller, because the frame is set to the same size as the first frame. 

Do you have any idea how I can make the axis somewhat dynamic, or do you see a better approach?

Here's the rest of the code to produce the effect. Its nearly hardcoded, because its still in experimental phase, sorry for that :P

The test does the following: Start with window -2 .. 1, -1.5 .. 1.5 and with linear interpolation zoom in to the window -0.82 .. -0.7, -0.2 .. -0.08.

 v is used to store the view explicitely, is explained afterwards

steps := 3:

d1 := (2-0.82)/steps:

d2 := (1+0.7)/steps:

d3 := (1.5-0.2)/steps:

d4 := (1.5+0.08)/steps:

l := Array(1 .. 4):

v := Array(1 .. 4):

for i from 0 to steps do

    v(i+1) := [d1*i-2 .. -d2*i+1, -1.5+d3*i .. 1.5-d4*i, 0 .. 0];

    l(i+1) := plot3d(0, d1*i-2 .. -d2*i+1, -1.5+d3*i .. 1.5-d4*i, orientation = [-90, 0], grid = [100, 100], style = patchnogrid, color = mandelbrot, lightmodel = None);

od:

display a single picture:

display(l(4));

create the animation:

ll := convert(l, list):

display(ll, insequence = true);

Another approach:

I tried to explicitely give the view as argument:

display(l(i), view = v(i)), v(i) stored in previous render loop. This works, but then the question to resolve is how to draw each frame into the same plot figure and not creating a new plot window for each one?

Now you should have all informations, any help is very welcome!

Thanks in advance,

Cheers,

geischtli

restart;
evalf(Int(cosh(t)/(cosh((17/15)*t)+cosh(t)), t = 0 .. infinity));

5.210624833

evalf(int(cosh(t)/(cosh((17/15)*t)+cosh(t)), t = 0 .. infinity));

2.720144439

I am trying to setup a general metric with the Physics package.  The metric is composed of the Minkowski metric plus the product of two null vectors. Here is the code:

retstart;

with(Physics);

Define(l[mu],eta[mu,nu]);

eta[mu,nu] := rhs(g_[]);

Setup(g[mu,nu]=eta[mu,nu]+l[mu]*l[nu]);

I get the following error:

Error, (in Physics:-Setup) wrong argument: g[mu, nu] = l[mu]*l[nu]+(Matrix(4, 4, {(1, 1) = -1, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (2, 2) = -1, (2, 3) = 0, (2, 4) = 0, (3, 3) = -1, (3, 4) = 0, (4, 4) = 1}, storage = triangular[upper], shape = [symmetric]))

My plan is to apply the rule g[~mu,~nu]*l[mu]*l[nu] = 0 and calculate the Christoffel symbols using the metric.

I am trying to use the Physics package because the DifferentialGeometry package seemed focused on Newman-Penrose.  I will not be using NP for the calculations, only a strict calculation of the Einstein field equations from the given metric.

Thank you.

Hi everybody,

In quantum mechanics, I want to draw the wave  of the electron in the hydrogen atom.  So in 3D, you will have a cricle and along the circle, there will be a cosine function that will meet after a full cercle (the standing wave that is possible).  The blue one and the red one at the same time.  But in 3D

The closest that I come to succes is with this commands:

But as you will see, the wave end at the same place as it started.  But instead of going in a cercle, is as a funny shape.

Thank you in advance for you trouble.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987

AOA... There are three question

1. I want to convert exp(Iota*theta) into ternometric function i.e., 

exp(Iota*theta) = cos(theta)+Iota*sin(theta)

Is there any comand pl help...

2. Also i want to rationalize the complex number...

3. I want to seprate real and imaginary parts of a comaplex numbers

I try to sort a polynomial using the graded reverse lexicograpic order. According to the Documentation this is achieved via tdeg.
So here is what i tried:

with(Groebner):
sort( x+y+z, order = tdeg(x,y,z));

or

sort(x+y+z, [x,y,z], tdeg);

In both cases maple returns "x+y+z" instead of the expected "z+y+x". What am i doing wrong?

With the package VectorCalculus we can study the speed and acceleration to their respective components. Considering the visualizaccion and algebraic calculations and to check with their respective commands. Both 2D and 3D.

Velocidad-Aceleració.mw     (in spanish)

Lenin Araujo Castillo

Physics Pure

Computer Science

Hello,

I have to simplify a piecewise function and Maple gets a more complicated solution than needed.




I don't know how to handle this kind of problems with Maple?
I don't understand why Maple doesn't see this?
Am I doing something wrong?

Thanks in advance for your help / advice.


# the code of my example
restart:
Mf(x):=piecewise(x<=L/2,1/2*x*F,x>1/2*L,1/2*x*F-F*(x-1/2*L));
# Make a dimensionless function:
# -    Mf(x):= Mf(xi)*F*L
# -    variable ξ  ( xi:=x/L )
eq[1]:=Mf(xi)*F*L=Mf(x);
Mf(xi):=solve(eq[1],Mf(xi));
Mf(xi):=subs(x=xi*L,Mf(xi));
# F is the Force and L is the Length of the beam:
Mf(xi):=simplify(Mf(xi)) assuming F>0,L>0;
print("When I simplify this function by hand it will be");
Mf(xi):=piecewise(xi<=1/2,1/2*xi,xi>1/2,-1/2*xi+1/2);

Hello, everyone. I faced some promblems on maple. Hope you can help me.

I want to set up the commutation rules like Pauli sigma matrix in Physics Packages,

but failed to get the correct one.

Here is my maple code:

restart;
with(Physics);
Setup(mathematicalnotation = true);

Setup(quantumop = {Q}, algebrarule = {%AntiCommutator(Q[j], Q[k]) = 2*KroneckerDelta[j, k], %Commutator(Q[j], Q[k]) = 2*I*(Sum(LeviCivita[i, k, l]*Q[l], l = 1 .. 3))});

Commutator(Q[1], Q[2]);
                              
AntiCommutator(Q[1], Q[2]);
                              
AntiCommutator(Q[1], Q[1]);
                              
Commutator(Q[1], Q[1]);

In fact, Commutator(Q[1], Q[2]) give a incorrect result 0 while the correct answer is  2*I*Q[3].

Do anyone know how to resolve this?

Hello,

can I remove tick marks from axis in plots ?

Hi Everybody.

Doing some calculation in quantum mechanics, I stuble on that integral:

I see immediately that the answer is 1/2.  But Maple 18 cannot give an answer other than a limit that he cannot evaluate.  I even try assumption that p and hbar are realcons.  I get infinity.

Any idea?

Thank you in advance for your help.

--------------------------------------
Mario Lemelin

Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 -  64 bits
messagerie : mario.lemelin@cgocable.ca
téléphone :  (819) 376-0987