Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi all. I have two questions about polynomials over Q.

i) Let f,g be polynomials over Q. How to show (with Maple) that the decomposition field of f is included in the decomposition field of g?

ii) Let f be a polynomial of degree n  over Q with Galois group C_n, the cyclic group of order n. Then, there is a polynomial P with coefficients in Q, of degree n-1, s.t. the iteration: u0=some root of f, u_{k+1}=P(u_k) gives all the roots of f; how to find P (with Maple)?

Thanks in advance.

Hello.

I have this problem when executing the entire worksheet or selected groups.
Also Maple can crash by itself, to its heart's content)
What I can do to solve this problem?
OS: W7 x64, Java is up to date

Thx.

(a) Plot the graph of  
                       sin(x)*exp( -x^2)
 for x in the interval [-2,2]. 
(b) Find to 10 decimal digits the maximum and minimum values of 
                         sin(x)*exp( -x^2)
 for x in [-2,2] AND find the corresponding values of x. [So if the maximum occurs at x=a, you should also compute sin(x)*exp( -x^2)   both to 10 digits. Similarly for the minimum. Using unapply to make the expression into a function will be useful here.]  

So far I have this for a

> j := exp(-x^2)*sin(x);

> plot(j, x = -2 .. 2);

 

There is a one-to-one correspondence between subsets of {1, 2, . . . , n} and binary lists of length n, that is, lists L = [x1, x2 , . . . , xn] where x1, x2, . . . , xn are elements of the set {0,1}.  The correspondence is given by associating to the set S the list L where xi = 1 if i is in S and 0 if not. For example, the set {1,3,5} corresponds to the list [1,0,1,0,1,0,0] if n = 7.

(a) Write a procedure list_to_set whose input is a binary list and whose output is the corresponding set. E. g., list_to_set([1,0,1,0,1]) will return the set {1,3,5}. Note that nops(L) is the length of a list.

(b) Write a procedure set_to_list whose input is a pair S,n where S is a subset of {1, 2, . . . , n} and n is a positive integer and whose output is the binary list of length n corresponding to the set S. E. g., if n = 5 then set_to_list({1,3,5},5) will return [1,0,1,0,1].

(c) Show by a few tests that each procedure works. Then apply set_to_list to each set in the powerset of {1, 2, 3, 4} to form all binary lists of length 4. Make a program to print out a table of the following form. (But the order need not be the same as that started below.)

   [0,0,0,0] <-->  {  }
   [1,0,0 0] <--> { 1 }
   [0,1,0,0] <--> { 2 } 
    ........
    etc

Some extra commas in the output is okay. You may obtain the power set of the set {1,2,...,n} by the command powerset(n); but you must first load the package combinat.

The expression X^n, for an interger power n and any X, can be computed using the following formulae, which represent a negative power in terms of a positive power (-n), and a positive power in terms of a smaller non-negative power (either n/2 or n-1), and use only multiplication and division: 
                                    
                          X^n = 1/X^(-n)                      if n < 0 
                          X ^n = I*d                             if n = 0
                          X ^n  = X^(n/2) * X^(n/2)      if n is even
                          X^n  = X*X^(n-1)                  if n is odd.
These formulae lead to an efficient recursive algorithm for computing integer powers using the minimal number of multiplications. 

(a) Write a procedure MatPow(X,n::integer) to implement this algorithm for computing powers of matrices. Test MatPow(<<1|2>,<3|4>>,12) and MatPow(<<1|2>,<3|4>>,-12).

(b) Write a procedure PolyPow(X,n::integer) to implement this algorithm for computing powers of numbers and polynomials. Your procedure needs to exapnd each product of polynomials in order to be effective. Test PolyPow(123,12), PolyPow(123,-12), PolyPow(x^2+1,12) and PolyPow(x^2+1,-12). 

Write a procedure using the variable args that will take an unspecified finite number of numbers, delete the smallest and the largest, and return the average of the rest as a floating point number. If there are fewer than 3 arguments have an error message say: "There are not enough arguments. There should be at least three." 

You should have no input parameters in the definition of the procedure. You may write it directly or you may use the Maple command sort as a part of your program. Do ?sort to see how sort works. Test your procedure with each of the following argument sequences:  

   50,40,40, 40, 40, 10

   1,2

   seq(100 - i, i = 1..100)

   seq(modp(n,111), n=1..1000).


