ON SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS

Size: px
Start display at page:

Download "ON SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS"

Transcription

1 MATHEMATICS OF COMPUTATION Voume 66, Number 27, January 997, Pages S (97) ON SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS C. FIEKER, A. JURK, AND M. POHST Abstract. Let Q E Fbeagebraic number fieds and M Fafree o E -modue. We prove a theorem which enabes us to determine whether a given reative norm equation of the form N F/E (η) = θ has any soutions η M at a and, if so, to compute a compete set of nonassociate soutions. Finay we formuate an agorithm using this theorem, consider its agebraic compexity and give some exampes.. Introduction Soving norm equations is a centra probem in the area of agebraic number theory. Athough there is an agorithm for soving absoute norm equations (e.g. see [] or [8, 5.3, 6.4]), none (except the absoute one) exists in the reative case. We outine a new agorithm to decide whether a reative norm equation has soutions at a and then, if there are soutions to compute a compete set of nonassociate (with respect to units of reative norm ) soutions. Finay we discuss the compexity of this agorithm and give some exampes. 2. Preiminaries First, we introduce severa definitions and notations. We consider the foowing situation: F = E(β) n E = Q(α) m Q Let α = α (),..., α (m) be the roots of the monic irreducibe poynomia f Q[t], E := Q(α). Furthermore, et β be a root of a monic irreducibe poynomia g E[t]and F:= E(β). We assume that α (),..., α (m) denote the rea roots of f and that α (m+) = α (m+m2+),...,α (m+m2) =α (m+2m2) are in C\R. For an arbitrary η E we define η (i) ( i m) as the image of η under the Q-isomorphism from E to E (i) := Q(α (i) )which maps α to α (i). The norm of an eement η of E is defined in the usua way: N(η) :=N E/Q (η):= m i= η(i). The definitions for F are essentiay the same, but here we have to be carefu about which fied we are using as base fied. In genera, the conjugates cannot be ordered in rea and pairs of compex ones and of course we get different conjugates Received by the editor August 30, 994 and, in revised form, March 27, 995 and Juy 20, Mathematics Subject Cassification. Primary Y40. Key words and phrases. Agebraic number theory, norm equations, reative norm equations, reative extensions. 399 c 997 American Mathematica Society License or copyright restrictions may appy to redistribution; see

2 400 C. FIEKER, A. JURK, AND M. POHST when we consider F as an extension of E rather than as an extension of Q. We wi discuss this in detai ater. The ring of integers of a fied K is denoted by o K in the seque. 3. Absoute norm equations In this section we give a short review on the Fincke Pohst agorithm to sove absoute norm equations. Given an arbitrary, but fixed k Z >0 and M EafreeZmodue of E, wewant to find a η M subject to (3.) or (3.2) N(η) =k N(η) = k, i.e., we want to determine if soutions exist and, if so, compute a of them. If M is of fu rank m, we fix a maxima system of independent units ɛ,..., ɛ r (r = m +m 2 ) of its ring of mutipiers. For practica computations it is advisabe to use LLL-reduction in the ogarithmic attice to produce units for which the absoute vaues of each of their conjugates are cose to. The next emma gives expicit bounds for a compete set of nonassociate soutions. Lemma 3.. Let η be a soution of (3.2). Then there exists a unit ɛ and a soution η = ηɛ subject to (3.3) where η(i) R m i k R i ( i m), R i := exp 2 m +m 2 j= og( ɛ (i) j ). Proof. See [, (6.3) Lemma] or [8, Theorem (4.2), Chapter 6]. If M is not of rank m, then it is not known how to compute reaistic bounds on the conjugates of the (finitey many) soutions of (3.2). Hence, in that case we need to stipuate bounds R i such that any soution η satisfies (3.3). This works we if we are ony interested in sma soutions, say with coefficients bounded by 0 6. Furthermore, we need the foowing, rather technica, emma from [, (6.23) Satz] or [8, Theorem (3.8), Chapter 5]: Lemma 3.2. For arbitrary γ, r R >0 define the functions h : R > R : t t t og t and g : R > R : t ( h(t))t h(t) + h(t)t h(t). Then there exists a unique zero λ = λ(γ,r) of g(t) ( + γ r )2/m. Using the previous two emmas, one can prove the foowing theorem, which wi aow us to sove absoute norm equations of the type (3.2) or (3.). Theorem 3.3. For arbitrary γ R >0 et λ = λ(γ,k) as in Lemma 3.2. Define constants L i := 2ogRi og λ, U i := 2ogRi og λ ( i m) with R i as in Lemma 3.. License or copyright restrictions may appy to redistribution; see

