Applied Numerical Methods with MATLAB for Engineers and Scientists Fourth Edition Steven C. Chapra Dr. download
Applied Numerical Methods with MATLAB for Engineers and Scientists Fourth Edition Steven C. Chapra Dr. download
https://textbookfull.com/product/applied-numerical-methods-with-
matlab-for-engineers-and-scientists-fourth-edition-steven-c-
chapra-dr/
https://textbookfull.com/product/numerical-methods-for-engineers-
and-scientists-using-matlab-second-edition-esfandiari/
https://textbookfull.com/product/numerical-methods-for-engineers-
and-scientists-using-matlab-ramin-s-esfandiari/
https://textbookfull.com/product/numerical-methods-for-
engineers-8ed-2021-8th-edition-chapra-s/
https://textbookfull.com/product/essential-matlab-for-engineers-
and-scientists-hahn/
Essential MATLAB for engineers and scientists Hahn
https://textbookfull.com/product/essential-matlab-for-engineers-
and-scientists-hahn-2/
https://textbookfull.com/product/essential-matlab-for-engineers-
and-scientists-hahn-3/
https://textbookfull.com/product/modern-physics-with-modern-
computational-methods-for-scientists-and-engineers-3rd-edition-
john-c-morrison/
https://textbookfull.com/product/variational-methods-for-
engineers-with-matlab-1st-edition-cursi/
https://textbookfull.com/product/essential-matlab-for-engineers-
and-scientists-sixth-edition-edition-hahn/
Applied Numerical Methods
with MATLAB® for Engineers and Scientists
Fourth Edition
Steven C. Chapra
Berger Chair in Computing and Engineering
Tufts University
APPLIED NUMERICAL METHODS WITH MATLAB® FOR ENGINEERS AND SCIENTISTS,
FOURTH EDITION
Published by McGraw-Hill Education, 2 Penn Plaza, New York, NY 10121. Copyright © 2018 by McGraw-Hill
Education. All rights reserved. Printed in the United States of America. Previous editions © 2012, 2008, and
2005. No part of this publication may be reproduced or distributed in any form or by any means, or stored in a
database or retrieval system, without the prior written consent of McGraw-Hill Education, including, but not
limited to, in any network or other electronic storage or transmission, or broadcast for distance learning.
Some ancillaries, including electronic and print components, may not be available to customers outside the
United States.
1 2 3 4 5 6 7 8 9 LCR 21 20 19 18 17
ISBN 978-0-07-339796-2
MHID 0-07-339796-2
All credits appearing on page or at the end of the book are considered to be an extension of the copyright page.
mheducation.com/highered
To
My brothers,
John and Bob Chapra
and
Steve Chapra teaches in the Civil and Environmental Engineering Department at Tufts
University, where he holds the Louis Berger Chair in Computing and Engineering. His other
books include Numerical Methods for Engineers and Surface Water-Quality Modeling.
Steve received engineering degrees from Manhattan College and the University of
Michigan. Before joining the faculty at Tufts, he worked for the Environmental Protection
Agency and the National Oceanic and Atmospheric Administration, and taught at Texas
A&M University and the University of Colorado. His general research interests focus on
surface water-quality modeling and advanced computer applications in environmental
engineering.
He has received a number of awards for his scholarly contributions, including the
Rudolph Hering Medal, the Meriam/Wiley Distinguished Author Award, and the Chandler-
Misener Award. He has also been recognized as the outstanding teacher at Texas A&M
University (1986 Tenneco Award), the University of Colorado (1992 Hutchinson Award),
and Tufts University (2011 Professor of the Year Award).
Steve was originally drawn to environmental engineering and science because of his
love of the outdoors. He is an avid fly fisherman and hiker. An unapologetic nerd, his love
affair with computing began when he was first introduced to Fortran programming as an
undergraduate in 1966. Today, he feels truly blessed to be able to meld his love of math-
ematics, science, and computing with his passion for the natural environment. In addition,
he gets the bonus of sharing it with others through his teaching and writing!
Beyond his professional interests, he enjoys art, music (especially classical music,
jazz, and bluegrass), and reading history. Despite unfounded rumors to the contrary, he
never has, and never will, voluntarily bungee jump or sky dive.
If you would like to contact Steve, or learn more about him, visit his home page at
http://engineering.tufts.edu/cee/people/chapra/ or e-mail him at [email protected].
iv
CONTENTS
Preface xiv
CHAPTER 1
Mathematical Modeling, Numerical Methods,
and Problem Solving 4
1.1 A Simple Mathematical Model 5
1.2 Conservation Laws in Engineering and Science 12
1.3 Numerical Methods Covered in This Book 13
1.4 Case Study: It’s a Real Drag 17
Problems 20
CHAPTER 2
MATLAB Fundamentals 27
2.1 The MATLAB Environment 28
2.2 Assignment 29
2.3 Mathematical Operations 36
2.4 Use of Built-In Functions 39
2.5 Graphics 42
2.6 Other Resources 46
2.7 Case Study: Exploratory Data Analysis 46
Problems 49
CHAPTER 3
Programming with MATLAB 53
3.1 M-Files 54
3.2 Input-Output 61
v
vi CONTENTS
CHAPTER 4
Roundoff and Truncation Errors 99
4.1 Errors 100
4.2 Roundoff Errors 106
4.3 Truncation Errors 114
4.4 Total Numerical Error 125
4.5 Blunders, Model Errors, and Data Uncertainty 130
Problems 131
CHAPTER 5
Roots: Bracketing Methods 138
5.1 Roots in Engineering and Science 139
5.2 Graphical Methods 140
5.3 Bracketing Methods and Initial Guesses 141
5.4 Bisection 146
5.5 False Position 152
5.6 Case Study: Greenhouse Gases and Rainwater 156
Problems 159
CHAPTER 6
Roots: Open Methods 164
6.1 Simple Fixed-Point Iteration 165
6.2 Newton-Raphson 169
6.3 Secant Methods 174
6.4 Brent’s Method 176
6.5 MATLAB Function: fzero 181
6.6 Polynomials 183
6.7 Case Study: Pipe Friction 186
Problems 191
CONTENTS vii
CHAPTER 7
Optimization 198
7.1 Introduction and Background 199
7.2 One-Dimensional Optimization 202
7.3 Multidimensional Optimization 211
7.4 Case Study: Equilibrium and Minimum Potential Energy 213
Problems 215
CHAPTER 8
Linear Algebraic Equations and Matrices 227
8.1 Matrix Algebra Overview 229
8.2 Solving Linear Algebraic Equations with MATLAB 238
8.3 Case Study: Currents and Voltages in Circuits 240
Problems 244
CHAPTER 9
Gauss Elimination 248
9.1 Solving Small Numbers of Equations 249
9.2 Naive Gauss Elimination 254
9.3 Pivoting 261
9.4 Tridiagonal Systems 264
9.5 Case Study: Model of a Heated Rod 266
Problems 270
CHAPTER 10
LU Factorization 274
10.1 Overview of LU Factorization 275
10.2 Gauss Elimination as LU Factorization 276
10.3 Cholesky Factorization 283
10.4 MATLAB Left Division 286
Problems 287
viii CONTENTS
CHAPTER 11
Matrix Inverse and Condition 288
11.1 The Matrix Inverse 288
11.2 Error Analysis and System Condition 292
11.3 Case Study: Indoor Air Pollution 297
Problems 300
CHAPTER 12
Iterative Methods 305
12.1 Linear Systems: Gauss-Seidel 305
12.2 Nonlinear Systems 312
12.3 Case Study: Chemical Reactions 320
Problems 323
CHAPTER 13
Eigenvalues 326
13.1 Mathematical Background 328
13.2 Physical Background 331
13.3 The Power Method 333
13.4 MATLAB Function: eig 336
13.5 Case Study: Eigenvalues and Earthquakes 337
Problems 340
CHAPTER 14
Linear Regression 346
14.1 Statistics Review 348
14.2 Random Numbers and Simulation 353
14.3 Linear Least-Squares Regression 358
14.4 Linearization of Nonlinear Relationships 366
14.5 Computer Applications 370
14.6 Case Study: Enzyme Kinetics 373
Problems 378
CONTENTS ix
CHAPTER 15
General Linear Least-Squares and Nonlinear Regression 385
15.1 Polynomial Regression 385
15.2 Multiple Linear Regression 389
15.3 General Linear Least Squares 391
15.4 QR Factorization and the Backslash Operator 394
15.5 Nonlinear Regression 395
15.6 Case Study: Fitting Experimental Data 397
Problems 399
CHAPTER 16
Fourier Analysis 404
16.1 Curve Fitting with Sinusoidal Functions 405
16.2 Continuous Fourier Series 411
16.3 Frequency and Time Domains 414
16.4 Fourier Integral and Transform 415
16.5 Discrete Fourier Transform (DFT) 418
16.6 The Power Spectrum 423
16.7 Case Study: Sunspots 425
Problems 426
CHAPTER 17
Polynomial Interpolation 429
17.1 Introduction to Interpolation 430
17.2 Newton Interpolating Polynomial 433
17.3 Lagrange Interpolating Polynomial 441
17.4 Inverse Interpolation 444
17.5 Extrapolation and Oscillations 445
Problems 449
CHAPTER 18
Splines and Piecewise Interpolation 453
18.1 Introduction to Splines 453
18.2 Linear Splines 455
18.3 Quadratic Splines 459
18.4 Cubic Splines 462
18.5 Piecewise Interpolation in MATLAB 468
18.6 Multidimensional Interpolation 473
18.7 Case Study: Heat Transfer 476
Problems 480
x CONTENTS
CHAPTER 19
Numerical Integration Formulas 488
19.1 Introduction and Background 489
19.2 Newton-Cotes Formulas 492
19.3 The Trapezoidal Rule 494
19.4 Simpson’s Rules 501
19.5 Higher-Order Newton-Cotes Formulas 507
19.6 Integration with Unequal Segments 508
19.7 Open Methods 512
19.8 Multiple Integrals 512
19.9 Case Study: Computing Work with Numerical Integration 515
Problems 518
CHAPTER 20
Numerical Integration of Functions 524
20.1 Introduction 524
20.2 Romberg Integration 525
20.3 Gauss Quadrature 530
20.4 Adaptive Quadrature 537
20.5 Case Study: Root-Mean-Square Current 540
Problems 544
CHAPTER 21
Numerical Differentiation 548
21.1 Introduction and Background 549
21.2 High-Accuracy Differentiation Formulas 552
21.3 Richardson Extrapolation 555
21.4 Derivatives of Unequally Spaced Data 557
21.5 Derivatives and Integrals for Data with Errors 558
21.6 Partial Derivatives 559
21.7 Numerical Differentiation with MATLAB 560
21.8 Case Study: Visualizing Fields 565
Problems 567
CONTENTS xi
CHAPTER 22
Initial-Value Problems 579
22.1 Overview 581
22.2 Euler’s Method 581
22.3 Improvements of Euler’s Method 587
22.4 Runge-Kutta Methods 593
22.5 Systems of Equations 598
22.6 Case Study: Predator-Prey Models and Chaos 604
Problems 609
CHAPTER 23
Adaptive Methods and Stiff Systems 615
23.1 Adaptive Runge-Kutta Methods 615
23.2 Multistep Methods 624
23.3 Stiffness 628
23.4 MATLAB Application: Bungee Jumper with Cord 634
23.5 Case Study: Pliny’s Intermittent Fountain 635
Problems 640
CHAPTER 24
Boundary-Value Problems 646
24.1 Introduction and Background 647
24.2 The Shooting Method 651
24.3 Finite-Difference Methods 658
24.4 MATLAB Function: bvp4c 665
Problems 668
McGraw-Hill Connect®
Learn Without Limits
Connect is a teaching and learning platform that is proven
to deliver better results for students and instructors.
Analytics
12.7%, and exam scores by 9.1%.
Connect Insight®
Connect Insight is Connect’s new one-of-a-kind
visual analytics dashboard—now available for both
instructors and students—that provides at-a-
glance information regarding student performance,
which is immediately actionable. By presenting
assignment, assessment, and topical performance
results together with a time metric that is easily visible for
aggregate or individual results, Connect Insight gives the user
the ability to take a just-in-time approach to teaching and
learning, which was never before available. Connect Insight presents
Students can view
data that empowers students and helps instructors improve class their results for any
performance in a way that is efficient and effective. Connect course.
Mobile
Connect’s new, intuitive mobile interface gives students
and instructors flexible and convenient, anytime–anywhere
access to all components of the Connect platform.
Adaptive
THE ADAPTIVE
READING EXPERIENCE
DESIGNED TO TRANSFORM
THE WAY STUDENTS READ
SmartBook®
Proven to help students improve grades and study more
efficiently, SmartBook contains the same content within
the print book, but actively tailors that content to the
needs of the individual. SmartBook’s adaptive technology
provides precise, personalized instruction on what the
student should do next, guiding the student to master
and remember key concepts, targeting gaps in knowledge
and offering customized feedback, and driving the student
toward comprehension and retention of the subject matter.
Available on tablets, SmartBook puts learning at the student’s
fingertips—anywhere, anytime.
www.mheducation.com
PREFACE
This book is designed to support a one-semester course in numerical methods. It has been
written for students who want to learn and apply numerical methods in order to solve prob-
lems in engineering and science. As such, the methods are motivated by problems rather
than by mathematics. That said, sufficient theory is provided so that students come away
with insight into the techniques and their shortcomings.
MATLAB® provides a great environment for such a course. Although other en-
vironments (e.g., Excel/VBA, Mathcad) or languages (e.g., Fortran 90, C++) could
have been chosen, MATLAB presently offers a nice combination of handy program-
ming features with powerful built-in numerical capabilities. On the one hand, its
M-file programming environment allows students to implement moderately compli-
cated algorithms in a structured and coherent fashion. On the other hand, its built-in,
numerical capabilities empower students to solve more difficult problems without try-
ing to “reinvent the wheel.”
The basic content, organization, and pedagogy of the third edition are essentially pre-
served in the fourth edition. In particular, the conversational writing style is intentionally
maintained in order to make the book easier to read. This book tries to speak directly to the
reader and is designed in part to be a tool for self-teaching.
That said, this edition differs from the past edition in three major ways: (1) new
material, (2) new and revised homework problems, and (3) an appendix introducing
Simulink.
1. New Content. I have included new and enhanced sections on a number of topics. The
primary additions include material on some MATLAB functions not included in previ-
ous editions (e.g., fsolve, integrate, bvp4c), some new applications of Monte Carlo
for problems such as integration and optimization, and MATLAB’s new way to pass
parameters to function functions.
2. New Homework Problems. Most of the end-of-chapter problems have been modified,
and a variety of new problems have been added. In particular, an effort has been made
to include several new problems for each chapter that are more challenging and dif-
ficult than the problems in the previous edition.
3. I have developed a short primer on Simulink which I have my students read prior to
covering that topic. Although I recognize that some professors may not choose to
cover Simulink, I included it as a teaching aid for those that do.
xiv
PREFACE xv
Aside from the new material and problems, the fourth edition is very similar to the
third. In particular, I have endeavored to maintain most of the features contributing to its
pedagogical effectiveness including extensive use of worked examples and engineering and
scientific applications. As with the previous edition, I have made a concerted effort to make
this book as “student-friendly” as possible. Thus, I’ve tried to keep my explanations straight-
forward and practical.
Although my primary intent is to empower students by providing them with a sound
introduction to numerical problem solving, I have the ancillary objective of making this
introduction exciting and pleasurable. I believe that motivated students who enjoy engi-
neering and science, problem solving, mathematics—and yes—programming, will ulti-
mately make better professionals. If my book fosters enthusiasm and appreciation for these
subjects, I will consider the effort a success.
Finally, I want to thank my family, and in particular my wife, Cynthia, for the love,
patience, and support they have provided through the time I’ve spent on this project.
Steven C. Chapra
Tufts University
Medford, Massachusetts
[email protected]
PEDAGOGICAL TOOLS
Theory Presented as It Informs Key Concepts. The text is intended for Numerical
Methods users, not developers. Therefore, theory is not included for “theory’s sake,” for ex-
ample no proofs. Theory is included as it informs key concepts such as the Taylor series, con-
vergence, condition, etc. Hence, the student is shown how the theory connects with practical
issues in problem solving.
Introductory MATLAB Material. The text includes two introductory chapters on how to
use MATLAB. Chapter 2 shows students how to perform computations and create graphs
in MATLAB’s standard command mode. Chapter 3 provides a primer on developing
numerical programs via M ATLAB M-file functions. Thus, the text provides students with
the means to develop their own numerical algorithms as well as to tap into MATLAB’s
powerful built-in routines.
Algorithms Presented Using MATLAB M-files. Instead of using pseudocode, this book
presents algorithms as well-structured MATLAB M-files. Aside from being useful com-
puter programs, these provide students with models for their own M-files that they will
develop as homework exercises.
Worked Examples and Case Studies. Extensive worked examples are laid out in detail
so that students can clearly follow the steps in each numerical computation. The case stud-
ies consist of engineering and science applications which are more complex and richer than
the worked examples. They are placed at the ends of selected chapters with the intention
of (1) illustrating the nuances of the methods and (2) showing more realistically how the
methods along with MATLAB are applied for problem solving.
Problem Sets. The text includes a wide variety of problems. Many are drawn from en-
gineering and scientific disciplines. Others are used to illustrate numerical techniques and
theoretical concepts. Problems include those that can be solved with a pocket calculator as
well as others that require computer solution with MATLAB.
1.1 MOTIVATION
What are numerical methods and why should you study them?
Numerical methods are techniques by which mathematical problems are formulated
so that they can be solved with arithmetic and logical operations. Because digital comput-
ers excel at performing such operations, numerical methods are sometimes referred to as
computer mathematics.
In the pre–computer era, the time and drudgery of implementing such calculations
seriously limited their practical use. However, with the advent of fast, inexpensive digital
computers, the role of numerical methods in engineering and scientific problem solving
has exploded. Because they figure so prominently in much of our work, I believe that
numerical methods should be a part of every engineer’s and scientist’s basic education.
Just as we all must have solid foundations in the other areas of mathematics and science,
we should also have a fundamental understanding of numerical methods. In particular, we
should have a solid appreciation of both
their capabilities and their limitations.
Beyond contributing to your overall
education, there are several additional
reasons why you should study numerical
methods:
1. N umerical methods greatly expand the
types of problems you can address.
They are capable of handling large sys-
tems of equations, nonlinearities, and
complicated geometries that are not
uncommon in engineering and science
and that are often impossible to solve
analytically with standard calculus. As
such, they greatly enhance your prob-
lem-solving skills.
2. Numerical methods allow you to use
“canned” software with insight. During
1
2 PART 1 Modeling, Computers, and Error Analysis
your career, you will invariably have occasion to use commercially available prepack-
aged computer programs that involve numerical methods. The intelligent use of these
programs is greatly enhanced by an understanding of the basic theory underlying the
methods. In the absence of such understanding, you will be left to treat such packages
as “black boxes” with little critical insight into their inner workings or the validity of
the results they produce.
3. Many problems cannot be approached using canned programs. If you are conversant
with numerical methods, and are adept at computer programming, you can design
your own programs to solve problems without having to buy or commission expensive
software.
4. Numerical methods are an efficient vehicle for learning to use computers. Because nu-
merical methods are expressly designed for computer implementation, they are ideal for
illustrating the computer’s powers and limitations. When you successfully implement
numerical methods on a computer, and then apply them to solve otherwise intractable
problems, you will be provided with a dramatic demonstration of how computers can
serve your professional development. At the same time, you will also learn to acknowl-
edge and control the errors of approximation that are part and parcel of large-scale
numerical calculations.
5. Numerical methods provide a vehicle for you to reinforce your understanding of math-
ematics. Because one function of numerical methods is to reduce higher mathematics
to basic arithmetic operations, they get at the “nuts and bolts” of some otherwise
obscure topics. Enhanced understanding and insight can result from this alternative
perspective.
With these reasons as motivation, we can now set out to understand how numerical
methods and digital computers work in tandem to generate reliable solutions to mathemati-
cal problems. The remainder of this book is devoted to this task.
This book is divided into six parts. The latter five parts focus on the major areas of nu-
merical methods. Although it might be tempting to jump right into this material, Part One
consists of four chapters dealing with essential background material.
Chapter 1 provides a concrete example of how a numerical method can be employed
to solve a real problem. To do this, we develop a mathematical model of a free-falling
bungee jumper. The model, which is based on Newton’s second law, results in an ordinary
differential equation. After first using calculus to develop a closed-form solution, we then
show how a comparable solution can be generated with a simple numerical method. We
end the chapter with an overview of the major areas of numerical methods that we cover in
Parts Two through Six.
Chapters 2 and 3 provide an introduction to the MATLAB® software environment.
Chapter 2 deals with the standard way of operating MATLAB by entering commands one
at a time in the so-called calculator, or command, mode. This interactive mode provides
a straightforward means to orient you to the environment and illustrates how it is used for
common operations such as performing calculations and creating plots.
1.2 PART ORGANIZATION 3
Chapter 3 shows how MATLAB’s programming mode provides a vehicle for assem-
bling individual commands into algorithms. Thus, our intent is to illustrate how MATLAB
serves as a convenient programming environment to develop your own software.
Chapter 4 deals with the important topic of error analysis, which must be understood
for the effective use of numerical methods. The first part of the chapter focuses on the
roundoff errors that result because digital computers cannot represent some quantities
exactly. The latter part addresses truncation errors that arise from using an approximation
in place of an exact mathematical procedure.
1
Mathematical Modeling,
Numerical Methods,
and Problem Solving
CHAPTER OBJECTIVES
The primary objective of this chapter is to provide you with a concrete idea of what
numerical methods are and how they relate to engineering and scientific problem
solving. Specific objectives and topics covered are
• Learning how mathematical models can be formulated on the basis of scientific
principles to simulate the behavior of a simple physical system.
• Understanding how numerical methods afford a means to generate solutions in a
manner that can be implemented on a digital computer.
• Understanding the different types of conservation laws that lie beneath the models
used in the various engineering disciplines and appreciating the difference
between steady-state and dynamic solutions of these models.
• Learning about the different types of numerical methods we will cover in this
book.
S
uppose that a bungee-jumping company hires you. You’re given the task of
predicting the velocity of a jumper (Fig. 1.1) as a function of time during the
free-fall part of the jump. This information will be used as part of a larger
analysis to determine the length and required strength of the bungee cord for jumpers
of different mass.
You know from your studies of physics that the acceleration should be equal to the ratio
of the force to the mass (Newton’s second law). Based on this insight and your knowledge
4
1.1 A Simple Mathematical Model 5
Upward force of physics and fluid mechanics, you develop the following mathematical model for the rate
due to air of change of velocity with respect to time,
resistance
c
dυ = g − ___
___ md υ2
dt
where υ = downward vertical velocity (m/s), t = time (s), g = the acceleration due to
gravity (≅ 9.81 m/s2), cd = a lumped drag coefficient (kg/m), and m = the jumper’s
mass (kg). The drag coefficient is called “lumped” because its magnitude depends on fac-
tors such as the jumper’s area and the fluid density (see Sec. 1.4).
Because this is a differential equation, you know that calculus might be used to obtain
an analytical or exact solution for υ as a function of t. However, in the following pages, we
will illustrate an alternative solution approach. This will involve developing a computer-
oriented numerical or approximate solution.
Aside from showing you how the computer can be used to solve this particular prob-
Downward lem, our more general objective will be to illustrate (a) what numerical methods are and
force due (b) how they figure in engineering and scientific problem solving. In so doing, we will also
to gravity show how mathematical models figure prominently in the way engineers and scientists use
numerical methods in their work.
FIGURE 1.1
Forces acting
on a free-falling
bungee jumper.
Dependent
variable
= f (
independent
variables
, parameters,
forcing
functions )
(1.1)
where the dependent variable is a characteristic that typically reflects the behavior or state
of the system; the independent variables are usually dimensions, such as time and space,
along which the system’s behavior is being determined; the parameters are reflective of
the system’s properties or composition; and the forcing functions are external influences
acting upon it.
The actual mathematical expression of Eq. (1.1) can range from a simple algebraic
relationship to large complicated sets of differential equations. For example, on the basis
of his observations, Newton formulated his second law of motion, which states that the
time rate of change of momentum of a body is equal to the resultant force acting on it. The
mathematical expression, or model, of the second law is the well-known equation
F = ma (1.2)
where F is the net force acting on the body (N, or kg m/s2), m is the mass of the object (kg),
and a is its acceleration (m/s2).
6 Mathematical Modeling, Numerical Methods, and Problem Solving
The second law can be recast in the format of Eq. (1.1) by merely dividing both sides
by m to give
a = __ F
m (1.3)
where a is the dependent variable reflecting the system’s behavior, F is the forcing func-
tion, and m is a parameter. Note that for this simple case there is no independent variable
because we are not yet predicting how acceleration varies in time or space.
Equation (1.3) has a number of characteristics that are typical of mathematical models
of the physical world.
• It describes a natural process or system in mathematical terms.
• It represents an idealization and simplification of reality. That is, the model ignores
negligible details of the natural process and focuses on its essential manifestations.
Thus, the second law does not include the effects of relativity that are of minimal
importance when applied to objects and forces that interact on or about the earth’s
surface at velocities and on scales visible to humans.
• Finally, it yields reproducible results and, consequently, can be used for predictive
purposes. For example, if the force on an object and its mass are known, Eq. (1.3) can
be used to compute acceleration.
Because of its simple algebraic form, the solution of Eq. (1.2) was obtained easily.
However, other mathematical models of physical phenomena may be much more complex,
and either cannot be solved exactly or require more sophisticated mathematical techniques
than simple algebra for their solution. To illustrate a more complex model of this kind,
Newton’s second law can be used to determine the terminal velocity of a free-falling body
near the earth’s surface. Our falling body will be a bungee jumper (Fig. 1.1). For this case,
a model can be derived by expressing the acceleration as the time rate of change of the
velocity (dυ/dt) and substituting it into Eq. (1.3) to yield
dυ = __
___ F (1.4)
dt m
where υ is velocity (in meters per second). Thus, the rate of change of the velocity is equal
to the net force acting on the body normalized to its mass. If the net force is positive, the
object will accelerate. If it is negative, the object will decelerate. If the net force is zero, the
object’s velocity will remain at a constant level.
Next, we will express the net force in terms of measurable variables and parameters.
For a body falling within the vicinity of the earth, the net force is composed of two
opposing forces: the downward pull of gravity FD and the upward force of air resistance FU
(Fig. 1.1):
F = FD + FU (1.5)
If force in the downward direction is assigned a positive sign, the second law can be
used to formulate the force due to gravity as
FD = mg (1.6)
2
where g is the acceleration due to gravity (9.81 m/s ).
1.1 A Simple Mathematical Model 7
Air resistance can be formulated in a variety of ways. Knowledge from the science
of fluid mechanics suggests that a good first approximation would be to assume that it is
proportional to the square of the velocity,
where cd is a proportionality constant called the lumped drag coefficient (kg/m). Thus, the
greater the fall velocity, the greater the upward force due to air resistance. The parameter
cd accounts for properties of the falling object, such as shape or surface roughness, that af-
fect air resistance. For the present case, cd might be a function of the type of clothing or the
orientation used by the jumper during free fall.
The net force is the difference between the downward and upward force. Therefore,
Eqs. (1.4) through (1.7) can be combined to yield
c
dυ = g − ___
___ md υ2 (1.8)
dt
Equation (1.8) is a model that relates the acceleration of a falling object to the forces
acting on it. It is a differential equation because it is written in terms of the differential rate
of change (dυ/dt) of the variable that we are interested in predicting. However, in contrast
to the solution of Newton’s second law in Eq. (1.3), the exact solution of Eq. (1.8) for the
velocity of the jumper cannot be obtained using simple algebraic manipulation. Rather,
more advanced techniques such as those of calculus must be applied to obtain an exact or
analytical solution. For example, if the jumper is initially at rest (υ = 0 at t = 0), calculus
can be used to solve Eq. (1.8) for
( )
___ ___
gm gc
υ(t) = ___
c tanh √___
√ md t (1.9)
d
where tanh is the hyperbolic tangent that can be either computed directly1 or via the more
elementary exponential function as in
x −x
ex − e−x
tanh x = _______ (1.10)
e +e
Note that Eq. (1.9) is cast in the general form of Eq. (1.1) where υ(t) is the dependent
variable, t is the independent variable, cd and m are parameters, and g is the forcing function.
1
MATLAB allows direct calculation of the hyperbolic tangent via the built-in function tanh(x).
8 Mathematical Modeling, Numerical Methods, and Problem Solving
(√ )
_________
9.81(68.1) 9.81(0.25)
√
υ(t) = _________
0.25
tanh _________
68.1
t = 51.6938 tanh(0.18977t)
t, s υ, m/s
0 0
2 18.7292
4 33.1118
6 42.0762
8 46.9575
10 49.4214
12 50.6175
∞ 51.6938
According to the model, the jumper accelerates rapidly (Fig. 1.2). A velocity of
49.4214 m/s (about 110 mi/hr) is attained after 10 s. Note also that after a sufficiently
FIGURE 1.2
The analytical solution for the bungee jumper problem as computed in Example 1.1. Velocity
increases with time and asymptotically approaches a terminal velocity.
60
Terminal velocity
40
υ, m/s
20
0
0 4 8 12
t, s
1.1 A Simple Mathematical Model 9
long time, a constant velocity, called the terminal velocity, of 51.6983 m/s (115.6 mi/hr)
is reached. This velocity is constant because, eventually, the force of gravity will be in
balance with the air resistance. Thus, the net force is zero and acceleration has ceased.
FIGURE 1.3
The use of a finite difference to approximate the first derivative of υ with respect to t.
υ(ti+1)
True slope
dυ/dt
∆υ
Approximate slope
υ(ti ) ∆υ = υ(ti+1) − υ(ti)
∆t ti+1 − ti
ti ti+1 t
∆t
10 Mathematical Modeling, Numerical Methods, and Problem Solving
υ(ti+1) − υ(ti )
___________ ___ cd 2
ti+1 − ti = g − m υ(ti )
[
υ(ti+1) = υ(ti ) + g − ___
c
]
md υ(ti )2 (ti+1 − ti ) (1.12)
Notice that the term in brackets is the right-hand side of the differential equation itself
[Eq. (1.8)]. That is, it provides a means to compute the rate of change or slope of υ. Thus,
the equation can be rewritten more concisely as
dυ
υi+1 = υi + ___
i Δt (1.13)
dt
where the nomenclature υi designates velocity at time ti, and Δt = ti+1 − ti.
We can now see that the differential equation has been transformed into an equation that
can be used to determine the velocity algebraically at ti+1 using the slope and previous values
of υ and t. If you are given an initial value for velocity at some time ti, you can easily compute
velocity at a later time ti+1. This new value of velocity at ti+1 can in turn be employed to extend
the computation to velocity at ti+2 and so on. Thus at any time along the way,
This approach is formally called Euler’s method. We’ll discuss it in more detail when we
turn to differential equations later in this book.
Solution. At the start of the computation (t0 = 0), the velocity of the jumper is zero.
Using this information and the parameter values from Example 1.1, Eq. (1.12) can be used
to compute velocity at t1 = 2 s:
[ ]
0.25 (0)2 × 2 = 19.62 m/s
υ = 0 + 9.81 − ____
68.1
For the next interval (from t = 2 to 4 s), the computation is repeated, with the result
[ ]
0.25 (19.62)2 × 2 = 36.4137 m/s
υ = 19.62 + 9.81 − ____
68.1
1.1 A Simple Mathematical Model 11
60
Terminal velocity
Approximate,
40 numerical solution
υ, m/s
Exact, analytical
solution
20
0
0 4 8 12
t, s
FIGURE 1.4
Comparison of the numerical and analytical solutions for the bungee jumper problem.
t, s υ, m/s
0 0
2 19.6200
4 36.4137
6 46.2983
8 50.1802
10 51.3123
12 51.6008
∞ 51.6938
The results are plotted in Fig. 1.4 along with the exact solution. We can see that the numeri-
cal method captures the essential features of the exact solution. However, because we have
employed straight-line segments to approximate a continuously curving function, there is
some discrepancy between the two results. One way to minimize such discrepancies is to
use a smaller step size. For example, applying Eq. (1.12) at 1-s intervals results in a smaller
error, as the straight-line segments track closer to the true solution. Using hand calcula-
tions, the effort associated with using smaller and smaller step sizes would make such
numerical solutions impractical. However, with the aid of the computer, large numbers of
calculations can be performed easily. Thus, you can accurately model the velocity of the
jumper without having to solve the differential equation exactly.
12 Mathematical Modeling, Numerical Methods, and Problem Solving
As in Example 1.2, a computational price must be paid for a more accurate numeri-
cal result. Each halving of the step size to attain more accuracy leads to a doubling of the
number of computations. Thus, we see that there is a trade-off between accuracy and com-
putational effort. Such trade-offs figure prominently in numerical methods and constitute
an important theme of this book.
Although Eqs. (1.14) and (1.15) might appear trivially simple, they embody the two funda-
mental ways that conservation laws are employed in engineering and science. As such, they
will form an important part of our efforts in subsequent chapters to illustrate the connection
between numerical methods and engineering and science.
Exploring the Variety of Random
Documents with Different Content
világon, és egy jó házból való leányt nem lehet csak úgy mir nichts,
dir nichts megszöktetni! Szeretném én látni azokat a bírákat, a kik
egy apát könnyekkel a szemében nem hallgatnak meg! És én őt meg
fogom mutatni neked becsukva!
Elza. – Egyébként, bánom is én, akármit csinálsz vele! Egy olyan
szamár, a ki elhiszi, hogy én az ő szép szeméért szököm meg
hazulról, és hogy titokban meg fogok vele esküdni, az ilyen szamár
nem érdemel kíméletet!
A báró. – Egy ilyen szamár! Át fogom adni a rendőrségnek! Hol
van?
Elza. – A negyedik emeleten.
A báró. – A negyedik emeleten?
Elza. – Igen. Csak nem fogok vele egy emeleten aludni?!
Iszonyuan hortyog!
A báró. – Isten, a ki az égben vagy!… Te tudod, hogy ő
hortyog?!…
Elza. – Nemcsak én tudom, az egész Berlin tudja. Ott már nem
akartak neki szobát adni. És igazuk volt, mert az egyszerűen
rémséges. Az első éjjel persze még semmit se sejtettem, és közös
appartement-t vettünk, mert akkor még volt pénzünk. Hát egyszerre
csak hallom a másik szobából ezt a horkolást! Rákopogtam és a zárt
ajtón átkiabáltam hozzá: »Hallja, barátom, ezt nem lehet kiállani!
Menjen le a portáshoz, és nyittasson magának szobát az épület túlsó
részében!« Úgy tett, és másnap a szomszédja iszonyu skandalumot
csinált a hôtelben. Mikor aztán az eset egy másik szomszéddal is
megismétlődött, ki akarták tenni a hôtelből. Na, mondhatom, ezt a
kirándulást nem fogja az élete legszebb napjai között emlegetni! Hát
mikor Berlinből a harmadik osztályon jöttünk Brémába, és hatvanan
ültünk egy waggonban! Én csak kaptam ülőhelyet, de neki majdnem
az egész úton állania kellett.
A báró. – Hallatlan! Hallatlan! A leányom megszökik hazulról és a
harmadik osztályon utazik Berlinből Brémába! Igazán nem tudom,
hogy ez a fej itten az én saját fejem-e, vagy egyszerre egy más
embernek a feje?!… És én ne csukassam be ezt a gazembert, a kit
fel fogok pofozni! De, szerencsétlen gyermek, hát mi jutott
eszedbe?! Remélem, ki fogod magadat magyarázni?!…
Elza. – Azonnal megmagyarázom az egészet. Persze nem fogod
helyeselni, a mit tettem, de legalább meg fogod érteni.
A báró. – Borzasztó! Az én leányom, az én szép, fiatal, kedves
leányom, egy millióval, a mi szintén nem utolsó dolog, arája az egyik
legjelesebb fiatal majoreskónak, egy grófnak, a kit közelebb
legfelsőbb helyről hitbizomány engedélyezésével fognak kitüntetni,
az én okos leányom! – megszökik egy szerkesztővel, a kinek a neve
nem jut eszembe, és a harmadik osztályon utazik Berlinből Brémába!
Ettől meg lehet bolondulni!
Elza. – Mondom, hogy a szerkesztő nem számít. Jobban
szerettem volna egy kocsist vagy egy czigányt, a legjobban egy
hordárt, egy öreg hordárt, de először: senki se hitte volna el, hogy
egy öreg hordárba szerettem bele, s másodszor: csak ez volt kéznél.
És ez legalább mindjárt ki is írhatja az esetet az ujságba. Ne
szörnyűlködjél; igen, én azt kerestem, hogy a botrány minél nagyobb
legyen. Hogy kit használok fel erre a czélra, az mindegy volt. A
szerkesztőről tehát nem érdemes beszélni; ne is említsük többet,
mert különben soha se magyarázhatom ki magamat. Ő itt mellékes;
ő itt nem szerepel.
A báró. – Hát ki szerepel, ha nem ő?!
Elza. – Az, a kit te az egyik legjelesebb fiatal majoreskónak
mondasz! A gróf, a kit legközelebb hitbizománynyal fognak
kitüntetni! A vőlegényem, a pimasz!
A báró. – A Miska! Nem értem. Történt netalán valami, mit jó
atyádnak nem mondtál el?!… Nem értem.
Elza. – Semmi se történt, csak észrevettem, hogy ez a csibész ott
akar hagyni a faképnél.
A báró. – Mi jut eszedbe?! Ezt te csak beképzelted magadnak.
Elza. – Rám bízhatod!… Ha egyszer mondom, hát úgy van! Na, és
felháborodtam. Mit? Ez akar engem a faképnél hagyni, ez akar
engem nevetség és szánalom tárgyává tenni?!… Ez a mitugrász, ez a
barázdabillegető, ez a zsúrhiéna, a ki végigette a várost, és tánczolt
érte, mert ő gróf?! Hiszen emlékszel rá, meddig környékezett,
mennyit ólálkodott körülöttem! Mennyire szeretett volna elhalászni,
és én rá se néztem! Emlékezhetel rá, hogyan nézett ki. Nem volt
egyebe, csak a neve, két rendbeli jó ruhája, tiszta volt az inge,
tudott francziául és mindig röhögött. Sohase bocsátom meg
magamnak azt az éktelen szamárságot, hogy a menyasszonya
lettem. Ha még szerettem volna egy kicsit, mint a Miczut! De nem.
Csak untam már, hogy annyi viszontagság után még mindig pártában
vagyunk; és mert ez a svihák mulatságosan tud káromkodni, meg,
mint az éhenkórászok rendesen, tud mindenféle
svarczkünsztlerséget, bebeszéltettem vele magamnak, hogy ő tetszik
nekem. Gondoltam magamban: »Fene bánja, legalább igazi grófné
leszek.« De mit tesz Isten? Egyszerre csak kiüt a fajtájában a kolera,
s az egész nemzetségéből egyes-egyedül ő marad meg irmagnak.
Másképpen ugyan nem juthatott volna pénzhez, az istenadta, egy
vörös rézgarashoz se, mert ennyit se tud megkeresni, és a mit eddig
megevett, megivott, elszítt, azt mind úgy kapta ajándékba, a
néniktől, a bácsiktól, a magas kormánytól, a Stefi gróftól, a kinek
fölsegítette a kabátját, a vak bárónétól, a kihez eljárt vizitába,
keresztény vagy nem keresztény házigazdáktól, a kiknek a leányai
után járt, no, meg a férjektől, a kiknek a felesége körül settenkedett.
Egy félországnak kellett meghalnia, hogy valamit örökölhessen, egy
sereg embernek kellett a csatatéren maradnia, a fűbe harapnia,
hogy ez a holló jóllakhasson, és a csoda megtörtént. Na, és
megjelent, kikefélkedve, kiszapulva, azzal az öntudatlan
leereszkedéssel, a melynél jobban semmi se sérti az embert, s azzal
a szemtelenül vihogó tekintettel, a mely megmagyarázza az egész
franczia forradalmat. A hernyóból pillangó lett, a smaroczerből
oligarcha. Máról holnapra megnőtt, megférfiasodott, megszépült.
Nem lehetett ráismerni. Mintha nem ő volna a tegnapi potyafráter;
mintha ő rohant volna ki Szigetvárból. A hogy belépett, láttam, hogy
a faképnél akar hagyni. És csakugyan már első intrádára czélozgatott
rá, hogy ez, meg az, késleltetni fogja a házasságunkat. Eleinte igen
komikusnak találtam, aztán iszonyuan szenvedett a hiuságom. Azt
mondtam magamban: nesze neked, nem szeretted, mégis hozzá
mentél volna, most derogálsz neki, és tessék, úgy otthagy, mint egy
elcsábított varróleányt!… nem kellesz neki, hoppon maradsz, ki
fognak nevetni, vagy szánakozni fognak rajtad, ez az Isten
büntetése!… Mert az Istent akkor találjuk meg a leghamarabb, a
mikor büntetni akar bennünket… Végül nem akartam hinni, hogy jól
látok; elhitettem magammal, hogy tévedek. Kelepczéket állítottam
neki, három kelepczét; mind a háromba belement. Először kértem,
hogy ezért, meg azért, távolabbi időpontot tűzzünk ki az esküvőre;
örömmel egyezett bele. Másodszor szerelmet szimuláltam és
sürgettem, hogy határozzuk meg a házasság napját. Elfanyalodott az
arcza és ezer kifogást talált. Harmadszor azt kivántam, ügyesen,
hogy meg ne lássa a vermet: halaszszuk el az esküvőt egy egész
esztendővel. Boldoggá tettem ezzel; nem is nagyon maszkirozta az
örömét. Most már bizonyos voltam. Nem akart, nem is igen tudott
egyszerre szakítani, mert nagyon angazsálta magát, s nem volt
mersze ahhoz, hogy czinikusan végezze a dolgát: félt egy pár
embernek az itéletétől. De föltette magában, hogy lassankint vissza
fog vonulni; leste az alkalmat a szakításra, s közben ügyesen
megkezdte a hátrálást. Talán nem is bánta, ha kitalálom a
szándékát. Ő jó volt nekem, mikor még csak ágrólszakadt grófocska
volt és éhesen szaglálta a pénzemet, de most, hogy urrá lett, én már
nem voltam jó neki. Most már neki is van pénze, nem szorult rám,
nem kellettem neki többé! A pimasz! (Kiköp a padlóra.)
A báró (egy pillanatig halluczinál: úgy tetszik neki, mintha egy rég
elmult században élne, és Szandeczben volna, egy nagyon szegény
lakásban, és hogy az, a kit maga előtt lát, egy ismeretlen, gyönyörű
zsidóleány… Már előre megbocsátott mindent; előre megbocsátotta
azt is, a mi nem történt meg. De most el van ragadtatva; és minden
igyekezete az, hogy elragadtatását palástolja:) – Te, te, te, kis
bolond!…
Elza. – Nem is vagyok annyira bolond! Így szóltam magamban:
»Urfi, te elszántad magadat, hogy így vagy úgy, alkalomadtán, vagy
a nélkül, megragadva valami később kinálkozó ürügyet, vagy
lassankint kopva ki a házból, de okvetetlenül a faképnél fogsz
hagyni, minden jó lélek részvétének adva át azt a szegény teremtést,
a ki, mióta igazi gróffá lettél, már nem méltó hozzád! Na hát, urfi,
abból nem eszel! A mit te akarsz, az nem létezik! Olyan nincs!
Ellenkezőleg, te leszel az, a kit a faképnél fognak hagyni, s ha
valakin nevetni fognak, ez a valaki te vagy!« Igen, de hogyan? Ha
visszaadom a gyürüjét, ha kiadom az útját, ha férjhez megyek
máshoz: ezzel, azonkívül, hogy csak neki szerzek vele örömet,
sohase mosom le magamról, hogy voltaképpen mégis én vagyok az,
a kit otthagytak. Csak egy mód kinálkozott. Csuffá teszem a
menyasszonyát, és ezzel csuffá teszem őt. Így elérem a czélomat.
Azt fogják mondani: a menyasszonya inkább megszökött egy
szegény ördöggel, a kit szeretett, mint hogy férjhezmenjen hozzá.
Megtettem. És ne hidd, hogy nem mértem a hiuságára érzékeny
csapást. Az oligarcha egy kicsit tüsszögni fog. És a maga körében, a
leányos házaknál, egy pár évig, a míg el nem felejtik ezt a históriát,
a kontesszek nem igen fognak veszekedni ezen az épouseur-ön, a ki
elől egy lipótvárosi leány a szégyenbe menekült valami obskúrus
fiatal emberrel! Az én emberem egy darabig nem fog henczegni!
Megboszultam magam!
A báró. – De milyen áron! Igy kompromittálni magadat, szegény
gyermekem!
Elza. – Hja, kényszerhelyzetben voltam! Csak két rossz között
választhattam.
A báró. – Mit fognak mondani az emberek? Nem fogják elhinni,
még csak nem is akarod elmondani nekik az igazat!…
Elza. – Az emberek így csak megszólnak, amúgy lesajnáltak
volna. Szivesebben választom ezt a megoldást. Talán nem bocsátják
meg, hogy megszöktem, de majd elfelejtik; ha hoppon maradok, azt
nem felejtették volna el.
A báró. – És nem sajnáltál bennünket? Ilyen szomoruságot, ilyen
aggodalmakat okozni nekünk!…
Elza. – De látod, milyen hamar előtelegrafáltalak!
A báró. – Szegény kis leányom! Milyen helyre juttatott a
meggondolatlanságod! És nélkülöznöd is kellett!
Elza. – Ó, az roppant mulatságos volt! Különösen tegnap és
tegnapelőtt látni a gavalléromat! Általában a kaland, az sokat ért!
Soha életemben nem mulattam olyan jól, mint Berlintől Brémáig a
harmadik osztályon. Különben hozhattam volna pénzt magammal, de
nem is gondoltam rá. A Rejtvényié pedig hamarabb elfogyott, mint
képzelte. Addig biztattam, hogy: »Csak magas borravalókat adjon!
Csak magas borravalókat adjon!« – mig egyszerre kiderült, hogy
nincsen miből.
A báró. – És ez a szamár, a kit nyakon fogok ütni, mert neki már
benőhetett volna a feje lágya, azt képzelte legyen, hogy te majd
titokban össze fogsz esküdni vele?
Elza. – Azt tökéletesen bolonddá tettem, holmi színésznőktől
levetett mosolyokkal, szemforgatásokkal és sóhajokkal. Még
Berlinben is azt képzelte, hogy imádom.
A báró. – Na, és látod, milyen súlyosak a meggondolatlanság
következései!
Elza. – Hogyhogy?
A báró. – Mert tessék, te nem szereted ezt az embert, azt a kit,
hiába beszélsz, mégis csak meg fogok fenyíteni, és a világ, hogy
reparálva legyen ez a szerencsétlen ügy, azt fogja követelni, hogy
nőül menj hozzá, egy ilyen nem tudom mihez, és nem tudom kihez,
mert a neve nem jut az eszembe.
Elza. – Kihez? Ehhez a tökfilkóhoz? Nem gárgyultam meg!
A báró. – És eljössz vele Brémába! De hát mit kerestél
Brémában?!
Elza. – Elhitettem vele, hogy csak Amerikában lehetünk
biztonságban, mert Európában ránk akadnátok és visszavinnétek.
Olyan messzire akartam elmenni hazulról, hogy maga a hely komoly
színben tüntesse fel a szökést. Ki hisz az olyan szerelemben, a mely
Harasztiba menekül?
A báró. – És pénz nélkül akartatok menni Amerikába?
Elza. – Azt mondtam neki, hogy ide várom a pénzt egy bizalmas
emberemtől, a kit még otthon avattam bele a titokba.
A báró. – Nekem úgy tetszik, ha én, megfelelő fenyiték után,
szubvenczionálom őt, otthon a pletykát még el lehetne némítani.
Elza. – Csakhogy én nem azért jöttem Brémába, hogy a pletyka
elnémuljon, hanem azért, hogy a pletyka beszéljen. Nagyon
csudálnám különben, ha a dolog titok maradna. Mert az oligarchának
megírtam.
A báró. – Csak legalább ezt ne tetted volna!
Elza. – Hisz ez volt a gyönyörüségek gyönyörüsége! Elképzelni,
hogy micsoda pofát vág hozzá!…
A báró. – De azt elképzelted-e, hogy szegény Blanka nővéredre
minő hatást tett ezen esemény?! Sírógörcsöt kapott, és többször
ismételte, természetesen a legnagyobb titokban, csakis kedves
anyád jelenlétében, a cselédség kizárásával, kiknek számára
különben egy rokonokhoz való hirtelen utazást adtunk ürügyül, hogy
most már többé nem fog férjhez mehetni.
É
Elza. – És mit szóltam volna én, mikor Blanka kompromittálta
magát és ezzel kompromittált engem is?! Azt, hogy, ha nem megyek
hamarosan férjhez, ez Blanka miatt van! Nos, most én is
kompromittáltam magam, s Blankával kvittek vagyunk! Legalább
mind a kettőnknek marad egy vigasztalásunk, ha végképpen
hajadonok maradunk: én vádolhatom őt, s ő vádolhat engem!
A báró. – Igazán pedig csakis azon egyéneket lehet vádolni, kik
kombináczióba kerültek; és kedves mamád betegségét, mert ő
bizony, szegény, nem mozoghat eléggé ezen czél érdekében, mint
némely anyák, a kik csakis ennek élnek…
Elza. – És a kik a legmesszebb menő engedményeket teszik a
vőlegényeknek, csakhogy leányaikat főkötő alá juttathassák! De ezt
nem is kívánhatjuk tőle, és én azt kezdem hinni, hogy akármilyen
egészséges ő, és akármennyit fáradozna értünk, az már lehetetlen,
hogy mi férjhez menjünk, mert a végzet nem akarja, s mert ez már
így van megírva a csillagokban!
A báró. – Ne légy ilyen babonás, gyermekem, sőt ellenkezőleg,
szépen kérlek, legyetek szivesek, legalább egyszer, csak egyetlen
egyszer, menjetek már férjhez!… ha a magatok kedvéért nem
akartok, hát az én kedvemért!… Annyi baj van veletek! És ez mind
nem volna, ha egyszer férjhez mentek! Mikor ez eszembe jut, már
azt se bánnám, ha ehhez az izéhez mennél férjhez, a mennyiben
szereted, és ebben az egy esetben, kivételesen nem ragaszkodnám a
megfenyítéséhez.
Elza. – Nem, ehhez nem! És ha ő volna a világ legtökéletesebb
embere, akkor se mennék hozzá, olyan rettenetesen horkol! De ha
szépen megkérsz, valaki máshoz majd csak férjhez megyek. Ha jól
elterjed a szökésem híre, ez divatba hoz, és a kérők nem fognak
hiányozni, majd meglátod! De ha én férjhez megyek a kedvedért,
szépen kérlek, te is tegyél meg nekem egy szivességet.
A báró. – Hogy érted ezt?
Elza. – Szépen kérlek, vigyél el valahová reggelizni, mert
rettenetesen éhes vagyok.
A báró (megrémülve). – Talán csak nem fogytatok ki annyira a
pénzből, hogy már a táplálkozásban is nélkülöznöd kellett volna?!…
Elza. – Még nem, de éppen jó, hogy ma megérkeztél!…
A báró. – Na hát meg kell fenyítenem azt az embert, határozottan
meg kell őt fenyítenem! Hol van? A negyedik emeleten? Hányadik
szám alatt? Pofon fogom őt vágni; egy apának joga van ehhez.
Elza. – Ne bántsd, szegényt! Inkább, ha elviszel reggelizni,
küldjünk neki egy kis köménymagos levest, mert tegnap már nem
evett meleget. De, apuskám, menjünk, mert szörnyen éhes vagyok!
XV.
EGY CSÖNDES ZUGBAN.
– A Berzsenyi-palota úgy ki van világítva, hogy az egész Nádor-utcza
fényben úszik. A Fürdő-utczától az országházig hármas kocsisor. A
portás ünnepi díszben: Dsingiszkán süvegével és jogarral. A fogatok
a kapun úgy járnak be és ki, mint a cselédek a krájzler-bolt
lépcsőjén. Benn nagy estély, czigánynyal. Az uj idők találmányaival
megkorrigált ezeregyéj. Majdnem az összes termek tárva. A
tánczteremben és a körülötte lévő helyiségekben akkora emberár
hömpölyög, a minőt csak nagy hazafiak temetésén látni. Az urak
egymás lábára hágnak és a hölgyek alkalmazkodnak a viszonyokhoz.
Csak a dohányzó mellett van egy csöndes zug: a könyvtárterem,
mely nem esik a nagy útvonalba. Ez a csöndes zug feltűnik
mindazoknak, a kik bizalmasabban óhajtanak beszélgetni. A csöndes
zugban egy órától kettőig a következő párbeszédek hangzanak el: –
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
textbookfull.com