If n people (numbered 1 to n) stand in a circle and someone starts going around the circle and eliminating every other person till only one person is left, the number J(n) of the person left at the end is given by 

    J(n) = 1                           if n = 1
    J(n) = 2*J(n/2) - 1          if n > 1 and n is even
    J(n) = 2*J((n-1)/2) + 1   if  n > 1 and n is odd

(i) Write a recursive procedure to compute J. [As a check the first 16 values (starting with 1) of J(n) are 1,1,3,1,3,5,7,1,3,5,7,9,11,13,15,1]. 
(ii)Compute the value of J(10000). 
(iii) Can you explain why this is so much faster than our recursive procedure to compute the n-th Fibonacci number?

Has anyone installed and run maple under ubuntu installed the Windows Subsystem for Linux?   We are having trouble doing this: specifically running programs with text input files, etc.

 

 

I have a say sum(F(k) ,k=1..n)  It fails unless n is an actual integer.  But sum(F(k), k=1..a) works. Then then eval(%,a=n) completes it. That is probably correct but not necessarily a valid assumption. Details in attached. Also it is hard to understand the error message.
 

restart

``

``

Sinta := (k^2+n^2-k)*n/((k^2+n^2-2*k+1)*(k^2+n^2))

(k^2+n^2-k)*n/((k^2+n^2-2*k+1)*(k^2+n^2))

(1)

Ai := sum(Sinta, k = 1 .. n)

-((1/2)*I)*(2*n^2+I*n)*Psi(n-I*n)/(4*n^2+1)+((1/2)*I)*(2*n^2-I*n)*Psi(n+I*n)/(4*n^2+1)+((1/2)*I)*(-2*n^2+I*n)*Psi(n+1-I*n)/(4*n^2+1)-((1/2)*I)*(-2*n^2-I*n)*Psi(n+1+I*n)/(4*n^2+1)+((1/2)*I)*(2*n^2+I*n)*Psi(-I*n)/(4*n^2+1)-((1/2)*I)*(2*n^2-I*n)*Psi(I*n)/(4*n^2+1)-((1/2)*I)*(-2*n^2+I*n)*Psi(1-I*n)/(4*n^2+1)+((1/2)*I)*(-2*n^2-I*n)*Psi(1+I*n)/(4*n^2+1)

(2)

Souta := n/(k^2+n^2-k)

n/(k^2+n^2-k)

(3)

NULL

sum(Souta, k = 1 .. n)

Error, (in assuming) when calling 'Dfnt_4'. Received: 'when calling 'Dfnt_4'. Received: 'when calling 'unknown'. Received: 'invalid input: Dfnt_4 expects its 3rd argument, fpts, to be of type Or(list, piecewise), but received 0'''

 

"(->)"

Error, invalid input: evalf[10] expects 1 argument, but received 0

 

 

 

This works

sum(Souta, k = 1 .. a)

