Novel Approach to Analysis of Nonlinear Recursions. 1 Department of Physics, Bar-Ilan University, Ramat-Gan, ISRAEL
|
|
- Hannah Page
- 5 years ago
- Views:
Transcription
1 Novel Approach to Analysis of Nonlinear Recursions G.Berkolaiko 1 2, S. Rabinovich 1,S.Havlin 1 1 Department of Physics, Bar-Ilan University, 529 Ramat-Gan, ISRAEL 2 Department of Mathematics, Voronezh University, Voronezh, RUSSIA (April 26, 1996) We propose a general method to map nonlinear recursions on a linear (but matrix) one. The solution of the recursion is represented as a product of matrices whose elements depend only on the form of the recursion and not on the initial conditions. First we consider the method for polynomial recursions of arbitrary order. Then we discuss the generalizations of the method for systems of polynomial recursions and arbitrary analytic recursions as well. The only restriction for these equations is to be solvable with respect to the highest order term (existence of a normal form). I. INTRODUCTION Recursions take a central place in various elds of science. Numerical solution of dierential equations and dierent models of evolution of a system involve, in general, recursions. By now, only linear recursions with constant coecients could be solved [1{3]. However, the nonlinearity changes the situation drastically. Solution of a rather simple recursion, the logistic map, y n+1 = y n (1 ; y n ) is far from being so simple as one might guess. The analysis of its behavior, while based on a roundabout approach, has revealed many amazing properties. In this paper we suggest a new approach to the solution of nonlinear recursions. It turns out, that the coecients of the i-th iteration of the polynomial depend linearly on the coecients of the (i ; 1)-th iteration. Using this fact we succeed to write down the general 1
2 solution of the polynomial recursion and to generalize it for any recursion y n+1 = f(y n ) with y y where f(x) is an analytic function. In the paper we study the conditions that the analytic function must satisfy to make this generalization possible. The gist of the method is the construction of a special transfer matrix, T. It allows to represent the solution in the form y n = het n yi (1) where he is the rst unit vector, yi is a vector of initial values and T n is the n-th power of the matrix T. II. FIRST-ORDER POLYNOMIAL RECURSION Here we consider a rst-order recursion equation in its normal form y n+1 = P (y n ) (2) where P (x) is a polynomial of degree m: P (x) = mx a k x k a m 6=: (3) Let y y be an initial value for the recursion (2). We denote by yi the column vector of powers of y yi = fy g 1 and the vector he is a row vector he =[ = 1] 1 : It should be = emphasized, that runs from, since in the general case a 6=. In this notation heyi is a scalar product that yields heyi = y: (4) Theorem 1. For any recursion of type Eq.(2) there exists a matrix T = ft k g 1 such that 2
3 y n = het n yi: (5) The matrix power T n exists for all n and all the operations in the right-hand side of Eq. (5) are associative. Proof. For n = the statement of the theorem is valid (see Eq. (4)). We introduce the n column vector y 1 i def = y 1 o 1 = where y 1 = P (y). Let T be a matrix such that y 1 i = Tyi : (6) The existence of this matrix will be proven later on. If such a matrix exists, then, analogically to Eq. (4), we have y 1 = hey 1 i = hetyi. Therefore, the statement of the theorem is true for n =1aswell. Assume, that Eq. (5) is valid for n = l for any initial value y. Theny l+1 can be represented as y l+1 = het l y 1 i, where y 1 = P (y) is considered as a new initial value of the recursion. Then, using Eq. (6) one gets y l+1 = het l y 1 i = hett l yi = het l+1 yi : Now we prove the existence of the matrix T. Onehashy 1 def =[P (y)] 1 =. In turn P (y) is the m-th degree polynomial P (y) = mx a i y i = mx T k y k (7) and we infer that T = ft k g 1 obeys Eq. (6). Note, that for and k satisfying k m we have T k, therefore, each row is nite (i.e., there is only a nite number of nonzero matrix elements in each row). This proves the existence of powers of T and associativity in Eq. (5). Thus, the proof is complete. III. SPECIAL CASES In some special cases the matrix T has the rather simple form. 3
4 (a) The binomial: P (x) =a p x p + a q x q. As one can see, in the general case elements of the matrix T have a form of rather complicated sums. However, they are degenerated to a fairly simple expression, when the polynomial (3) has only two terms. In this case one has P (y) =(a p y p + a q y q ) = X i a ;i p a i q yp(;i)+qi : Denoting k = p( ; i)+qi i = l(k) =(q ; p) ;1 (k ; p) we have P (y) = qx k=p y k l(k) a ;l(k) p a l(k) q : Thus, the matrix elements T k are T k = a ;l(k) p a l(k) q : (8) l(k) Example 1. To demonstrate this let us consider the following recursion: y n+1 = y n (1 ; y n ) with y y: (9) Plugging p =1,q =2,a p = ;a q = into Eq. (8) one immediately obtains T k =(;1) k; k ; : (1) Thus, we recover the result of Rabinovich et al [4] for the recursion equation known as the logistic mapping. form. (b) The trinomial P (x) =a + a p x p + a q x q and a 6=. Then the matrix T has a special Lemma 1. Let P (x) =a + a p x p + a q x q. Then, the following decomposition exists T = AT where T is a matrix corresponding to the polynomial P (x) = a p x p + a q x q and A is a triangular matrix. 4
5 Proof. Consider P (x) =a p x p + a q x q and the corresponding matrix T. It yields n o T yi = y i def 1 = P 1 (y) = : For the matrix T one gets Tyi = y 1 i def = fp (y) g 1 = P (y) = X i a ;i (a p y p + a q y q ) i : Denoting in the last line A i a ;i i one obtains y 1 i = Ay1i = AT yi = Tyi and T = AT. The lemma isproven. IV. POLYNOMIAL RECURSIONS WITH NON-CONSTANT COEFFICIENTS As shown in [4] a generalization of Eq. (9), y n+1 = n y n (1 ; y n ) with y y (11) can be solved using a similar approach. The solution is y n = het n T 2 T 1 yi (12) where the matrix elements of T i are now i-dependent: (T i ) k =(;1) k; ( i ) : (13) k ; The same argument isvalid for the general case. Theorem 2. Let the polynomial in the Eq. (2) depends on n. Then the solution Eq. (5) takes the form of Eq. (12) with the obvious changes (a ::: a m become i-dependent functions) in corresponding matrix elements. To illustrate this Theorem we consider Example 2. We are going to apply the previous material to the Riccati equation. Usually, people use this name for the equation y n+1 y n + a n y n+1 + b n y n + c n =: 5
6 However, by appropriate variable change [1,2] this equation can be reduced to a linear one and then treated by conventional techniques. Here we shall be dealing with the following recursion: y n+1 = a n + b n y n + c n y 2 n with y y: This is another possible (asymmetric) discrete analog of the Riccati dierential equation [5{7]. It is well-known that the latter cannot be solved in quadratures. The general results of the two previous sections can be employed to write down the solution of this recursion. Namely, the solution reads y n = het n T 2 T 1 yi where the matrix T i is a product of two matrices T i = A i S i : The matrix elements are given by (A i ) k = k a ;k i and (S i ) k = k ; b 2;k i c k; i : V. SYSTEM OF FIRST-ORDER NONLINEAR EQUATIONS Actually, very little is known about systems of nonlinear equations [8]. We now extend our method of Sect. II to handle systems of nonlinear equations. Let us demonstrate it on the following example: 8 >< >: u n+1 = u n (1 ; v n ) with u u v n+1 = v n (1 ; u n ) with v v: (14) Proceeding here as in Sect. II, we are checking the transformation of a product u v k : X [u(1 ; v)] [v(1 ; u)] k = u (;1) r v r k v k (;1) s k u s r s r s X = u p v q k (;1) (p;)+(q;k) k : (15) p q q ; k p ; 6
7 One can proceed with the aid of multidimensional matrices [11], but here we prefer to use more traditional two-dimensional matrices. Introduce now avector uvi, which consists of powers u v k of variables u and v, arranged in some order. Namely,introduce a biection ( ): N 2 N, wheren is the set of the natural numbers and zero. Then, the monomial u v k is the ( k)-th component of the vector uvi. For example, one can use ( k) =k + 1 ( + k)( + k +1): 2 Introducing a matrix T with the elements T ( k)(p q) =(;1) (p;)+(q;k) k q ; k p ; k we basically return to the familiar transfer-matrix construction. Indeed, we have u n = he (1 ) T n uvi v n = he ( 1) T n uvi where he ( k) = f ( k) i g 1 : Theorem 3. Consider a system of m rst-order nonlinear equations: x (i) n+1 = P i (x (1) ::: x(m) ) i =1 ::: m: (16) n n Then the matrix T exists such that x (i) n = he (i) T n x (1) x (m) i where the vectors he (i) and x (1) x (m) i are dened asabove, with the aid of some biection ( ::: ): N m N. Below, we present the scheme of the proof of a transfer-matrix representation for arbitrary systems of rst-order nonlinear equations. 7
8 As before we are checking the product P 1 1 P m m P (1 ::: m): The polynomial P (1 ::: m) depends on m variables x (1) ::: x (m) and therefore can be represented as P (1 ::: m)(x (1) ::: x (m) )=ht (1 ::: m)x (1) x (m) i where ht (1 ::: m) is the constant vector of coecients of the polynomial P (1 ::: m). This vector ht (1 ::: m) is the ( 1 ::: m )-th row of the transfer matrix T. The above procedure can be easily generalized on polynomials with variable coecients, P i (n x (1) ::: x (m) ), in complete analog to the Theorem 2. VI. RECURSION SOLUTION FOR ARBITRARY ANALYTIC FUNCTIONS IN THE RHS Let the series converge absolutely in [;r r]. Then the series a k x k f(x) (17) ~a k x k ~ f(x) (18) where ~a k = a k, also converge in [;r r]. We construct the following absolutely convergent in [;r r] [9,1] series " 1 X a k x k# = a + 1 " X 1 ~a k x k# =~a + 1 a 1 a ;1 x + ~a 1 ~a ;1 x + " 1 " 1 a 2 a ;1 + 2 ~a 2 ~a ;1 + 2 a 2 1 a;2 ~a 2 1 ~a;2 # # x 2 + :::= x 2 + :::= T k x k (19) ~T k x k : (2) Denition 1. The matrix T = ft k g 1 is said to be the transfer matrix of the analytic function f(x) if P 1 T kx k = f (x). 8
9 Using Eq. (19) one can estimate the coecients T k for a xed k: T k + k k a ;k max k a i C(a + ) (21) ik for arbitrary >, a constant C = C(), and >k. Theorem 4. Let the series b k x k g(x) (22) converge absolutely in [;R R]. The matrix S be the transfer matrix of the function g(x), and ~r r be such that the function ~ f(x), dened by(18), satises ~ f(x): [;~r ~r] [ R]. Then the matrix product ST exists and is equal to the transfer matrix of the composition of the functions f(x) and g(x), g f(x), i.e., the -th row of the matrix ST is the vector of the coecients of Maclaurin's decomposition of the function [g f] (x). Proof. First we prove the existence of the matrix ST. One can estimate the ( k)-th element of the matrix ST using Eq. (21): (ST) k S i T ik C 1 + S i T ik C 1 + C i=k+1 i=k+1 S i (a + ) i : Since min f(x) ~ = f() ~ = a <Rwe canchoose > suciently small, that a + <R, therefore the series in the last expression converges. Thus, the matrix ST exists. Now we should study the series Since the series converge absolutely: S i T ik x k (ST) k x k = S i ~T ik x k we can change the order of summation [1]: (ST) k x k = S i T ik x k = S i T ik x k : (23) h i S i f(x) ~ i < 1: (24) S i [f(x)] i =[g(f(x))] =[g f] (x): The theorem is proved. Restricting ourselves to the rst row of the matrix S only we get the following 9
10 Corollary 1. Let the functions f(x) and g(x) satisfy the conditions of Theorem 4 the bra-vector hb is the vector of coecients of the series (22). Then the vector hc = hbt consists of the coecients of the series c k x k = g f(x): (25) It is easy to see, that calculating the scalar product of the vector hc and the vector xi = fx i g 1 we obtain hbtxi = g f(x): (26) In this way we can write down the solution of the recursion y n+1 = f(y n ) with y y: (27) Theorem 5. Let the series (17) converge absolutely on [;r r], the matrix T be the transfer matrix of the function f(x) and ~r <rbe such that (see Eq. (18)) ~ f(x): [;~r ~r] [ ~r]. Then for any initial conditions y 2 [;~r ~r] y n = het n yi (28) where the bra-vector he =[ i 1 ] 1 and the ket-vector yi = fyi g 1. In order to solve a generalization of the recursion (27), y n+1 = f(n y n ), we should consider n-dependent matrices T n instead of the matrix T. One can use the same technique to solve the more general case of multivariable recursions, y n = f(y n;1 ) f: R n R n f(x) =ff i (x)g n i=1 x 2 Rn Here we should consider series decomposition of the functions Q n i=1 [f i (x)] a i, instead of the series (19), and then proceed by the analogy to multivariable polynomial recursions. However, since the notation becomes very cumbersome, we shall not write down this formalism. The theorems of this section can be easily generalized on complex mappings. Namely, the statement of the Theorem 4 for the complex variable functions reads: 1
11 Theorem 6. Let the series (17), (18) and (22) be the series of complex variable z with the complex coecients. Series (17) and (18) converge in the circle z r and series (22) converge in z R. The matrices S and T be the transfer matrices of the functions g(z) and f(z) correspondingly and ~r r ~r 2 R be such that f(~r) ~ R. Then the matrix product ST exists and equals to the transfer matrix of the composition of the functions f(z) and g(z), g f(z), i.e., the -th row of the matrix ST is the vector of the coecients of Maclaurin's decomposition of the function [g f] (z) convergent in z ~r. VII. EXAMPLES In order to show that the condition ~ f(x): [;~r ~r] [ R] is essential (see Theorem 4) we present the following Example 3. Let g(x) =x 4 ; 1 4 x x8 ; 1 16 x1 + ::: i.e. a 2i = (;1)i (i;1) 2 a 2i+1 = i 2 and radius of convergence is R =1 letf(x) bedened by f(x) =2x(1 ; x 2 ) f:[;1 1] [;1 1]. The transfer matrix T of the function f(x) according to Eq. (8) is T k =(;1) (k;)=2 (k;)=2 2 for even ; k. The sign in the sequences fs 1 2i g 1 i=1 and ft 2i 2kg 1 i=1 alternates. Therefore the sign in the sequence fs 1 2iT 2i 2k g 1 i=1 is constant and for the 2k-th component of the rst row of the matrix ST one has c 1 2k = S 1 2i T 2i 2k = =4 k;2 1 (k ; 3) 2 (2k ; 4)(2k ; 5) 2 S 1 2i T 2i 2k S 1 2(k;2) T 2(k;2) 2k 4 k;2 for any k > 3. Thus, the components of the the rst row cannot be the coecients of Maclaurin's decomposition of the function g f because of its divergence for x > 1=4. The following example presents a family of function whose transfer matrices form is invariant under the multiplication. 11
12 Example 4. Let the function f(x) = ax=(b + x). Then the components of the transfer matrix T are T k = k ; 1 ; 1 a b ;k for k > and T k = k for =ork = (29) where n m = 8 >< >: n m(n;m) integer m n m n any other m n and nm is the Kronecker symbol. If g(x) =cx=(d + x) ands is the transfer matrix of g(x), then the components of the matrix ST are (ST) k = = k ; 1 ; 1 S i T ik = ac a + b i ; 1 c d ;i k ; 1 ; 1 i ; 1 ;k bd : a + b a i b ;k = k ; 1 k; X c b ;k (a=d) l+ ; 1 l= k ; l One can see, that the matrix elements (ST) k have the same form as in Eq. (29). VIII. SUMMARY In this paper wehave presented a new method to obtain the solution of arbitrary polynomial recursions. The method has been generalized to the systems of multivariable recursions and analytical recursions. We found a class of analytical function recursions to which the method can be applied. Particularly, this class contains functions which are analytical over all space. Generally, the solution is obtained in the form of a matrix power, applied to the vectors of initial values. We have presented a way to construct such a matrix. Famous and important examples, such asthelogistic map and the Riccati recursion, have been considered and the corresponding matrices have been written down explicitly. The following generalizations are also can be done: a) Multivariable analytic recursions. This can be done, for example, as in the Section V. 12
13 b) System of higher-order nonlinear equations. The scheme is quite obvious: introduction of new variables to bring each equation to the rst-order structure and, then, construction of a transfer matrix, as in the previous case. IX. ACKNOWLEDGMENT We thank Dany Ben-Avraham for critical reading. [1] R. P. Agarwal: Dierence Equations and Inequalities, Marcel Dekker, New York, 1992 [2] V. Lakshmikantham, D. Trigiante: Theory of Dierence Equations: Numerical Methods and Applications, Academic Press, New York, 1988 [3] F. B. Hildebrand: Finite-Dierence Equations and Simulations, Prentice-Hall, Englewood Clits, NJ, 1968 [4] S. Rabinovich, G. Berkolaiko, S. Buldyrev, A. Shehter and S. Havlin: Physica A218, 457 (1995) [5] D. Zwillinger: Handbook of Dierential Equations, 2nd ed., Academic Press, Boston, 1992 [6] W. T. Reid: Riccati Dierential Equation, Academic Press, New York, 1972 [7] S. Bittanti, A. J. Laub, J. C. Willems, eds: The Riccati Equation, Springer, Berlin, 1991 [8] V. L. Kocic and G. Ladas: Global Behavior of Nonlinear Dierence Equations of Higher Order with Applications, Kluwer, Dordrecht, 1993 [9] R. Courant: Dierential and Integral Calculus, Blackie & Son, London, v.1, 1963 [1] E. Goursat: A Course in Mathematical Analysis, Dover, New York, v.1, 1956 [11] S. Rabinovich, G. Berkolaiko, S. Havlin: preprint, to be published 13
Analysis of Carleman Representation of Analytical Recursions
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 224, 8190 1998 ARTICLE NO. AY985986 Analysis of Carleman Representation of Analytical Recursions G. Berolaio Department of Physics, Bar-Ilan Uniersity,
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationVector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)
Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational
More informationLinear Algebra, 4th day, Thursday 7/1/04 REU Info:
Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector
More informationGeneralization of Hensel lemma: nding of roots of p-adic Lipschitz functions
Generalization of Hensel lemma: nding of roots of p-adic Lipschitz functions (joint talk with Andrei Khrennikov) Dr. Ekaterina Yurova Axelsson Linnaeus University, Sweden September 8, 2015 Outline Denitions
More informationNumerical Integration exact integration is not needed to achieve the optimal convergence rate of nite element solutions ([, 9, 11], and Chapter 7). In
Chapter 6 Numerical Integration 6.1 Introduction After transformation to a canonical element,typical integrals in the element stiness or mass matrices (cf. (5.5.8)) have the forms Q = T ( )N s Nt det(j
More informationInstitute for Advanced Computer Studies. Department of Computer Science. On the Convergence of. Multipoint Iterations. G. W. Stewart y.
University of Maryland Institute for Advanced Computer Studies Department of Computer Science College Park TR{93{10 TR{3030 On the Convergence of Multipoint Iterations G. W. Stewart y February, 1993 Reviseed,
More informationContents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces
Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and
More information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationHomework 7 Solutions to Selected Problems
Homework 7 Solutions to Selected Prolems May 9, 01 1 Chapter 16, Prolem 17 Let D e an integral domain and f(x) = a n x n +... + a 0 and g(x) = m x m +... + 0 e polynomials with coecients in D, where a
More information1. Introduction Let the least value of an objective function F (x), x2r n, be required, where F (x) can be calculated for any vector of variables x2r
DAMTP 2002/NA08 Least Frobenius norm updating of quadratic models that satisfy interpolation conditions 1 M.J.D. Powell Abstract: Quadratic models of objective functions are highly useful in many optimization
More informationOn the classication of algebras
Technische Universität Carolo-Wilhelmina Braunschweig Institut Computational Mathematics On the classication of algebras Morten Wesche September 19, 2016 Introduction Higman (1950) published the papers
More informationOctober 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable
International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität
More informationARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (1996), 13 { 27. ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION
ARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (996), 3 { 27 ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION P. J. Y. Wong and R. P. Agarwal Abstract. We oer sucient conditions for the
More informationNumerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
Available online at http://ijim.srbiau.ac.ir Int. J. Industrial Mathematics Vol. 1, No. 2 (2009)147-161 Numerical Solution of Hybrid Fuzzy Dierential Equation (IVP) by Improved Predictor-Corrector Method
More information2 ALGEBRA II. Contents
ALGEBRA II 1 2 ALGEBRA II Contents 1. Results from elementary number theory 3 2. Groups 4 2.1. Denition, Subgroup, Order of an element 4 2.2. Equivalence relation, Lagrange's theorem, Cyclic group 9 2.3.
More informationEXAMPLES OF PROOFS BY INDUCTION
EXAMPLES OF PROOFS BY INDUCTION KEITH CONRAD 1. Introduction In this handout we illustrate proofs by induction from several areas of mathematics: linear algebra, polynomial algebra, and calculus. Becoming
More informationFinding Low Degree Annihilators for a Boolean Function Using Polynomial Algorithms
Finding Low Degree Annihilators for a Boolean Function Using Polynomial Algorithms Vladimir Bayev Abstract. Low degree annihilators for Boolean functions are of great interest in cryptology because of
More informationFraction-free Row Reduction of Matrices of Skew Polynomials
Fraction-free Row Reduction of Matrices of Skew Polynomials Bernhard Beckermann Laboratoire d Analyse Numérique et d Optimisation Université des Sciences et Technologies de Lille France bbecker@ano.univ-lille1.fr
More informationTHE FUNDAMENTAL THEOREM OF ALGEBRA VIA PROPER MAPS
THE FUNDAMENTAL THEOREM OF ALGEBRA VIA PROPER MAPS KEITH CONRAD 1. Introduction The Fundamental Theorem of Algebra says every nonconstant polynomial with complex coefficients can be factored into linear
More informationAutomatica, 33(9): , September 1997.
A Parallel Algorithm for Principal nth Roots of Matrices C. K. Koc and M. _ Inceoglu Abstract An iterative algorithm for computing the principal nth root of a positive denite matrix is presented. The algorithm
More informationThe Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov
More informationRough Sets, Rough Relations and Rough Functions. Zdzislaw Pawlak. Warsaw University of Technology. ul. Nowowiejska 15/19, Warsaw, Poland.
Rough Sets, Rough Relations and Rough Functions Zdzislaw Pawlak Institute of Computer Science Warsaw University of Technology ul. Nowowiejska 15/19, 00 665 Warsaw, Poland and Institute of Theoretical and
More informationA BASIC FAMILY OF ITERATION FUNCTIONS FOR POLYNOMIAL ROOT FINDING AND ITS CHARACTERIZATIONS. Bahman Kalantari 3. Department of Computer Science
A BASIC FAMILY OF ITERATION FUNCTIONS FOR POLYNOMIAL ROOT FINDING AND ITS CHARACTERIZATIONS Bahman Kalantari 3 Department of Computer Science Rutgers University, New Brunswick, NJ 08903 Iraj Kalantari
More informationProperties of Arithmetical Functions
Properties of Arithmetical Functions Zack Clark Math 336, Spring 206 Introduction Arithmetical functions, as dened by Delany [2], are the functions f(n) that take positive integers n to complex numbers.
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationA note on continuous behavior homomorphisms
Available online at www.sciencedirect.com Systems & Control Letters 49 (2003) 359 363 www.elsevier.com/locate/sysconle A note on continuous behavior homomorphisms P.A. Fuhrmann 1 Department of Mathematics,
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More information7. Baker-Campbell-Hausdorff formula
7. Baker-Campbell-Hausdorff formula 7.1. Formulation. Let G GL(n,R) be a matrix Lie group and let g = Lie(G). The exponential map is an analytic diffeomorphim of a neighborhood of 0 in g with a neighborhood
More informationAbstract. 1 Introduction. 2 Idempotent filters. Idempotent Nonlinear Filters
Idempotent Nonlinear Filters R.K. Pearson and M. Gabbouj Tampere International Center for Signal Processing Tampere University of Technology P.O. Box 553, FIN-33101 Tampere, Finland E-mail: {pearson,gabbouj}@cs.tut.fi
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
More informationVector Spaces ปร ภ ม เวกเตอร
Vector Spaces ปร ภ ม เวกเตอร 5.1 Real Vector Spaces ปร ภ ม เวกเตอร ของจ านวนจร ง Vector Space Axioms (1/2) Let V be an arbitrary nonempty set of objects on which two operations are defined, addition and
More information4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial
Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationTaylor series based nite dierence approximations of higher-degree derivatives
Journal of Computational and Applied Mathematics 54 (3) 5 4 www.elsevier.com/locate/cam Taylor series based nite dierence approximations of higher-degree derivatives Ishtiaq Rasool Khan a;b;, Ryoji Ohba
More information2 EBERHARD BECKER ET AL. has a real root. Thus our problem can be reduced to the problem of deciding whether or not a polynomial in one more variable
Deciding positivity of real polynomials Eberhard Becker, Victoria Powers, and Thorsten Wormann Abstract. We describe an algorithm for deciding whether or not a real polynomial is positive semidenite. The
More informationPROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS
PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix
More informationDiscrete Halanay-type inequalities and applications
Nonlinear Analysis 55 (2003) 669 678 www.elsevier.com/locate/na Discrete Halanay-type inequalities and applications Eduardo Liz a;, Anatoli Ivanov b, Juan Bosco Ferreiro c a Departamento de Matematica
More informationMORE CONSEQUENCES OF CAUCHY S THEOREM
MOE CONSEQUENCES OF CAUCHY S THEOEM Contents. The Mean Value Property and the Maximum-Modulus Principle 2. Morera s Theorem and some applications 3 3. The Schwarz eflection Principle 6 We have stated Cauchy
More information1 Groups Examples of Groups Things that are not groups Properties of Groups Rings and Fields Examples...
Contents 1 Groups 2 1.1 Examples of Groups... 3 1.2 Things that are not groups....................... 4 1.3 Properties of Groups... 5 2 Rings and Fields 6 2.1 Examples... 8 2.2 Some Finite Fields... 10
More informationNP-hardness of the stable matrix in unit interval family problem in discrete time
Systems & Control Letters 42 21 261 265 www.elsevier.com/locate/sysconle NP-hardness of the stable matrix in unit interval family problem in discrete time Alejandra Mercado, K.J. Ray Liu Electrical and
More informationChapter 11 - Sequences and Series
Calculus and Analytic Geometry II Chapter - Sequences and Series. Sequences Definition. A sequence is a list of numbers written in a definite order, We call a n the general term of the sequence. {a, a
More informationThe Pacic Institute for the Mathematical Sciences http://www.pims.math.ca pims@pims.math.ca Surprise Maximization D. Borwein Department of Mathematics University of Western Ontario London, Ontario, Canada
More informationRoots of Unity, Cyclotomic Polynomials and Applications
Swiss Mathematical Olympiad smo osm Roots of Unity, Cyclotomic Polynomials and Applications The task to be done here is to give an introduction to the topics in the title. This paper is neither complete
More informationTaylor and Maclaurin Series. Copyright Cengage Learning. All rights reserved.
11.10 Taylor and Maclaurin Series Copyright Cengage Learning. All rights reserved. We start by supposing that f is any function that can be represented by a power series f(x)= c 0 +c 1 (x a)+c 2 (x a)
More informationON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay
ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationX. Hu, R. Shonkwiler, and M.C. Spruill. School of Mathematics. Georgia Institute of Technology. Atlanta, GA 30332
Approximate Speedup by Independent Identical Processing. Hu, R. Shonkwiler, and M.C. Spruill School of Mathematics Georgia Institute of Technology Atlanta, GA 30332 Running head: Parallel iip Methods Mail
More informationLEBESGUE INTEGRATION. Introduction
LEBESGUE INTEGATION EYE SJAMAA Supplementary notes Math 414, Spring 25 Introduction The following heuristic argument is at the basis of the denition of the Lebesgue integral. This argument will be imprecise,
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationMATH32031: Coding Theory Part 15: Summary
MATH32031: Coding Theory Part 15: Summary 1 The initial problem The main goal of coding theory is to develop techniques which permit the detection of errors in the transmission of information and, if necessary,
More informationAbstract Minimal degree interpolation spaces with respect to a nite set of
Numerische Mathematik Manuscript-Nr. (will be inserted by hand later) Polynomial interpolation of minimal degree Thomas Sauer Mathematical Institute, University Erlangen{Nuremberg, Bismarckstr. 1 1, 90537
More informationform, but that fails as soon as one has an object greater than every natural number. Induction in the < form frequently goes under the fancy name \tra
Transnite Ordinals and Their Notations: For The Uninitiated Version 1.1 Hilbert Levitz Department of Computer Science Florida State University levitz@cs.fsu.edu Intoduction This is supposed to be a primer
More informationNew Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems
Systems & Control Letters 43 (21 39 319 www.elsevier.com/locate/sysconle New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems E. Fridman Department of Electrical
More informationMath 361: Homework 1 Solutions
January 3, 4 Math 36: Homework Solutions. We say that two norms and on a vector space V are equivalent or comparable if the topology they define on V are the same, i.e., for any sequence of vectors {x
More informationOn reaching head-to-tail ratios for balanced and unbalanced coins
Journal of Statistical Planning and Inference 0 (00) 0 0 www.elsevier.com/locate/jspi On reaching head-to-tail ratios for balanced and unbalanced coins Tamas Lengyel Department of Mathematics, Occidental
More informationNotes on Iterated Expectations Stephen Morris February 2002
Notes on Iterated Expectations Stephen Morris February 2002 1. Introduction Consider the following sequence of numbers. Individual 1's expectation of random variable X; individual 2's expectation of individual
More information4. Higher Order Linear DEs
4. Higher Order Linear DEs Department of Mathematics & Statistics ASU Outline of Chapter 4 1 General Theory of nth Order Linear Equations 2 Homogeneous Equations with Constant Coecients 3 The Method of
More informationPade approximants and noise: rational functions
Journal of Computational and Applied Mathematics 105 (1999) 285 297 Pade approximants and noise: rational functions Jacek Gilewicz a; a; b;1, Maciej Pindor a Centre de Physique Theorique, Unite Propre
More informationDoc. Math. J. DMV 321 Chern Classes of Fibered Products of Surfaces Mina Teicher Received: March 16, 1998 Revised: September 18, 1998 Communicated by
Doc. Math. J. DMV 31 Chern Classes of Fibered Products of Surfaces Mina Teicher Received: March 16, 1998 Revised: September 18, 1998 Communicated by Thomas Peternell Abstract. In this paper we introduce
More information1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary
Prexes and the Entropy Rate for Long-Range Sources Ioannis Kontoyiannis Information Systems Laboratory, Electrical Engineering, Stanford University. Yurii M. Suhov Statistical Laboratory, Pure Math. &
More informationTHE LARGEST INTERSECTION LATTICE OF A CHRISTOS A. ATHANASIADIS. Abstract. We prove a conjecture of Bayer and Brandt [J. Alg. Combin.
THE LARGEST INTERSECTION LATTICE OF A DISCRIMINANTAL ARRANGEMENT CHRISTOS A. ATHANASIADIS Abstract. We prove a conjecture of Bayer and Brandt [J. Alg. Combin. 6 (1997), 229{246] about the \largest" intersection
More informationContents. 6 Systems of First-Order Linear Dierential Equations. 6.1 General Theory of (First-Order) Linear Systems
Dierential Equations (part 3): Systems of First-Order Dierential Equations (by Evan Dummit, 26, v 2) Contents 6 Systems of First-Order Linear Dierential Equations 6 General Theory of (First-Order) Linear
More informationAlgebraic Expressions and Identities
9 Algebraic Epressions and Identities introduction In previous classes, you have studied the fundamental concepts of algebra, algebraic epressions and their addition and subtraction. In this chapter, we
More informationCalculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm
Calculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm Hartmut Führ fuehr@matha.rwth-aachen.de Lehrstuhl A für Mathematik, RWTH Aachen
More informationChern characters via connections up to homotopy. Marius Crainic. Department of Mathematics, Utrecht University, The Netherlands
Chern characters via connections up to homotopy Marius Crainic Department of Mathematics, Utrecht University, The Netherlands 1 Introduction: The aim of this note is to point out that Chern characters
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More information8 Singular Integral Operators and L p -Regularity Theory
8 Singular Integral Operators and L p -Regularity Theory 8. Motivation See hand-written notes! 8.2 Mikhlin Multiplier Theorem Recall that the Fourier transformation F and the inverse Fourier transformation
More informationNew concepts: Span of a vector set, matrix column space (range) Linearly dependent set of vectors Matrix null space
Lesson 6: Linear independence, matrix column space and null space New concepts: Span of a vector set, matrix column space (range) Linearly dependent set of vectors Matrix null space Two linear systems:
More informationLinear Algebra, 3rd day, Wednesday 6/30/04 REU Info:
Linear Algebra, 3rd day, Wednesday 6/30/04 REU 2004. Info: http://people.cs.uchicago.edu/laci/reu04. Instructor: Laszlo Babai Scribe: Richard Cudney Rank Let V be a vector space. Denition 3.. Let S V,
More informationOn the Third Order Rational Difference Equation
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 27, 32-334 On the Third Order Rational Difference Equation x n x n 2 x (a+x n x n 2 ) T. F. Irahim Department of Mathematics, Faculty of Science Mansoura
More informationLecture Notes: Mathematical/Discrete Structures
Lecture Notes: Mathematical/Discrete Structures Rashid Bin Muhammad, PhD. This booklet includes lecture notes, homework problems, and exam problems from discrete structures course I taught in Fall 2006
More informationSCTT The pqr-method august 2016
SCTT The pqr-method august 2016 A. Doledenok, M. Fadin, À. Menshchikov, A. Semchankau Almost all inequalities considered in our project are symmetric. Hence if plugging (a 0, b 0, c 0 ) into our inequality
More information2 IM Preprint series A, No. 1/2003 If vector x(k) = (x 1 (k) x 2 (k) ::: x n (k)) denotes the time instants in which all jobs started for the k th tim
Eigenvectors of interval matrices over max-plus algebra Katarna Cechlarova Institute of mathematics, Faculty ofscience, P.J. Safarik University, Jesenna 5, 041 54 Kosice, Slovakia, e-mail: cechlarova@science.upjs.sk
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationAsymptotic expansion of multivariate conservative linear operators
Journal of Computational and Applied Mathematics 150 (2003) 219 251 www.elsevier.com/locate/cam Asymptotic expansion of multivariate conservative linear operators A.-J. Lopez-Moreno, F.-J. Muñoz-Delgado
More information1 Introduction Future high speed digital networks aim to serve integrated trac, such as voice, video, fax, and so forth. To control interaction among
On Deterministic Trac Regulation and Service Guarantees: A Systematic Approach by Filtering Cheng-Shang Chang Dept. of Electrical Engineering National Tsing Hua University Hsinchu 30043 Taiwan, R.O.C.
More informationusing the Hamiltonian constellations from the packing theory, i.e., the optimal sphere packing points. However, in [11] it is shown that the upper bou
Some 2 2 Unitary Space-Time Codes from Sphere Packing Theory with Optimal Diversity Product of Code Size 6 Haiquan Wang Genyuan Wang Xiang-Gen Xia Abstract In this correspondence, we propose some new designs
More informationA version of for which ZFC can not predict a single bit Robert M. Solovay May 16, Introduction In [2], Chaitin introd
CDMTCS Research Report Series A Version of for which ZFC can not Predict a Single Bit Robert M. Solovay University of California at Berkeley CDMTCS-104 May 1999 Centre for Discrete Mathematics and Theoretical
More informationGENERAL SOLUTION OF A FIBONACCI-LIKE RECURSION RELATION AND APPLICATIONS
# GENERAL SOLUTION OF A FIBONACCI-LIKE RECURSION RELATION AND APPLICATIONS ALAIN J. PHARES Villanova University, Villanova, PA 19085 (Submitted June 1982) 1. INTRODUCTION: THE COMBINATORICS FUNCTION
More informationAlternative theorems for nonlinear projection equations and applications to generalized complementarity problems
Nonlinear Analysis 46 (001) 853 868 www.elsevier.com/locate/na Alternative theorems for nonlinear projection equations and applications to generalized complementarity problems Yunbin Zhao a;, Defeng Sun
More informationContents. 2 Partial Derivatives. 2.1 Limits and Continuity. Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v. 2.
Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v 260) Contents 2 Partial Derivatives 1 21 Limits and Continuity 1 22 Partial Derivatives 5 23 Directional Derivatives and the Gradient
More informationINRIA Rocquencourt, Le Chesnay Cedex (France) y Dept. of Mathematics, North Carolina State University, Raleigh NC USA
Nonlinear Observer Design using Implicit System Descriptions D. von Wissel, R. Nikoukhah, S. L. Campbell y and F. Delebecque INRIA Rocquencourt, 78 Le Chesnay Cedex (France) y Dept. of Mathematics, North
More informationof Classical Constants Philippe Flajolet and Ilan Vardi February 24, 1996 Many mathematical constants are expressed as slowly convergent sums
Zeta Function Expansions of Classical Constants Philippe Flajolet and Ilan Vardi February 24, 996 Many mathematical constants are expressed as slowly convergent sums of the form C = f( ) () n n2a for some
More informationGenerating Functions of Partitions
CHAPTER B Generating Functions of Partitions For a complex sequence {α n n 0,, 2, }, its generating function with a complex variable q is defined by A(q) : α n q n α n [q n ] A(q). When the sequence has
More informationREMARKS ON THE TIME-OPTIMAL CONTROL OF A CLASS OF HAMILTONIAN SYSTEMS. Eduardo D. Sontag. SYCON - Rutgers Center for Systems and Control
REMARKS ON THE TIME-OPTIMAL CONTROL OF A CLASS OF HAMILTONIAN SYSTEMS Eduardo D. Sontag SYCON - Rutgers Center for Systems and Control Department of Mathematics, Rutgers University, New Brunswick, NJ 08903
More informationsatisfying ( i ; j ) = ij Here ij = if i = j and 0 otherwise The idea to use lattices is the following Suppose we are given a lattice L and a point ~x
Dual Vectors and Lower Bounds for the Nearest Lattice Point Problem Johan Hastad* MIT Abstract: We prove that given a point ~z outside a given lattice L then there is a dual vector which gives a fairly
More informationCentro de Processamento de Dados, Universidade Federal do Rio Grande do Sul,
A COMPARISON OF ACCELERATION TECHNIQUES APPLIED TO THE METHOD RUDNEI DIAS DA CUNHA Computing Laboratory, University of Kent at Canterbury, U.K. Centro de Processamento de Dados, Universidade Federal do
More informationVector Spaces ปร ภ ม เวกเตอร
Vector Spaces ปร ภ ม เวกเตอร 1 5.1 Real Vector Spaces ปร ภ ม เวกเตอร ของจ านวนจร ง Vector Space Axioms (1/2) Let V be an arbitrary nonempty set of objects on which two operations are defined, addition
More informationPOLYNOMIALS OF AN ELLIPTIC CURVE. non-torsion point M such that M mod p is non-zero and non-singular, and it is asymptotically
EXPLICIT VALUATIONS OF DIVISION POLYNOMIALS OF AN ELLIPTIC CURVE y J. Cheon and z S. Hahn Abstract. In this paper, we estimate valuations of division polynomials and compute them explicitely at singular
More informationBalance properties of multi-dimensional words
Theoretical Computer Science 273 (2002) 197 224 www.elsevier.com/locate/tcs Balance properties of multi-dimensional words Valerie Berthe a;, Robert Tijdeman b a Institut de Mathematiques de Luminy, CNRS-UPR
More informationMethods for a Class of Convex. Functions. Stephen M. Robinson WP April 1996
Working Paper Linear Convergence of Epsilon-Subgradient Descent Methods for a Class of Convex Functions Stephen M. Robinson WP-96-041 April 1996 IIASA International Institute for Applied Systems Analysis
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig Uniformly distributed measures in Euclidean spaces by Bernd Kirchheim and David Preiss Preprint-Nr.: 37 1998 Uniformly Distributed
More informationISSUED BY KENDRIYA VIDYALAYA - DOWNLOADED FROM Chapter - 2. (Polynomials)
Chapter - 2 (Polynomials) Key Concepts Constants : A symbol having a fixed numerical value is called a constant. Example : 7, 3, -2, 3/7, etc. are all constants. Variables : A symbol which may be assigned
More informationCoecient bounds for certain subclasses of analytic functions of complex order
Hacettepe Journal of Mathematics and Statistics Volume 45 (4) (2016), 1015 1022 Coecient bounds for certain subclasses of analytic functions of complex order Serap Bulut Abstract In this paper, we introduce
More information