MaplePrimes Questions

In the link below is my code which generates numerical data of a Fourier series of an impulse function (not to be confused with the DIRAC DELTA IMPULSE).  The issue is as k increases the harmonic frequencies become faster than the time step which causes resolution issues at higher k values.  Is there a way to iteratively change the time step for each iteration of k value?  However, the range for t would have to correspondingly decrease or the file size would balloon by the corresponding change in k value.

So lets say my time step for k = 1 is 0.001 and my time range is 0..1, now for k = 1000 my time step would be 0.000001, but my time range would shrink from 0..1 to 0..0.001.  I hope this is not too confusing.  Maybe experimenting with my code will aid in understanding my request.

Thanks for any assistance

MATRIX_loop.mw


 

`` 

``

sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2 = 0, diff(diff(diff(g(eta), eta), eta), eta)+f(eta)*(diff(diff(g(eta), eta), eta))+2*g(eta)*(diff(diff(f(eta), eta), eta))-3*(diff(g(eta), eta))*(diff(f(eta), eta)) = 0, (diff(diff(g(eta), eta), eta))*(diff(diff(g(eta), eta), eta))-g(eta)*(diff(diff(diff(g(eta), eta), eta), eta))*f(0) = 1, (D(f))(0) = 1, (D(f))(5) = 0, g(0) = 1, (D(g))(0) = -1, (D(g))(5) = 0], numeric, method = bvp); plots[odeplot](sol1, color = red)

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

 

``

``


 In this link i attache the code , which is working on this Download Three_euqtions.mw and i got an error which is i couldnot understand. So, kindly can you help on this ??

Thats the solution

a)   -(2x+10)+8=4*(-x+1)

b)   5x+27=-9*(x-10)+7

c)   12x+2*(x-3)=7*(x-5)+1

d)    -4*(x+8) +14=7*(x-2)+29

 

 

Using a Maple procedure is a very common task. For instance, here, F is a differential operator applied to any expression f

restart;
F := proc (f) options operator, arrow; x*(diff(f, y))-y*(diff(f, x)) end proc;
                            / d   \     / d   \
                F := f -> x |--- f| - y |--- f|
                            \ dy  /     \ dx  /
F(sin(x*y));
                    2             2         
                   x  cos(x y) - y  cos(x y)

It is also possible to index a procedure. There is an appropriate syntax. At the end, one should not forget to assign the procedure to a table. Here, F simply returns its indices.

restart;
`index/F` := proc (idx) return idx end proc;
F := table(F);
                       F := TABLE(F, [])
F[i, j, k, 1, 2, 3];
                       [i, j, k, 1, 2, 3]

My question is how to define a procedure that uses both indexed and non-indexed arguments? That is, a procedure that returns a result according inputs of the type:

F[i, j, k,...](x, y, z,...)

I need such a syntax to define a differential operator (like the first example) that is further defined as a tensor within the Physics package (Physics:-Define). This differential operator thus needs to be indexed.

 

Hi, there! I came across a problem that needs your help. When I added the following codes in the "Click action" (I'm not sure if it is this name, My Maple version is not English), it returns the error: "too many levels of recursion". Please help me out. The codes are

if elenum = 1 then
    faultcom := theDiagnosis[1][1][1][1];
    faultloc := theDiagnosis[1][1][1][2];
else
    for numi from 1 to elenum do
        faultcom := faultcom union theDiagnosis[1][1][i][1];
        faultloc := faultloc union theDiagnosis[1][1][i][2];
    end do;
end if;

Plot the two curves defined by

−x* y^2+4 *x=5

and

(1/3)*x^3+y^2=1

and find the number of intersections of these curves in the square

 

Help me out.

Cheers.

 

This is what I typed in the Maple and only got 1 answer...but the answer should be 4..

eqn1:=−x* y^2+4 *x=5
eqn2:=(1/3)*x^3+y^2=1;

sol:=fsolve({eqn1,eqn2},{x=-10..10,y=-10..10});
              {x = 1.324605526, y = -0.4746486082}

To answer this question you need to create a Maple function using Maple's arrow (->) notation.

Your function should take a Maple list of complex numbers as its input and return the largest modulus from that list.

Enter your function in the box below.

 

Really need help 

Cheers.

Hello,

I'm new as to write programs in Maple.

My problem is: How can i write program code into an empty worksheet?

How can i insert new commands into an existing program/procedure?

What do i have to do?

Can you help me please for my first steps, i didn't find any helpful information about this.

Thanks in advance

Volker

Hi, I am new to Maple 2017. I would appreciate advice on how to plot the following integral over the range a = -10 to a = 10

I(a):=int(exp(-x^2)*sin(a*x),x=0..infinity)

 

Thank you.

Aidan

Where I can found inverted pendulum on maple sim or good book?thanks

Enter the matrix

To prevent typing errors you may copy and paste following Maple command for entering A.

 

A := Matrix([[4, 7, 5, 2, 6, 4, 9, 7, 5, 8, 5], [7, 7, 3, 9, 4, 5, 2, 3, 5, 6, 8], [5, 6, 5, 6, 6, 3, 7, 4, 9, 1, 9], [9, 1, 4, 8, 9, 4, 4, 2, 1, 3, 3], [5, 6, 8, 2, 2, 1, 2, 9, 5, 4, 1], [4, 7, 9, 6, 9, 2, 3, 5, 6, 3, 5], [7, 2, 5, 7, 3, 4, 9, 8, 1, 8, 9], [5, 7, 1, 2, 7, 9, 9, 2, 2, 4, 7], [1, 6, 5, 3, 9, 9, 1, 5, 7, 4, 2], [3, 7, 8, 7, 8, 3, 6, 4, 6, 9, 2]]);

 

[Note that Maple does not display the full matrix when it has more than 10 rows or 10 columns.  However, you can still work with the matrix using Maple commands.]

Use Maple to create the vector b that is column 6 from A and the matrix C that is made from columns 1 to 5 and 7 to 11 of A (in the same order as the columns of A).

Now solve the matrix equation

x = b

and enter the 6th component of the unique vector solution for x in the box below.  (Your answer should be an exact fraction, not a decimal)

The answer is 63695593/46022164

Can someone help me out?

Cheers.

 

The question is to find principle argument of the 5 roots of the polynomial z^5-5*z^4-z^2-2 

I used the Maple command:

Solve(z^5-5*z^4-z^2-2);

evalf(%);

argument(%);

But the Maple told me that expecting 1 argument but got 5.

I want to find all the argument and then use Max command to find the principle argument...it's doesn't work

Another similar thing is to find the moduli of a 5 roots of the polynomial z^5-5*z^4-z^2-2 

I got the same error as previous one ..

Can someone help me out ?

Cheers.

The question is to fin limit (n^1.5)*sum of 2*k/(2*k+3) k=1 to n as n goes to infinity 

I got the answer is float infinity...

How can I get a correct answer ?

Can someone help me out .

Cheers.

 

The question is to find integral of [e^(-x)*sin(x^2/)]/(2+x) from x=1 to infinity.

I typed in Maple and used evalf(%) 

But I got a float(infinity) on the upper limits instead of a correct numeric answer...

Can someone help me ?

Cheers.

1. Where do we report potential bugs

 

2. My finding:

https://imgur.com/a/w688r

The image uploading on this forum is not letting me and I find it to be less user friendly. Please view the sequence of images in the given link above.

This is not a problem, just my finiding.

 

First 820 821 822 823 824 825 826 Last Page 822 of 2322