3 SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS 40 If (3.2) is sovabe, then there exists r =(r,...,r m ) Z m and a soution η of (3.2) subject to (3.3), and m (3.4) λ rj η (j) 2 m(k + γ) 2/m. j= The coordinates of r satisfy m j= r j =0,L j r j U j ( j m) and r m+j = r m+m 2+j for a j m 2 with at most one exception, where we have r m+j + = r m+m 2+j. Forrank M = m, any soution of (3.2) is associate with one satisfying (3.4). Remark 3.4. Representing η in a basis of M, the eft hand side of (3.4) becomes a positive definite quadratic form over M, so that the soutions of the inequaity are attice points inside an eipsoid. The basic idea in the proof of the theorem above is an observation due to M. Pohst which guarantees (under certain conditions) the existence of some vector λ R m >0 with m j= λ j η (j) 2 = mk for every soution η of (3.2). The main task in the proof then is to obtain a finite set of candidates for the λ s. This was achieved by U. Fincke [] who transformed our probem into a discrete one by considering ony vectors of the form (λ r,...,λ rm ), r j Z. In the next step, we obtain bounds for the r j s from the bounds for the soutions. Finay, we reduce the set of the admissibe exponent vectors r =(r,...,r m ). To use this theorem for soving (3.2), one has to cacuate the set of a admissibe r Z m, for each such r to compute the set of a η M subject to (3.4) and, finay, to determine the soutions of (3.2) among them. The foowing statement concerns the agebraic compexity when using the Fincke-Pohst method to enumerate the points of the eipsoids. Theorem 3.5. Let k Z and the signature (m,m 2 ) be fixed. Then there exists a v > depending ony on the conjugates of a Z basis of M and a sufficienty arge U max m j=i U j, U j as in Theorem 3.3, such that the number of arithmetic operations (addition, mutipication, division, cacuation of square roots) used to sove (3.2) with the agorithm described above is bounded by (m 2 +)(3 4ogk mog λ + 4ogU og λ )m+m2 (2Um 2/3 vλ /2 (+ λ k )/m ) m/2 (2m 3 +24m 2 ). Proof. See [, (7.5) Satz]. 4. Reative norm equations Before we sketch the theory and the agorithm, we make some further definitions. Let M F be a free o E -modue. As above, we assume the rank of M to be n in order to produce bounds. For an arbitrary θ o E we want to decide if there are any η M satisfying (4.) N F/E (η) =θ or (4.2) N F/E (η) = θ and, if so, we wi compute a compete set of nonassociate soutions. Athough the agorithm described in the previous section coud, in principe, be used to accompish this (certainy one can sove N(η) = N E/Q (N F/E (η)) = License or copyright restrictions may appy to redistribution; see

4 402 C. FIEKER, A. JURK, AND M. POHST F ()... F (n)... F ((m )n+)... F (nm) E ()... Q E (m) Figure. The ordering of the conjugate fieds N E/Q (θ) and then test the soutions), it turns out to be far too expensive in terms of computation time. Hence, the approach for the absoute case needs to be changed appropriatey. To dea with the difficuties concerning the conjugates of F arising from the fact that we can consider F both as an extension of E and of Q, we fix the ordering of the conjugate fieds as described in Figure. To simpify notation, we define the foowing abbreviation: N (j) F/E (γ) :=N (γ ((j )n) )= n F ((j )n+) /E (j) i= γ((j )n+i). We require that the conjugate fieds E (),..., E (m) of E are ordered as described in 2. For the rea conjugate fieds E (i) ( i m ) we stipuate furthermore that the conjugate fieds of F, which are extensions of these, are ordered in a simiar manner. For i m we require that F ((i )n+j) R ( j s i )and F ((i )n+si+j) = F ((i )n+si+ti+j) R ( j t i ), s i +2t i =n,where(s i,t i ) denotes the reative signature of F ((i )n+) over E (i) [4, Lemmas ]. As in the absoute case, the first task is to obtain bounds for a set containing a non-associate soutions of (4.2) or (4.). In the previous section we used units of the ring of mutipiers of M to transform arbitrary soutions to those contained in a bounded region. In order to preserve (4.), one has to use units ɛ where N F/E (ɛ) =, whereas for (4.2), it suffices to require that N F/E (ɛ) is a torsion unit of E. Here we wi restrict ourseves to units of the form ɛ n /N F/E (ɛ). Ceary, we have N F/E (ɛ n /N F/E (ɛ)) = N F/E (ɛ n )/N F/E (ɛ) n =. Letɛ,..., ɛ r be a maxima system of independent units of the ring of mutipiers of M. We have the foowing emma, which corresponds to Lemma 3.: Lemma 4.. Let η be a soution of (4.2) or (4.). Then there exists a unit ɛ and a soution η = ηɛ of (4.2) and (4.) subject to η((j )n+i) (4.3) R R n (j )n+i ( i m, j n), (j )n+i θ (j) where R (j )n+i =exp ( n 2 r = og( ɛ ((j )n+i) ) + 2 ) r og( N (j) F/E (ɛ ) ). = License or copyright restrictions may appy to redistribution; see

