Elementary mathematical and computational tools

for electrical and computer engineers using

MATLAB 2nd ed Edition Manassah pdf download

https://ebookgate.com/product/elementary-mathematical-and-
computational-tools-for-electrical-and-computer-engineers-using-
matlab-2nd-ed-edition-manassah/

Get Instant Ebook Downloads – Browse at https://ebookgate.com

Instant digital products (PDF, ePub, MOBI) available
Download now and explore formats that suit you...

MATLAB for Engineers 2nd Edition Holly Moore

https://ebookgate.com/product/matlab-for-engineers-2nd-edition-holly-
moore/

ebookgate.com

Analysis and Design of Control Systems Using MATLAB 2nd ed

Edition Dukkipati

https://ebookgate.com/product/analysis-and-design-of-control-systems-
using-matlab-2nd-ed-edition-dukkipati/

ebookgate.com

Introduction to Computational Modeling Using C and Open

Source Tools Garrido

https://ebookgate.com/product/introduction-to-computational-modeling-
using-c-and-open-source-tools-garrido/

ebookgate.com

Computational Colour Science using MATLAB 1st Edition

Professor Stephen Westland

https://ebookgate.com/product/computational-colour-science-using-
matlab-1st-edition-professor-stephen-westland/

ebookgate.com
MATLAB for Control Engineers Katsuhiko Ogata

https://ebookgate.com/product/matlab-for-control-engineers-katsuhiko-
ogata/

ebookgate.com

Essential MATLAB for Scientists and Engineers Second

Edition Brian Hahn

https://ebookgate.com/product/essential-matlab-for-scientists-and-
engineers-second-edition-brian-hahn/

ebookgate.com

Mathematical tools 1st ed Edition Franklin Coyle

https://ebookgate.com/product/mathematical-tools-1st-ed-edition-
franklin-coyle/

ebookgate.com

Mathematical Physics Applied Mathematics for Scientists

and Engineers 2nd Edition Bruce R. Kusse

https://ebookgate.com/product/mathematical-physics-applied-
mathematics-for-scientists-and-engineers-2nd-edition-bruce-r-kusse/

ebookgate.com

MATLAB for control system engineers Second Edition

Dukkipati

https://ebookgate.com/product/matlab-for-control-system-engineers-
second-edition-dukkipati/

ebookgate.com
7425_C000.qxd 6/9/06 16:53 Page i

ELEMENTARY
MATHEMATICAL and
COMPUTATIONAL TOOLS
for ELECTRICAL and
COMPUTER ENGINEERS
USING MATLAB®
Second Edition
7425_C000.qxd 6/9/06 16:53 Page ii
7425_C000.qxd 6/9/06 16:53 Page iii

ELEMENTARY
MATHEMATICAL and
COMPUTATIONAL TOOLS
for ELECTRICAL and
COMPUTER ENGINEERS
USING MATLAB®
Second Edition

Jamal T. Manassah

Boca Raton London New York

CRC is an imprint of the Taylor & Francis Group,

an informa business
7425_C000.qxd 6/9/06 16:53 Page iv

MATLAB® is a trademark of The MathWorks, Inc. and is used with permission. The MathWorks does not
warrant the accuracy of the text or exercises in this book. This book’s use or discussion of MATLAB®
software or related products does not constitute endorsement or sponsorship by The MathWorks of a
particular pedagogical approach or particular use of the MATLAB® software.

CRC Press
Taylor & Francis Group
6000 Broken Sound Parkway NW, Suite 300
Boca Raton, FL 33487-2742

© 2007 by Taylor & Francis Group, LLC

CRC Press is an imprint of Taylor & Francis Group, an Informa business

No claim to original U.S. Government works

Printed in the United States of America on acid-free paper
10 9 8 7 6 5 4 3 2 1

International Standard Book Number-10: 0-8493-7425-1 (Hardcover)

International Standard Book Number-13: 978-0-8493-7425-8 (Hardcover)

This book contains information obtained from authentic and highly regarded sources. Reprinted material
is quoted with permission, and sources are indicated. A wide variety of references are listed. Reasonable
efforts have been made to publish reliable data and information, but the author and the publisher cannot
assume responsibility for the validity of all materials or for the consequences of their use.

No part of this book may be reprinted, reproduced, transmitted, or utilized in any form by any electronic,
mechanical, or other means, now known or hereafter invented, including photocopying, microfilming,
and recording, or in any information storage or retrieval system, without written permission from the
publishers.

