Supplementary Material: Accelerating Adaptive Online Learning by Matrix Approximation

Size: px
Start display at page:

Download "Supplementary Material: Accelerating Adaptive Online Learning by Matrix Approximation"

Transcription

1 Supplemenary Maerial: Acceleraing Adapive Online Learning by Marix Approximaion Yuanyu Wan and Lijun Zhang Naional Key Laboraory for Novel Sofware Technology Nanjing Universiy, Nanjing 1003, China A Theoreical Analysis In his supplemenary maerial, we provide proof of Theorems 1 and. A.1 Supporing Resuls The following resuls are used hroughou our analysis. Lemma 1. (Proposiion 3 of [1]). Le sequence {β } be generaed by ADA-RP. We have R(T ) 1 η [ BΨ+1(β, β + B Ψ (β, β + ] + 1 η B Ψ 1 (β, β + η f (β ) Ψ. Lemma. Le X = i=1 x ix i and A denoe he pseudo-inverse of A, hen x, (X 1/ ) x x, (X 1/ T ) x = r(x 1/ T ). Lemma can be proved in he same way as Lemma 10 of [1]. Theorem 3. (Theorem.3 of [11]). Le 0 <, δ < 1 and S = 1 k R R k n where he enries R i,j of R are independen sandard normal random variables. Then if k = Θ( d+log(1/δ) ), hen for any fixed n d marix A, wih probabiliy 1 δ, simulaneously for all x R d, (1 ) Ax SAx (1 + ) Ax. Based on he above heorem, we derive he following corollary. Corollary 1. Le 0 <, δ < 1 and each enry of r R τ is a Gaussian random variable drawn from N (0, 1/ τ) independenly. Then, if τ = Ω( r+log(t/δ) ), wih probabiliy 1 δ, simulaneously for all = 1,..., T, (1 )C C S S (1 + )C C.

2 Theorem. (Theorem 10 of [0]). Le R be a Gaussian random marix of size p n. Le C = diag(c 1,..., c p ) and S = diag(s 1,..., s p ) be p p diagonal marices, where c i 0 and c i + s i = 1 for all i. Le M = C + 1 n SRR S and r = i s i. P r(λ 1 (M) 1 + ) q exp ( cn ), max i (s i )r ) P r(λ p (M) 1 ) q exp ( cn max i (s i )r, where he consan c is a leas 1/3, and q is he rank of S. Based on he above heorem, we derive he following corollary. σ i, r α+σi = max 1 T σ 1. Le K = αi d + C C, K = αi d + S S, and Corollary. Le c 1/3, α > 0, σ i = λ i(c C ), r = i max r and σ 1 = k T Ĩ = K / K K / r. If τ σ 1 c (α+σ 1 1 δ, simulaneously for all = 1,, T, dt log ) δ (1 )I d Ĩ (1 + )I d., hen wih probabiliy a leas A. Proof of Theorem 1 Le X denoe S S. Firs, we consider bounding he firs erm in he upper bound of Lemma 1. Wih probabiliy 1 δ, we have B Ψ+1 (β, β + B Ψ (β, β + = 1 β β +1, ( X 1/ 1/ +1 X )(β β + 1 β β +1, 1 + X 1/ +1 (β β + 1 β β +1, 1 X 1/ (β β + 1 β β +1, (X 1/ +1 X1/ )(β β + + β β +1, (X 1/ +1 + X1/ )(β β + 1 β β +1 (X 1/ +1 X1/ ) + β β +1, (X 1/ +1 + X1/ )(β β + 1 β β +1 r(x 1/ +1 X1/ ) + β β +1, (X 1/ +1 + X1/ )(β β +

3 where he firs inequaliy is due o Corollary 1. Thus, we can ge 1 [ BΨ+1(β, β + B Ψ (β, β + ] β β +1 r(x 1/ +1 X1/ ) + β β +1, (X 1/ +1 + X1/ )(β β + 1 max T β β r(x 1/ T ) 1 β β 1 r(x 1/ + max T β β X 1/ β β 1 r(x 1/. Noe ha β 1 = 0, hen B Ψ1 (β, β = 1 β 1/, (σi d + X β 1 σ β + + () β r(x 1/ where he inequaliy is due o Corollary 1. Then, we consider he upper bound of T f (β ) Ψ. Wih probabiliy 1 δ, we have 1 f (β ) Ψ = 1/ g, (σi d + X ) g 1 g, (X ) 1/ g = l (β x ) x, (X ) 1/ x 1 1 (3) where he inequaliy is due o Corollary 1. According o Lemma, we have f (β ) Ψ T l (β x ) x, (X ) 1/ x 1 max T l (β x ) T x, (X ) 1/ x 1 1 max T l (β x ) r(x 1/ T ). We complee he proof by subsiuing (3), (), and (5) ino Lemma 1. (5) A.3 Proof of Theorem Inspired by he proof of Theorem 1, we can derive Theorem by respecively bounding each erm in he upper bound of Lemma 1. Before ha, we need o derive he lower and upper bounds of (S S ) 1/ based on Corollary.

4 Le he SVD of C be C = UΣV where U R d d, Σ R d d, V R d. According o Corollary, wih probabiliy a leas 1 δ, simulaneously for all = 1,..., T, S S = K αi d = K 1/ Ĩ K 1/ αi d and (1 + )K αi d = (1 + )C C + αi d = U ( (1 + )ΣΣ + αi d ) U S S + αi d = K αi d + αi d = K 1/ Ĩ K 1/ αi d + αi d (1 )K αi d + αi d = (1 )C C. Then simulaneously for all = 1,..., T, we have and (S S ) 1/ 1 + U(ΣΣ) 1/ U + αui d U = 1 + X 1/ + αi d (6) (S S ) 1/ = (S S ) 1/ + αi d αi d ( (S S ) + αi d ) 1/ αid 1 X 1/ αi d. Then we consider bounding he firs erm in he upper bound of Lemma 1. Le X denoe S S. Simulaneously for all = 1,..., T, we have B Ψ+1 (β, β + B Ψ (β, β + = 1 β β +1, ( X 1/ 1/ +1 X )(β β + 1 β β +1, 1 + X 1/ +1 (β β + 1 β β +1, 1 X 1/ )(β β β β +1, αi d (β β + = 1 β β +1, 1 + X 1/ +1 (β β + 1 β β +1, 1 X 1/ )(β β + + α (β β + 1 β β +1 r(x 1/ +1 X1/ ) + β β +1, (X 1/ +1 + X1/ )(β β + + α (β β + (7)

5 where he firs inequaliy is due o (6), (7) and he las inequaliy has been proved in he proof of Theorem 1. Thus, we can ge [ BΨ+1 (β, β + B Ψ (β, β + ] Noe ha β 1 = 0, hen 1 max T β β r(x 1/ T ) 1 β β 1 r(x 1/ + max T β β β β 1 r(x 1/ X 1/ + α(t 1) max T β β. B Ψ1 (β, β = 1 β 1/, (σi d + X β 1 σ β + + β r(x 1/ + 1 α β. (8) (9) Before considering he upper bound of T f (β ) Ψ, we need o derive he upper bound of H. Le he SVD of S be S = UΣV where U R d d, Σ R d d, V R d. We also have, for all = 1,..., T, H = σi d + (S S ) 1/ = U ( σi d + (ΣΣ) 1/) U U ( αi d + (ΣΣ) ) 1/ U = (αi d + S S ) 1/ due o σ α λ i (S S ) + α λ i (S S ) for all i = 1,, d. Then according o Corollary, wih probabiliy a leas 1 δ, simulaneously for all = 1,..., T, H ( (αi d + S S ) 1/) ( 1/ = (K Ĩ K 1/ 1 (K ) 1/ 1 ( = (αid + X ) ) 1/. 1 1 ) ) 1/ Thus, we can ge f (β ) Ψ = g, H g g, ( (αi d + X ) ) 1/ g 1 = l (β x ) 1 x, (X ) 1/ x.

6 According o Lemma, we have f (β ) Ψ 1 max T l (β x ) x, (X ) 1/ x (10) max 1 T l (β x ) r(x 1/ T ). We complee he proof by subsiuing (8), (9), and (10) ino Lemma 1. A. Proof of Corollary 1 Le C = UΣV be he singular value decomposiion of C. Noice ha U R r, ΣV R r d. According o Theorem 3, we have if τ = Θ( r+log(1/δ) ), hen simulaneously x R r, wih probabiliy 1 δ, (1 ) Ux R Ux (1 + ) Ux Le y R d be arbirary vecor, hen C y = UΣV y = Ux where x = ΣV y R r. Then we have y S S y = y C R R C y = R Ux (1 + ) Ux = (1 + )y C C y and y S S y = y C R R C y = R Ux (1 ) Ux = (1 )y C C y. Then, we have (1 )C C S S (1+)C C wih probabiliy 1 δ, provided τ = Ω( r+log(1/δ) ). Using he union bound, we have if τ = Ω( r+log(t/δ) ), wih probabiliy 1 δ, simulaneously for all = 1,..., T, (1 )C C S S (1 + )C C. A.5 Proof of Corollary Le he SVD of C be C = UΣV where U R d d, Σ R d d, V R d. Then we have K = U(αI d + ΣΣ )U and Ĩ = K / K K / = K / (αi d + C R R C )K / ( = U αi d (αi d + ΣΣ) + (αi p + ΣΣ) / ΣV R R VΣ(αI d + ΣΣ ) /) U ( = U αi d (αi d + ΣΣ) + (αi p + ΣΣ) / ΣRR Σ(αI d + ΣΣ ) /) U where R = V R R d τ is a Gaussian random marix due o ha V is an orhogonal marix and R is a Gaussian random marix.

7 α α+σ i Le c i = a leas 1 δ, and s i = r σ 1 σ i α+σi. Then according o Theorem, wih probabiliy (1 )I d Ĩ (1 + )I d provided τ c (α+σ1 δ where he consan c is a leas 1/3. Using he union bound, we complee he proof. ) log d

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Exercises: Similarity Transformation

Exercises: Similarity Transformation Exercises: Similariy Transformaion Problem. Diagonalize he following marix: A [ 2 4 Soluion. Marix A has wo eigenvalues λ 3 and λ 2 2. Since (i) A is a 2 2 marix and (ii) i has 2 disinc eigenvalues, we

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3

More information

Math 315: Linear Algebra Solutions to Assignment 6

Math 315: Linear Algebra Solutions to Assignment 6 Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen

More information

Tracking Adversarial Targets

Tracking Adversarial Targets A. Proofs Proof of Lemma 3. Consider he Bellman equaion λ + V π,l x, a lx, a + V π,l Ax + Ba, πax + Ba. We prove he lemma by showing ha he given quadraic form is he unique soluion of he Bellman equaion.

More information

The General Linear Test in the Ridge Regression

The General Linear Test in the Ridge Regression ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge

More information

Approximation Algorithms for Unique Games via Orthogonal Separators

Approximation Algorithms for Unique Games via Orthogonal Separators Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

Machine Learning 4771

Machine Learning 4771 ony Jebara, Columbia Universiy achine Learning 4771 Insrucor: ony Jebara ony Jebara, Columbia Universiy opic 20 Hs wih Evidence H Collec H Evaluae H Disribue H Decode H Parameer Learning via JA & E ony

More information

On a Fractional Stochastic Landau-Ginzburg Equation

On a Fractional Stochastic Landau-Ginzburg Equation Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn

More information

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004 ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform

More information

U( θ, θ), U(θ 1/2, θ + 1/2) and Cauchy (θ) are not exponential families. (The proofs are not easy and require measure theory. See the references.

U( θ, θ), U(θ 1/2, θ + 1/2) and Cauchy (θ) are not exponential families. (The proofs are not easy and require measure theory. See the references. Lecure 5 Exponenial Families Exponenial families, also called Koopman-Darmois families, include a quie number of well known disribuions. Many nice properies enjoyed by exponenial families allow us o provide

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

2 Some Property of Exponential Map of Matrix

2 Some Property of Exponential Map of Matrix Soluion Se for Exercise Session No8 Course: Mahemaical Aspecs of Symmeries in Physics, ICFP Maser Program for M 22nd, January 205, a Room 235A Lecure by Amir-Kian Kashani-Poor email: kashani@lpensfr Exercise

More information

Empirical Process Theory

Empirical Process Theory Empirical Process heory 4.384 ime Series Analysis, Fall 27 Reciaion by Paul Schrimpf Supplemenary o lecures given by Anna Mikusheva Ocober 7, 28 Reciaion 7 Empirical Process heory Le x be a real-valued

More information

Appendix to Online l 1 -Dictionary Learning with Application to Novel Document Detection

Appendix to Online l 1 -Dictionary Learning with Application to Novel Document Detection Appendix o Online l -Dicionary Learning wih Applicaion o Novel Documen Deecion Shiva Prasad Kasiviswanahan Huahua Wang Arindam Banerjee Prem Melville A Background abou ADMM In his secion, we give a brief

More information

Optimal Paired Choice Block Designs. Supplementary Material

Optimal Paired Choice Block Designs. Supplementary Material Saisica Sinica: Supplemen Opimal Paired Choice Block Designs Rakhi Singh 1, Ashish Das 2 and Feng-Shun Chai 3 1 IITB-Monash Research Academy, Mumbai, India 2 Indian Insiue of Technology Bombay, Mumbai,

More information

DISCRETE GRONWALL LEMMA AND APPLICATIONS

DISCRETE GRONWALL LEMMA AND APPLICATIONS DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is

More information

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990), SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

GRADIENT ESTIMATES FOR A SIMPLE PARABOLIC LICHNEROWICZ EQUATION. Osaka Journal of Mathematics. 51(1) P.245-P.256

GRADIENT ESTIMATES FOR A SIMPLE PARABOLIC LICHNEROWICZ EQUATION. Osaka Journal of Mathematics. 51(1) P.245-P.256 Tile Auhor(s) GRADIENT ESTIMATES FOR A SIMPLE PARABOLIC LICHNEROWICZ EQUATION Zhao, Liang Ciaion Osaka Journal of Mahemaics. 51(1) P.45-P.56 Issue Dae 014-01 Tex Version publisher URL hps://doi.org/10.18910/9195

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

V L. DT s D T s t. Figure 1: Buck-boost converter: inductor current i(t) in the continuous conduction mode.

V L. DT s D T s t. Figure 1: Buck-boost converter: inductor current i(t) in the continuous conduction mode. ECE 445 Analysis and Design of Power Elecronic Circuis Problem Se 7 Soluions Problem PS7.1 Erickson, Problem 5.1 Soluion (a) Firs, recall he operaion of he buck-boos converer in he coninuous conducion

More information

ψ(t) = V x (0)V x (t)

ψ(t) = V x (0)V x (t) .93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

More information

Double system parts optimization: static and dynamic model

Double system parts optimization: static and dynamic model Double sysem pars opmizaon: sac and dynamic model 1 Inroducon Jan Pelikán 1, Jiří Henzler 2 Absrac. A proposed opmizaon model deals wih he problem of reserves for he funconal componens-pars of mechanism

More information

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type EMS SCM join meeing Barcelona, May 28-30, 2015 On sochasic parial differenial equaions of parabolic ype Isván Gyöngy School of Mahemaics and Maxwell Insiue Edinburgh Universiy 1 I. Filering problem II.

More information

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations

Research Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations Hindawi Publishing Corporaion Boundary Value Problems Volume 11, Aricle ID 19156, 11 pages doi:1.1155/11/19156 Research Aricle Exisence and Uniqueness of Periodic Soluion for Nonlinear Second-Order Ordinary

More information

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given

More information

Math-Net.Ru All Russian mathematical portal

Math-Net.Ru All Russian mathematical portal Mah-NeRu All Russian mahemaical poral Roman Popovych, On elemens of high order in general finie fields, Algebra Discree Mah, 204, Volume 8, Issue 2, 295 300 Use of he all-russian mahemaical poral Mah-NeRu

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Online Convex Optimization Example And Follow-The-Leader

Online Convex Optimization Example And Follow-The-Leader CSE599s, Spring 2014, Online Learning Lecure 2-04/03/2014 Online Convex Opimizaion Example And Follow-The-Leader Lecurer: Brendan McMahan Scribe: Sephen Joe Jonany 1 Review of Online Convex Opimizaion

More information

Robert Kollmann. 6 September 2017

Robert Kollmann. 6 September 2017 Appendix: Supplemenary maerial for Tracable Likelihood-Based Esimaion of Non- Linear DSGE Models Economics Leers (available online 6 Sepember 207) hp://dx.doi.org/0.06/j.econle.207.08.027 Rober Kollmann

More information

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.

Homework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross. Homework (Sas 6, Winer 7 Due Tuesday April 8, in class Quesions are derived from problems in Sochasic Processes by S. Ross.. A sochasic process {X(, } is said o be saionary if X(,..., X( n has he same

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

Chapter Three Systems of Linear Differential Equations

Chapter Three Systems of Linear Differential Equations Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

Lie Derivatives operator vector field flow push back Lie derivative of

Lie Derivatives operator vector field flow push back Lie derivative of Lie Derivaives The Lie derivaive is a mehod of compuing he direcional derivaive of a vecor field wih respec o anoher vecor field We already know how o make sense of a direcional derivaive of real valued

More information

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!

Finish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details! MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his

More information

Model Reduction for Dynamical Systems Lecture 6

Model Reduction for Dynamical Systems Lecture 6 Oo-von-Guericke Universiä Magdeburg Faculy of Mahemaics Summer erm 07 Model Reducion for Dynamical Sysems ecure 6 v eer enner and ihong Feng Max lanck Insiue for Dynamics of Complex echnical Sysems Compuaional

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

Asymptotic Equipartition Property - Seminar 3, part 1

Asymptotic Equipartition Property - Seminar 3, part 1 Asympoic Equipariion Propery - Seminar 3, par 1 Ocober 22, 2013 Problem 1 (Calculaion of ypical se) To clarify he noion of a ypical se A (n) ε and he smalles se of high probabiliy B (n), we will calculae

More information

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology

in Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology Risk and Saey in Engineering Pro. Dr. Michael Havbro Faber ETH Zurich, Swizerland Conens o Today's Lecure Inroducion o ime varian reliabiliy analysis The Poisson process The ormal process Assessmen o he

More information

Supplementary Document

Supplementary Document Saisica Sinica (2013): Preprin 1 Supplemenary Documen for Funcional Linear Model wih Zero-value Coefficien Funcion a Sub-regions Jianhui Zhou, Nae-Yuh Wang, and Naisyin Wang Universiy of Virginia, Johns

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

More information

Lecture 4: Processes with independent increments

Lecture 4: Processes with independent increments Lecure 4: Processes wih independen incremens 1. A Wienner process 1.1 Definiion of a Wienner process 1.2 Reflecion principle 1.3 Exponenial Brownian moion 1.4 Exchange of measure (Girsanov heorem) 1.5

More information

ψ ( t) = c n ( t ) n

ψ ( t) = c n ( t ) n p. 31 PERTURBATION THEORY Given a Hamilonian H ( ) = H + V( ) where we know he eigenkes for H H n = En n we ofen wan o calculae changes in he ampliudes of n induced by V( ) : where ψ ( ) = c n ( ) n n

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

be a Householder matrix. Then prove the followings H = I 2 uut Hu = (I 2 uu u T u )u = u 2 uut u

be a Householder matrix. Then prove the followings H = I 2 uut Hu = (I 2 uu u T u )u = u 2 uut u MATH 434/534 Theoretical Assignment 7 Solution Chapter 7 (71) Let H = I 2uuT Hu = u (ii) Hv = v if = 0 be a Householder matrix Then prove the followings H = I 2 uut Hu = (I 2 uu )u = u 2 uut u = u 2u =

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

Singular Value Decomposition

Singular Value Decomposition Singular Value Decomposition Motivatation The diagonalization theorem play a part in many interesting applications. Unfortunately not all matrices can be factored as A = PDP However a factorization A =

More information

Fréchet derivatives and Gâteaux derivatives

Fréchet derivatives and Gâteaux derivatives Fréche derivaives and Gâeaux derivaives Jordan Bell jordan.bell@gmail.com Deparmen of Mahemaics, Universiy of Torono April 3, 2014 1 Inroducion In his noe all vecor spaces are real. If X and Y are normed

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Lecture 12: Multiple Hypothesis Testing

Lecture 12: Multiple Hypothesis Testing ECE 830 Fall 00 Saisical Signal Processing insrucor: R. Nowak, scribe: Xinjue Yu Lecure : Muliple Hypohesis Tesing Inroducion In many applicaions we consider muliple hypohesis es a he same ime. Example

More information

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems. Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need

More information

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.

More information

Operators related to the Jacobi setting, for all admissible parameter values

Operators related to the Jacobi setting, for all admissible parameter values Operaors relaed o he Jacobi seing, for all admissible parameer values Peer Sjögren Universiy of Gohenburg Join work wih A. Nowak and T. Szarek Alba, June 2013 () 1 / 18 Le Pn α,β be he classical Jacobi

More information

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS

BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS BOUNDEDNESS OF MAXIMAL FUNCTIONS ON NON-DOUBLING MANIFOLDS WITH ENDS XUAN THINH DUONG, JI LI, AND ADAM SIKORA Absrac Le M be a manifold wih ends consruced in [2] and be he Laplace-Belrami operaor on M

More information

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k) Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be

More information

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Multi-component Levi Hierarchy and Its Multi-component Integrable Coupling System

Multi-component Levi Hierarchy and Its Multi-component Integrable Coupling System Commun. Theor. Phys. (Beijing, China) 44 (2005) pp. 990 996 c Inernaional Academic Publishers Vol. 44, No. 6, December 5, 2005 uli-componen Levi Hierarchy and Is uli-componen Inegrable Coupling Sysem XIA

More information

Lecture 2 October ε-approximation of 2-player zero-sum games

Lecture 2 October ε-approximation of 2-player zero-sum games Opimizaion II Winer 009/10 Lecurer: Khaled Elbassioni Lecure Ocober 19 1 ε-approximaion of -player zero-sum games In his lecure we give a randomized ficiious play algorihm for obaining an approximae soluion

More information

Supplementary Material

Supplementary Material Dynamic Global Games of Regime Change: Learning, Mulipliciy and iming of Aacks Supplemenary Maerial George-Marios Angeleos MI and NBER Chrisian Hellwig UCLA Alessandro Pavan Norhwesern Universiy Ocober

More information

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial

More information

Stochastic Modelling in Finance - Solutions to sheet 8

Stochastic Modelling in Finance - Solutions to sheet 8 Sochasic Modelling in Finance - Soluions o shee 8 8.1 The price of a defaulable asse can be modeled as ds S = µ d + σ dw dn where µ, σ are consans, (W ) is a sandard Brownian moion and (N ) is a one jump

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

Monotonic Solutions of a Class of Quadratic Singular Integral Equations of Volterra type

Monotonic Solutions of a Class of Quadratic Singular Integral Equations of Volterra type In. J. Conemp. Mah. Sci., Vol. 2, 27, no. 2, 89-2 Monoonic Soluions of a Class of Quadraic Singular Inegral Equaions of Volerra ype Mahmoud M. El Borai Deparmen of Mahemaics, Faculy of Science, Alexandria

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

MATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.

MATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices. MATH 2050 Assignmen 9 Winer 206 Do no need o hand in Noe ha he final exam also covers maerial afer HW8, including, for insance, calculaing deerminan by row operaions, eigenvalues and eigenvecors, similariy

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Transform Techniques. Moment Generating Function

Transform Techniques. Moment Generating Function Transform Techniques A convenien way of finding he momens of a random variable is he momen generaing funcion (MGF). Oher ransform echniques are characerisic funcion, z-ransform, and Laplace ransform. Momen

More information

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive. LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis

More information

6. Stochastic calculus with jump processes

6. Stochastic calculus with jump processes A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio

More information

SINGULAR VALUE DECOMPOSITION

SINGULAR VALUE DECOMPOSITION DEPARMEN OF MAHEMAICS ECHNICAL REPOR SINGULAR VALUE DECOMPOSIION Andrew Lounsbury September 8 No 8- ENNESSEE ECHNOLOGICAL UNIVERSIY Cookeville, N 38 Singular Value Decomposition Andrew Lounsbury Department

More information

Backward stochastic dynamics on a filtered probability space

Backward stochastic dynamics on a filtered probability space Backward sochasic dynamics on a filered probabiliy space Gechun Liang Oxford-Man Insiue, Universiy of Oxford based on join work wih Terry Lyons and Zhongmin Qian Page 1 of 15 gliang@oxford-man.ox.ac.uk

More information

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn

ψ ( t) = c n ( t) t n ( )ψ( ) t ku t,t 0 ψ I V kn MIT Deparmen of Chemisry 5.74, Spring 4: Inroducory Quanum Mechanics II p. 33 Insrucor: Prof. Andrei Tokmakoff PERTURBATION THEORY Given a Hamilonian H ( ) = H + V ( ) where we know he eigenkes for H H

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Representation of Stochastic Process by Means of Stochastic Integrals

Representation of Stochastic Process by Means of Stochastic Integrals Inernaional Journal of Mahemaics Research. ISSN 0976-5840 Volume 5, Number 4 (2013), pp. 385-397 Inernaional Research Publicaion House hp://www.irphouse.com Represenaion of Sochasic Process by Means of

More information

Basic Entropies for Positive-Definite Matrices

Basic Entropies for Positive-Definite Matrices Journal of Mahemaics and Sysem Science 5 (05 3-3 doi: 0.65/59-59/05.04.00 D DAVID PUBLISHING Basic nropies for Posiive-Definie Marices Jun Ichi Fuii Deparmen of Ars and Sciences (Informaion Science, Osaa

More information

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial

More information

Online Learning with Partial Feedback. 1 Online Mirror Descent with Estimated Gradient

Online Learning with Partial Feedback. 1 Online Mirror Descent with Estimated Gradient Avance Course in Machine Learning Spring 2010 Online Learning wih Parial Feeback Hanous are joinly prepare by Shie Mannor an Shai Shalev-Shwarz In previous lecures we alke abou he general framework of

More information

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems

On Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes

More information

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind

Orthogonal Rational Functions, Associated Rational Functions And Functions Of The Second Kind Proceedings of he World Congress on Engineering 2008 Vol II Orhogonal Raional Funcions, Associaed Raional Funcions And Funcions Of The Second Kind Karl Deckers and Adhemar Bulheel Absrac Consider he sequence

More information

The Singular Value Decomposition and Least Squares Problems

The Singular Value Decomposition and Least Squares Problems The Singular Value Decomposition and Least Squares Problems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo September 27, 2009 Applications of SVD solving

More information

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering

More information