5 SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS 403 Proof. The proof is essentiay the same as in the absoute case. The bounds are worse than those in Lemma 3. because of the restricted set of units used for the transformation [4, Lemma 4.]. Assuming that we have bounds for the conjugates of soutions from a different source, we can again generaize the method to modues M of rank ess than n, as in the ast section. We note that the bounds here are not a generaization of those obtained in Lemma 3.. Having estabished this emma, and using Lemma 3.2 of the preceding section, we obtain the foowing theorem: Theorem 4.2. For arbitrary γ j R >0 et λ j = λ(γ j, θ (j) )( j m)as in Lemma 3.2. Define constants L (j) i := 2ogR (j )n+i og λ j, U (j) i := 2ogR (j )n+i og λ j ( i n) with R (j )n+i as in Lemma 4.. Then for every soution ν of (4.) there exists a unit ɛ with N F/E (ɛ) =,r (j) =(r (j),...,r(j) n ) Zn ( j m) and a soution η = ɛν of (4.) subject to (4.3) such that () n i i= λr(j) (2) m i= r(j) i (3) L (j) i j η ((j )n+i) 2 n( θ (j) + γ j ) 2/n, j m; =0, j m; U (j) i, i n, j m; r (j) i (4) For j m we have in addition r (j) s j+i = r(j) s j+t j+i for a i t j with at most one exception, where we have r (j) s j+i +=r(j) s. j+t j+i Proof. Let ɛ and η = ν be as in Lemma 4.. For j m and i n define y (j) i := θ (j) 2/n η ((j )n+i) 2. Ceary, m i= y(j) i =. Sincey (j) i and λ j are positive, we can find r (j) i Z, ẽ (j) i [0, ) with (4.4) y (j) i = λ r(j) i j and 0 = n = r(j) = n = ẽ(j). By performing some engthy computations, we can change the ẽ s and the r s to fufi conditions (2) and (4) and then verify (3) using (4.4). As in the absoute case, the vaidity of () is a consequence of Lemma 3.2 [4, Satz 4.4]. A straightforward computation gives N (j) F/E (η) θ(j) +γ j for a η satisfying () of Theorem 4.2. Anaogous to the description given in the paragraph after Remark 3.4, we get an agorithm for soving (4.2). We note that the inequaities () describe eipsoids defined via the positive definite quadratic forms [4, Lemma 4.6] x =(x,...,x n ) n i= +ẽ (j) i λ r(j) i j x j 2 ( j m). Their attice points can be cacuated with a modified Fincke-Pohst enumeration agorithm. 5. A modified Fincke-Pohst agorithm A main part in soving norm equations is the enumeration of a points in suitabe eipsoids. In this section we present a modification of the Fincke-Pohst method License or copyright restrictions may appy to redistribution; see