n*Psi(a+1/2-(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)-n*Psi(a+1/2+(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)-n*Psi(1/2-(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)+n*Psi(1/2+(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)

(4)

Ao := eval(%, a = n)

n*Psi(n+1/2-(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)-n*Psi(n+1/2+(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)-n*Psi(1/2-(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)+n*Psi(1/2+(1/2)*(-4*n^2+1)^(1/2))/(-4*n^2+1)^(1/2)

(5)

"(->)"

n*Psi(n+.5000000000-.5000000000*(-4.*n^2+1.)^(1/2))/(-4.*n^2+1.)^(1/2)-1.*n*Psi(n+.5000000000+.5000000000*(-4.*n^2+1.)^(1/2))/(-4.*n^2+1.)^(1/2)-1.*n*Psi(.5000000000-.5000000000*(-4.*n^2+1.)^(1/2))/(-4.*n^2+1.)^(1/2)+n*Psi(.5000000000+.5000000000*(-4.*n^2+1.)^(1/2))/(-4.*n^2+1.)^(1/2)

(6)

````

``


 

Download Summation_problem.mw

Dear Users!

Hope you would be fine. I want to construct system of equations by comparing the likes powers of x^i*y^j*t^k1*exp(k2*eta) for an expression H1 present in attached file. Please see the fix my problem. I shall be very thankful for your kind help. 

Help.mw

From  time to tiime Malel emits a sort of "Boiing" noise when I make a keyboard error. I am not sure exaclty when it does this, but it is loud. Is it possible to turn it off? I have looked quite extensively and cannot find out anything.

In the following problem at two example are given. For Z=2 the sum is converging whereas at Z=4 it is not converging. Thank you

 

PROBLEM.mw

I have a problem writing a program for the numerical solution of nonlinear volterra integral equation using the method of reproducing kernel space. I have my algorithm as well as the program I tried to write, though they are full of error messages. Please could anyone give me a clue on how to go about my challenges. The algorithm is as follows:

Step 1. Fix π‘Ž ≤ π‘₯ and 𝑑 ≤ 𝑏.
If 𝑑 ≤ π‘₯, set 𝑅π‘₯(𝑑) = 1 − π‘Ž + 𝑑.
Else set 𝑅π‘₯(𝑑) = 1 − π‘Ž + π‘₯.
Step 2. For 𝑖 = 1, 2, . . . , π‘š set π‘₯i = (𝑖 − 1)/(π‘š − 1).

Set πœ“i(π‘₯) = 𝐿t𝑅π‘₯(𝑑)|𝑑=π‘₯i .
Step 3. Set 𝑒0(π‘₯1) = 𝑒(π‘₯1).
Step 4. For 𝑖 = 1, 2, . . . , π‘š set 𝛾ij = [πœ“-1]ij.
Step 5. 𝑛 = 1.
Step 6. Set Sn = Σ𝑛
π‘˜=1 𝛾nk𝑒k-1(π‘₯k).
Step 7. Set 𝑒n(π‘₯) = Σ𝑛
𝑖=1 Siπœ“i(π‘₯).
Step 8. If 𝑛 < π‘šthen set 𝑛 = 𝑛 + 1 and go to step 6.
Else stop.

So I have this system of equations with which I am not sure if the result is the same or not using "series" and "limit" or what is going on here.

I hope it is clear what I mean.


 

restart; with(MathematicalFunctions); Assume(k__2H2O > 0, `k__HA+OH` > 0, `k__A+H2O` > 0, `k__H3O+OH` > 0, `k__HA+H2O` > 0, `k__H3O+A` > 0, HA__0 > 0, H2O > 0); sys := k__2H2O*H2O^2+`k__A+H2O`*H2O*(HA__0-HA)-(H3O*`k__H3O+OH`+HA*`k__HA+OH`)*OH = 0, k__2H2O*H2O^2+`k__HA+H2O`*H2O*HA-(`k__H3O+A`*(HA__0-HA)+`k__H3O+OH`*OH)*H3O = 0, (H2O*`k__HA+H2O`+OH*`k__HA+OH`)*HA-(H2O*`k__A+H2O`+H3O*`k__H3O+A`)*(HA__0-HA) = 0; sys := `~`[simplify]([eval(eval(sys, HA = HA__0+OH-H3O), HA__0 = x__HA0*H2O)]); sol := solve(sys, [OH, H3O]); sol := sol[1]; OH__sol := simplify(rhs(sol[1])); H3O__sol := simplify(rhs(sol[2])); simplify(OH__sol*H3O__sol); OHH3O := simplify(limit(%, `k__HA+OH` = 0)); series(OHH3O, x__HA0 = 0, 2); collect(convert(%, polynom), x__HA0, simplify, factor); r1 := limit(%, x__HA0 = 0); r2 := radnormal(limit(OHH3O, x__HA0 = 0)); simplify(r1-r2)

[`&Intersect`, `&Minus`, `&Union`, Assume, Coulditbe, Evalf, Get, Is, SearchFunction, Sequences, Series]

 

{H2O::(RealRange(Open(0), infinity))}, {HA__0::(RealRange(Open(0), infinity))}, {k__2H2O::(RealRange(Open(0), infinity))}, {`k__A+H2O`::(RealRange(Open(0), infinity))}, {`k__H3O+A`::(RealRange(Open(0), infinity))}, {`k__H3O+OH`::(RealRange(Open(0), infinity))}, {`k__HA+H2O`::(RealRange(Open(0), infinity))}, {`k__HA+OH`::(RealRange(Open(0), infinity))}

 

k__2H2O*H2O^2+`k__A+H2O`*H2O*(HA__0-HA)-(H3O*`k__H3O+OH`+HA*`k__HA+OH`)*OH = 0, k__2H2O*H2O^2+`k__HA+H2O`*H2O*HA-(`k__H3O+A`*(HA__0-HA)+`k__H3O+OH`*OH)*H3O = 0, (H2O*`k__HA+H2O`+OH*`k__HA+OH`)*HA-(H2O*`k__A+H2O`+H3O*`k__H3O+A`)*(HA__0-HA) = 0

 

[-OH^2*`k__HA+OH`+((-x__HA0*`k__HA+OH`-`k__A+H2O`)*H2O+H3O*(`k__HA+OH`-`k__H3O+OH`))*OH+k__2H2O*H2O^2+`k__A+H2O`*H2O*H3O = 0, (x__HA0*`k__HA+H2O`+k__2H2O)*H2O^2+`k__HA+H2O`*(OH-H3O)*H2O+(-`k__H3O+A`*H3O+OH*(`k__H3O+A`-`k__H3O+OH`))*H3O = 0, H2O^2*x__HA0*`k__HA+H2O`+((x__HA0*`k__HA+OH`+`k__A+H2O`+`k__HA+H2O`)*OH-H3O*(`k__A+H2O`+`k__HA+H2O`))*H2O+(OH-H3O)*(H3O*`k__H3O+A`+OH*`k__HA+OH`) = 0]

 

-RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)*H2O^2*(-`k__A+H2O`*RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)+k__2H2O)/(-`k__H3O+OH`*RootOf(-x__HA0*`k__A+H2O`^2*`k__HA+H2O`+k__2H2O^2*`k__H3O+A`-k__2H2O*`k__A+H2O`^2-k__2H2O*`k__A+H2O`*`k__HA+H2O`+(2*x__HA0*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-k__2H2O*`k__A+H2O`*`k__H3O+A`+k__2H2O*`k__A+H2O`*`k__H3O+OH`+k__2H2O*`k__H3O+OH`*`k__HA+H2O`)*_Z+(-x__HA0*`k__H3O+OH`^2*`k__HA+H2O`-k__2H2O*`k__H3O+A`*`k__H3O+OH`+`k__A+H2O`^2*`k__H3O+OH`+`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`)*_Z^2+(`k__A+H2O`*`k__H3O+A`*`k__H3O+OH`-`k__A+H2O`*`k__H3O+OH`^2-`k__H3O+OH`^2*`k__HA+H2O`)*_Z^3)+`k__A+H2O`)

 

-(k__2H2O*`k__H3O+A`^2-2*`k__A+H2O`^2*`k__H3O+A`+`k__A+H2O`^2*`k__H3O+OH`-2*`k__A+H2O`*`k__H3O+A`*`k__HA+H2O`+2*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`+`k__H3O+OH`*`k__HA+H2O`^2)*`k__A+H2O`*`k__HA+H2O`*H2O^2*x__HA0/((`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`)*(k__2H2O*`k__H3O+A`^2-`k__A+H2O`^2*`k__H3O+OH`-2*`k__A+H2O`*`k__H3O+OH`*`k__HA+H2O`-`k__H3O+OH`*`k__HA+H2O`^2))-(k__2H2O*`k__H3O+A`-`k__A+H2O`^2-`k__A+H2O`*`k__HA+H2O`)*H2O^2*(`k__A+H2O`+`k__HA+H2O`)/(`k__H3O+A`*(`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`))

 

-(k__2H2O*`k__H3O+A`-`k__A+H2O`^2-`k__A+H2O`*`k__HA+H2O`)*H2O^2*(`k__A+H2O`+`k__HA+H2O`)/(`k__H3O+A`*(`k__A+H2O`*`k__H3O+A`-`k__A+H2O`*`k__H3O+OH`-`k__H3O+OH`*`k__HA+H2O`))

 

k__2H2O*H2O^2/`k__H3O+OH`

 

-`k__A+H2O`*(-(`k__A+H2O`+`k__HA+H2O`)^2*`k__H3O+OH`+k__2H2O*`k__H3O+A`^2)*H2O^2/(`k__H3O+OH`*`k__H3O+A`*((-`k__A+H2O`-`k__HA+H2O`)*`k__H3O+OH`+`k__A+H2O`*`k__H3O+A`))

(1)

``


 

Download Mapleprimes_-_Ionproduct.mw

is there example data that can verify maxwell equations?

First 123 124 125 126 127 128 129 Last Page 125 of 2238
ο»Ώ