For permission to photocopy or use material electronically from this work, please access www.copyright.
com (http://www.copyright.com/) or contact the Copyright Clearance Center, Inc. (CCC) 222 Rosewood
Drive, Danvers, MA 01923, 978-750-8400. CCC is a not-for-profit organization that provides licenses and
registration for a variety of users. For organizations that have been granted a photocopy license by the
CCC, a separate system of payment has been arranged.

Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are
used only for identification and explanation without intent to infringe.

Library of Congress Cataloging-in-Publication Data

Manassah, Jamal T.
Elementary mathematical and computational tools for electrical and computer
engineers using MATLAB / Jamal T. Manassah. -- 2nd ed.
p. cm.
Includes bibliographical references and index.
ISBN-13: 978-0-8493-7425-8 (alk. paper)
ISBN-10: 0-8493-7425-1 (alk. paper)
1. Electric engineering--Mathematics. 2. Computer science--Mathematics. 3.
MATLAB. I. Title.

TK153.M362 2007
510.2’46213--dc22 2 006018608

Visit the Taylor & Francis Web site at

http://www.taylorandfrancis.com

and the CRC Press Web site at

http://www.crcpress.com
7425_C000.qxd 6/9/06 16:53 Page v

Author

Jamal T. Manassah has been a professor of electrical engineering at the City

College of New York since 1981. He earned his B.Sc. degree in physics from
the American University of Beirut, and his M.A. and Ph.D. in theoretical
physics from Columbia University. Dr. Manassah was a member of the
Institute for Advanced Study. His current research interests are in theoretical
and computational quantum and nonlinear optics, and in photonics.
7425_C000.qxd 6/9/06 16:53 Page vi
7425_C000.qxd 6/9/06 16:53 Page vii

Preface to the First Edition

This book is mostly based on a series of notes for a primer course in electrical
and computer engineering that I taught at the City College of New York
School of Engineering. Each week, the class met for an hour of lecture and a
three-hour computer laboratory session where students were divided into
small groups of 12 to 15 students each. The students met in an informal learn-
ing community setting, a computer laboratory, where each student had the
exclusive use of a PC. The small size of the groups permitted a great deal of
individualized instruction, which was a key ingredient to cater successfully
to the needs of students with heterogeneous high school backgrounds.
A student usually takes this course in the second semester of his or her
freshman year. Typically, the student would have completed one semester of
college calculus, and would be enrolled in the second course of the college
calculus sequence and in the first course of the physics sequence for students
in the physical sciences and engineering.
My purpose in developing this book is to help bring the beginner engi-
neering student’s analytical and computational skills to a level of competency
that would permit him or her to participate, enjoy, and succeed in subsequent
electrical and computer engineering courses. My experience indicates that the
lack of mastery of fundamental quantitative tools is the main impediment to
a student’s progress in engineering studies.
The specific goals of this book are:

1. To make you more comfortable applying the mathematics and

physics that you learned in high school or in college courses,
through interactive activities.
2. To introduce you, through examples, to many new practical tools of
mathematics, including discrete variables material that is essential to
your success in future electrical engineering courses.
3. To instruct you in the use of a powerful computer program,
MATLAB®, which was designed to be simultaneously user-friendly
and powerful in tackling efficiently the most demanding problems
of engineering and sciences.
4. To give you, through the applications and examples covered,
glimpses of some of the fascinating problems that an electrical or
computer engineer solves in the course of completing many of his
or her design projects.

My experience indicates that you can achieve the above goals through the
following work habits that I usually recommend to my own students:

● Read carefully the material from this book that is assigned to you
by your instructor for the upcoming week, and make sure to solve
7425_C000.qxd 6/9/06 16:53 Page viii

the suggested preparatory exercises in advance of the weekly

lecture.
● Attend the lecture and follow closely the material presented, in par-
ticular the solutions to the more difficult preparatory exercises and
the demonstrations.
● Following the lecture, make a list of questions on the preparatory
material to which you still seek answers, and ask your instructor
for help and clarification on these questions, preferably in the first
30 minutes of your computer lab session.
● Complete the in-class exercises during the computer lab session. If
you have not finished solving all in-class exercises, make sure you
complete them on your own, when the lab is open, or at home if you
own a computer, and certainly before the next class session, along
with the problems designated in the book as homework problems
and assigned to you by your instructor.

In managing this course, I found it helpful for both students and instructors
to require each student to solve all problems in a bound notebook. The advan-
tage to the student is to have easy access to his or her previous work, personal
notes, and reminders that he or she made as the course progressed. The advan-
tage to the instructor is to enhance his or her ability to assess, more easily and
readily, an individual student’s progress as the semester progresses.
This book may be used for self-study by readers with perhaps a little more
mathematical maturity acquired through a second semester of college calcu-
lus. The advanced reader of this book who is familiar with numerical meth-
ods will note that, in some instances, I did not follow the canonical order for
the sequence of presentation of certain algorithms, thus sacrificing some opti-
mality in the structure of some of the elementary programs included. This
was necessitated by the goal I set for this book, which is to introduce both
analytical and computational tools simultaneously.
The sections of this book that are marked with asterisks include material
that I assigned as projects to students with either strong theoretical interest
or more mathematical maturity than a typical second semester freshman
student. Although incorporated in the text, they can be skipped in a first
reading. I hope that, by their inclusion, I will facilitate to the interested
reader a smooth transition to some new mathematical concepts and compu-
tational tools that are of particular interest to electrical engineers.
This text greatly benefited from course material previously prepared by
my colleagues in the departments of electrical engineering and computer
science at City College of the City University of New York, in particular,
P. Combettes, I. Gladkova, B. Gross, and F. Thau. They provided either the start-
ing point for my subsequent efforts in this course, or the peer critique for the
early versions of this manuscript. I owe them many thanks and, of course, do
not hold them responsible for any of the remaining imperfections in the text.
The preparation of this book also owes a lot to my students. Their ques-
tions and interest in the material contributed to many modifications in the
7425_C000.qxd 6/9/06 16:53 Page ix

order and in the presentation of the different chapters. Their desire for work-
ing out more applications led me to expand the scope of the examples and
exercises included in the text. I am grateful to all of them. I am also grateful
to Erwin Cohen, who introduced me to the fine team at CRC Press, and to
Jerry Papke whose stewardship of the project from start to end at CRC Press
was most supportive and pleasant. The editorial and production teams at
CRC in particular, Samar Haddad, the project editor, deserve credit for the
quality of the final product rendering. Naomi Fernandes and her colleagues
at The MathWorks Inc. kindly provided me with a copy of the new release
of MATLAB for which I am grateful.
I dedicate this book to Azza, Tala, and Nigh whose support and love
always made difficult tasks a lot easier.
Jamal T. Manassah
New York, January 2001
7425_C000.qxd 6/9/06 16:53 Page x
7425_C000.qxd 6/9/06 16:53 Page xi

Preface to the Second Edition

Comments received from readers and my own students invariably sug-

gested the desirability for expanding Chapter 1 in a manner that will reduce
the need to refer as often to the MATLAB® online help or to a MATLAB user
manual. It was also suggested that the reference value of this book might be
enhanced if a chapter on symbolic manipulation with MATLAB was added.
This second edition is responsive to both of these requests.
The MATLAB syntax was also updated to conform to the current version
of MATLAB. I have also added a number of special sections designated to
the more advanced students or for those students seeking additional chal-
lenges. Typically, these sections were assigned to outstanding students who
supplemented the regular course homework with special assignments. The
continuous need to challenge these gifted students is specially rewarding
when teaching in an urban university where many promising future scien-
tific talents are discovered only after these students are enrolled in college-
level courses.
Now that my modest experiment with introducing an integrated, empiri-
cal, and interdisciplinary style to teaching beginner engineering students
applied mathematics is approximately 10 years in progress, it may be useful
to reflect on the successes and failures of this particular effort:

• Almost everybody in my academic circle now recognizes the use-

fulness of supplementing and enhancing the regular rigorous stan-
dard mathematical curriculum for engineering students with a
hands-on elementary applied course.
• Those of my colleagues who still prefer to have the different mate-
rial taught in a more compartmentalized manner and at a later
stage of a student career argue that students may not be willing or
able to devote the blocks of time required to integrate the skills that
both the analytical and numerical techniques demand within the
same course. While I accept the premise that my adopted approach
certainly requires a higher level of motivation and time commit-
ment by students, it is also true that students who were willing to
invest the required time in a course based on this text ended up
ready to plunge in more advanced engineering courses much ear-
lier than others. I am particularly gratified by the performance of
those of my students who, following the successful completion of
this text went on to advanced engineering courses and were able to
concentrate and excel in the new engineering material covered
7425_C000.qxd 6/9/06 16:53 Page xii

there while many of their classmates were bogged down acquiring

the computational tools which those selected had already mastered.
• There is still no consensus in my own academic circle as to whether
it is advisable to let students explore through computer experi-
ments new mathematical results prior (or even much prior) to the
formal introduction of the underlying theoretical principles in a
more structured environment. My own observations are that a stu-
dent’s intuitive skills are strongly enhanced by accumulated math-
ematical empirical discovery and experience; however, it will be
probably up to psychologists working in learning theories to come
up with a more definitive answer to the value of an empirical
approach in teaching applied mathematics.
• Teachers who adopted this text successfully informed me that they
incorporated a certain amount of flexibility in different students’
assignments to reflect the different levels of skills present in a nor-
mal student population distribution. My own experience was not
different.

My advice to beginner students of engineering using this text is to take

their time and finish each chapter with its solved examples and suggested
problems before moving on to the next chapter. I do not recommend to the
beginner to merely skim through the core material. The reader will quickly
discover that many computational skills are acquired and internalized best
when developed slowly and cumulatively.
The use of this text as reference material by a number of practicing scien-
tists and engineers is often brought to my attention by some pleasant com-
ments in a faculty meeting, professional conference, or social function. This
encouragement by peers is most appreciated.
This new edition benefited from comments, emails, and letters that I
received from many students, readers, and colleagues. I am grateful for all
those who took the time to make suggestions, point out misprints, and/or
argued the pedagogical approach that I adopted in presenting certain mate-
rial. I owe special thanks to Herman Cummins and Robert Brenart for their
valuable comments on this second edition’s book proposal, and to Allison
Taub, my editor at Taylor & Francis who took the initiative in actively and
enthusiastically championing the publication of a second edition of this text,
and in supporting its preparation at this time. Courtney Esposito and her
colleagues at The MathWorks Inc. continue to provide me with copies of the
most current releases of MATLAB and their supporting documentation for
which I am most grateful. Julie Spadaro, project editor at Taylor & Francis,
assured a quality final product, successfully coordinating production teams
on different continents.
My wife Azza and daughters Tala and Nigh urged me to prepare a second
edition of this book because partially, I suspect, they liked the fact that the
mail folder that I bring home with me every day now had many communi-
cations (even fan letters!) from other than the dozen or so colleagues, friends,
7425_C000.qxd 6/9/06 16:53 Page xiii

and editors who regularly read, comment, and respond to my research

papers. But above all, I know that deep in their hearts they supported me
unconditionally in what they knew was to me a labor of love. To them I
rededicate this book.
Jamal T. Manassah
New York, September 2006
7425_C000.qxd 6/9/06 16:53 Page xiv
7425_C000.qxd 6/9/06 16:53 Page xv

Contents

1. Starting with MATLAB and Exploring Its Graphics

Capabilities ....................................................................................................1
1.1 First Steps ..........................................................................................1
1.2 Basic Algebraic Operations and Functions ..................................2
1.3 Plotting Points ..................................................................................6
1.3.1 Axes Commands ..............................................................8
1.3.2 Labeling a Graph ............................................................9
1.3.3 Plotting a Point in 3-D ....................................................9
1.4 M-files ..............................................................................................10
1.5 MATLAB Simple Programming ..................................................11
1.5.1 The for Iterative Loop ................................................11
1.5.2 if-else-end Structures ............................................12
1.5.3 The while Loop ............................................................15
1.6 Arrays ..............................................................................................16
1.6.1 Array Relational Operations ........................................17
1.6.2 Array Algebraic Operations ........................................18
1.6.3 Combining Arrays Relational and Algebraic
Operations: Alternative Syntax to the if Statement ..19
1.6.4 Plotting Arrays ..............................................................20
1.7 Data Analysis ..................................................................................22
1.7.1 Manipulation of Data ....................................................22
1.7.2 Displaying Data ............................................................24
1.7.3 Normal Distribution ......................................................25
1.8 Parametric Equations ....................................................................27
1.8.1 Definition ........................................................................27
1.8.2 More Examples ..............................................................28
1.8.3 Oscilloscope Graphics ..................................................30
1.9 Polar Plots ......................................................................................32
1.10 3-D Plotting ....................................................................................34
1.10.1 Straight-Edge Geometric Figures ................................34
1.10.2 Parametric Equations for a 3-D Curve ......................38
1.10.3 Plotting a 3-D Surface ..................................................40
1.10.4 Contour: A Powerful Tool for Exploring 2-D
Geometries ......................................................................43
1.11 Animation ......................................................................................48
1.12 Specialized Plots: Velocity, Gradient, etc. ..................................49
1.12.1 Velocity Plots from Parametric Equations ................49
1.12.2 Gradient of a Potential ..................................................50
7425_C000.qxd 6/9/06 16:53 Page xvi

1.13 Printing and Saving Work in MATLAB ....................................52

1.13.1 Printing a Figure ............................................................52
1.13.2 Printing a Program File (Script M-file) ........................54
1.13.3 Converting a MATLAB Graphic into an
Image File ........................................................................54
1.13.4 Saving Values of Variables ............................................54
1.13.5 Loading the Value of a Variable ..................................54
1.13.6 Saving a MATLAB Session ..........................................54
1.13.7 Loading workspace.mat ..........................................54
1.14 MATLAB Commands Review ....................................................55

2. Difference Equations ................................................................................59

2.1 Simple Linear Forms ....................................................................59
2.2 Amortization ..................................................................................61
2.3 An Iterative Geometric Construct: The Koch Curve ................64
2.4 Solution of Linear Constant Coefficients
Difference Equations ....................................................................67
2.4.1 Homogeneous Solution ................................................68
2.4.2 Particular Solution ........................................................70
2.4.3 General Solution ............................................................70
2.5 Convolution-Summation of a First-Order System with
Constant Coefficients ....................................................................74
2.6 General First-Order Linear Difference Equations ....................76
2.7 Nonlinear Difference Equations ..................................................78
2.7.1 Computing Irrational Numbers ..................................78
2.7.2 The Logistic Equation ..................................................79
2.8 Fractals and Computer Art ..........................................................81
2.8.1 Mira’s Model ..................................................................81
2.8.2 Hénon’s Model ..............................................................84
2.9 Generation of Special Functions from Their Recursion
Relations ..........................................................................................86

3. Elementary Functions and Some of Their Uses ..................................89

3.1 Function Files ................................................................................89
3.2 Examples with Affine Functions ................................................91
3.2.1 Application to a Simple Circuit ..................................91
3.2.2 Further Consideration of Figure 3.1 ............................94
3.3 Examples with Quadratic Functions ..........................................95
3.4 Examples with Polynomial Functions ........................................96
3.5 Examples with the Trigonometric Functions ............................99
3.6 Examples with the Logarithmic Function ................................100
3.6.1 Ideal Coaxial Capacitor ..............................................100
3.6.2 The Decibel Scale ........................................................101
3.6.3 Entropy ..........................................................................101
3.7 Examples with the Exponential Function ................................103
3.7.1 Application to a Simple RC Circuit ..........................103
7425_C000.qxd 6/9/06 16:53 Page xvii

3.8 Examples with the Hyperbolic Functions and

Their Inverses ..............................................................................105
3.8.1 Application to the Capacitance of
Two Parallel Wires ......................................................105
3.9 Commonly Used Signal Processing Functions ........................105
3.10 Animation of a Moving Rectangular Pulse ..............................111
3.11 Use of the Function Handle ........................................................112
3.12 MATLAB Commands Review ....................................................114

4. Differentiation, Integration, and Solutions of Ordinary

Differential Equations ............................................................................115
4.1 Limits of Indeterminate Forms ..................................................115
4.2 Derivative of a Function ............................................................117
4.3 Infinite Sums ................................................................................119
4.4 Numerical Integration ................................................................121
4.5 A Better Numerical Differentiator ............................................124
4.5.1 Application ..................................................................127
4.6 A Better Numerical Integrator: Simpson’s Rule ......................128
4.7 Numerical Solutions of Ordinary Differential Equations ......135
4.7.1 First-Order Iterator ......................................................137
4.7.2 Higher-Order Iterators: The Runge–Kutta
Method ..........................................................................142
4.7.3 MATLAB ODE Solvers ................................................146
4.8 Integral Equations ........................................................................152
4.9 MATLAB Commands Review ..................................................158

5. Root Solving and Optimization Methods ..........................................159

5.1 Finding the Real Roots of a Function of One Variable ..........159
5.1.1 Graphical Method ........................................................159
5.1.2 Numerical Methods ....................................................160
5.1.3 MATLAB fzero Built-in Function ..........................167
5.1.4 Application: Zeros of the Zero-Order
Bessel Function ............................................................168
5.2 Roots of a Polynomial ................................................................170
5.3 Optimization Methods for Functions of One Variable ..........171
5.3.1 Graphical Method ........................................................172
5.3.2 Numerical Method: The Golden Section Method ..173
5.3.3 MATLAB fminbnd Built-in Function ......................175
5.4 The Zeros and the Minima of Functions in Two Variables ..178
5.4.1 The MATLAB fsolve Built-in Command ..............178
5.4.2 The MATLAB fminsearch Built-in Command ....181
5.5 Finding the Minima of Functions with
Constraints Present ......................................................................185
5.5.1 Lagrange Multipliers ..................................................186
5.5.2 MATLAB fmincon Built-in Function ......................187
5.6 MATLAB Commands Review ..................................................189
7425_C000.qxd 6/9/06 16:53 Page xviii

6. Complex Numbers....................................................................................191
6.1 Introduction ..................................................................................191
6.2 The Basics ......................................................................................191
6.2.1 Addition ........................................................................192
6.2.2 Multiplication by a Real or Imaginary
Number ........................................................................193
6.2.3 Multiplication of Two Complex Numbers ..............193
6.3 Complex Conjugation and Division ........................................195
6.3.1 Application to Division ..............................................196
6.4 Polar Form of Complex Numbers ............................................197
6.4.1 New Insights into Multiplication and Division
of Complex Numbers ..................................................198
6.4.2 Roots of Complex Numbers ......................................200
6.4.3 The Function y ⫽ e j␪ ....................................................201
6.5 Analytical Solutions of Constant Coefficients ODE ..............204
6.5.1 Transient Solutions ......................................................205
6.5.2 Solution in the Presence of a Source:
Green Function Technique ..........................................207
6.5.3 Steady-State Solutions ................................................210
6.5.4 Applications to Circuit Analysis ................................211
6.6 Phasors ..........................................................................................213
6.6.1 Phasor of Two Added Signals ....................................213
6.6.2 Total Phasor of Many Signals ....................................214
6.7 Interference and Diffraction of Electromagnetic Waves ........216
6.7.1 The Electromagnetic Wave ........................................216
6.7.2 Addition of Two Electromagnetic Waves ................217
6.7.3 Generalization to N-waves ........................................218
6.8 Solving ac Circuits with Phasors: The Impedance
Method ..........................................................................................220
6.8.1 RLC Circuit Phasor Analysis ......................................221
6.8.2 The Infinite LC Ladder ................................................222
6.9 Transfer Function for a Difference Equation with
Constant Coefficients ..................................................................224
6.10 MATLAB Commands Review ..................................................234

7. Vectors ........................................................................................................235
7.1 Vectors in Two Dimensions ........................................................235
7.1.1 Addition ........................................................................235
7.1.2 Multiplication of a Vector by a Real Number ..........235
7.1.3 Cartesian Representation ............................................236
7.1.4 MATLAB Representation of Vectors ........................238
7.2 Dot (or Scalar) Product ..............................................................239
7.2.1 MATLAB Representation of the Dot Product ..........241
7.3 Components, Direction Cosines, and Projections ..................243
7.3.1 Components ..................................................................243
7.3.2 Direction Cosines ........................................................243
7.3.3 Projections ....................................................................244
7425_C000.qxd 6/9/06 16:53 Page xix

7.4 The Dirac Notation and Some General Theorems ..................244

7.4.1 Cauchy–Schwartz Inequality ....................................246
7.4.2 Triangle Inequality ......................................................248
7.5 Cross Product and Scalar Triple Product ................................249
7.5.1 Cross Product ..............................................................249
7.5.2 Geometric Interpretation of the Cross Product ......250
7.5.3 Scalar Triple Product ..................................................250
7.6 Tangent, Normal, and Curvature ..............................................253
7.7 Velocity and Acceleration Vectors in Polar Coordinates ........255
7.8 Line Integral ..................................................................................261
7.9 Infinite-Dimensional Vector Spaces ..........................................263
7.10 MATLAB Commands Review ..................................................272

8. Matrices ......................................................................................................273
8.1 Setting Up Matrices ....................................................................273
8.1.1 Creating a Matrix by Keying in the Elements ........273
8.1.2 Retrieving Special Matrices from the
MATLAB Library ........................................................275
8.1.3 Functional Construction of Matrices ........................277
8.2 Adding Matrices ..........................................................................280
8.3 Multiplying a Matrix by a Scalar ..............................................280
8.4 Multiplying Matrices ..................................................................281
8.5 Inverse of a Matrix ......................................................................282
8.6 Solving a System of Linear Equations ......................................285
8.7 Application of Matrix Methods ................................................288
8.7.1 dc Circuit Analysis ......................................................288
8.7.2 dc Circuit Design ........................................................289
8.7.3 ac Circuit Analysis ......................................................290
8.7.4 Accuracy of a Truncated Taylor Series ....................291
8.7.5 Reconstructing a Function from Its Fourier
Components ..................................................................294
8.7.6 Interpolating the Coefficients of an (n ⫺ 1)-
Degree Polynomial from n Points ............................296
8.7.7 Least-Squares Fit of Data ............................................297
8.7.8 Numerical Solution of Fredholm Equations ............298
8.8 Eigenvalues and Eigenvectors ..................................................299
8.8.1 Finding the Eigenvalues of a Matrix ........................299
8.8.2 Finding the Eigenvalues and Eigenvectors
Using MATLAB ............................................................302
8.9 The Cayley–Hamilton and Other Analytical Techniques ......305
8.9.1 Cayley–Hamilton Theorem ........................................305
dX
8.9.2 Solution of Equations of the Form ⫽ AX ..........306
dt
8.9.3 Solution of Equations of the Form
dX
= AX ⫹ B(t)..............................................................308
dt
7425_C000.qxd 6/9/06 16:53 Page xx

8.9.4 Pauli Spinors ................................................................312

8.10 Special Classes of Matrices ........................................................319
8.10.1 Hermitian Matrices ......................................................319
8.10.2 Unitary Matrices ..........................................................323
8.10.3 Unimodular Matrices ..................................................324
8.11 Transfer Matrices ..........................................................................330
8.12 Covariance Matrices ....................................................................334
8.12.1 Parametric Estimation ................................................336
8.12.2 Karhunen–Loeve Transform ......................................339
8.13 MATLAB Commands Review ..................................................341

9. Transformations ........................................................................................343
9.1 Two-Dimensional Geometric Transformations ......................343
9.1.1 Construction of Polygonal Figures ............................343
9.1.2 Inversion about the Origin and Reflection
about the Coordinate Axes ........................................344
9.1.3 Rotation around the Origin ........................................345
9.1.4 Scaling ............................................................................346
9.1.5 Translation ....................................................................347
9.2 Homogeneous Coordinates ........................................................347
9.3 Manipulation of 2-D Images ......................................................351
9.3.1 Geometrical Manipulation of Images ......................351
9.3.2 Digital Image Processing ............................................352
9.3.3 Encrypting an Image ..................................................353
9.4 Lorentz Transformation ..............................................................355
9.4.1 Space–Time Coordinates ............................................355
9.4.2 Addition Theorem for Velocities ..............................357
9.5 Iterative Constructs ....................................................................358
9.5.1 The Koch Curve ..........................................................358
9.5.2 The Serpenski Curve ..................................................361
9.6 MATLAB Commands Review ..................................................363

10. A Taste of Probability Theory ................................................................365

10.1 Introduction ..................................................................................365
10.2 Basics ..............................................................................................366
10.3 Addition Laws for Probabilities ................................................371
10.4 Conditional Probability ..............................................................375
10.4.1 Total Probability Theorem ..........................................377
10.4.2 Bayes Theorem ............................................................377
10.5 Repeated Trials ............................................................................380
10.6 Generalization of Bernoulli Trials ............................................382
10.7 The Poisson and the Normal Distributions ............................382
10.7.1 The Poisson Distribution ............................................383
10.7.2 The Normal Distribution ............................................385
7425_C000.qxd 6/9/06 16:53 Page xxi

Appendix A: Review of Elementary Functions............................................389

A.1 Affine Functions ..........................................................................389
A.2 Quadratic Functions ....................................................................390
A.3 Polynomial Functions ..................................................................395
A.4 Trigonometric Functions ............................................................395
A.5 Inverse Trigonometric Functions ..............................................397
A.6 The Natural Logarithmic Function ..........................................397
A.7 The Exponential Function ..........................................................398
A.8 The Hyperbolic Functions ..........................................................399
A.9 The Inverse Hyperbolic Functions ............................................401

Appendix B: Determinants ..............................................................................403

Appendix C: Symbolic Calculations with MATLAB ..................................407

C.1 Symbolic Manipulation ..............................................................408
C.1.1 Creating Symbolic Expressions ................................408
C.1.2 Algebraic Manipulation ..............................................411
C.1.3 Plotting Symbolic Expressions ..................................414
C.2 Symbolic Solution of Algebraic and Transcendental
Equations ......................................................................................417
C.3 Symbolic Calculus ......................................................................418
C.4 Symbolic Linear Algebra ............................................................423
C.5 z-Transform and Laplace Transform ........................................428
C.5.1 z-Transform ..................................................................428
C.5.2 Solving Constant Coefficients Linear Difference
Equations Using z-Transform ....................................431
C.5.3 Laplace Transform ......................................................432
C.5.4 Solving Constant Coefficients Linear ODE
Using Laplace Transform ............................................436

Appendix D: Some Useful Formulae ............................................................439

Appendix E: Text Formatting ..........................................................................443

Selected References ..........................................................................................445

Index......................................................................................................................447
7425_C000.qxd 6/9/06 16:53 Page xxii
7425_C001.qxd 9/2/06 12:03 PM Page 1

1
Starting with MATLAB and Exploring Its
Graphics Capabilities

MATLAB can be thought of as a library of programs that will prove very use-
ful in solving many electrical engineering computational problems. MATLAB
is an ideal tool for numerically assisting you in obtaining answers, which is a
major goal of engineering analysis and design. This program is very useful in
research on circuit analysis, device design, signal processing, filter design,
control system analysis, antenna design, microwave engineering, photonics
engineering, computer engineering, and every other subfields of electrical
engineering. It is also a powerful graphic and visualization tool of beneficial
use for both quantitative courses and in your future career. In this chapter, we
start our exploration of MATLAB; in subsequent chapters, we will add to our
knowledge about this powerful tool of the modern engineer.

1.1 First Steps

The first step in using MATLAB is to know how to call it. It is important to
remember that although the front-end and the interfacing for machines with
different operating systems are sometimes different, once you are inside
MATLAB, all programs and routines are written in the same manner. Only
those few commands that are for file management and for interfacing with
external devices such as printers may be different for different operating sys-
tems. After entering MATLAB, usually by double clicking on the MATLAB
program icon, you should see the prompt >>, which means the program
interpreter is waiting for you to enter instructions. The interpreter goes to
work when you press the Return key. In case you wish to go to a new line
for the purpose of continuing your instructions, you should enter ellipses
( … ) before you hit the Return key.
In case the program is already opened, and to make sure that nothing is
saved from a previous session, type and enter:

>>clear all
1
7425_C001.qxd 9/2/06 12:03 PM Page 2

2 Elementary Mathematical and Computational Tools Using MATLAB®

This command removes all variables from the workspace. Also type and
enter (we shall henceforth say just enter):

>>clf

This command creates a graph window (if one does not already exist) or
clears an existing graph window. The command figure creates a new fig-
ure window, leaving the existing ones unaltered.
The command quit stops MATLAB.
Because it is not the purpose of this text to explain the function of every
MATLAB command, how would you get information on a particular syntax?
The MATLAB program has extensive help documentation available with
simple commands. For example, if you wanted help on a function called
roots (we will use this function often), you would type help roots.
Note that the help facility cross-references other functions that may have
related uses. This requires that you know the function name. If you want an
idea of the available help files in MATLAB, type help. This gives you a list
of topics included in MATLAB.
To get help on a particular topic such as elementary matrices and matrix
manipulation, type matlab\elmat. This gives you a list of all relevant
functions pertaining to that area. Now you may type help for any function
listed there. For example, try help size.

1.2 Basic Algebraic Operations and Functions

The MATLAB environment can be used, on the most elementary level, as a
tool to perform simple algebraic manipulations and function evaluations.

EXAMPLE 1.1

Exploring the calculator functions of MATLAB: The purpose of this example

is to show how to manually enter data and how to use the basic algebraic
operations of MATLAB. Note that the statements will be executed immedi-
ately after they are typed and entered (no equal sign is required).
Enter the text that follows the >> prompt to find out the MATLAB
responses to the following:

>> format short

>> 2.54376+2.32e1 %2.32e1=2.32*(10^1)
ans =
25.7438
>> 5.45^2
ans =
7425_C001.qxd 9/2/06 12:03 PM Page 3

Starting with MATLAB® and Exploring Its Graphics Capabilities 3

29.7025
>> ceil(7.58)
ans =
8
>> 2*sin(pi/4)
ans =
1.4142

The % symbol is used so that one can type comments in a program. (Comments
following the % symbol are ignored by the MATLAB interpreter.) Then enter:

>> format short e

>> 2.54376+2.32e1
ans =
2.5744e+001
>> 5.45^2
ans =
2.9703e+001
>> ceil(7.58)
ans =
8
>> 2*sin(pi/4)
ans =
1.4142e+000 ⵧ

The last command gave twice the sine of
/4. Note that the argument of the
function was enclosed in parentheses directly following the name of the
function. Therefore, if you wanted to find sin3(x), the proper MATLAB syn-
tax would be

>>sin(pi/4)^3

The numeric functions of MATLAB are:

ceil(x) rounds the number to nearest integer toward 
fix(x) rounds the number to nearest integer toward 0
floor(x) rounds the number to nearest integer toward 
round(x) rounds the number to nearest integer

The numeric formats of MATLAB are:

format short (default value) gives the number with four decimal
digits
format long gives the number with 16 digits
format short e gives the number with four decimals plus an exponent
format long e gives the number with 15 digits plus an exponent
format bank gives the number with 2 digits after the decimal (use-
ful for financial calculations)
7425_C001.qxd 9/2/06 12:03 PM Page 4

4 Elementary Mathematical and Computational Tools Using MATLAB®

To facilitate its widespread use, MATLAB has all the standard elementary
mathematical functions as built-in functions. Type help elfun to obtain a
listing of these functions. Remember that this is just a subset of the available
functions in the MATLAB library.

>>help elfun

The response to the last command will give you a list of these elementary
functions, some of which may be new to you, but all of which will be used
in your future engineering studies. Many of these functions will also be used
in later chapters of this book.
In the following table, the MATLAB expressions for the most common
mathematical functions are given:

MATLAB form Name Mathematical form

sqrt(x) squareroot 兹x
exp(x) exponential exp(x)  ex
log(x) natural logarithm ln(x)
log10(x) base10 logarithm log(x)  (ln(x)/ln(10))
sin(x) sine sin(x)
cos(x) cosine cos(x)
sec(x) secant sec(x)
tan(x) tangent tan(x)  (sin(x)/cos(x))
cot(x) cotangent cot(x)  (1/tan(x))
asin(x) inverse sine sin1(x)  arcsin(x)
acos(x) inverse cosine cos1(x)
asec(x) inverse secant sec1(x)
atan(x) inverse tangent tan1(x)
acot(x) inverse cotangent cot1(x)
sinh(x) hyperbolic sine sinh(x)  (ex  ex)/2
cosh(x) hyperbolic cosine cosh(x)  (ex  ex)/2
sech(x) hyperbolic secant sech(x)  1/cosh(x)
tanh(x) hyperbolic tangent tanh(x)  sinh(x)/cosh(x)
coth(x) hyperbolic cotangent coth(x)  1/tanh(x)
asinh(x) inverse hyperbolic sine (
sinh1 ( x)  ln x  x 2  1 )
where    x  

acosh(x) inverse hyperbolic cosine (

cosh1 ( x)  ln x  x 2  1 )
where x 1
1 1 
asech(x) inverse hyperbolic secant sech1 ( x)  ln   2  1 
x x 
where 0  x  1
1  1 x 
atanh(x) inverse hyperbolic tangent tanh1 ( x)  ln  
2  1 x 
where 1  x  1

1  x 1 
acoth(x) inverse hyperbolic cotangent coth1 ( x) 
ln  
2  x 1 
where x  1 or x  1
7425_C001.qxd 9/2/06 12:03 PM Page 5

Starting with MATLAB® and Exploring Its Graphics Capabilities 5

EXAMPLE 1.2

Assigning and calling values of parameters: In addition to inputting data

directly to the screen, you can assign a symbolic constant or constants to rep-
resent data and perform manipulations on them.
For example, enter and note the answer to each of the following:

>> format short

>> a=2.65;
>> b=3.765;
>> c=a*(a+b)
c =
16.9997
>> d=a*b
d =
9.9772
>> e=a/b
e =
0.7039
>> f=a^3/b^2
f =
1.3128
>> g=a+3*b^2
g =
45.1757

(The spacing shown above is obtained by selecting Compact for Numerical

Display in the Command Window Preferences in the File pull-down menu.)
ⵧ

If we desire at any point to know the list of the current variables in a MATLAB
session, we enter:

>>who

Note on variable names: Variable names must begin with a letter, contains
letters, digits, and underscore characters, and must contain less than 32 char-
acters. MATLAB is case-sensitive.

Question: From the above, can you deduce the order in which MATLAB
performs the basic algebraic operations?

Answer: The order of precedence is as follows:

1. Parentheses, starting with the innermost
2. Exponentiation from left to right
7425_C001.qxd 9/2/06 12:03 PM Page 6

6 Elementary Mathematical and Computational Tools Using MATLAB®

3. Multiplication and division with equal precedence, also evaluated

from left to right
4. Addition and subtraction, with equal precedence, also evaluated
from left to right

IN-CLASS EXERCISE

Pb. 1.1 Using the values, a  0.875 and b  1.5786, find the values of the fol-
lowing expressions (give your answers in the engineering short format):

h  sin(a) sin3/2(b)
j  a1/3 b3/7
k  (sin1(a/b))2/3
l  sinh(a) ab

1.3 Plotting Points

In this section, you will learn how to use some simple MATLAB graphics
commands to plot points. We use these graphics commands later in the text
for plotting arrays and for visualizing their properties. To view all the func-
tions connected with two-dimensional (2-D) graphics, type

>>help plot

All graphics functions connected with three-dimensional curves graphics

can be looked up by typing

>>help plot3

A point P in the x–y plane is specified by two coordinates: the x-coordinate

measures the horizontal distance of the point from the y-axis, while the
y-coordinate measures the vertical distance above or below the x-axis. These
coordinates are called Cartesian coordinates, and any point in the plane can
be described in this manner. We write for the point, P(x, y).
Other representations can also be used to locate a point with respect to a
particular set of axes. For example, in the polar representation, the point is
specified by an r-coordinate that measures the distance of the point from the
origin, while the -coordinate measures the angle which the line passing
through the origin and this point makes with the positive x-axis. The angle
is measured anticlockwise from the positive x-axis.
The purpose of the following two examples is to learn how to represent
points in a Cartesian plane and to plot them using MATLAB.
7425_C001.qxd 9/2/06 12:03 PM Page 7

Starting with MATLAB® and Exploring Its Graphics Capabilities 7

EXAMPLE 1.3

Plot the point P(3, 4). Mark this point with a red asterisk.
Solution: Enter the following:
>> x1=3;
y1=4;
plot(x1,y1,'*r') ⵧ
Note that the semicolon is used in the above commands to suppress the
echoing of the values of the inputs. The '*r' is used to mark the location and
the color of the point that we are plotting.
The symbols for the most common markers and colors for display in
MATLAB plot-functions are:

Markers Symbol Color Symbol

Circle o Blue b
Cross x Cyan c
Diamond d Green g
Dot ⴢ Magenta m
Pentagram p Yellow y
Plus sign ⴙ Red r
Square s White w
Star * Black k

EXAMPLE 1.4

Plot the second point, R(2.5, 4), on the graph while keeping point P of the
previous example on the same graph. Mark the new point with a small green
circle.

Solution: If we went ahead, defined the coordinates of R, and attempted to

plot the point R through the following commands:
>>x2=2.5;
y2=4;
plot(x2,y2,'og')
We would find that the last plot command erases the previous plot output.
Thus, what should we do if we want both points plotted on the same graph?
The answer is to use the hold on command after the first plot.
The following illustrates the instructions that you should have entered
instead of entering the above:
>>hold on
x2=2.5;
y2=4;
plot(x2,y2,'og')
hold off
The hold off turns off the hold on feature. ⵧ
7425_C001.qxd 9/2/06 12:03 PM Page 8

8 Elementary Mathematical and Computational Tools Using MATLAB®

NOTES
1. There is no limit to the number of plot commands you can type
before the hold is turned off.
2. An alternative method for viewing multiple points on the same
graph is available: we may instead, following the entering of the
values of x1, y1, x2, y2, enter:

>>plot(x1,y1,'*r',x2,y2,'og')

This has the advantage, in MATLAB, if no color is specified, that the

program assigns automatically a different color to each point.

1.3.1 Axes Commands

You may have noticed that MATLAB automatically adjusts the scale on a
graph to accommodate the coordinates of the points being plotted. The axis
scaling can be manually enforced by using the command

>>axis([xmin xmax ymin ymax])

Make sure that the minimum axis value is less than the maximum axis value
or an error will result.
In addition to being able to adjust the scale of a graph, you can also change
the aspect ratio of the graphics window. This is useful when you wish to see
the correct x to y scaling. For example, without this command, a circle will look
more like an ellipse, even if we have chosen ymax  ymin  xmax  xmin.

EXAMPLE 1.5

Plot the vertices of the square, formed by the points (1, 1), (1, 1),
(1, 1), (1, 1), keeping the same aspect ratio as in a graph paper, i.e., 1.

Solution: Enter the following:

>>x1=-1;y1=-1;
>>x2=1;y2=-1;
>>x3=-1;y3=1;
>>x4=1;y4=1;
>>plot(x1,y1,'o',x2,y2,'o',x3,y3,'o',x4,y4,'o')
>>axis([-2 2 -2 2])
>>axis square %makes the aspect ratio 1 ⵧ

Note that prior to the axis square command, the square looked like a rect-
angle. If you want to go back to the default aspect ratio, type axis normal.
In specific instances, you may need to choose other than the default
tick marks chosen by MATLAB, you can achieve this by adding the
7425_C001.qxd 9/2/06 12:03 PM Page 9

Starting with MATLAB® and Exploring Its Graphics Capabilities 9

following command:

set(gca,'Xtick',[xmin:dx:xmax],'Ytick',[ymin:dy:ymax])

Here gca stands for get current axes, and dx and dy are respectively the spac-
ings between tick marks that you desire in the x- and y-directions.
You can also make the graph easier to read, by using the command:

>>grid

This command displays gridlines at the tick marks.

If we repeat Example 1.5 with these new features added, the set of instruc-
tions will now read:

>>x1=-1;y1=-1;x2=1;y2=-1;x3=-1;y3=1;x4=1;y4=1;
>>plot(x1,y1,'o',x2,y2,'o',x3,y3,'o',x4,y4,'o')
>>axis([-2 2 -2 2])
>>axis square
>>set(gca,'Xtick',[-2:0.5:2],'Ytick',[-2:0.5:2])
>>grid

The tick marks were chosen in this instance to be separated by –12 units in both
the x- and y-directions.

1.3.2 Labeling a Graph

To add labels to your graph, the functions xlabel, ylabel, and title can
be used as follows:

>>xlabel('x-axis')
>>ylabel('y-axis')
>>title('points in a plane')

If you desire to also add a caption anywhere in the graph, you can use the
MATLAB command

>>gtext('caption')

and place it at the location on the graph of your choice, by clicking the mouse
on the desired location when the crosshair is properly centered there.
Different fonts and symbols can be used in labeling MATLAB graphs.
LaTeX symbols are used by MATLAB for this purpose (see Appendix E).

1.3.3 Plotting a Point in 3-D

In addition to being able to plot points on a plane (2-D space), MATLAB is
also able to plot points in a three-dimensional space (3-D space). For this, we
use the plot3 function.
7425_C001.qxd 9/2/06 12:03 PM Page 10

10 Elementary Mathematical and Computational Tools Using MATLAB®

EXAMPLE 1.6

Plot the points P(1, 1, 1), Q(4, 5, 6) and R(2, 5, 3). Show a grid in your graph.
Solution: Enter the following commands:

>>x1=1;x2=4;x3=2;
>>y1=1;y2=5;y3=5;
>>z1=1;z2=6;z3=3;
>>plot3(x1,y1,z1,'o',x2,y2,z2,'*',x3,y3,z3,'d')
>>axis([0 6 0 6 0 6])
>>grid ⵧ

NOTE You can also plot multiple points in a 3-D space in exactly the same
way as you did on a plane. Axis adjustment can still be used, but the vector
input into the axis command must now have six entries, as follows:

axis([xmin xmax ymin ymax zmin zmax])

You can similarly label your 3-D figure using xlabel, ylabel, zlabel, and
title.
The grid command is also valid in conjunction with the plot3 command.

1.4 M-files
In the last section, we found that to complete a figure with a caption, we had
to enter several commands one by one in the command window. Typing
errors will be time-consuming to fix because if you are working in the com-
mand window, and if you make an error, you will need to retype all or part
of the program. Even if you do not make any mistakes (!), all of your work
may be lost if you inadvertently quit MATLAB and have not taken the nec-
essary steps to save the contents of the important program that you just fin-
ished developing. This will be time-consuming. Especially, if you are
simulating a process, and all that you want to change in successive runs are
the parameters of the problem. To preserve large sets of commands, you can
store them in a special type of file called an M-file.
MATLAB supports two types of M-files: script and function M-files. To hold
a large collection of commands, we use a script M-file. The function M-file is
discussed in Chapter 3. To make a script M-file, you need to open a file using
the built-in MATLAB Menu. First select New from the File menu. Then select
the M-file entry from the pull-down menu. After typing the M-file contents,
you need to save the file. For this purpose, use the save as command from the
File window. A field will pop up in which you can type in the name you have
chosen for this file. In the pop-up lower window, indicate that it is an .m file.
7425_C001.qxd 9/2/06 12:03 PM Page 11

Starting with MATLAB® and Exploring Its Graphics Capabilities 11

NOTES
1. Avoid naming a file by a mathematical abbreviation, the name of a
mathematical function, a MATLAB command, or a number.
2. To run your script M-file, just enter in the command window the
filename omitting the .m extension in the file name at its end at the
MATLAB prompt. To be able to access this file from the command
window, make sure that the folder in which you saved this file is in
the MATLAB Path. If this is not the case use the Set Path in the File
pull-down menu and follow the screen prompts.

EXAMPLE 1.7

For practice, go to the File pull-down menu, select New and M-file to create
the following file that you will name myfile_17.
Enter in the editor window, the following set of instructions:

x1=1;y1=.5;x2=2;y2=1.5;x3=3;y3=2;
plot(x1,y1,'o',x2,y2,'+',x3,y3,'*')
axis([0 4 0 4])
xlabel('xaxis')
ylabel('yaxis')
title('3points in a plane')

After creating and saving this file as myfile_17 in your work folder, go to the
MATLAB command window and enter myfile_17. MATLAB will execute the
instructions in the same sequence as the statements stored in that file. ⵧ

1.5 MATLAB Simple Programming

In this section, we introduce the programming flow control commands for,
if, and while.

1.5.1 The for Iterative Loop

The power of computers lies in their ability to perform a large number of
repetitive calculations. To do this without entering the value of a parameter
or variable each time that these are changed, all computer languages have
control structures that allow commands to be performed and controlled by
counter variables, and MATLAB is no different. For example, the MATLAB
for loop allows a statement or a group of statements to be repeated a pre-
scribed number (i.e., positive integer) of times.
7425_C001.qxd 9/2/06 12:03 PM Page 12

12 Elementary Mathematical and Computational Tools Using MATLAB®

EXAMPLE 1.8

Generate the square of the first ten integers.

Solution: Edit and execute the following script M-file:

for m=1:10
x(m)=m^2;
end ⵧ

In this case, the number of repetitions is controlled by the index variable m,

which takes on the values m=1 to 10 in intervals of 1. Therefore, 10 assign-
ments are made. What the above loop is doing is sequentially assigning the
different values of m^2 (i.e., m2) for each element of the “x-array.” An array
is just a data structure that can hold multiple entries. An array can be 1-D
such as in a vector or 2-D such as in a matrix. More will be said about vec-
tors and matrices in subsequent chapters. At this time, assume 1-D and 2-D
arrays as pigeonholes with numbers or ordered pair of numbers respectively
assigned to them.
To find the value of a particular slot of the array, such as slot 3, enter:

>>x(3)

To read all the values stored in the array, enter:

>>x

Question: What do you get if you enter m, following the execution of the
above program?
Answer: The value 10, which is the last value of m read by the counter.

Homework Problem
Pb. 1.2 A couple establishes a college savings account for their newly
born baby girl, deposit $2000 in this account, and commit themselves to
deposit an equal amount at each future birthday of the new born. They
chose to invest this money with an investment firm that guarantees a
minimum annual return of 5%. What will be the minimum balance of the
account when the young lady turns 18? (Include in your calculation the
amount deposited on her 18th birthday.)

1.5.2 if-else-end Structures

If a sequence of commands must be conditionally evaluated based on a rela-
tional test, the programming of this logical relationship is executed with
some variation of an if-else-end structure.
7425_C001.qxd 9/2/06 12:03 PM Page 13

Starting with MATLAB® and Exploring Its Graphics Capabilities 13

A. The simplest form of this structure is:

if expression
Commands evaluated if expression is True
else
Commands evaluated if expression is False
end

NOTES
1. The commands between the if and else statements are evaluated
for all elements in the expression which are true.
2. The conditional expression uses the Boolean logical symbols &
(AND), | (OR), and ⬃ (NOT) to connect different propositions.

EXAMPLE 1.9

Find for integer a, 0  a  10, the values of C, defined as follows:

 ab for a  5

C 3
 ab for a  5
2

and b  15.

Solution: Edit and execute the following script M-file:

>>b=15;
for a=1:10
if a>5
C(a)=a*b;
else
C(a)=(a*b)*(3/2);
end
end

Check that the values of C that you will obtain by typing C are:

22.5 45 67.5 90 112.50 90 105 120 135 150 ⵧ

B. When there are three or more alternatives for the conditional, the
if-else-end structure takes the form:
if expression 1
Commands 1 evaluated if expression 1 is True
elseif expression 2
Commands 2 evaluated if expression 2 is True
7425_C001.qxd 9/2/06 12:03 PM Page 14

14 Elementary Mathematical and Computational Tools Using MATLAB®

elseif expression 3
Commands 3 evaluated if expression 3 is True
...

....
else
Commands evaluated if no other expression is True
end

In this form, only the commands associated with the first True expression
encountered are evaluated, ensuing relational expressions are not tested.

EXAMPLE 1.10

Find for integers a, 1  a  15, the value of D defined by

 a3 for a  5

D   2a for 5  a  10

 a  7 for 10  a  15

Solution: Edit and execute the following script M-file:

for a=1:15
if a<=5
D(a)=a^3;
elseif a>5 & a<10
D(a)=2*a;
else
D(a)=7+a;
end
end
D ⵧ

Homework Problem
Pb. 1.3 For the values of integer a going from 1 to 10, find the values of
C such that

 a2 if a is even
C
 a
3
if a is odd

Use the stem command to graphically show C. (Hint: Look up in the help
file the function mod.)
7425_C001.qxd 9/2/06 12:03 PM Page 15

Starting with MATLAB® and Exploring Its Graphics Capabilities 15

1.5.3 The while Loop

The while loop is used when the iteration process has to terminate when a
priori specified condition is satisfied. The difference with the for command
is that there the number of iterations is specified in advance.
The syntax for the while loop is as follows:

Initial assignment
while relational condition
expression
end

NOTE Special care should be exercised when using the while loop, as in
some cases the looping may never stop.

EXAMPLE 1.11

The unknown quantity x has the form 2n, where n is an integer. We desire to
find the largest such number subject to the condition that x  75,345.
Solution: Using the above syntax of the while loop, we enter:

>>x=2;
while 2*x<75345
x=2*x;
end
x

The unknown variable is initially assigned the value 2, it has this value until
it encounters the statement x=2*x, which should be read that x(n  1) 
2*x(n), the value of x is then changed to xnew. Before each pass through the
loop, the value of x is checked to see whether x  75,345. If this condition is
satisfied, the next iteration is carried through; otherwise the looping is
stopped. The final result to our problem is the value of x before any more
loop is executed.
If the statement of the problem is changed, and we had asked instead the
question to list all x’s that satisfied the condition x  75,345. We would have
entered instead of the above the following instructions:

>>x=2;
while x<75345
disp(x)
x=2*x;
end

The disp (display) command displays the results of all iterations satisfying
the desired condition. ⵧ
7425_C001.qxd 9/2/06 12:03 PM Page 16

16 Elementary Mathematical and Computational Tools Using MATLAB®

Homework Problem
Pb. 1.4 Using while, find the largest number x which can be written in
the form

x  3n  2n  3

and such that n is an integer and x  35,724.

1.6 Arrays
As mentioned earlier, think of an array as a string of ordered numbers, i.e.,
to determine an array, we need to specify the value and position of each ele-
ment in this string.
As pointed out earlier, an element in the array is addressed by writing the
name of the array followed by the position of the element within parenthe-
ses. For example, x(3) will be the value of the third element in the array x.
The basic array functions are:

length(x) finds the length of the array x (i.e., number of elements)

find(x) computes a new array where the successive elements indicate
the positions where the components of the array x are nonzero
max(x) gives the value of the largest element in the array x
min(x) gives the value of the smallest element in the array x
sum(x) gives the sum of all elements in the array x

EXAMPLE 1.12

Enter the array x  [0 3 6 2 11 0 7 9], and find its third element, length,
nonzero elements, value of its largest and smallest elements, and sum of all
its elements.

Solution:

>> x=[0 3 6 -2 11 0 7 9];

>> x(3)
ans =
6
>> length(x)
ans =
8
>> find(x)
7425_C001.qxd 9/2/06 12:03 PM Page 17

Starting with MATLAB® and Exploring Its Graphics Capabilities 17

ans =
2 3 4 5 7 8
>> max(x)
ans =
11
>> min(x)
ans =
-2
>> sum(x)
ans =
34 ⵧ

1.6.1 Array Relational Operations

MATLAB has six relational operations that compare the elements of two
arrays of the same length.
The relational operators are:

= = equal to
~ = not equal to
< smaller than
<= smaller or equal to
> larger than
>= larger than or equal

If we combine within parentheses the symbols for two arrays with a rela-
tional operator, the result will be a new array consisting of 0’s and 1’s, where
the 1’s are in the positions where the relation between the two elements in
the same position from the two arrays is satisfied, and 0’s otherwise.

EXAMPLE 1.13

Study the following printout:

x = [1 4 7 5 11 2];
y = [3 -2 7 2 12 -4];
>> z1=(x==y)
z1 =
0 0 1 0 0 0
>> z2=(x~=y)
z2 =
1 1 0 1 1 1
>> z3=(x<y)
z3 =
1 0 0 0 1 0
>> z4=(x<=y)
7425_C001.qxd 9/2/06 12:03 PM Page 18

18 Elementary Mathematical and Computational Tools Using MATLAB®

z4 =
1 0 1 0 1 0
>> z5=(x>y)
z5 =
0 1 0 1 0 1
>> z6=(x>=y)
z6 =
0 1 1 1 0 1
>> z7=(z1==z2)
z7 =
0 0 0 0 0 0
>> z8=(z3==z6)
z8 =
0 0 0 0 0 0
>> z9=(z4==z5)
z9 =
0 0 0 0 0 0

Comments: The values of z7, z8, and z9 are identically zeros. This means that
no corresponding elements of z1 and z2, z3 and z6, z4 and z5, respectively are
ever equal, i.e., each pair represents mutually exclusive conditions. ⵧ

1.6.2 Array Algebraic Operations

In Section 1.5, we used for loops repeatedly. However, this kind of loop-
programming is very inefficient and must be avoided as much as possible in
MATLAB. In fact, ideally, a good MATLAB program will always minimize
the use of loops because MATLAB is an interpreted language — not a com-
piled one. As a result, any looping process is very inefficient. Nevertheless,
at times we use the for loops, when necessitated by pedagogical reasons.
To understand array operations more clearly, consider the following:
a=1:3 % a starts at 1, goes to 3 in increments of 1.
If the increment is not 1, you must specify the increment; for example,

b=2:0.2:6 % b starts at 2 ends at 6 in steps of 0.2

To distinguish arrays operations from either operations on scalars or on

matrices, the symbol for multiplication becomes .*, that of division ./,
and that of exponentiation .^. Thus, for example,

c=a.*b % takes every element of a and multiplies

% it by the element of b in the same array location

Similarly, for exponentiation and division:

d=a.^b
e=a./b
7425_C001.qxd 9/2/06 12:03 PM Page 19

Starting with MATLAB® and Exploring Its Graphics Capabilities 19

If you try to use the regular scalar operations symbols, you will get an error
message.
Note that array operations such as the above require the two arrays to
have the same length (i.e., the same number of elements). To verify that two
arrays have the same number of elements (dimension), use the length com-
mand. The exception to the equal length arrays rule is when one of the
arrays is a scalar constant.

NOTE The expression x=linspace(0,10,200) is also a generator for an

x-array. Its first element equal to 0, its last element equal to 10, and it has 200
equally spaced points between 0 and 10. Here, the number of points rather
than the increment is specified; that is, length(x)=200.

1.6.3 Combining Arrays Relational and Algebraic Operations: Alternative

Syntax to the if Statement
As an alternative to the if syntax, we can use a combination of arrays rela-
tional and algebraic operations to generate a complicated array from a sim-
ple one.

EXAMPLE 1.14

Using array operations, find for integers a, 1  a  15, the value of D defined by

 a4 for a  5

D   2a for 5  a  10

 a  7 for 10  a  15

Solution:
>>a=1:15;
D=(a.^4).*(a<=5)+2*a.*(5<a).*(a<10)+(a+7).*(10<=a)...
.*(a<=15) ⵧ
NOTE We have used in the above the property that the logical AND link-
ing two separate relational operations can be represented by the array prod-
uct of these relations.

Homework Problem
Pb. 1.5 The array n is all integers from 1 to 110. The array y is defined by

 n
1/ 2
if n divisible by 3, 5, or 7
y=
 0 otherwise

Find the value of the sum of the elements of y.

7425_C001.qxd 9/2/06 12:03 PM Page 20

20 Elementary Mathematical and Computational Tools Using MATLAB®

1.6.4 Plotting Arrays

If we are given two arrays x and y, we can plot the y-array as function of the
x-array. What MATLAB is actually doing is displaying all the points
P(xi, yi) and connecting each two consecutive points by a straight line.
The same options for markers and color that were available for plotting
points and described earlier are also valid in this case. Additionally, we can
specify here the style of the line that connects the different points.

Style MATLAB Symbol

Solid line –
Dashed line ––
Dash-dot line –.
Dotted line :

If the points are close together, the points will look connected by a continu-
ous line making a smooth curve; we say that the program graphically inter-
polated the discrete points into a continuous curve.
The commands for labeling, axis and tick marks used in plotting points are
also valid here. We can zoom to a particular region of the graph by using the
command zoom, And we can read off the coordinates of particular point(s)
in the graph, by using

[x,y]=ginput(n)

This command serves to read the values of the coordinates of points off a dis-
played graph (n is the number of points). The command pops out the
crosshair. The operator manually zeros on the points of interest and clicks. As
a result, MATLAB prints the coordinates of the points in the same order as
those of the clicks of the mouse.

EXAMPLE 1.15

The three arrays x, y, and z are defined as follows:

x  0:0.1:10; y  x 2  6x  20; z  3x  5

Plot y as function of x, and z as function of x. Show the traces as red dotted

line and blue solid line, respectively, label the x- and y-axis, and label graph
as Figure of Example 1.15.

Solution:
>>x=0:0.1:10;
y=x.^2-6*x-20;
z=2*x-5;
plot(x,y,'r:',x,z,'b-')
xlabel('x axis')
7425_C001.qxd 9/2/06 12:03 PM Page 21

Starting with MATLAB® and Exploring Its Graphics Capabilities 21

ylabel('y-z axis')
title('Figure of Example 1.15')
grid

In the above solution, we plotted both curves on the same graph. Often, it may
be required to plot the two curves in separate graphs but in the same figure
window. This can be achieved using the subplot command. The arguments
of the subplot command subplot(m,n,p) specify m the number of rows par-
titioning the graph, n the number of columns, and p the number of the partic-
ular subgraph chosen (enumerated through the left to right, top to bottom
convention). The program to plot the above arrays in two graphs in a column
would read

>>x=0:0.1:10;
y=x.^2-6*x-20;
z=2*x-5;
subplot(2,1,1)
plot(x,y,'r:')
xlabel('x axis')
ylabel('y axis')
title('Figure of Example 1.14-a')
subplot(2,1,2)
plot(x,z,'b-')
xlabel('x axis')
ylabel('z axis')
title('Figure of Example 1.14-b') ⵧ

All the plotting that we did so far has been using the linear scale in both the
horizontal and vertical axes. Sometimes when the elements of an array can
change values by several decades over the range of interest, it is well advised
to use the different kinds of the log-plots. As the need arises we can, instead
of the plot command, use one of the following alternatives:

semilogy creates a semilog plot with logarithmic scale on the y-axis

semilogx creates a semilog plot with logarithmic scale on the x-axis
loglog creates a log–log graph

Homework Problem
Pb. 1.6 The arrays y and z are generated from the array of positive inte-
gers x  1:20, such that:
y  2x and z  51/x
Plot the array z as function of the array y. (You should be able to read off
the value of the coordinates of any point on the plot over the entire range
of the graph.)
7425_C001.qxd 9/2/06 12:03 PM Page 22

22 Elementary Mathematical and Computational Tools Using MATLAB®

1.7 Data Analysis

1.7.1 Manipulation of Data
The most convenient representation for data collected from experiments is in
the form of histograms. Typically, you collect data and want to sort it out in
different bins. But prior to getting to this point, let us review and introduce
some array-related commands that are useful in data manipulation.
Let {yn} be a data set. It can be represented in MATLAB by an array y.

length(y) gives the length of the array (i.e., the number of data points)
y(i) gives the value of the element of y in the ith position
y(i:j) gives a new array, consisting of all elements of y between and
including the locations i and j
find(y) computes a new array where the successive elements indicate
the positions where the components of the array y are nonzero
max(y) gives the value of the largest element in the array
min(y) gives the value of the smallest element in the array
sort(y) sorts the elements of y in ascending order and gives a new
array of the same size as y
sum(y) sums the elements of y. The answer is a scalar
mean(y) gives the mean value of the elements of y
std(y) computes the standard deviation of the elements of y

The definitions of the mean and of the standard deviation are, respectively,
given by

N
1
y
N
∑ y(i)
i1

N
( y(i)  y )2
 ∑
i1 ( N  1)

where N is the length of the array.

EXAMPLE 1.16

An old recorder at Bill’s backyard in Hamilton Heights in New York down-

loaded the following raw data for the temperature readings on July 25
between 12:01 p.m. and 1:00 p.m.:

T=[95 95.5 96 95.8 97 96.5 96 95 55 96 95 97 96 ...

96.5 96.2 125 96 95 95.5 96]
Other documents randomly have
different content
l’autre côté du petit lac, le Palais lumineux flamboie, allumant
comme des lanternes japonaises dans l’interstice des feuillages. Du
jardin, des valses tziganes montent, dansent et pleurent, raclées par
un orchestre invisible. A une des fenêtres ouvertes, un des musiciens
se penche avec des contorsions de torse et des roulements de
prunelles qui semblent accompagner son rythme. «Vovos Elek, me
chuchote-t-on à l’oreille, le plus demandé des tziganes de Budapest
et un peu mieux que Rigo, n’est-ce pas?»
Et pas une femme à table pour nous glisser à l’oreille: «Donne
donc cinq louis au tzigane.» Il n’y a que des hommes: Dulong,
l’oseur de tant de jolies architectures du Champ de Mars; Soulié,
autre architecte, autre artiste; Lafitte, l’homme de tous les sports;
Chincholle, qui vient de porter un toast; Octave Uzanne, etc.
Onze heures. «—Branle-bas sur le pont: c’est l’heure où
Madeleine Lemaire s’habille en porte Binet pour son bal. —En porte
Binet! Non, cette chose monstrueuse et grotesque... —Oui, elle en
paraîtra coiffée, comme jadis du buste de Dumas, feu madame
Aubernon. —Madeleine Lemaire coiffée de la porte Binet! Ses pires
ennemis n’auraient pas trouvé cela!»

Mercredi 30 mai.—A l’Opéra-Comique, dix heures. «Hænsel et

Gretel». Le charme frais, l’émotion attendrie de l’éternelle légende
des enfants perdus dans la forêt, le plus fin joyau d’art peut-être de
toute cette année que cette représentation —dans quel cadre exquis
et voulu!— du conte théâtral de madame Adélaïde Wette et de M.
Humperdinck.
Et, dans des décors d’une réalité légendaire, intérieur de
chaumière et forêt de sapins, déjà rencontrés en Allemagne, aux
bords des lacs du Tyrol bavarois, ce sont les chansons, les danses,
les querelles, les extases, les angoisses et les épouvantes d’Hænsel
et de Gretel, fuyant la correction maternelle pour tomber aux mains
de l’horrible fée Grignotte, l’ancestrale ogresse de toutes les histoires
et de tous les contes de fées, l’abracadabrante «baba Yaga» des
légendes russes, si merveilleusement rendue avec ses gestes
hésitants, ses mines voraces et sa démarche cauteleuse par Mme
Delna, Delna, qui vient de trouver là une création en tous points
digne de son talent comique et l’unanime succès une fois déjà
rencontré dans «Falstaff».
Autour de sa silhouette fantômatique et grotesque de «maman
Lèchefrite», ce sont les ébats et les mines effarées des deux plus
jolies poupées de Nuremberg qu’aient jamais rêvées buveur de bière
et fumeur d’opium: Rioton-Gretel et de Craponne-Hænsel, toutes
deux si gosses, si démantibulées de gestes et d’une terreur si
touchante dans l’ombre crépusculaire et hantée de la forêt; et c’est,
silhouette exquise de rêve, Mlle Mastiana dans l’«Homme au sable»,
à côté de —vision matinale trop vite évanouie— Mlle Daffetye dans
l’«Homme à la rosée». Et alors des trouvailles de mise en scène où
se reconnaît le génie de Carré: l’escalier d’or où s’étage en deux
rangs d’ailes harmonieuses la descente des anges penchés sur le
sommeil des enfants; les mirages de l’étang, où l’horreur des
souches grimaçantes s’aggrave de phosphorescentes prunelles de
hiboux et la chevauchée à travers les sapins de l’ogresse, le fatidique
balai du sabbat entre les jambes.
A la partition colorée et chantante de M. Humperdinck, d’une
orchestration si substantielle et si savante, d’un charme populaire de
vieux lieds et de rondes célèbres en pays de Rhin, la traduction de
M. Catulle Mendès prête une poésie ailée et puérile qui restitue à
toute l’œuvre son atmosphère wagnérienne, son parfum de petite
fleur légendaire tombée de la couronne d’Elisabeth dans les «Maîtres
chanteurs».

Vendredi 1er juin.—Le Grand Bazar. Ce qu’ils en disent, ce qu’ils

en pensent, ce qu’il faut avoir vu: «—Les tapisseries tissées de fils
d’or au pavillon de l’Espagne. —Les Watteau de l’empereur au
pavillon de l’Allemagne. —Les Espagnols à la feria. —Les ciboires, les
ornements d’église et les vêtements d’apparat bossués d’émaux du
pavillon de la Hongrie. —Le mobilier de l’Angleterre. —Les faïences
du Danemark. —De l’extase pour la Finlande est bien portée.
Remarquez l’atmosphère de calme et d’honnêteté dégagée par les
peuples du Nord; venir se reposer dans l’ombre fraîche de leur vertu
de la luxure brutale de l’Orient, sujet de conversations et de
chroniques. —Divagations tout indiquées sur les clochetons et les
toits ajourés de la Suède et de la même Finlande; les opposer aux
moucharabiehs et minarets plâtreux de l’Algérie, vomir sur les
soukhs tunisiens, dithyramber encore, au besoin, sur le phare de
l’Allemagne, conspuer le style italien. —Au petit palais, n’admettre
que le dix-huitième siècle, reconnaître à première vue un Riesener
d’un Crescent, se pâmer sur l’enseigne de Watteau, la fameuse
enseigne, éviter de parler des «Trois Grâces», la trop clamée et
réclamée pendule de M. de Camondo (le sujet commence à être
rebattu), découvrir si possible d’autres Falconnet, affecter le plus
profond mépris pour les œuvres entassées dans le grand palais,
dédaigner la peinture contemporaine, être très documenté sur
l’exposition des voitures, carrosses de sacre, berlines de l’émigré,
coucous, pataches et diligences, connaître par le menu les noms des
propriétaires des anciennes vinaigrettes: pour éviter les gaffes, relire
au besoin le beau livre d’Octave Uzanne. —Ne pas se vanter de dîner
tous les soirs à l’Exposition: on vous croirait sans famille ou au ban
de vos relations; éviter de dire: «Nous en sommes à notre trentième
restaurant», pour ne pas s’entendre répondre: «Ça fait l’éloge de
votre estomac.»
Bal Madeleine Lemaire. Paraître très informé sur la genèse de
chaque costume, admirer sans restriction le fabuleux rajah bleu
turquoise de M. de Montesquiou, savoir que les bijoux étaient prêtés
par Sarah. —S’étonner de la mise en loge des Altesses et de la
complaisance des invités consentant à défiler la parade devant la
princesse Hélène et les grands-ducs, lancer d’un ton indigné: «Moi,
jamais je n’aurais voulu me donner en spectacle», pour risquer
aussitôt cette restriction perfide: «Bah! il y avait tous les théâtres».
Scandale du jour, fausse nouvelle «—Et ce mariage Pougy-
Lorrain? Qu’y a-t-il de vrai dans tout cela? —Le plaisir que Liane a à
le démentir. —Et Lorrain à le laisser croire. —Pourquoi? —Par
rosserie pour rappeler des mariages précédents —identiques— et
embêter le ménage... —Ne le nommez pas.»

Samedi 2 juin.—Le Grand Bazar. La merveille du Trocadéro et

peut-être de toute l’Exposition, le Pnom et la pagode des Bouddhas,
à l’Indo-Chine, derrière les pavillons de la Martinique et de la
Réunion.
Très haut dans les arbres, au-dessus des longues feuilles des
bananiers, les toits recourbés d’une pagode luisent, éperonnés d’or
et lambrissés, comme d’autant d’écailles, de petites tuiles de
nuances délicates. Ils luisent très haut dans le ciel, ces toits on dirait
laqués, gardés le long d’une vaste terrasse par vingt chimères, vingt
monstres hilares, moitié dogues et moitié requins, écartant tous une
croupe rebondie d’où s’élance une rigide queue verticale.
C’est la pagode des Bouddhas. Une admirable frise représentant
une guirlande de bayadères court le long de la terrasse; des idoles
mitrées veillent quatre par quatre, enclavées dans un pylône, aux
bords des parapets. Ce sont des déesses accroupies, la fleur de lotus
à la main, et leur solennité tranquille encadre, aux quatre coins du
temple, le rire immense et répété des monstres. Une énorme
coupole aux enroulements de turban hindou se renfle et s’effile
derrière les toits; elle couronne un vaste cube de pierre, percé de
trois ouvertures; quatre géants flanquent l’entablement du dôme et
le gardent, appuyés sur des massues. Architecture religieuse et
barbare, un peu menaçante et pourtant harmonieuse, tant elle est
voulue. Chaque angle, chaque ligne, chaque statue y renferme un
symbole, et tous les détails d’ornementation concourent, on le sent,
comme dans les cathédrales gothiques, à l’affirmation d’un dogme
ou d’une foi. C’est le lent et magnifique épanouissement d’un mythe
inscrit dans la pierre et le métal, un «credo» d’architecture où
chaque marche d’escalier, chaque statue a un sens mystique. Et une
grande admiration vous prend pour les peuples disparus (car la race
Khmer, dont ce monument atteste la puissance, est depuis
longtemps abolie), et une grande admiration, dis-je, vous prend pour
ces peuples lointains chez qui l’idée religieuse fut si forte qu’à travers
les siècles et les espaces la reconstitution rapetissée d’un de leurs
temples nous impose, à nous autres modernes, le respect d’une
ferveur et le regret d’une foi.
Un raide escalier y conduit par quarante marches, gardées par
dix monstres, les mêmes dogues hilares, dentés comme des requins:
ils sont là musclés et trapus, toutes les cinq marches, escortant de
leur grimace immobile la lente montée des visiteurs... autrefois, des
pèlerins. Deux géants appuyés sur des massues font sentinelle à la
porte du temple, à niveau de la terrasse... et, le long des lourds
parapets, à la pointe recourbée des toits, comme au faîte de grands
mâts plantés çà et là parallèles aux arbres, des clochettes
tintinnabulent avec des cliquetis de métal, mélancoliques et
mystiques sonneries qui, là-bas, dans les forêts de l’Inde, dénoncent
l’approche des lieux saints à la dévotion des fidèles comme à l’effroi
des parias.
De chaque côté de l’escalier, à niveau des soubassements du
temple, deux petits dômes enturbannés se bombent et
tirebouchonnent au-dessus de deux cubes de pierre: l’entrée des
souterrains.

Dimanche 3 juin.—Le Grand Bazar, neuf heures du soir.

Marchera-t-il? Ne marchera-t-il pas? C’est du château qu’il s’agit, le
château d’eau dont les frises illuminées devaient éclairer tout le
Champ de Mars pendant que, dans le cintre, des fontaines
jaillissantes, des retombées d’eau liquide et des cascades d’écume
devaient offrir à l’ébahissement des foules une apothéose
hydraulique incendiée de toutes les lueurs et de toutes les nuances
du prisme.
Des hauteurs du Trocadéro jusqu’au fond du Champ de Mars,
une marée humaine, un océan de têtes curieuses se tient figé, enlizé
par la masse même des groupes, dans l’attente du spectacle promis.
Cette foule! On n’y jetterait pas une aiguille! Les trains de plaisir de
la Pentecôte y ont versé, depuis le matin, de véritables caravanes.
Les journaux du lendemain publieront la statistique des entrées avec
les recettes fabuleuses du trottoir roulant: plus de six cent cinquante
mille visiteurs auront fait réaliser à la Compagnie du trottoir-clou de
l’Exposition, pour ce seul dimanche, quatre-vingt-dix mille francs de
recette. On est revenu aux plus beaux jours de 1889. Mais que de
papiers gras, que de débris de charcuterie et que de litres vides sur
les degrés des palais, les gazons des rares pelouses et les assises de
la tour Eiffel!
Oui, tous les trains de plaisir ont donné, tous les arrivages de
province et de banlieue et tous les faubourgs ouvriers aussi. Public
de kermesse, tout ce monde a dîné dehors, sur des bancs ou sur des
chaises, les autres assis par terre, à la bonne franquette, comme en
plein bois de Vincennes ou de Boulogne; personne d’entre ces
pèlerins de M. Picard n’est rentré dîner au logis ou à l’hôtel: tous ont
emporté le panier de provisions, avec le rond de saucisson et le
cervelas à l’ail obligatoires.
Aussi la cohue est-elle énorme et mal odorante. «Ça fouette»,
selon l’expression consacrée dans l’argot imagé de l’atelier. C’est,
mêlée à la senteur âcre de la poussière, l’odeur de paquebot mal
tenu spéciale à la foule, car tout ce monde a beaucoup marché
depuis l’aube. Mais, dans l’excitation de la joie de voir enfin
fonctionner le château d’eau et l’électricité, tous surmontent leur
lassitude, tous attendent, heureux, les yeux fixes et la bouche bée.
Marchera-t-il ou ne marchera-t-il pas? Tous les journaux, depuis
deux jours, annoncent pour ce soir... enfin, le fonctionnement
complet de cette huitième merveille. Des cris d’admiration anticipée
courent à travers les groupes. Quand ça va-t-il commencer? Tout
autour, l’embrasement des palais éclairés «a giorno» fait plus obscur
et plus noir le grand trou d’ombre où dort encore l’hypothétique
apothéose du château d’eau.
Eh bien, le château d’eau marche mal ou plutôt ne marche pas.
Les frises s’allument bien, rouges et vertes, comme envahies dans
leurs détails d’ornementation par des serpents versicolores; mais les
jaillissements d’eau et les cascades des bassins étagés du bas
s’obstinent à demeurer dans l’ombre. Cela ne s’éclaire que sous les
projections intermittentes de la tour Eiffel.
Parfois, comme des lueurs vertes sourdent au ras de terre, des
lividités vaporeuses s’échevèlent peut-être à un mètre au-dessus de
la foule; mais le motif central demeure noir. Et, tout à coup, toutes
les frises s’éteignent. Quand elles se rallument, les lueurs vertes des
jets d’eau s’évanouissent à leur tour: impossible d’obtenir un
ensemble. C’est une déception pour tous et pour toutes que ces
crissements d’eau et ces vagues blancheurs écumantes devinées
dans le clair-obscur.

Jeudi 7 juin.—Le Grand Bazar, au Petit Palais.

A droite en entrant, dans la claire et haute rotonde en lanterne
qui termine chaque galerie, parmi la lumière fine et pure des
grandes fenêtres Louis XVI ouvertes sur la travée du fleuve et la
verdure des jardins, trois bibelots merveilleux requièrent et
retiennent au seuil des longues salles, où triomphe l’art français du
dix-huitième siècle: un traîneau et deux chaises à porteurs.
Le traîneau: Une immense tortue d’eau s’écartèle, éperdue, sur
deux montants dorés et peints, à soixante centimètres du sol, une
énorme tortue brune à la tête et aux pattes dardées hors de sa
carapace dans la tension d’un violent effort. Elle soutient une conque
dorée du plus pur style Louis XV, où s’enchâssent un siège et des
coussins de vieux velours vert. La sellette du cocher s’érige en
arrière sur la queue dressée d’un dauphin; des branches de chêne
en or moulu courent le long de la coquille.
Aujourd’hui immobile et remisé dans une galerie de musée, ce
traîneau, qui dut jadis emporter sur la pièce d’eau des Suisses les
terreurs amusées de quelque favorite, toute de zibeline et de velours
sous le loup de satin qui protégeait du froid, raconte d’anciennes
splendeurs entre deux chaises à porteurs du temps, l’une de
cardinal, dont le panneau principal porte encore les insignes entre
des nudités de déesses envolées sur fond d’or; l’autre, bibelot choyé
de quelque petite maîtresse, encadre de guirlandes et de fins
rinceaux du plus pur style rocaille des marines on dirait de Claude
Lorrain, tant elles sont pompeuses et décoratives dans le vert de
leurs eaux et l’ambre de leurs ciels.
Mais ne cherchez pas à connaître la provenance de ces pièces
précieuses. Une déception cruelle attend au Petit Palais les
amoureux nostalgiques du passé: aucune autre indication n’existe
que le petit carton écrit à la main, mentionnant, au-dessus de
quelques-uns, qu’ils appartiennent à des marchands.
Le catalogue du Petit Palais n’existe pas... en français, mais tout
le monde peut lire en anglais la nomenclature détaillée des trésors
entassés là. Il fait bon d’être d’outre-Manche. «Les Anglais chez
eux», telle pourrait être la devise de ce temple de notre art
rétrospectif. Son catalogue en français ne paraîtra que dans
quelques jours. Or nous sommes le 7 juin; l’Exposition est ouverte le
15 avril.

Samedi 9 juin.—Pour Jean de Bonnefon. Le Saint-Siège en

coquetterie avec le ministère et le gouvernement.
A un des derniers dîners du comte d’Haussonville, en plein
orléanisme intransigeant et militant du Faubourg, le nonce se
répandait en éloges et en aménités sur le ministère actuel, se
félicitant des bons procédés de nos Excellences vis-à-vis Rome et la
papauté, vantant l’urbanité de Pierre et de Jacques et ne se
plaignant que d’une chose: le quartier affecté à la résidence
apostolique, ce quartier de Monceau, un peu trop mondain, un peu
trop élégant pour le caractère ecclésiastique d’un mandataire papal.
C’est le faubourg Saint-Germain qu’il eût fallu à la nonciature pour y
trouver un logement plus en rapport avec son caractère. Une
première mise de fonds de cinquante mille francs eût fait le reste.
C’était à la jeunesse catholique, à la noblesse française de se
remuer, de se saigner un peu pour assurer au représentant du Saint-
Père une résidence digne de lui. Et, comme l’assistance, un peu
refroidie par les précédents éloges du monsignor, accueillait sans
enthousiasme l’insinuation domiciliaire et son appel de fonds, notre
Italien, sans se démonter: «On dit M. de Rothschild très aimable,
très généreux et d’une cordialité parfaite quand on veut bien
s’adresser à lui; il est on ne peut mieux disposé pour nos œuvres et
l’intérêt de l’Eglise prime tout, n’est-il pas vrai? Comme la première
vertu chrétienne est l’esprit de sacrifice, peut-être devrais-je lui
rendre visite, aller trouver M. le baron?»
A quoi un des convives, un peu impatienté: «Vous y êtes déjà
allé, monseigneur!»

Dimanche 17 juin.—Le Turf à domicile. Il y a huit jours, c’était la

grande poussière, la grande chaleur et la grande cohue du Grand
Prix de Paris, six cent mille curieux de province et d’ailleurs rués à
Longchamps et mijotant au soleil sur le gazon de la pelouse; il y a
huit jours, c’étaient les ovations au roi de Suède, le triomphe de
Semendria et le jeu de massacre de la tribune officielle.
Aujourd’hui, dans l’ombre fraîche des persiennes closes, c’est la
contemplation muette, secouée de petits rires, de tout le pesage
évoquée par Sem dans son album le «Turf». Ils défilent tous, les
gros propriétaires et les entraîneurs, princes du paddock et rois de la
cote, et c’est, silhouettés avec une verve bon-enfant et à peine
caricaturale, les carrures et les maigreurs, les dos voûtés ou les
sveltesses connues des hautes personnalités des courses; leur tenue
habituelle a même été saisie, le chiffonnage de leur cravate et la
façon de porter en arrière ou en avant la jumelle ou la sacoche. Et
c’est la barbe blanche et carrée de M. Aumont à côté du petit
pardessus mastic de M. Edmond Blanc. Voici la canne de M. Abeille
et le melon-cap de Morand, le pince-nez de Veil-Picard, la moustache
et le pantalon tirebouchonné du baron Finot, la rose rouge de M.
Schickler, les sourcils méphistophéliques de Rochefort, et les favoris
en nageoires du baron de Rothschild, et c’est Pratt, et c’est,
démantibulé et ressemblant à crier, M. Ledat. La lorgnette et la
canne de Brémond voisinent avec l’ombrelle et la cravate à pois de
Saint-Albin, et c’est Deschamps, voûté et pensif, et c’est la baronne
de Rothschild et, merveilleux entre tous, les trois joyaux comiques
de cette série: le prince Troubetskoï indiqué par trois traits de
maître, le comte Boni de Castellane, pincé, sanglé, cambré, d’une
jolie prétention précautionneuse et proprette, l’air d’un chat de
marquise déguisé en clubman, et alors une Polaire écrasée sous
l’auréole d’un immense chapeau rouge, l’air d’une goule d’Egypte
hilare et lubrique avec son long sourire et ses yeux plus longs
encore, sur une taille de guêpe exaspérée, douloureuse de minceur.

Lundi 18 juin.—Le Grand Bazar. Coins d’Exposition. La douleur

d’Aïscha. A l’«Andalousie», dans la petite courette des Ouled; huit
heures du soir, avant la représentation.
Les Ouled viennent de dîner, et, avant de descendre dans le patio
dallé de la cour des Lions, l’amusante reconstitution du Généralife où
les visiteurs peuvent les admirer tous les soirs, pittoresquement
groupées autour de la fontaine, affalées en tas dans les oripeaux
barbares du désert et des poses abandonnées de fauves au repos,
elles «farnientent» et prennent le frais sur la petite terrasse de terre
battue, fumant, qui une cigarette, prenant, qui une tasse de kaoua,
disséminées au hasard des chaises, toutes bruissantes de joyaux et
de soie, le coude aux tables, le geste las et la pensée ailleurs.
C’est l’heure où les Ouled songent à la patrie absente. Oh! le bois
de palmiers de l’oasis de Biskra!
Dans un petit groupe formé de trois Ouled, l’une tient une lettre
déployée à la main et furtivement s’essuie les yeux: elle pleure. Une
curiosité me prend et aussi une pitié. La lettre est écrite en arabe;
j’ai longtemps habité l’Algérie, je m’approche de l’Ouled: «Vous vous
ennuyez ici? vous regrettez l’Afrique? —Oh! oui, ici, jamais sortir,
jamais libre. Si j’avais su, serais jamais venue... Je compte les
jours... Encore cinq mois! C’est mes parents qui ont voulu! Les
reverrai-je jamais?» Et elle me montre la lettre.
 Aïscha (elle s’appelle Aïscha) n’est pas une Ouled, c’est une
danseuse.
Elle est née à Constantine et danse à Biskra dans le café-concert
des spahis, auprès du quartier de cavalerie; ses parents viennent la
voir deux fois par an dans l’oasis. Aïscha n’est pas une vulgaire
Ouled, une fille du désert vouée par sa naissance même à la
prostitution; elle n’est pas de la tribu décriée: c’est sa famille qui lui
a fait embrasser la carrière de danseuse, comme elle l’a contrainte à
suivre un barnum à Paris pour l’Exposition.
Paris! Elle est venue aussi un peu par curiosité de ce Paris qu’elle
ne visitera jamais, ce Paris qu’elle sent gronder, rire et haleter
derrière les murailles de l’Andalousie, et où on ne la mène jamais!
Les danses commencent à deux heures de l’après-midi et finissent à
onze heures et demie du soir. Aïscha et ses compagnes partiront fin
octobre sans connaître la grande ville monstrueuse et sonore dont le
mirage les a sûrement attirées du fond du Sahara, au delà des
immensités bleues de la Méditerranée. Et une langueur de
prisonnière accable le front et les yeux d’Aïscha. Une lettre reçue, ce
soir même, de Constantine a réveillé sa peine et la danseuse arabe
pleure.
D’autres danseuses, toutes sonnaillantes d’anneaux, cuirassées
de pierreries et luisantes de soie, se sont approchées de nous,
curieuses. Leurs grands yeux noirs gouachés de kohl nous
observent; leurs regards cherchent à comprendre, à saisir.
Là-bas, au théâtre espagnol, les contorsions de Pepe le gitane
commencent à attirer la foule, qu’appellent les «olle» et les
battements de main de sa troupe... et j’abandonne la courette des
Ouled et la douleur d’Aïscha.

Mardi 19 juin.—Le Grand Bazar. M. de Max, directeur de théâtre;

un drame de M. le comte de Pesquidoux au théâtre égyptien.
 C’est dans ce cadre que M. de Max va monter et incarner le
«Ramsès» de M. de Pesquidoux, le poète mondain, qui donna, il y a
un an, une «Salomé» au Nouveau-Théâtre.
Ramsès au théâtre égyptien, les Pharaons dans le temple de
Dandour, sous la colonnade de Philæ et parmi les Osiris et les dieux
à tête de chien des Memphis et des Thèbes aux cent portes de
l’Egypte légendaire, voilà de la couleur locale ou je ne m’y connais
pas.
Et l’on parle d’un décor prestigieux, d’une Memphis verte tout en
bronze et en marbre vert, avec toute une figuration vêtue de vert,
un immense et mouvant joyau d’émail vert, ciels de turquoise
malade, costumes mordorés et glauques, où, parmi des lueurs de
métal et des frissonnements de palmes, M. de Max, drapé d’étoffes
changeantes et vertes, paré de colliers de jade et casqué de
l’épervier d’or, apparaîtra comme un grand scarabée humain à côté
de madame Eugénie Nau en Juive, l’unique et délicate blancheur de
ce drame... de clair-obscur et de cauchemar...

Mercredi 20 juin.—Le Grand Bazar. Au Trocadéro. Le Stéréorama

mouvant. Une des plus belles choses de l’Exposition, le coin hanté,
visité par les artistes et les peintres, une vision d’art et de réalité
comme n’en ont jamais encore donné les panoramas et maréoramas
des précédentes Expositions, une joie et une émotion, toute
l’envolée et la sensation du départ, de la vie libre des traversées
dans la mélancolie et la gaieté des ciels changeant d’heure en heure,
sous la large et remuante caresse de la mer.
Et c’est Bône apparue dans les gris de lin et le rose émoi d’une
aurore en mer; des barques de pêcheurs sortent du port; les vagues
remuent des paillettes de nacre. Puis voici le bleu profond et la
dureté d’émail de la Méditerranée au large. Des lames courent,
éperdues, dans des festons d’écume; des lueurs tombent d’un ciel
orageux qui les étament; elles déroulent un remous de plomb en
fusion sous un horizon blanc de vapeurs... Le cap Carbon s’y dresse;
ses pentes d’un gris qui se violace, montent, grandissent et pointent
en éperon sur vous..., un promontoire de roches abruptes d’une
douceur à l’œil de soie et de pétales de fleurs: ce sont des mauves
et des roses atténués et fondus où les montagnes ont l’air de
grandes arabesques... O lumière de l’Algérie! Puis voici les fumées et
les hautes coques de l’escadre. Alger apparaît comme une carrière
entre les hauteurs vertes de Mustapha et les rochers de Saint-
Eugène. Voici sa kasbah soleilleuse et blanche et son aspect de
vieille falaise. Qui reconnaîtrait une ville dans cette brèche de pierre
calcaire? C’est Alger pourtant!
Et puis voilà la mer encore, une mer d’huile caressée par une
lumière jeune et calme et enfin, dans des zébrures d’or rose et des
floconnements de braise éparse sur du vert, voici Oran et Mers-el-
Kébir dans l’ardeur enflammée d’un coucher de soleil!
Le bruit de la machine à vapeur, qui fait mouvoir les cylindres
invisibles où se développent toutes ces toiles, ajoute encore à
l’illusion: c’est bien le halètement d’une machine de paquebot,
l’effort continu, trépidant d’une hélice. On est à bord d’un steamer.
Auprès de moi, dans l’ombre, une dame chancelle, se cramponne
éperdument au bras de l’homme qui l’accompagne. On la conduit au
banc qui règne au fond de la salle, on la déplace, on la fait asseoir.
L’illusion est trop forte: la spectatrice a le mal de mer.
Une réflexion: «Epatant! murmure un peintre. Il n’y a plus de
peinture: tous les tableaux f... le camp à côté de cela!»
Et l’auteur de cet émerveillement est un M. Gadan, un nom hier
inconnu et que maintenant il faut retenir.

Vendredi 22 juin.—La «ville en feu». Dix heures du soir, sur le

pont de la Concorde. La ville en feu, c’est l’autre ville, la provisoire,
celle dont les palais, les restaurants, les mosquées, les buvettes, les
pagodes et les beuglants, les tréteaux et les cathédrales flambent et
rougeoient, tous les dimanches, pour la foule et, tous les vendredis,
pour les privilégiés à cinq tickets, du pont de la Concorde à la
passerelle d’Iéna: le Paris de l’Exposition silhouetté, ce soir, dans un
reflet d’incendie sous la voûte nocturne, apparue plus profonde et
comme reculée devant toutes ces lueurs.
Le château d’eau marche enfin et dans un braisillement de vitrail
encadre sous son fronton versicolore des jaillissements et des
cascades de saphirs et de rubis liquides, puis de topazes et de
sardoines («Orgeat, absinthe, limonade, bière!» comme
goguenardaient déjà les curieux devant les fontaines lumineuses de
1889). Le Trocadéro, faufilé de gros points d’or dans tous les détails
de son architecture, échafaude dans la nuit un Orient de ramadan et
de bazar: minarets auréolés de lampions, dômes et terrasses
illuminés a giorno pour la grande joie des yeux, amusés par ces jeux
de l’électricité et des ténèbres. Mais le grand spectacle est dans la
travée sombre et miroitante du fleuve, dans la Seine brusquement
étranglée par les palais de la rue des Nations et les serres de la rue
de Paris et charriant dans ses eaux des reflets et des flammes, la
Seine changée en une coulée de lave incandescente entre les pierres
des quais et les piliers des ponts.
Oh! la magie de la nuit, de la nuit multiforme et changeante! La
porte Binet et ses grotesques pylônes devenus d’émail translucide y
prennent une certaine grandeur. Symbole de ce temps, c’est une
porte de Byzance, le dôme de je ne sais quelle entrée de Tiflis ou de
Samarcande que domine la «Parisienne». Paris, ville d’Orient, Paris
conquis par le Levantin, voilà ce qu’enseigne à la foule, s’écrasant
sur les ponts et dans les Champs-Elysées, le quotidien incendie de la
monstrueuse foire qui bat, là-bas, son plein.
Dans les groupes. «—Ça devient presque beau, cette architecture
salamandre. —Oui, quand on ne la voit pas! —La nuit arrange tout.
Et puis c’est de l’architecture provisoire: ça ne restera pas. —On a
voulu nous éviter des regrets... attention délicate des architectes de
M. Picard!... —Aussi, ce que Paris me paraît beau depuis l’Exposition!
J’y découvre tous les jours d’anciens coins qui me ravissent. —En
dehors des guichets? —Naturellement! —Ainsi, du pont Alexandre
(et celui-là, je vous l’accorde, il est d’une belle audace), le grandiose
du Louvre et du pont Neuf apparus après les verdures des Tuileries.
—Et les fonds de Billancourt et des coteaux de Meudon, le vaporeux
même des usines de Javel, commandées par la grande «Liberté» du
pont de Grenelle, comme cela se compose, le soir, vers les six
heures, vu du pont d’Iéna! —Et les dômes du Sacré-Cœur, dont on a
tant médit, comme ils couronnent bien Montmartre et en font une
jolie ville italienne, on dirait d’un primitif, vus de l’avenue Montaigne!
—En effet, ils s’échelonnent bien, ces dômes romans, et vous
consolent quand on sort de la rue de Paris. —Ils en ont pourtant fait
couler, de la bonne encore! L’a-t-on assez insultée et vilipendée,
cette architecture du Sacré-Cœur! —Mais c’était avant le Manoir à
l’envers. —Et la fresque des «Bonshommes Guillaume». —Et l’île
Saint-Louis? Jamais je ne l’ai tant aimée que maintenant. —Et le
chevet de Notre-Dame! —Vous voyez que l’Exposition a du bon.»

Dimanche 24 juin.—Au Val-Meudon, chez Rodin, une heure après

midi.
«Ce qui entend le plus de bêtises, c’est un tableau un jour de
vernissage», est-il écrit dans le journal des Goncourt. Ce qui a fait
couler le plus d’encre et de toutes les nuances, c’est bien ce pauvre
grand Rodin, que les uns ont voulu voir satanique, les autres mage,
astrologue et même thaumaturge, quand il est simplement un grand
artiste amoureux et fervent de la beauté sous toutes ses formes, la
beauté multiple et diverse, dont il a su saisir et fixer les aspects sous
son pouce de sculpteur génial.
Il est là, présidant la table, amusant à regarder avec sa longue
barbe d’anachorète, ses petits yeux embusqués sous son grand front
de penseur, toute la sensualité de sa nature éprise de beauté
indiquée dans la vibration des narines et la mobilité du long nez
bougeur: tête fruste, intelligente et madrée de faune devenu ermite
et qui, sous la bure du moine, guetterait encore la dryade dans le
mystère des soirs. Il y a là Thaulow et sa femme, tous deux bien
norvégiens avec leurs tailles de géant et leurs regards limpides; l’air
d’un roi de légende en vérité, Thaulow, avec sa belle barbe ondée et
ses grands yeux rieurs, mais d’un roi Gambrinus, monarque
débonnaire et grand buveur de bière. Il y a là Escudier et sa face
nerveuse Paul Escudier, le promoteur du pavillon Rodin et de
l’exposition de l’avenue de Montaigne; il y a là de Braisnes et
d’autres encore.
On parle du congrès féministe, du roi de Suède et de la soirée de
la veille. Rodin en a rapporté le programme. L’Institut ne s’est pas
mis en frais: on dirait un prospectus de parfumeur. Thaulow, qui en
sa qualité de Norvégien goûte juste le roi de Suède, explique sa
popularité dans une phrase définitive: «Il descend si bien le perron
du château!» Escudier, qui, l’avant-veille, a reçu le congrès féministe
à l’Hôtel de Ville, est forcé de constater, à quelques exceptions près,
l’inélégance des femmes progressistes et leur indéniable laideur. «—
Mais que veulent-elles, en somme, toutes ces féministes? —Quand
une femme se préoccupe tant de la question sociale, c’est que les
hommes ne s’occupent plus d’elle. —Le rêve de la féministe, je vais
vous le dire, moi: avoir un homme pour elle seule, tandis que
l’homme rêve de se partager toutes les femmes des autres. D’où le
désaccord... —... éternel. —Question sociale, question sexuelle. —En
somme, ce qu’elles veulent, c’est autre chose. —Et ce que c’est
femme!»
Par les fenêtres ouvertes, c’est le plus merveilleux panorama
peut-être des environs de Paris, toute la vallée de la Seine dominée
et commandée depuis le pont de Billancourt jusqu’à celui de Saint-
Cloud, toute la travée du fleuve profondément encaissée entre deux
coteaux de verdure, les hauts ombrages de Bellevue et de Meudon
moutonnant à l’infini jusqu’aux cimes lointaines du parc de la
Manufacture et, sur la coulée d’eau luisante de la rivière, la douce et
noble grisaille de pierre du vieux pont de Sèvres: un fond de
paysage qui fait songer à ceux de Watteau et que Rodin, ce
voluptueux et ce caressant amant de la beauté, contemple avec des
yeux si clairs qu’on les dirait lavés de larmes.
 Lundi 25 juin.—Le Grand Bazar. Dix heures du soir, dans le
grouillement de la rue de Paris. —Tous ceux et toutes celles qui
veulent être vus sont là; manteaux du soir mousseux, mousselines
de soie, ruches et franfreluches, raglans flottants et cravates
blanches. Boniments sur les tréteaux, parades et lazzi. Etalé sur
deux rangs de chaises, c’est Paris qui regarde dédaigneusement se
bousculer la province; vautrée sur les mêmes chaises, c’est la
province qui regarde curieusement défiler tout Paris.
Dans les feuillages, de grosses oranges lumineuses, qui sont
autant de lanternes vénitiennes, prolongent de l’Alma aux Invalides
une chimérique foire de Neuilly.
Ils et Elles causent. «—Comme vous arrivez tard! —Ne m’en
parlez pas! Nous sortons de «Ramsès». Ça a commencé à neuf
heures moins le quart. —Ça vaut la peine? —Il y a un beau décor. —
Et ça a marché? —L’électricité, mal. —Au Trocadéro, parbleu! —Et la
pièce? —Elle doit être de Reinach. —? —Comment? —Oui: il y a là
dedans un bon Juif, un invraisemblable Juif de la captivité d’Egypte
qui, de désespoir de voir sa fille aimée par le Pharaon, lui plonge un
poignard dans le sein. —Virginius d’Israël. Le fait est rare—... et peu
biblique, quand on songe que Mardochée fit mariner six mois dans
les aromates et les fards sa nièce Esther, qu’il destinait à Assuérus,
et que les grands-rabbins de Béthulie dépêchèrent Judith à
Holopherne pour délivrer la ville assiégée. Les bons Juifs d’alors ne
répugnaient pas à se servir des femmes pour dénouer et précipiter
les événements. —Oui, le père juif de M. de Pesquidoux est un
unique exemple... —... qu’ils n’ont pas suivi d’ailleurs, car vous vous
souvenez, au moment de l’Affaire...»
Et le groupe s’enfonce dans les groupes.
La façade de la Loïe Fuller rutile et flamboie dans un savant et
subtil éclairage des longues draperies qui la décorent; l’innombrable
figuration des Auteurs-Gais —paillasses, paillettes et paillons— se
cambre, s’agite et chatoie et se ploie.
 Devant une des baraques, un pitre lugubre ressasse ce boniment
de mauvais présage: «Entrez, mesdames et messieurs, entrez!
Demain, il sera peut-être trop tard. Vous pouvez recevoir sur la tête
une passerelle: nous sommes à l’Exposition (sic).»

Mercredi 27 juin.—Au Trocadéro, trois heures. L’union chorale des

étudiants d’Upsal.
 Entends-nous, Svea, notre mère à tous!—Fais-nous lutter pour ton bien jusqu’à
la mort!—Jamais nous ne te trahirons,—reçois-en notre serment, toujours
inébranlable!—A outrance nous défendrons—le pays libre qui encore est le nôtre,—
chaque parcelle de l’héritage—que tu laissas dans nos sagas et dans nos champs.
—Mais, si, par la ruse, la félonie, par la discorde et la violence, tu es menacée—
nous nous confions en l’Eternel—comme jadis nos Pères.
Il est beau alors, il est beau—d’être vainqueur dans le combat,—mais plus
beau encore,—ô mère, de mourir pour toi!

Et les voix sonores et pures, unies dans un merveilleux accord,

mieux qu’unies, mêlées et fondues s’enflent comme une mer,
montent comme une sève et s’épanouissent en une espèce de
floraison d’âmes et d’harmonies qui est l’âme même de leur race, à
ces Suédois blonds et guerriers qui déchaînèrent Charles XII en
ouragan sur l’Europe et gardent encore, eux, le solide amour du
pays, le culte des traditions et leur sang intact de vieux Northmans.
Ils sont là, groupés, tassés sur l’immense scène du Trocadéro,
leur casquette blanche à la main, tous étudiants de cette université
d’Upsal où, mêlés à leurs sujets, les rois de Suède étudient deux
ans, dociles à une vieille tradition du royaume; et, parmi les
chanteurs assemblés là, il y a des médecins, des avocats, des
magistrats, tous anciens étudiants, demeurés solidaires de
l’université et venus avec les jeunes révéler et affirmer à Paris la
patrie suédoise et l’union d’Upsal.
Et c’est une réconfortante chose que de voir et d’entendre
combien ces braves gens à la voix si pure, si vierge, pour ainsi dire,
aiment passionnément et fièrement leur Suède, dont ils célèbrent
dans leurs chants et les joies populaires et les vieilles légendes.
Et c’est le chant de Suomi, tout retentissant du murmure des
sapins et du mugissement des torrents, et c’est le poème d’Olaf
Trygvason, la tragique et l’épique ballade sur la mort du vieux roi de
Norvège trahi et coulé sur son navire par le roi de Danemark et
pleuré par son fidèle Erling, et c’est, sautillant et léger comme le
lièvre dans les bruyères, alerte comme le vent du matin, le «chant
d’Ingrid», tout scintillant de rosée, tout brillant de soleil.
Holà! houp! fais-tu si hauts tes sauts
Sur la bruyère?

Puis, mélancolique, d’une nostalgie d’exil et de regret, voici la

«chanson du Neck», entendue déjà sur les lèvres de la Nilsson dans
l’«Hamlet» d’Ambroise Thomas, la légendaire et populaire mélodie
que l’auteur du «Caïd» a mise dans la bouche d’Ophélie, rythmée
par Barbier:
Pâle et blonde
Dort sous l’onde
La willis au regard de feu.

Et alors toute la joie à la fois simple et un peu brutale des

paysans de là-bas, mise en musique par le génie de Sœderman dans
la suite intitulée «Noces de paysans suédois».
Oh! la gaieté de la valse chantée du «cortège nuptial», le
caractère allègre de la «chanson des souhaits». la joie un peu
lourde, mais si fortement cadencée de la danse à la ferme et du
chœur des buveurs.
De la bière forte, nous en boirons
Jusqu’à ce que nous ayons le hoquet;
Le tonneau, nous le viderons
Rien ne doit être gardé!
Du pied marquez la mesure, garçons!
Voilà comme il faut danser!

Mais, entre toutes ces mélodies, la perle, le joyau est l’émotion

naïve, la simplicité touchante du chœur intitulé «A l’Eglise», sur ces
délicieuses paroles, en l’honneur des mariés:
«Sur le chemin qu’ils vont suivre tous deux.
Dirige-les, Seigneur, par ta sagesse,
Et que, l’un sur l’autre s’appuyant,
Ils soient soutenus par ta grâce!»
 Le couple s’arrête devant les degrés de l’autel
Où le pasteur les attend le livre à la main.
Il veut bénir cette union
Et nouer pour jamais ces liens sacrés.
Les vieilles sanglotent, et les vieux restent debout,
Les larmes aux yeux, les cœurs joyeux;
Et les filles d’honneur pensent certainement toutes:
«La prochaine fois, ce sera pour moi qu’on jouera la marche nuptiale!»

En l’honneur de ces braves gens qui nous mettent le cœur en

fête, crions allègrement: «Vive la Suède!» comme ils crient eux-
mêmes de tous leurs poumons: «Vive la France!»

Samedi 30 juin.—A l’Opéra-Comique. Mme Rose Caron dans

«Iphigénie». —Il n’y a pas à dire, le plus beau rôle de toute la
longue carrière de Brunehilde, de Salammbô et d’Elvire. Si
persistante que nous soit demeurée à tous la vision de Mme Caron
dans «Sigurd», quand, svelte et blanche dans un rai de lune et
couronnée de verveine, elle effeuillait —de quelle voix délicate et
pure!— et la sauge pourprée et les aveux de son âme dans le
courant du torrent, comme maîtrise de style, comme silhouette
héroïque et comme harmonie de gestes, on ne peut pas aller plus
loin que Mme Rose Caron dans cette dernière création d’«Iphigénie».
Malgré l’usure indéniable de la voix, elle trouve, au second acte
surtout, des accents de tendresse d’une mélancolie si touchante
qu’ils en effacent jusqu’au souvenir de Mme Raunay.
Mme Raunay, à la Renaissance! De quel dithyrambe ne l’avions-
nous pas accueillie quand, droite et svelte sous les longs voiles de la
prêtresse de Diane, elle nous apparut, l’automne dernier, dans la
noble et attendrissante partition de Gluck? Mais, si belles qu’aient
été les attitudes de l’Iphigénie de la rue de Bondy, rien ne peut lutter
avec la grâce contenue, le charme de tristesse et de résignation, le
parfum de pitié et de ferveur qui s’émanent, comme une atmosphère
de beauté psychique, de la bouche, de la physionomie, du port de
tête, de la démarche et du moindre geste de l’Iphigénie de la rue
Favart, Iphigénie parfaite qui, pareille à une admirable statue
Welcome to Our Bookstore - The Ultimate Destination for Book Lovers
Are you passionate about books and eager to explore new worlds of
knowledge? At our website, we offer a vast collection of books that
cater to every interest and age group. From classic literature to
specialized publications, self-help books, and children’s stories, we
have it all! Each book is a gateway to new adventures, helping you
expand your knowledge and nourish your soul
Experience Convenient and Enjoyable Book Shopping Our website is more
than just an online bookstore—it’s a bridge connecting readers to the
timeless values of culture and wisdom. With a sleek and user-friendly
interface and a smart search system, you can find your favorite books
quickly and easily. Enjoy special promotions, fast home delivery, and
a seamless shopping experience that saves you time and enhances your
love for reading.
Let us accompany you on the journey of exploring knowledge and
personal growth!

ebookgate.com