6 404 C. FIEKER, A. JURK, AND M. POHST adopted to reative norm equations. We consider the foowing situation: Let A (i) C n n be positive definite, C (i) > 0 subject to A (i) = A (i+m2), C (i) = C (i+m2) (m <i m +m 2 ). We present an agorithm to cacuate a x M subject to x tr (i) A(i) x (i) C (i), where x (i) =(x ((i )n+),...,x ((i )n+n) ) tr. Let ν,...,ν n be an o E -basis of M. For x M we have x= n = z ν with z i o E and x ((i )n+j) = n = z(i) ν ((i )n+j). Defining V i := (ν ((i )n+j),weset ) j n n Q i (x):=x tr (i) A(i) x (i) =(z,...,z n )V tr i A (i) V i (z,...,z n ) tr. Let B i := V tr i A(i) V i. Using the agorithm for quadratic suppement (see [3, (2.3)]), we obtain z tr B i z = n = q(i), z(i) + n s=+ q(i),s z(i) s 2. (Note that we have q (i), > 0 since B i is positive definite.) This yieds the foowing straightforward agorithm: Agorithm 5.. (Reative enumeration) Input: q (i),s, C(i) as above Output: A x M subject to Q i (x) C (i) Init: Step : := m, T (i) := C (i),u (i) S := {z o E z (i) U (i) := 0 ( i m). 2 T(i) i m}. Step 2: If S goto Step 4. Step 3: Set := +,if>mthen terminate ese goto Step 2. Step 4: Choose z in S arbitrary and deete z from S. Step 5: If =printx:= n s= z sν s goto Step 2. Step 6: :=, T (i) := T (i) + q(i) +,+ (z(i) + + U (i) (i) + ), U := n s=+ q(i),s z(i) s, ( i m) goto Step. It remains to give an agorithm to compute the set S in Step. We wi do this using a modified Fincke-Pohst agorithm. Instead of computing S, we compute a set S := {z o E m i= z(i) U (i) 2 C } with C := m T (i) i=.letµ q (i),...,µ m, be a Z-basis of o E and define µ ()... µ () m U () W :=... Cm (m+). µ (m)... µ (m) m q (i), U (m) Now we have m i= x(i) U (i) 2 =(z,...,z m,)w tr W (z,...,z m,) tr for every x = m s= z sµ s o E, z s Z. Ceary, W tr W is positive semidefinite and rea, so the S can be cacuated with the cassica Fincke-Pohst agorithm if we fix the ast coordinate to be. License or copyright restrictions may appy to redistribution; see

7 SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS The agorithm We are now abe to present a compete agorithm for soving (4.2). Agorithm 6.. (Computing a compete set of nonassociate soutions of (4.2)) Input: θ o E, γ R m >0, M and a compete set ɛ,...,ɛ r of independent units of the ring of mutipiers of M. Output: A compete set of nonassociate soutions of (4.2). Init: Compute λ j = λ(γ j, θ (j) ) as in Lemma 3.2, R (j )m+i as in Lemma 4. and L (j) i and U (j) i as in Theorem 4.2. Let S :=. Step : Compute I := {(r (j) i ) i n Z n m (2) (4) of Theorem 4.2 hod }. j m Step 2: Whie (I ) do Init enumeration: Choose r I arbitrary, and deete r from I. Enumerate: Determine a η M satisfying () of Theorem 4.2 (e.g., using Agorithm 5.). Check: For each soution of the previous step, check (4.2). If we have a soution, check whether it is associated with one aready in S. If we have a new soution, store it in S. Step 3: If S = return No soution ese return S. In order to sove (4.) instead of (4.2), one simpy has to change the condition in Check. To estimate the agebraic compexity of this agorithm, we introduce the T 2 -norm: T 2 : E R:η m i= η(i) 2 and the so caed reative T 2 -norm T F/E 2,j : F R:η n i= η((j )n+i) 2 ( j m). Let ν,..., ν n be an o E -basis of M and ν,..., νn the corresponding dua basis subject to Tr F/E (ν i νj )=δ i,j. Define d j := max T F/E 2,j (νi )( j m), L j := max{λ r(j) i j i n m C := d j L j (n θ (j) + γ j ) 2/n. j= r I, i n}, Furthermore, et µ,...,µ m be a Z-basis of o E and µ,..., µ m its dua basis subject to Tr E/Q (µ i µ j )=δ i,j. Letd:= max m T 2 (µ ). We begin with the compexity of Agorithm 5.. Lemma 6.2. () The ( agorithm for computing one set S i in Step of Agorithm 5. cacuates O (2 ( ( 4dC +m ) )) dc +) 4dC + points. For each point it needs O(m 2 ) operations, ( and additiona ( O(m 3 ) operations for precomputing. (2) Agorithm 5. needs O (mn 2 + m 3 ) (2 ( ( 4dC +m ) )) n ) dc +) 4dC + operations. Proof. In the same way as in [3, Proof of 3.5] (cf. [4, Lemma 3.7] or [4, Lemma 4.3]) we normaize the task to q (i),. Now C is an upper bound for a S i. () Since we fix the ast coordinate, we effectivey enumerate an m dimensiona eipsoid, so the statements foow at once from [3]. License or copyright restrictions may appy to redistribution; see

