I am trying to animate on maple. I have two solutions 𝑢1(𝑥, 𝑎) = 𝑐𝑜𝑠(𝑥 − 2𝑎) and 𝑢2(𝑥, 𝑎) = 𝑐𝑜𝑠(𝑥 + 2𝑎), for 𝑎 ∈ [1,10], how do I animate them?

S := Sum(1/i(i + 1)(i + 2)(i + 3), i = 1 .. infinity);
evalf(S);
                        -Float(infinity)
S1 := applyop(convert, [1, 2], S, parfrac, x);
                    infinity                      
                     -----                        
                      \                           
                       )              1           
              S1 :=   /     ----------------------
                     -----  i(i + 1)(i + 2)(i + 3)
                     i = 1                        
S2 := applyop(expand, 1, S1);
                    infinity                      
                     -----                        
                      \                           
                       )              1           
              S2 :=   /     ----------------------
                     -----  i(i + 1)(i + 2)(i + 3)
                     i = 1                        
S3 := map(op(0, S2), op(S2));
                                            /infinity    \
                                            | -----      |
                                            |  \         |
                                            |   )        |
                                            |  /     (-1)|
                                            | -----      |
                                            \ i = 1      /
            /infinity                      \              
            | -----                        |              
            |  \                           |              
            |   )                          |              
      S3 := |  /     i(i + 1)(i + 2)(i + 3)|              
            | -----                        |              
            \ i = 1                        /              
I am trying to find the limit approaching inf of the sum of 1/( 𝑖(𝑖 + 1)(𝑖 + 2)(𝑖 + 3)). I dont know what to do after this. 

lim 𝑛→+∞∑ 1 /(𝑖(𝑖 + 1)(𝑖 + 2)(𝑖 + 3))

Let 𝑦 and 𝑧 and be variables, let 𝑦 =( 𝑎+𝑏𝑧+𝑐𝑧^2)/( 𝑑+𝑒𝑧+𝑓𝑧 ^2) . Use the function solve in Maple to solve for 𝑧 in terms of 𝑦, that is, find a function 𝑧 = 𝑓(𝑦) (in fact, you should be able to find two). Let 𝑔(𝑧) = 𝑎+𝑏𝑧+𝑐𝑧^2)/( 𝑑+𝑒𝑧+𝑓𝑧 ^2)  and compute 𝑔(𝑓(𝑦)) for both of these functions and show that if you simplify the expressions you get a predictable answer. 

I ave already solved the values, however I am unable to simplify these expressions.

I am working with the Physics package in maple. I have a rank 6 tensor, I want to assign the elements of a set (with equal number of components as the tensor) as the corresponding components of a tensor. Entering manually all the components using 'RedefineTensorComponents' is not an option, as it has 4096 components. How can I do this?

Hello everyone!
I wrote proc for Spectral density estimation using Welch's method.

And I want to understand wich fuction is better for computation of Cross-correlation: SignalProcessing:-CrossCorrelation or Statistics:-CrossCorrelation?

Here is my programm with simple signal exmaples:

P.S. (7 hours after initial question) I've just found several mistakes in my PSDw-proc in overlapping.

So, I uploaded new version of my programm. But now I'm not sure in unit of measurement  of power spectral density, but algorithm works correctly.

Spectral_density-test-correct.mw

When plotting a sequence S, plot(S) a series of horizontal lines are plotted. Is there a way to plot a "time series" without having to reconstruct S, e.g., "seq([i,S[i]],i=1..numelems(S))" and plot that? Seems like a lot of work just to make a simple plot.

I am interested in finding the lower Riemann sum for a partition of unequal width.

The points of the partition are P = { -1, -1/4, 1/4, 3/4, 1} , and the function is f(x) = x^2.

My attempt:

restart: with(Student[Calculus1]):
RiemannSum(x^2,x={[ -1,-1/4],[-1/4, 1/4],[1/4,3/4],[3/4,1]},method=lower)

It said error range must be specified. I looked at the help page but I didnt find a specific command for entering partition points manually.

can anyone help me to calculate the exact  value of the eigenvalues of this matrix:
 

Download mat.mw

 Evaluate, to 10 significant figures

 integration of (e^(-x) *sin(x^2/2))/(3+x) with the limits of infinity and 1.

