2.9 Incomplete Elliptic Integrals
|
|
- Holly Melanie Todd
- 6 years ago
- Views:
Transcription
1 .9 Inomplete Ellipti Integrals A. Purpose An integral of the form R t, P t) 1/) dt 1) in whih Pt) is a polynomial of the third or fourth degree that has no multiple roots, and R is a rational funtion of t and Pt) 1/, is either elementary, or is an ellipti integral. It is always possible to express integrals of the form of Eq. 1) linearly in terms of elementary funtions and three ellipti integrals of anonial form. These funtions are desribed more ompletely in [1] and []. Several anonial forms have been proposed, but the most widely used are due to Jaobi, Legendre and Carlson. In eah of Eqs. ) 4) we present first Jaobi s and then Legendre s form of the anonial ellipti integrals: F ϕ, k) = = Eϕ, k) = = y ϕ y ϕ 1 t ) 1/ 1 k t ) 1/ dt 1 k sin θ ) 1/ dθ ) 1 t ) 1/ 1 k t ) 1/ dt 1 k sin θ ) 1/ dθ 3) Πϕ, α, k) y = 1 α t ) 1 1 t ) 1/ 1 k t ) 1/ dt = ϕ 1 α sin θ ) 1 1 k sin θ ) 1/ dθ 4) in whih y = sin ϕ. If ϕ is equal to π/, the integrals are said to be omplete, otherwise they are inomplete. Carlson s forms of the anonial ellipti integrals are R D a, b, ) = 3 t + a) 1/ t + b) 1/ t + ) 3/ dt 5) in whih a and b are nonnegative suh that a + b > and is positive; if either a or b is zero, the integral is omplete, otherwise it is inomplete, R F a, b, ) = 1 t + a) 1/ t + b) 1/ t + ) 1/ dt 6) 1997 Calif. Inst. of Tehnology, 15 Math à la Carte, In. in whih a, b and are nonnegative and at most one of them is zero; if one of a, b or is zero, the integral is omplete, otherwise it is inomplete, and R J a, b,, r) = 3 t + r) 1 t + a) 1/ t + b) 1/ t + ) 1/ dt 7) in whih a, b and are nonnegative, and at most one of them is zero, and r is nonzero; if one of a, b or is zero, the integral is omplete, otherwise it is inomplete. Notie that R D a, b, ) = R J a, b,, ). But the neessity to ompute R D a, b, ) arises frequently in pratie, and a proedure espeially tailored to ompute R D a, b, ) is more effiient than omputing R J a, b,, ). The funtion R C a, b) = R F a, b, b) is elementary, but also appears frequently. A proedure is provided to ompute R C a, b). Identify a, b and suh that a b, and assume a <. Then 3/ 3 R D a, b, ) = k sin 3 [F ϕ, k) Eϕ, k)] 8) ϕ 1/ F ϕ, k) R F a, b, ) = 9) sin ϕ 3/ 3 R J a, b,, r) = α sin 3 [Πϕ, α, k) F ϕ, k)] ϕ 1) where os ϕ = a/, k = b)/ a) and α = r)/ a). The subprograms desribed in this hapter evaluate the anonial forms of inomplete ellipti integrals, using either the Legendre or the Carlson parameterization. B. Usage B.1 Program Prototype, Single Preision, Legendre s Form, E and F PHI, K, F, E Assign values to PHI and K. B.1.a CALL SELEFI PHI, K, F, E, ) PHI [in] Argument, ϕ, of the ellipti integral. Require PHI π/. K [in] Modulus, k. Require K 1.. July 11, 15 Inomplete Ellipti Integrals.9 1
2 F [out] Fϕ, k) with ϕ given by PHI and k given by K. E [out] Eϕ, k) with ϕ given by PHI and k given by K. [out] status indiator: = no errors 1 = Magnitude of argument too large, PHI > π/. = Magnitude of Modulus too large, K > = PHI = π/ and K = 1, F is infinite. B. Program Prototype, Single Preision, Legendre s Form, Π PHI, K, ALPHA, PI Assign values to PHI, K and ALPHA. CALL SELPII PHI, K, ALPHA, PI, ) B..a PHI [in] Argument, ϕ, of the ellipti integral. Require PHI π/. K [in] Square of the modulus, k. Require k sin ϕ 1.. See Setion E. ALPHA [in] Charateristi, α. Require α sin ϕ 1.. See Setion E. PI [out] Πϕ, α, k), with ϕ given by PHI, α given by ALPHA and k given by K. [out] status indiator. If =, there were no errors. Other values are produed by proedures SRFVAL and SRJVAL see Setions B.5 and B.6) whih are used in omputing Πφ, α, k). B.3 Program Prototype, Single Preision, Carlson s Form, R C X, Y, RC Assign values to X and Y. B.3.a CALL SRCVAL X, Y, RC, ) X, Y [in] Arguments of the ellipti integral. Require X, Y. See Setion E. RC [out] The omputed value of R C X, Y). [out] Status indiator: = no errors 1 = X <. or Y =.. = X + Y too small See Setion E). 3 = X or Y or X + Y too large See Setion E). 4 = Y < and Y too large and X too small See Setion E). B.4 Program Prototype, Single Preision, Carlson s Form, R D X, Y, Z, RD Assign values to X, Y and Z. B.4.a CALL SRDVAL X, Y, Z, RD, ) X, Y, Z [in] Arguments of the ellipti integral. Require X, Y, X + Y >, Z >. See Setion E. RD [out] The omputed value of R D X, Y, Z). [out] Status indiator: = no errors 1 = X <. or Y <. or Z <.. = X + Y too small or Z too small See Setion E). 3 = X or Y or Z too large See Setion E). B.5 Program Prototype, Single Preision, Carlson s Form, R F X, Y, Z, RF Assign values to X, Y and Z. B.5.a CALL SRFVAL X, Y, Z, RF, ) X, Y, Z [in] Arguments of the ellipti integral. Require X, Y, Z, at most one of X, Y or Z equal zero. See Setion E. RF [out] The omputed value of R F X, Y, Z). [out] Status indiator: = no errors 1 = X <. or Y <. or Z <.. = X + Y or X + Z or Y + Z too small See Setion E). 3 = X or Y or Z too large See Setion E). B.6 Program Prototype, Single Preision, Carlson s Form, R J X, Y, Z, R, RJ Assign values to X, Y, Z and R. CALL SRJVAL X, Y, Z, R, RJ, ).9 Inomplete Ellipti Integrals July 11, 15
3 B.6.a X, Y, Z, R [in] Arguments of the ellipti integral. Require X, Y, Z, at most one of X, Y or Z equal zero, R. See Setion E. RJ [out] The omputed value of R J X, Y, Z, R). [out] Status indiator: = no errors 1 = X <. or Y <. or Z <. or R =.. = X + Y or X + Z or Y + Z or R too small See Setion E). 3 = X or Y or Z or R too large See Setion E). B.7 Modifiations for Double Preision For double preision usage, hange the type statements to DOUBLE PRECISION and hange the subprogram names SELEFI, SELPII, SRCVAL, SRDVAL, SRFVAL and SRJVAL to DELEFI, DELPII, DRCVAL, DRDVAL, DRFVAL and DRJVAL, respetively. π Λ α, β) = Λα, β, π/) [ = sin β R F, os α, 1) 1 ] 3 sin αr D, os α, 1) R F os β, 1 os α sin β, 1 ) 1 3 os α sin 3 β R F, os α, 1) R D os β, 1 os α sin β, 1) 1) The variants of Legendre s integrals used by Bulirsh in [4] and [5] are el1x, k ) = xr F x, 1 + x ), 13) elx, k, a, b) = axr F x, 1 + x ) b a)x3 R D x, 1 + x ) 14) ele3x, k, p) = xr F x, 1 + x ) p)x3 R J x, 1 + x, 1 + px ) 15) C. Examples and Remarks C.1 Related Funtions Logarithms, inverse irular funtions and inverse hyperboli funtions an be expressed in terms of R C, see [9, pp. 163, 186]: ln x)/x 1) = R C x), x), x > ; sin 1 x)/x = R C 1 x, 1), 1 x 1; sinh 1 x)/x = R C 1 + x, 1), < x < ; os 1 x)/1 x ) 1 = R C x, 1), x 1; osh 1 x)/x 1) 1 = R C x, 1), x 1; tan 1 x)/x = R C 1, 1 + x ), < x < ; tanh 1 x)/x = R C 1, 1 x ), 1 < x < 1; ot 1 x = R C x, x + 1), x < ; oth 1 x = R C x, x 1), x > 1. The first seven of these allow omputing nearly indeterminate forms with more auray than would be possible using the naïve formulation. Heuman s lambda funtion [3] is a variant of Legendre s third integral: elk, p, a, b) = ar F, k, 1) = 1 3 b pa)r J, k, 1, p ) C. Whih Proedure Should Be Used? 16) Several fators influene the hoie of proedure. If one needs to write a simple program and use it one, one should probably hoose the proedure that evaluates the funtions in the form most similar to the way the problem is posed. If one needs to write a program that will have substantial use, one should usually prefer SE- LEFI to SRDVAL and SRFVAL, as the former is up to 3 times faster than the latter two. An exeption to this rule ours if one needs to ompute R D a, b, ) with < maxa, b), in whih ase the parameters for SELEFI will be out of range. If auray is an issue but speed is not, one may prefer SRDVAL and SRFVAL to SELEFI, at least for omputing Fϕ, k). See testing in Setion D below). SRCVAL is somewhat slower, on an IBM PC/AT with a 1 os α sin β ) numeri data proessor, than using the equivalent Fortran intrinsi funtions. This is no surprise, as most 1/ os Λα, β, ϕ) of the intrinsi funtions are implemented by hardware. α sin β os β But the inverse hyperboli funtions are not. SRCVAL = sin ϕ R F os ϕ, 1 sin α sin ϕ, 1) is roughly the same speed as the proedures in Chapter.1. As mentioned above, it may be advantageous to R J os ϕ, use SRCVAL to ompute nearly indeterminate forms. sin α sin 3 ϕ os α sin β ) 1 sin α sin ϕ, 1, 1 sin α sin ϕ ) SELPII is implemented by using SRJVAL and SRFVAL 1 os α sin 11) see Eqs. 9) and 1) above). Thus, there is no speial β July 11, 15 Inomplete Ellipti Integrals.9 3
4 advantage in speed or auray to one or the other. The sole riterion is how losely the forms of the funtions evaluated diretly by the proedures math the forms of the funtions the user needs to evaluate. D. Funtional Desription D.1 Properties of the Funtions The first form given in Eqs. ) 4) is the Jaobi or algebrai form. When expressed in this form Eq. ) is finite for all real and omplex y, inluding, has a simple pole of order 1 for y =, and is logarithmially infinite for y = 1/α. D. Method of Computation The proedure SELEFI is based upon a proedure ELLPI developed by Allan V. Hershey and modified by Alfred H. Morris, desribed in [6]. The proedure uses series expansions due to DiDonato and Hershey, desribed in [7]. The proedure SELPII is based upon a proedure EPI developed by Alfred H. Morris, desribed in [5]. It omputes Πϕ, k, α ) using Eqs. 9) and 1), as omputed by SRFVAL and SRJVAL. The proedures SRCVAL, SRDVAL, SRFVAL and SRJVAL are based on proedures developed by B. C. Carlson and Elaine M. Notis, desribed in [8] and [9]. All of the referened proedures were revised to be onsistent with low level modules and naming onventions of MATH77. D.3 Testing The single preision programs for Eϕ, k), Fϕ, k), R D a, b, ) and R F a, b, ) were tested on an IBM PC/AT using IEEE arithmeti) by omparison to double preision results, as desribed below. The relative preision of IEEE single preision arithmeti is ρ = The auray of proedure SELEFI was assessed by omparing its results to double preision results obtained by applying Eqs. 8) and 9), with R D a, b, ) and R F a, b, ) evaluated by DRDVAL and DRFVAL, respetively. The auray of proedures SRDVAL and SRFVAL was assessed by omparing their results to double preision results obtained by applying Eqs. 8) and 9), with Eϕ, k) and Fϕ, k) evaluated by DELEFI. To test SELEFI, the retangular region ϕ π/ k 1 of the ϕ k plane was divided into regions, and a point was randomly seleted in eah region. To test SRDVAL and SRFVAL, the argument was set to 1., the retangular region a < 1 b < 1 of the a b plane was divided into regions, and a point was randomly seleted in eah region. The maximum relative and absolute errors are summarized in the following table. Max. Rel. Max. Abs. Funtion Error Error Eϕ, k).8ρ.98ρ Fϕ, k) 5.4ρ 15.91ρ R D a, b, 1) 1.ρ 3.1ρ R F a, b, 1) 1.35ρ.55ρ Errors in Fϕ, k) inrease as the arguments approah the infinite singularity at ϕ = π/ and k = 1. Referenes 1. Milton Abramowitz and Irene A. Stegun, Handbook of Mathematial Funtions, Applied Mathematis Series 55, National Bureau of Standards 1966) Chapter 17, Paul F. Byrd and Morris D. Friedman, Handbook of Ellipti Integrals for Engineers and Sientists, Springer Verlag, Berlin 1971). 3. H. Kuki, Tables of omplete ellipti integrals, J. Math. and Physis 1941) Roland Bulirsh, Numerial alulation of ellipti integrals and ellipti funtions, Numerishe Mathematik ) Roland Bulirsh, Numerial alulation of ellipti integrals and ellipti funtions, Numerishe Mathematik ) Alfred. H. Morris, Jr., NSWC Library of Mathematis Subroutines. Tehnial Report NSWCDD/TR-9/45, Naval Surfae Warfare Center, Dahlgren, VA USA Jan. 1993) Armido R. DiDonato and Allan V. Hershey, New formulas for omputing inomplete ellipti integrals of the first and seond kind, J. ACM ) B. C. Carlson, Computing ellipti integrals by dupliation, Numerishe Mathematik ) B. C. Carlson and Elaine M. Notis, Algorithm 577: Algorithms for inomplete ellipti integrals [S1], ACM Trans. on Math. Software 7, 3 Sept. 1981) B. C. Carlson, Speial Funtions of Applied Mathematis, Aademi Press, New York 1977). E. Error Proedures and Restritions The proedure SELEFI requires ϕ π/, and k 1. Proedure SELPII omputes Πϕ, k, α ) from R J a, b,, r) and R F a, b, ) using Eqs. 9) and 1). The initial values for the arguments are a = os ϕ, b = 1 k sin ϕ, r = 1 α sin ϕ, and = maxa, b, r). Then a, b and r are replaed by a, b and r, respetively. SELPII requires ϕ π/. Restritions on k and α are enfored indiretly by restritions on.9 4 Inomplete Ellipti Integrals July 11, 15
5 a, b and imposed by SRFVAL and SRJVAL, desribed below. The ranges for ϕ and k an be extended using formulae 113.1, 113., 114.1, 115.1, 115., 16., 161. and 16. from [], or formulae through from [1]. Denote the largest representable magnitude by Ω, and the smallest nonzero representable magnitude by ω. General restritions on the arguments to proedures SR- CVAL, SRDVAL, SRFVAL and SRJVAL were desribed above in Setion B. SRCVAL requires X + Y 5ω, X Ω/5, Y Ω/5, and, if Y <.36/ ω it requires X ωω) /5. Denote the mahine round-off level by ρ, that is, ρ is the smallest positive number suh that the representation of 1 + ρ is different from 1. Let ε be the solution of the equation ρ = 3ε 6 1 ε) 3/, Ω D = Ω /3 and ω D = εω /3 /1. SRDVAL requires X + Y ω D and Z ω D, X Ω D, Y Ω D and Z Ω D. SRFVAL requires X + Y 5ω, X + Z 5ω, Y + Z 5ω, X Ω/5, Y Ω/5 and Z Ω/5. Let Ω J = Ω/5) 1/3 /5 and ω J = 5ω) 1/3. SRJVAL requires X + Y ω J, Y + Z ω J, X + Z ω J, R ω J, X Ω J, Y Ω J, Z Ω J and R Ω J. The aessible ranges of the arguments may be extended beyond the ranges admissible in the proedures by using the homogeneity of the funtions: R F ka, kb, k) = k 1/ R F a, b, ), R J ka, kb, k, kr) = k 3/ R J a, b,, r). and If any of the restritions above is violated, all proedures return an error indiator in the argument named, and invoke the error message proessor see Chapter 19.) with LEVEL =. The proedure ERMSET see Chapter 19.) may be used to affet the default error proessing ation. F. Supporting Information The soure language for these subroutines is ANSI Fortran 77. The proedures SELEFI and SELPII were written by W. V. Snyder in Deember 199, based on earlier proedures desribed by Alfred H. Morris, Naval Surfae Warfare Center, Dahlgren, VA in [5]. The proedures SRCVAL, SRDVAL, SRFVAL and SRJVAL were written by W. V. Snyder in Deember 199, based on earlier proedures desribed by Carlson and Notis in [9]. Entry Required Files DELEFI AMACH, DELEFI, DERM1, DERV1, DLNREL, ERFIN, ERMSG DELPII AMACH, DELPII, DERM1, DERV1, DRCVAL, DRFVAL, DRJVAL, ERFIN, ERMSG DRCVAL AMACH, DERM1, DERV1, DRCVAL, ERFIN, ERMSG DRDVAL AMACH, DERM1, DERV1, DRDVAL, ERFIN, ERMSG DRFVAL AMACH, DERM1, DERV1, DRFVAL, ERFIN, ERMSG DRJVAL AMACH, DERM1, DERV1, DRCVAL, DRFVAL, DRJVAL, ERFIN, ERMSG SELEFI AMACH, ERFIN, ERMSG, SELEFI, SERM1, SERV1, SLNREL SELPII AMACH, ERFIN, ERMSG, SELPII, SERM1, SERV1, SRCVAL, SRFVAL, SRJVAL SRCVAL AMACH, ERFIN, ERMSG, SERM1, SERV1, SRCVAL SRDVAL AMACH, ERFIN, ERMSG, SERM1, SERV1, SRDVAL SRFVAL AMACH, ERFIN, ERMSG, SERM1, SERV1, SRFVAL SRJVAL AMACH, ERFIN, ERMSG, SERM1, SERV1, SRCVAL, SRFVAL, SRJVAL July 11, 15 Inomplete Ellipti Integrals.9 5
6 DRSELI program DRSELI >> DRSELI Krogh Changes t o use M77CON >> DRSELI WV Snyder Create s e p a r a t e s i n g l e and d o u b l e demos. >> DRSELI WV Snyder JPL O r i g i n a l ode. S r e p l a e s? : DR?ELI,?RCVAL,? ELEFI,? ELPII,?RDVAL,?RFVAL,?RJVAL Demonstration d r i v e r f o r inomplete e l l i p t i i n t e g r a l proedures. real ALPHA, E, F, K, K, PHI, PI, R, RC, RD, RF, RJ real SINPHI, T, U, X, Y, Z integer Compute ar s i n e x using ASIN and RC, f o r x =.5 print, I d e n t i t i e s from write up : x =. 5 e a l l s r v a l 1. e x x, 1. e, r, i e r r ) i f i e r r. eq. ) then t = asin x ) x r print ASIN. 5 ). 5 RC1.5,1) =, g15. 8 ), t else print SRCVAL r e t u r n s e r r o r s i g n a l, i 1 ), i e r r Evaluate i d e n t i t i e s g i v e n by e q u a t i o n s 8 1) in t h e write up with k = 1/, s i n phi ) = 1/4, alpha = 1/, = 1. From t h i s, we have a = 3/4, b = r = 7/8. alpha =. 5 e k = sqrt. 5 e ) k =. 5 e s i n p h i =. 5 e phi = asin s i n p h i ) r =.875 e x =. 7 5 e y =.875 e z = 1. e a l l s e l e f i phi, k, f, e, i e r r ) i f i e r r. ne. ) then print SELEFI r e t u r n s e r r o r s i g n a l, i 1 ), i e r r a l l s e l p i i phi, k, alpha, pi, i e r r ) i f i e r r. ne. ) then print SELPII r e t u r n s e r r o r s i g n a l, i 1 ), i e r r a l l s r d v a l x, y, z, rd, i e r r ) i f i e r r. ne. ) then print SRDVAL r e t u r n s e r r o r s i g n a l, i 1 ), i e r r a l l s r f v a l x, y, z, r f, i e r r ) i f i e r r. ne. ) then print SRFVAL r e t u r n s e r r o r s i g n a l, i 1 ), i e r r.9 6 Inomplete Ellipti Integrals July 11, 15
7 a l l s r j v a l x, y, z, r, r j, i e r r ) i f i e r r. ne. ) then print SRJVAL r e t u r n s e r r o r s i g n a l, i 1 ), i e r r u = sqrt z 3) rd t = 3. e / k s i n p h i 3) f e ) r = u t ) / u print Equation 8 ), LHS RHS)/LHS =, g15. 8 ), r u = sqrt z ) r f t = f / s i n p h i r = u t ) / u print Equation 9 ), LHS RHS)/LHS =, g15. 8 ), r u = sqrt z 3) r j t = 3 / alpha s i n p h i 3) p i f ) r = u t ) / u print Equation 1 ), LHS RHS)/LHS =, g15. 8 ), r 99 stop end ODSELI I d e n t i t i e s from write up : ASIN. 5 ). 5 RC1.5,1) = E 7 Equation 8 ), LHS RHS)/LHS = E 5 Equation 9 ), LHS RHS)/LHS =. Equation 1 ), LHS RHS)/LHS = E 5 July 11, 15 Inomplete Ellipti Integrals.9 7
7.1 Roots of a Polynomial
7.1 Roots of a Polynomial A. Purpose Given the oeffiients a i of a polynomial of degree n = NDEG > 0, a 1 z n + a 2 z n 1 +... + a n z + a n+1 with a 1 0, this subroutine omputes the NDEG roots of the
More information11.1 Polynomial Least-Squares Curve Fit
11.1 Polynomial Least-Squares Curve Fit A. Purpose This subroutine determines a univariate polynomial that fits a given disrete set of data in the sense of minimizing the weighted sum of squares of residuals.
More information3.2 Gaussian (Normal) Random Numbers and Vectors
3.2 Gaussian (Normal) Random Numbers and Vetors A. Purpose Generate pseudorandom numbers or vetors from the Gaussian (normal) distribution. B. Usage B.1 Generating Gaussian (normal) pseudorandom numbers
More information4.3 Singular Value Decomposition and Analysis
4.3 Singular Value Deomposition and Analysis A. Purpose Any M N matrix, A, has a Singular Value Deomposition (SVD) of the form A = USV t where U is an M M orthogonal matrix, V is an N N orthogonal matrix,
More informationChapter IX. Roots of Equations
14 Chapter IX Roots of Equations Root solving typially onsists of finding values of x whih satisfy a relationship suh as: 3 3 x + x = 1 + e x The problem an be restated as: Find the root of 3 f ( x) =
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationarxiv:math/ v1 [math.ca] 6 Sep 1994
NUMERICAL COMPUTATION OF REAL OR COMPLEX arxiv:math/909227v1 [math.ca] 6 Sep 199 ELLIPTIC INTEGRALS B. C. Carlson Ames Laboratory and Department of Mathematics, Iowa State University, Ames, Iowa 50011-020,
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More information6.3 Basic Linear Algebra Subprograms (BLAS1)
6.3 Basi Linear Algebra Subprograms (BLAS1) A. Purpose This is a set of subroutine and funtion subprograms for basi mathematial operations on a single vetor or a pair of vetors. The operations provided
More informationA EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM.
A EUCLIDEAN ALTERNATIVE TO MINKOWSKI SPACETIME DIAGRAM. S. Kanagaraj Eulidean Relativity s.kana.raj@gmail.om (1 August 009) Abstrat By re-interpreting the speial relativity (SR) postulates based on Eulidean
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationChapter 2: Solution of First order ODE
0 Chapter : Solution of irst order ODE Se. Separable Equations The differential equation of the form that is is alled separable if f = h g; In order to solve it perform the following steps: Rewrite the
More informationLOGISTIC REGRESSION IN DEPRESSION CLASSIFICATION
LOGISIC REGRESSIO I DEPRESSIO CLASSIFICAIO J. Kual,. V. ran, M. Bareš KSE, FJFI, CVU v Praze PCP, CS, 3LF UK v Praze Abstrat Well nown logisti regression and the other binary response models an be used
More informationRIEMANN S FIRST PROOF OF THE ANALYTIC CONTINUATION OF ζ(s) AND L(s, χ)
RIEMANN S FIRST PROOF OF THE ANALYTIC CONTINUATION OF ζ(s AND L(s, χ FELIX RUBIN SEMINAR ON MODULAR FORMS, WINTER TERM 6 Abstrat. In this hapter, we will see a proof of the analyti ontinuation of the Riemann
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationREVISION SHEET FP2 (Edx) CALCULUS. x x 0.5. x x 1.5. π π. Standard Calculus of Inverse Trig and Hyperbolic Trig Functions = + = + arcsin x = +
the Further Mathematis netwk www.fmnetwk.g.uk V 07 REVISION SHEET FP (Ed) CLCULUS The main ideas are: Calulus using inverse trig funtions & hperboli trig funtions and their inverses. Malaurin series Differentiating
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationA Numerical Method For Constructing Geo-Location Isograms
A Numerial Method For Construting Geo-Loation Isograms Mike Grabbe The Johns Hopkins University Applied Physis Laboratory Laurel, MD Memo Number GVW--U- June 9, 2 Introdution Geo-loation is often performed
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationMATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II
Journal of Applied Mathematis and Computational Mehanis 2014, 13(2), 141-147 MATHEMATICA AND NUMERICA BAI OF BINARY AOY OIDIFICATION MODE WITH UBTITUTE THERMA CAPACITY. PART II Ewa Węgrzyn-krzypzak 1,
More informationDifferential Equations 8/24/2010
Differential Equations A Differential i Equation (DE) is an equation ontaining one or more derivatives of an unknown dependant d variable with respet to (wrt) one or more independent variables. Solution
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationExplicit formulae for integrals of s and p type GTFs
ELECTRONIC JOURNAL OF THEORETICAL CHEMISTRY, VOL., 66 70 (997) Expliit formulae for integrals of s and p type GTFs CAROL A. BAXTER AND DAVID B. COOK The Department of Chemistry, The University of Sheffield,
More informationPlanck unit theory: Fine structure constant alpha and sqrt of Planck momentum
Plank unit theory: Fine struture onstant alpha and sqrt of Plank momentum Malolm Maleod e-mail: mail4malolm@gmx.de The primary onstants; G,, h, e, α, k B, m e... range in preision from low G (4-digits)
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationA Functional Representation of Fuzzy Preferences
Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationOn Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations
On Certain Singular Integral Equations Arising in the Analysis of Wellbore Reharge in Anisotropi Formations C. Atkinson a, E. Sarris b, E. Gravanis b, P. Papanastasiou a Department of Mathematis, Imperial
More informationA Combination of BIM and BEM for efficiently analyzing optical. elements
A Combination of BIM and BEM for effiiently analyzing optial elements SUN Fang, LIU Juan, ZHANG Guo-ting, HU Chuan-fei, SU Xiaoxing Institute of Lightwave Tehnology, hool of eletronis and information engineering,
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationKAMILLA OLIVER AND HELMUT PRODINGER
THE CONTINUED RACTION EXPANSION O GAUSS HYPERGEOMETRIC UNCTION AND A NEW APPLICATION TO THE TANGENT UNCTION KAMILLA OLIVER AND HELMUT PRODINGER Abstrat Starting from a formula for tan(nx in the elebrated
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationPerturbation Analyses for the Cholesky Factorization with Backward Rounding Errors
Perturbation Analyses for the holesky Fatorization with Bakward Rounding Errors Xiao-Wen hang Shool of omputer Siene, MGill University, Montreal, Quebe, anada, H3A A7 Abstrat. This paper gives perturbation
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.286: The Early Universe December 21, 2013 Prof. Alan Guth QUIZ 3 SOLUTIONS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.286: The Early Universe Deember 2, 203 Prof. Alan Guth QUIZ 3 SOLUTIONS Quiz Date: Deember 5, 203 PROBLEM : DID YOU DO THE READING? (35
More informationRelativistic Addition of Velocities *
OpenStax-CNX module: m42540 1 Relativisti Addition of Veloities * OpenStax This work is produed by OpenStax-CNX and liensed under the Creative Commons Attribution Liense 3.0 Abstrat Calulate relativisti
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More information5.1 Composite Functions
SECTION. Composite Funtions 7. Composite Funtions PREPARING FOR THIS SECTION Before getting started, review the following: Find the Value of a Funtion (Setion., pp. 9 ) Domain of a Funtion (Setion., pp.
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More informationBrazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle
Brazilian Journal of Physis, vol. 9, no. 3, September, 1999 51 Classial and Quantum Mehanis of a Charged Partile in Osillating Eletri and Magneti Fields V.L.B. de Jesus, A.P. Guimar~aes, and I.S. Oliveira
More informationThe Second Postulate of Euclid and the Hyperbolic Geometry
1 The Seond Postulate of Eulid and the Hyperboli Geometry Yuriy N. Zayko Department of Applied Informatis, Faulty of Publi Administration, Russian Presidential Aademy of National Eonomy and Publi Administration,
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationSome facts you should know that would be convenient when evaluating a limit:
Some fats you should know that would be onvenient when evaluating a it: When evaluating a it of fration of two funtions, f(x) x a g(x) If f and g are both ontinuous inside an open interval that ontains
More informationThe universal model of error of active power measuring channel
7 th Symposium EKO TC 4 3 rd Symposium EKO TC 9 and 5 th WADC Workshop nstrumentation for the CT Era Sept. 8-2 Kosie Slovakia The universal model of error of ative power measuring hannel Boris Stogny Evgeny
More informationInternational Journal of Electronics and Computer Science Engineering 817. Available Online at ISSN
International Journal of Eletronis and Computer Siene Engineering 817 Available Online at www.ijese.org ISSN- 2277-1956 A Duly Synhronized, Straightforward Approah For Realizing the General Charateristis
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationA NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL
Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS
More informationTheory. Coupled Rooms
Theory of Coupled Rooms For: nternal only Report No.: R/50/TCR Prepared by:. N. taey B.., MO Otober 00 .00 Objet.. The objet of this doument is present the theory alulations to estimate the reverberant
More informationThe Electromagnetic Radiation and Gravity
International Journal of Theoretial and Mathematial Physis 016, 6(3): 93-98 DOI: 10.593/j.ijtmp.0160603.01 The Eletromagneti Radiation and Gravity Bratianu Daniel Str. Teiului Nr. 16, Ploiesti, Romania
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationAcoustic Waves in a Duct
Aousti Waves in a Dut 1 One-Dimensional Waves The one-dimensional wave approximation is valid when the wavelength λ is muh larger than the diameter of the dut D, λ D. The aousti pressure disturbane p is
More informationCritical Reflections on the Hafele and Keating Experiment
Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As
More information(q) -convergence. Comenius University, Bratislava, Slovakia
Annales Mathematiae et Informatiae 38 (2011) pp. 27 36 http://ami.ektf.hu On I (q) -onvergene J. Gogola a, M. Mačaj b, T. Visnyai b a University of Eonomis, Bratislava, Slovakia e-mail: gogola@euba.sk
More informationMATHEMATICS MHT-CET TRIUMPH. Salient Features MULTIPLE CHOICE
Written in aordane with the latest MHT-CET Paper Pattern whih inludes topis based on Std. XII S. and relevant hapters of Std. XI S. (Maharashtra State Board) MHT-CET TRIUMPH MATHEMATICS MULTIPLE CHOICE
More informationA Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil
IOSR Journal of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 PP 11-15 www.iosrjen.org A Heuristi Approah for Design and Calulation of Pressure Distribution over Naa 4 Digit Airfoil G.
More informationThe transition between quasi-static and fully dynamic for interfaces
Physia D 198 (24) 136 147 The transition between quasi-stati and fully dynami for interfaes G. Caginalp, H. Merdan Department of Mathematis, University of Pittsburgh, Pittsburgh, PA 1526, USA Reeived 6
More information10.2 The Occurrence of Critical Flow; Controls
10. The Ourrene of Critial Flow; Controls In addition to the type of problem in whih both q and E are initially presribed; there is a problem whih is of pratial interest: Given a value of q, what fators
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationReview of Force, Stress, and Strain Tensors
Review of Fore, Stress, and Strain Tensors.1 The Fore Vetor Fores an be grouped into two broad ategories: surfae fores and body fores. Surfae fores are those that at over a surfae (as the name implies),
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More informationStrauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u
Strauss PDEs e: Setion 3.4 - Exerise 3 Page 1 of 13 Exerise 3 Solve u tt = u xx + os x, u(x, ) = sin x, u t (x, ) = 1 + x. Solution Solution by Operator Fatorization Bring u xx to the other side. Write
More informationDevelopment of the Numerical Schemes and Iteration Procedures Nielsen, Peter Vilhelm
Aalborg Universitet Development of the Numerial Shemes and Iteration Proedures Nielsen, Peter Vilhelm Published in: Euroaademy on Ventilation and Indoor Climate Publiation date: 2008 Doument Version Publisher's
More informationChapter 2: One-dimensional Steady State Conduction
1 Chapter : One-imensional Steay State Conution.1 Eamples of One-imensional Conution Eample.1: Plate with Energy Generation an Variable Conutivity Sine k is variable it must remain insie the ifferentiation
More informationOn maximal inequalities via comparison principle
Makasu Journal of Inequalities and Appliations (2015 2015:348 DOI 10.1186/s13660-015-0873-3 R E S E A R C H Open Aess On maximal inequalities via omparison priniple Cloud Makasu * * Correspondene: makasu@uw.a.za
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationMAC Calculus II Summer All you need to know on partial fractions and more
MC -75-Calulus II Summer 00 ll you need to know on partial frations and more What are partial frations? following forms:.... where, α are onstants. Partial frations are frations of one of the + α, ( +
More informationBuckling loads of columns of regular polygon cross-section with constant volume and clamped ends
76 Bukling loads of olumns of regular polygon ross-setion with onstant volume and lamped ends Byoung Koo Lee Dept. of Civil Engineering, Wonkwang University, Iksan, Junuk, 7-79, Korea Email: kleest@wonkwang.a.kr
More informationThe tanh - coth Method for Soliton and Exact Solutions of the Sawada - Kotera Equation
Volume 117 No. 13 2017, 19-27 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu The tanh - oth Method for Soliton and Exat Solutions of the Sawada -
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More informationAlgorithm 736: Hyperelliptic Integrals and the Surface Measure of Ellipsoids
Algorithm 736: Hyperellipti Integrals and the Surfae Measure of Ellipsoids CHARLES F. DUNKL and DONALD E. RAMIREZ University of Virginia The algorithm for omputing a lass of hyperellipti integrals and
More informationNAVAL UNDERWATER SYSTEMS CENTER NEW LONDON LABORATORY NEW LONDON, CONNECTICUT Technical Memorandum
NAVAL UNDERWATER SYSTEMS CENTER NEW LONDON LABORATORY NEW LONDON, CONNECTICUT 06320 Tehnial Memorandum A FORTRAN COMPUTER PROGRAM TO EVALUATE AND PLOT THE STATISTICS OF THE MAGNITUDE-SQUARED COHERENCE
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More information6.4 Dividing Polynomials: Long Division and Synthetic Division
6 CHAPTER 6 Rational Epressions 6. Whih of the following are equivalent to? y a., b. # y. y, y 6. Whih of the following are equivalent to 5? a a. 5, b. a 5, 5. # a a 6. In your own words, eplain one method
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More informationJoint simultaneous inversion of PP and PS angle gathers
Inversion of and angle gathers Joint simultaneous inversion of and angle gathers Brian H. ussell, aniel. Hampson, Keith Hirshe, and Janusz eron ABTACT e present a new approah to the joint simultaneous
More informationThe gravitational phenomena without the curved spacetime
The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,
More informationClassical Trajectories in Rindler Space and Restricted Structure of Phase Space with PT-Symmetric Hamiltonian. Abstract
Classial Trajetories in Rindler Spae and Restrited Struture of Phase Spae with PT-Symmetri Hamiltonian Soma Mitra 1 and Somenath Chakrabarty 2 Department of Physis, Visva-Bharati, Santiniketan 731 235,
More informationGeneral Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines
General Closed-form Analytial Expressions of Air-gap Indutanes for Surfaemounted Permanent Magnet and Indution Mahines Ronghai Qu, Member, IEEE Eletroni & Photoni Systems Tehnologies General Eletri Company
More informationQuasi-Monte Carlo Algorithms for unbounded, weighted integration problems
Quasi-Monte Carlo Algorithms for unbounded, weighted integration problems Jürgen Hartinger Reinhold F. Kainhofer Robert F. Tihy Department of Mathematis, Graz University of Tehnology, Steyrergasse 30,
More informationPseudo Spheres. A Sample of Electronic Lecture Notes in Mathematics. Eberhard Malkowsky.
Pseudo Spheres A Sample of Eletroni Leture Notes in Mathematis Eberhard Malkowsky Mathematishes Institut Justus Liebig Universität Gießen Arndtstraße D-3539 Gießen Germany /o Shool of Informatis Computing
More informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationCALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS
International Journal of Modern Physis A Vol. 24, No. 5 (2009) 974 986 World Sientifi Publishing Company CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS PAVEL SNOPOK, MARTIN
More informationLong time stability of regularized PML wave equations
Long time stability of regularized PML wave equations Dojin Kim Email:kimdojin@knu.a.kr Yonghyeon Jeon Email:dydgus@knu.a.kr Philsu Kim Email:kimps@knu.a.kr Abstrat In this paper, we onsider two dimensional
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More informationThe Impact of Information on the Performance of an M/M/1 Queueing System
The Impat of Information on the Performane of an M/M/1 Queueing System by Mojgan Nasiri A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master
More informationCasimir self-energy of a free electron
Casimir self-energy of a free eletron Allan Rosenwaig* Arist Instruments, In. Fremont, CA 94538 Abstrat We derive the eletromagneti self-energy and the radiative orretion to the gyromagneti ratio of a
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs e: Setion.1 - Exerise 1 Page 1 of 6 Exerise 1 Solve u tt = u xx, u(x, 0) = e x, u t (x, 0) = sin x. Solution Solution by Operator Fatorization By fatoring the wave equation and making a substitution,
More informationChapter 2 Linear Elastic Fracture Mechanics
Chapter 2 Linear Elasti Frature Mehanis 2.1 Introdution Beginning with the fabriation of stone-age axes, instint and experiene about the strength of various materials (as well as appearane, ost, availability
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationSampler-A. Secondary Mathematics Assessment. Sampler 521-A
Sampler-A Seondary Mathematis Assessment Sampler 521-A Instrutions for Students Desription This sample test inludes 14 Seleted Response and 4 Construted Response questions. Eah Seleted Response has a
More informationA Longitudinal Aerodynamic Data Repeatability Study for a Commercial Transport Model Test in the National Transonic Facility
NASA Tehnial Paper 3522 A Longitudinal Aerodynami Data Repeatability Study for a Commerial Transport Model Test in the National Transoni Faility R. A. Wahls and J. B. Adok Langley Researh Center Hampton,
More information