8 406 C. FIEKER, A. JURK, AND M. POHST (( (2) Ceary, we consider at most O (2 dc +) and for each point we need O(mn 2 + m 3 )operations. ( ( 4dC +m 4dC ) )) n ) + points, We note that we need O(mn 3 ) operations for cacuating the matrices B i and the q (i),s. Theorem 6.3. Define D j := 2 max i n og R (j )n+i, R (j )n+i as in Lemma 4., j m, ( Q := mn 3 +(mn 2 + m 3 ) (2 ( ) n 4dC + m dc +) +). 4dC The set I := {(r (j) i ) i n j m Z n m (2) (4) of Theorem 4.2 are vaid } contains at most m m 2 +m 2 ((t j +)( D j +3) sj+tj 2 ) ( D j +3) n og λ j og λ j j= j=m + eements. For each (r (j) i ) i n I (i.e., for each quadratic form considered) the j m agorithm needs O(Q) arithmetic operations to enumerate a η with () of Theorem 4.2. Hence Agorithm 6. requires arithmetic operations. O(( I)Q) Proof. See [4, Lemmas 4.2 and 4.3]. We note that the estimate for I is a straightforward computation using ony (2) (4) of Theorem 4.2 and Lemma 4.. Let E := Q(α), α a zero of 7. Exampe f(x) :=x 3 +2x 2 x+2. We consider the extension F := E(β) of absoute degree 9, where β is a zero of g(x) :=x 3 +x 2 +(+2α+α 2 )x+( 2+α 2 ). The signature of E is m = m 2 =. SinceFhas ony one rea embedding, the unit rank of F is 4. A system of fundamenta units of o F is ɛ := 2 (α + α2 ), ɛ 2 := 2 (2 + 2α +(2 α α2 )β), License or copyright restrictions may appy to redistribution; see

9 SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS 407 ɛ 3 := ( 3+α 2 +(2 2α α 2 )β+(+α)β), ɛ 4 := 2 ( α +5α2 +(6+29α+9α 2 )β+(4 9α 9α 2 )β 2 ). We consider the foowing probems: () Sove N F/Q (x) = k for x o F, k Z >0, using the absoute method. From Lemma 3. and Theorem 3.3 we obtain bounds (after an LLL-reduction of the unit attice): i L i = U i This requires us to consider eipsoids if a soutions of the norm equation are to be determined. (2) Sove N F/E (x) = θ for x o F and θ o E by (a) the absoute method: Lemma 4. and Theorem 3.3 yied bounds: i L i U i A tota of eipsoids need to be enumerated. (b) the reative method: Over E () we have the reative signature t = s =. By Lemma 4. and Theorem 4.2 we obtain: j 2 3 i L (j) i = U (j) i that is a tota of ony 4289 eipsoids. We note that we actuay computed a these eipsoids, so that the numbers L i, U i given above are exact. 8. Tabes A computations were carried out on an HP9000/735s with 96MB memory using software deveoped under KANT V4 [5]. The operating system on the machine was HP-UX 9.0. F. Grunewad in Düssedorf asked us the foowing question: Let E = Q(α), α {, 2, 3},β= a+bα for 3 a, b 3 such that F = E(β) isof absoute degree 4. How many nonassociate soutions do exist for (4.2) with sma θ, where nonassociate means moduo units of reative norm? In Tabes and 2 we can ony present a part of our resuts because of imited space. Another probem of W. Pesken in Aachen and H. Brückner in Hamburg was to sove severa norm equations in cycotomic fieds (see Tabe 3 for detais). We present soutions of (4.), meaning that no soutions exist, and ζ k denoting a primitive kth root of unity. License or copyright restrictions may appy to redistribution; see

10 408 C. FIEKER, A. JURK, AND M. POHST Tabe. Exampes E F θ soution re. norm Q( 2) E( ) 2 ( ) (2 2( + )) 2 2 (2 + 2( + )) 2 ± 2 3 ± (2 ± ( ) 2 2 ) ± 2 2 Q( 3) E( ) (± 3) 2 ( 3) 2 ( 2+ (± 3)) 2 ( 3) 2 ( + 3 ( + 3)) 2 ( (± 3)) (2 3+ (± 3)) 2 ( 7 3) 2 ( ) 2 (5 3 3) 2 (4 + ( +3 3)) 2 ( 5 3 3) 2 ( ( 2+2 3)) 2 (7 + 3) ( ) 2 (7 3) ( + 3) 2 3 License or copyright restrictions may appy to redistribution; see