Dear all

I solve the first-order PDE with a boundary condition contains a parameter  s

When I run the code there is no solution displayed using pdsolve

Many thanks for your help

PDEBCS.mw

with(Maplets[Elements]);

Button("Deduct", Evaluate(f = 'work(3)'))

Here I run fucntion 'work()' when click the button. Is there a way to run several functions when I click the same button?

I have a matrix as follows.

How do we equalize this matrix to zero matrix and solve it?
 

A:=Matrix(2, 2, [[-0.0001633261895*z[1, 2]^2 + 0.0002805135275*z[1, 2]*z[2, 2] - 0.0001200583046*z[2, 2]^2 + 0.0006934805795*z[1, 1]^2 - 0.001190280265*z[1, 1]*z[2, 1] + 0.00007689977894*z[1, 1]*z[1, 2] - 0.00009937418547*z[1, 1]*z[2, 2] + 0.0005090615773*z[2, 1]^2 - 0.00003303758400*z[2, 1]*z[1, 2] + 0.00005683264925*z[2, 1]*z[2, 2] + 0.7021232886*z[1, 1] - 0.3171553245*z[1, 2] - 0.08291569324*z[2, 1] + 0.04647270631*z[2, 2] - 0.1436869545, 0.0002939068385*z[2, 1]^2 + 0.4237544799*z[1, 1] - 0.03129537402*z[1, 2] - 0.06276282411*z[2, 1] + 0.02529757039*z[2, 2] + 0.0004003811990*z[1, 1]^2 + 0.0002177682527*z[1, 1]*z[1, 2] - 0.0006872086309*z[1, 1]*z[2, 1] - 0.0001976167183*z[1, 1]*z[2, 2] - 0.0001764013184*z[2, 1]*z[1, 2] + 0.0001600777394*z[2, 1]*z[2, 2] - 0.1237363898], [0.00006763201108*z[2, 1]*z[1, 2] - 0.0001020812322*z[1, 2]*z[2, 2] - 0.00001554113990*z[2, 1]*z[2, 2] - 0.00003577693711*z[1, 1]*z[1, 2] + 0.0004330743651*z[1, 1]*z[2, 1] - 0.00001941220415*z[1, 1]*z[2, 2] - 0.01736180925 + 0.5623450996*z[2, 1] - 0.2353707048*z[2, 2] - 0.0003226356619*z[1, 1]^2 + 0.00007598605473*z[1, 2]^2 - 0.0001392051452*z[2, 1]^2 + 0.00003283047567*z[2, 2]^2 + 0.04653058230*z[1, 1] - 0.03026711709*z[1, 2], -0.00008037012799*z[2, 1]^2 + 0.03994641178*z[1, 1] - 0.02291248064*z[1, 2] + 0.3140461555*z[2, 1] + 0.01853659924*z[2, 2] - 0.0001862737861*z[1, 1]^2 - 0.0001013147396*z[1, 1]*z[1, 2] + 0.0002500356011*z[1, 1]*z[2, 1] + 0.00005403916772*z[1, 1]*z[2, 2] + 0.00008206914192*z[2, 1]*z[1, 2] - 0.00004377396755*z[2, 1]*z[2, 2] - 0.01370765196]])

And then, I want to select the roots in which 0<z[i,j]<1. 

Finally, I want to create matrix Z from these roots as follows   

Z:=Matrix(2, 2, [z[1, 1], z[1, 2],z[2, 1], z[2, 2]]);

2.22 *10^(-87)

2.55 *10^(-90)

would like to get 2.22 and 2.55 only by ignoring any power

how to do this?

How can I plot the contours with label for the following?

The Code is as follows:

z := -y + sech(x - 3*t);

w := 10*sech(x - 3*t);

with(plots);

P1 := plot(eval(w, t = 0), x = -10 .. 10):
P2 := contourplot(eval(z, t = 0), x = -10 .. 10, y = -eval(w, t = 0) .. eval(w, t = 0), contours = 5, grid = [101, 101]):
display(P1, P2);

Thanks

I assume I'm missing the verb here, but how does one display the elements of a vector and an array?

I can only do it with browsing or changing the "number format". 
 

Download Vector_Display_Problem.mw