11 SOLVING RELATIVE NORM EQUATIONS IN ALGEBRAIC NUMBER FIELDS 409 Tabe 2. Exampes E F β soution re. norm Q( ) E( 2) + 2 2(± ) ± 2 (2 + 2( + )) ± 2 (±2+ 2(± )) ± 2 ± (± 3 ) 4 ± 3 2( 3 ± ) 4 ± 3 ±2 3 4 ± p E( ) p p + p + + p p + E( p + ) p + + p + + p + p + p p + License or copyright restrictions may appy to redistribution; see

12 40 C. FIEKER, A. JURK, AND M. POHST Tabe 3. Exampes E F β soution Q(ζ 5 ) R Q(ζ 5 ) 2 ζ 4 ζ 5 ζ 5 Q(ζ 7 ) R Q(ζ 7 ) 2 ζ 7 ζ 7 ζ 7 Q(ζ 8 ) R Q(ζ 8 ) 2 ζ 8 ζ 8 ζ 8 Q(ζ 9 ) R Q(ζ 9 ) 2 ζ 9 ζ 9 ζ 9 Q(ζ 0 ) R Q(ζ 0 ) 2 ζ 0 ζ 0 ζ 0 Q(ζ ) R Q(ζ ) 2 ζ ζ ζ Q(ζ 2 ) R Q(ζ 2 ) 2 ζ 2 ζ 2 ζ 2 Q(ζ 6 ) R Q(ζ 6 ) Q(ζ 32 ) R Q(ζ 32 ) Q(ζ 64 ) R Q(ζ 64 ) References. U. Fincke, Ein Eipsoidverfahren zur Lösung von Normgeichungen in agebraischen Zahkörpern, Thesis, Düssedorf U. Fincke, M. Pohst, A Procedure for Determining Agebraic Integers of Given Norm, Proceedings EUROCAL 83, Springer Lecture Notes in Computer Science No. 62, Springer, 983. MR 86k: U. Fincke, M. Pohst, Improved methods for Cacuating Vectors of Short Length in a Lattice, Incuding a Compexity Anaysis, Math. Comp. 44, (985). MR 86e: A. Jurk, Über die Berechnung von Lösungen reativer Normgeichungen in agebraischen Zahkörpern, Thesis, Düssedorf The Kant-Group, KANT V4, to appear in J. Symb. Comput. 6. K. Maher, Inequaities for Idea Bases in Agebraic Number Fieds, J. Austra. Math. Soc. 4, (964). MR 3: W. Narkiewicz, Eementary and Anaytic Theory of Agebraic Numbers, 2nd ed., Springer, 990. MR 9h:07 8. M. Pohst, H. Zassenhaus, Agorithmic Agebraic Number Theory, Encycopaedia of Mathematics and its Appications, Cambridge University Press, 989. MR 92b:074 Fachbereich 3 Mathematik, Sekretariat MA 8, Technische Universität Berin, Straße des 7. Juni 36, D-0623 Berin, Germany E-mai address: fieker@math.tu-berin.de Desdorfer Weg 5, 508 Bedburg, Germany Fachbereich 3 Mathematik, Sekretariat MA 8, Technische Universität Berin, Straße des 7. Juni 36, D-0623 Berin, Germany E-mai address: pohst@math.tu-berin.de License or copyright restrictions may appy to redistribution; see

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

Algorithms to solve massively under-defined systems of multivariate quadratic equations

Algorithms to solve massively under-defined systems of multivariate quadratic equations Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische

More information

A proposed nonparametric mixture density estimation using B-spline functions

A proposed nonparametric mixture density estimation using B-spline functions A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),

More information

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric

More information

A lower regulator bound for number fields

A lower regulator bound for number fields A lower regulator bound for number fields Claus Fieker and Michael E. Pohst Abstract 1 Introduction In [] a regulator bound for algebraic number fields F was derived which makes use of specific data of

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Small generators of function fields

Small generators of function fields Journa de Théorie des Nombres de Bordeaux 00 (XXXX), 000 000 Sma generators of function fieds par Martin Widmer Résumé. Soit K/k une extension finie d un corps goba, donc K contient un éément primitif

More information

Efficient Algorithms for Pairing-Based Cryptosystems

Efficient Algorithms for Pairing-Based Cryptosystems CS548 Advanced Information Security Efficient Agorithms for Pairing-Based Cryptosystems Pauo S. L. M. Barreto, HaeY. Kim, Ben Lynn, and Michae Scott Proceedings of Crypto, 2002 2010. 04. 22. Kanghoon Lee,

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

An Extension of Almost Sure Central Limit Theorem for Order Statistics

An Extension of Almost Sure Central Limit Theorem for Order Statistics An Extension of Amost Sure Centra Limit Theorem for Order Statistics T. Bin, P. Zuoxiang & S. Nadarajah First version: 6 December 2007 Research Report No. 9, 2007, Probabiity Statistics Group Schoo of

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Continued fractions with low complexity: Transcendence measures and quadratic approximation. Yann BUGEAUD

Continued fractions with low complexity: Transcendence measures and quadratic approximation. Yann BUGEAUD Continued fractions with ow compexity: Transcendence measures and quadratic approximation Yann BUGEAUD Abstract. We estabish measures of non-quadraticity and transcendence measures for rea numbers whose

More information

A UNIVERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIVE ALGEBRAIC MANIFOLDS

A UNIVERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIVE ALGEBRAIC MANIFOLDS A UNIERSAL METRIC FOR THE CANONICAL BUNDLE OF A HOLOMORPHIC FAMILY OF PROJECTIE ALGEBRAIC MANIFOLDS DROR AROLIN Dedicated to M Saah Baouendi on the occasion of his 60th birthday 1 Introduction In his ceebrated

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Volume 13, MAIN ARTICLES

Volume 13, MAIN ARTICLES Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi

More information

CONGRUENCES. 1. History

CONGRUENCES. 1. History CONGRUENCES HAO BILLY LEE Abstract. These are notes I created for a seminar tak, foowing the papers of On the -adic Representations and Congruences for Coefficients of Moduar Forms by Swinnerton-Dyer and

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

Selmer groups and Euler systems

Selmer groups and Euler systems Semer groups and Euer systems S. M.-C. 21 February 2018 1 Introduction Semer groups are a construction in Gaois cohomoogy that are cosey reated to many objects of arithmetic importance, such as cass groups

More information

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

MIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI

MIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI MIXING AUTOMORPHISMS OF COMPACT GROUPS AND A THEOREM OF SCHLICKEWEI KLAUS SCHMIDT AND TOM WARD Abstract. We prove that every mixing Z d -action by automorphisms of a compact, connected, abeian group is

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

THE PARTITION FUNCTION AND HECKE OPERATORS

THE PARTITION FUNCTION AND HECKE OPERATORS THE PARTITION FUNCTION AND HECKE OPERATORS KEN ONO Abstract. The theory of congruences for the partition function p(n depends heaviy on the properties of haf-integra weight Hecke operators. The subject

More information

On generalized quantum Turing machine and its language classes

On generalized quantum Turing machine and its language classes Proceedings of the 11th WSEAS Internationa onference on APPLIED MATHEMATIS, Daas, Texas, USA, March -4, 007 51 On generaized quantum Turing machine and its anguage casses SATOSHI IRIYAMA Toyo University

More information

Unconditional security of differential phase shift quantum key distribution

Unconditional security of differential phase shift quantum key distribution Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice

More information

MONTE CARLO SIMULATIONS

MONTE CARLO SIMULATIONS MONTE CARLO SIMULATIONS Current physics research 1) Theoretica 2) Experimenta 3) Computationa Monte Caro (MC) Method (1953) used to study 1) Discrete spin systems 2) Fuids 3) Poymers, membranes, soft matter

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Generalized Bell polynomials and the combinatorics of Poisson central moments

Generalized Bell polynomials and the combinatorics of Poisson central moments Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various

More information

arxiv: v1 [math.fa] 23 Aug 2018

arxiv: v1 [math.fa] 23 Aug 2018 An Exact Upper Bound on the L p Lebesgue Constant and The -Rényi Entropy Power Inequaity for Integer Vaued Random Variabes arxiv:808.0773v [math.fa] 3 Aug 08 Peng Xu, Mokshay Madiman, James Mebourne Abstract

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

Transcendence of stammering continued fractions. Yann BUGEAUD

Transcendence of stammering continued fractions. Yann BUGEAUD Transcendence of stammering continued fractions Yann BUGEAUD To the memory of Af van der Poorten Abstract. Let θ = [0; a 1, a 2,...] be an agebraic number of degree at east three. Recenty, we have estabished

More information

JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF SEVERAL VARIABLES

JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF SEVERAL VARIABLES PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Voume 128, Number 7, Pages 2075 2084 S 0002-99390005371-5 Artice eectronicay pubished on February 16, 2000 JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF

More information

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

Mixed Volume Computation, A Revisit

Mixed Volume Computation, A Revisit Mixed Voume Computation, A Revisit Tsung-Lin Lee, Tien-Yien Li October 31, 2007 Abstract The superiority of the dynamic enumeration of a mixed ces suggested by T Mizutani et a for the mixed voume computation

More information

Course 2BA1, Section 11: Periodic Functions and Fourier Series

Course 2BA1, Section 11: Periodic Functions and Fourier Series Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........

More information

Fitting affine and orthogonal transformations between two sets of points

Fitting affine and orthogonal transformations between two sets of points Mathematica Communications 9(2004), 27-34 27 Fitting affine and orthogona transformations between two sets of points Hemuth Späth Abstract. Let two point sets P and Q be given in R n. We determine a transation

More information

Smoothers for ecient multigrid methods in IGA

Smoothers for ecient multigrid methods in IGA Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third

More information

On prime divisors of remarkable sequences

On prime divisors of remarkable sequences Annaes Mathematicae et Informaticae 33 (2006 pp. 45 56 http://www.ektf.hu/tanszek/matematika/ami On prime divisors of remarkabe sequences Ferdinánd Fiip a, Kámán Liptai b1, János T. Tóth c2 a Department

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

Powers of Ideals: Primary Decompositions, Artin-Rees Lemma and Regularity

Powers of Ideals: Primary Decompositions, Artin-Rees Lemma and Regularity Powers of Ideas: Primary Decompositions, Artin-Rees Lemma and Reguarity Irena Swanson Department of Mathematica Sciences, New Mexico State University, Las Cruces, NM 88003-8001 (e-mai: iswanson@nmsu.edu)

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

Homogeneity properties of subadditive functions

Homogeneity properties of subadditive functions Annaes Mathematicae et Informaticae 32 2005 pp. 89 20. Homogeneity properties of subadditive functions Pá Burai and Árpád Száz Institute of Mathematics, University of Debrecen e-mai: buraip@math.kte.hu

More information

b n n=1 a n cos nx (3) n=1

b n n=1 a n cos nx (3) n=1 Fourier Anaysis The Fourier series First some terminoogy: a function f(x) is periodic if f(x ) = f(x) for a x for some, if is the smaest such number, it is caed the period of f(x). It is even if f( x)

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SIAM J. NUMER. ANAL. Vo. 0, No. 0, pp. 000 000 c 200X Society for Industria and Appied Mathematics VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW

More information

Minimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch

Minimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch The Eighth Internationa Symposium on Operations Research and Its Appications (ISORA 09) Zhangiaie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp. 402 408 Minimizing Tota Weighted Competion

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

Factorization of Cyclotomic Polynomials with Quadratic Radicals in the Coefficients

Factorization of Cyclotomic Polynomials with Quadratic Radicals in the Coefficients Advances in Pure Mathematics, 07, 7, 47-506 http://www.scirp.org/journa/apm ISSN Onine: 60-0384 ISSN Print: 60-0368 Factoriation of Cycotomic Poynomias with Quadratic Radicas in the Coefficients Afred

More information

PRIME TWISTS OF ELLIPTIC CURVES

PRIME TWISTS OF ELLIPTIC CURVES PRIME TWISTS OF ELLIPTIC CURVES DANIEL KRIZ AND CHAO LI Abstract. For certain eiptic curves E/Q with E(Q)[2] = Z/2Z, we prove a criterion for prime twists of E to have anaytic rank 0 or 1, based on a mod

More information

Restricted weak type on maximal linear and multilinear integral maps.

Restricted weak type on maximal linear and multilinear integral maps. Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0 Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

On norm equivalence between the displacement and velocity vectors for free linear dynamical systems

On norm equivalence between the displacement and velocity vectors for free linear dynamical systems Kohaupt Cogent Mathematics (5 : 95699 http://dxdoiorg/8/33835595699 COMPUTATIONAL SCIENCE RESEARCH ARTICLE On norm equivaence between the dispacement and veocity vectors for free inear dynamica systems

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Lower Bounds for the Relative Greedy Algorithm for Approximating Steiner Trees

Lower Bounds for the Relative Greedy Algorithm for Approximating Steiner Trees This paper appeared in: Networks 47:2 (2006), -5 Lower Bounds for the Reative Greed Agorithm for Approimating Steiner Trees Stefan Hougard Stefan Kirchner Humbodt-Universität zu Berin Institut für Informatik

More information

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence

More information

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27 Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that

More information

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY

School of Electrical Engineering, University of Bath, Claverton Down, Bath BA2 7AY The ogic of Booean matrices C. R. Edwards Schoo of Eectrica Engineering, Universit of Bath, Caverton Down, Bath BA2 7AY A Booean matrix agebra is described which enabes man ogica functions to be manipuated

More information

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry

More information