SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
|
|
- Chastity Stafford
- 5 years ago
- Views:
Transcription
1 SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, Abstract The idea of using a symbolic manipulator as the engine for analyzing chemical reactions is developed. With the reactions as the input, it is shown how to symbolically derive the governing differential equations for the concentration of each species, the conservation laws for the scheme, and certain steady states. It is also shown how to use the information derived symbolically to construct, and then execute, a numerical scheme for solving the problem. This is done by linking the symbolic manipulator with a X-window display program that produces a graphical interface relating the data, the numerical solution, and the plotting routines. The routines are demonstrated using the Oregonator model that arises in the study of excitable systems.
2 Introduction A major difficulty in chemical kinetics is having to deal with large systems of nonlinear ordinary differential equations. The size of the systems encountered can be significant, for example, the reaction scheme for the formation of water involves 30 reversible reactions [Yablonskii, 1991 #23], the scheme proposed for methanol pyrolysis consists of 66 reversible reactions [Norton, 1990 #13], and the scheme for ethylene oxidation and pyrolysis has 307 reactions [Dagaut, 1990 #3]. The reduction of such systems using the conservation laws, or an analysis to determine the conditions that would lead some species to an equilibrium state, is therefore limited. In fact, the only viable approach for dealing with such problems is to solve them numerically. The limitation with this is that a numerical solution is confined to specific cases and it is therefore not conductive to a general analysis that determines how the solution depends on the parameters. An alternative approach is to use a symbolic program, like Maple, to study the problem ([Heck, 1993 #25];[Ellis, 1992 #26]). An advantage in doing this is that general, or exact, expressions can be derived. So, for example, the conservation laws of very large kinetic schemes are easy to determine. With such information, when the solution needs to be evaluated, the system that actually has to be solved can be significantly smaller. It is worth pointing out that computing languages like Maple are a compartively recent development and general purpose programs have been commercially available for only about ten years. Consequently earlier studies, such as the one by ([Szamosi, 1984 #28]), rely on more conventional floating point calculations. Stated in general terms, the objective of this paper is to demonstrate how to use a symbolic manipulator as an interface between the mathematical formulation of a physical problem and the determination of the quantitative characteristics of the solution. Of particular interest here is how to use a symbolic manipulator as the engine for analyzing chemical reactions. For the algorithms that are described, all that is needed are the reactions. For example, the Oregonator model used in the study of excitable systems consists of the five reactions [Field, 1974 #27] A + Y X + P, X + Y 2P, A + X 2X + 2Z, (1) 2X A + P, Z 1 2 Y. 2 9/29/04
3 Given this information the program can symbolically derive the governing differential equations for the concentration of each species, the conservation laws for the scheme, and certain steady states. It can also solve the equations numerically and then display the results in a variety of formats (e.g., as a function of time, in a two-dimensional phase plane, in three-dimensional phase space, etc.). Because of the ubiquitous nature of the laws of mass action, such capabilities are useful both for teaching as well as research. It is important to point out that some of the problems discussed here can be solved by other means using commercially available software, e.g., DASAL and CRAMS ([Marsili, 1990 #29]). Also, symbolic programs can be somewhat slower than routines that are run using programming languages like FORTRAN or C. However, what is gained by using a symbolic system is an exact result so the calculations need only be done once (unlike a FORTRAN or C program that must be rerun every time the parameters are changed). This is particularly important when one wants to determine the rate parameters from experiment. Also, it is easier to understand the physics since one obtains an answer that depends on the material parameters in a general way. For complex problems, carrying out a symbolic analysis may not produce the solution but it can result in a significant reduction in the system. This is valuable as it can simplify the problem that then must be solved numerically. Consequently, symbolic manipulators provide an important tool for analyzing chemical kinetics problems that can be used in conjunction with classical scientific computing methods. In the next section the theory from chemical kinetics, on which the algorithms are based, is discussed. In conjunction with this, the basic constituents of the symbolic program are described. To help demonstrate what is involved, certain of the commands used in the program are written out for the Oregonator model in (1). After this the X- window graphical interface used to numerically solve the differential equations is described. In the last section, as an example, the complete program is demonstrated using the Oregonator model. Chemical Kinetics and Symbolic Computing For the general form of the schemes considered here we assume there are n reactions involving m distinct species X 1, X 2,..., X m. Therefore, the scheme is composed of individual reactions of the form m α ij X j j=1 m β ij X j for i = 1,..., n. (2) j=1 3 9/29/04
4 It is not assumed that these are elementary reactions, but we do assume that the stoichiometric coefficients (i.e., the α ij 's and β ij 's ) are non-negative. Given the reaction scheme the first step in the symbolic algorithm is the determination of the species involved. For example, the input for the scheme in (1) is a list of reactions such as the following > R := [ A + Y > X + P, X + Y > 2*P, A + X > 2*X + 2*Z, 2*X > A + P, Z > 0.5*Y ]; (3) Identifying the species that are present is simply a matter of searching the input list, R, to find the indeterminants. In Maple this is accomplished by issuing the command > species := indets(r); species := {X, P, A, Y, Z} To determine the number of species (i.e., the value of m in (2)) the following command is used > m := nops(species); 5 With these results, m = 5 and species[1] = X, species[2] = P,..., species[m] = Z. Once the species are known then stoichiometric coefficient product and reactant matrices, M p and M r, are constructed that, when taken together, are algebraically equivalent to the reaction scheme. Assuming there are n reactions and m species, as in (2), then the nxm matrix M p has i,j-entry β ij and the nxm matrix M r has i,j-entry α ij. Determining these coefficients is a sequential process that involves running through the list of reactions, R, and identifying the coefficients of the reactants and products. As an example, for the Oregonator model given in (1) the coefficient matrices are X P A Y Z X P A Y Z M r = M p = /29/04
5 To help identify where the entries of the matrices originate, the species associated with each column are included in the above expressions for M r and M p. In general, the ith row of each matrix corresponds to the ith reaction and the jth column is connected with the jth species. What is important here is that given the two stoichiometric matrices then it is possible to uniquely reconstruct the reaction scheme. This is not possible if one uses only the stoichiometric matrix M s, which is defined as M s = M p M r. (4) For example, the stoichiometric matrix M s = ( 1 1 1) applies to the scheme consisting of the single reaction A + 2B B + C. However, it also applies to the scheme consisting of the single reaction A + B C. Once M p and M r are known it is then an easy matter to determine the rate r i of the ith reaction. To state the formula obtained using the laws of mass action, consider the system given in (2). If the ith row of M r contains only zeros then r i = 0, otherwise m r i = k i j=1 x j M r [i,j]. (5) In this expression, k i is the rate constant for the ith reaction and M r [i,j] is the i,j-th element of M r. Also, x j designates the concentration of species X j. The Maple commands used to implement the formula in (5) are the following > for i from 1 to n do > r.i := k.i: > for j from 1 to m do > if (Mr[i,j]<>0) then r.i := r.i*(x[j]^mr[i,j]) fi: > od: > od: In the above statements, r.i, k.i, Mr[i,j] and x[j] designate r i, k i, M r [i,j] and x j, respectively. Once the rates have been determined then the kinetic equation for the time evolution of the concentration of the jth species is n d dt x j = i=1 (M p [i,j] M r [i,j])r i. (6) 5 9/29/04
6 The construction of the reaction rates and subsequent determination of the differential equations is simply an iteration process that exhausts all the possibilities in the scheme using the formulas in (5) and (6). The coefficient matrices can also be used to determine the conservation laws for the system [Holmes, 1991 #24]. These can be found by finding the linearly independent solutions of the equation M s c = 0, (7) where M s is the stoichiometric matrix given in (4). In other words, if c = (c 1,..., c m ) T is a nonzero solution of (7) then a conservation law for the system is m c j x j = constant. (8) j=1 The determination of the linearly independent solutions of (7) is equivalent to finding a basis for the kernel, or null space, of the stoichiometric matrix defined in (4). Using Maple, the command to find a basis for the kernel of M s is > kernel( Ms, 'nlaws' ); The output from this command is a listing of the linearly independent solutions of (7) and the value of nlaws, which is the number of solutions found. The rest of the procedure for finding the conservation laws consists of a few lines of code for writing out the laws as given in (8). As a final comment, in physical chemistry conservation laws are sometimes restricted to cases of when the c i 's are non-negative [Érdi, 1989 #20], although, this condition is not universal [Othmer, 1985 #21]. The above algorithm allows for negative coefficients in the conservation laws. The determination of equilibrium states, and obtaining analytic expressions that involve the rate parameters, are tasks for which symbolic computing is ideally suited. There is more information in these analytic expressions, if they can be obtained, than in a numerical approximation. A case for this has already been made by Kreye, et al. (1988), where they use SMP to find rate constants for linear reaction schemes [Kreye, 1988 #7], and by Liddell, et al. (1990), who find rate constants for certain types of schemes using equilibrium states [Liddell, 1990 #10]. Just how successful the program is in this endeavor, however, depends on the particular problem since the reaction schemes can produce highly nonlinear systems that may be unsolvable in closed form. The Oregonator 6 9/29/04
7 model in (1) is such an example. Even so it still possible to reduce the system considerably. To illustrate how, suppose the system is known to approach an equilibrium state. In this case, x j 0 as t, irrespective of the initial concentrations, whenever the scheme contains an equation of the form d dt x i = kx j α, (9) where α is positive and k is a rate constant. It is a simple matter to find all such reactions using the stoichiometric coefficient matrices. This is done by first finding which columns of M s have a single nonzero entry and then checking the appropriate row of M r to determine whether the rate depends on a single species. The algorithm does this and then sets the respective concentration to zero. It then goes back to see if any new species with a zero equilibrium state can be found. This process is continued until no new species are obtained. To summarize the above discussion, the symbolic procedures in the kinetics program, and their functions, are as follows: 1) procedure "odes" a) determine the species and number of reactions from an input list containing the reactions; b) loop through the list and construct the stoichiometric coefficient matrices; c) determine the rate of each reaction; c) construct the kinetic equations; 2) procedure "laws" a) find the system's linearly independent conservation laws; 3) procedure "steady" a) find the zero equilibrium states The program, which is written for Maple, is approximately 550 lines long (this includes the help pages as well as the graphical interface described below). It is available through anonymous ftp at Internet site ftp.rpi.edu, and the file to download is pub/math/kinetics.tar.z. Interface Between Symbolic and Numerical Computing Even though symbolic manipulators, such as Maple and Mathematica, have the capability of finding exact solutions to certain types of differential equations, the nonlinear nature of most reactions makes this a remote possibility. This brings up the issue of how to employ the information derived symbolically to construct, and then execute, a numerical 7 9/29/04
8 scheme. The easiest approach is to use the numerical solver available within the symbolic program. In the case of Maple this involves the Runge-Kutta-Fehlberg (RKF45) method, although other methods, such as a seventh-eighth order continuous Runge-Kutta method, are also available ([Burden, 1993 #30]). One can use other algorithms but this aspect of the procedure is not of concern here as numerical methods for initial value problems are well studied [Gear, 1971 #22]. The real difficulty is keeping track of the information needed to compute the solution, particularly for large reaction schemes. Our approach is to link the symbolic manipulator with system files that contain the necessary data to compute the solution (e.g., rate constants, initial concentrations, etc.). To facilitate the handling of the data, and to make displaying the results as simple as possible, a Motif X-window display program has been written. This program is launched from within Maple and it produces a graphical interface relating the data, the numerical solution, and the plotting routines. An example of the window that comes up is shown in Fig. 1. The data window has entries for all the parameters needed to compute the solution, such as the initial concentration for each species and the rate constants (the latter are labeled k1, k2, k3,... in the program). The user can also specify the time interval over which the problem is solved and the number of points used in computing the solution. On the left side of the window are toggle buttons that enable the user to specify how the solution is displayed. The possibilities include: i) a 2-D plot showing one or more species as a function of time, ii) a 2-D plot showing one or more species as a function of one of the other species, and iii) a 3-D plot where X and Y are specified species (or time) and there are one or more Z selections. The user can also specify the color of each solution curve using a pop-up menu. Finally, there are four buttons for saving, plotting, quitting, and help. The display program waits until the user enters the required data and then presses the plot button. At this point all of the data is collected from the field and buttons, and the interface program writes a Maple procedure using this information. Within the procedure the system of differential equations, and associated initial conditions, are defined and numerically solved using the "dsolve" command. After this the procedure creates a PLOT or PLOT3D Maple structure based on the solution and the plotting data provided by the interface. Once the procedure is complete the display program writes the code to a file that Maple reads and then executes. The X-window commands used to display the data window make use of the Motif X windows libraries and a variety of widgets for the different fields in the data window. 8 9/29/04
9 The interface was written and tested with the GNU C compiler and should be widely portable across most machines running UNIX with the Motif X window libraries. Example To demonstrate the program, a sample session is given below using the Oregonator model in (1). The session was carried out using Maple V, Release 3, on a SPARCstation 2 with 16 MBytes of RAM and running X-windows. The Maple commands begin with a > and each command is preceded by a comment line that indicates what the command does (the comment lines begin with a # ). Also, where appropriate, the response from Maple follows the command. # Read in the program... > read kinetics; # Enter the Oregonator reaction scheme... > R := [ A + Y > X + P, X + Y > 2*P, A + X > 2*X + 2*Z, 2*X > A + P, Z > 0.5*Y ]; R := [X + P < A + Y, 2 P < X + Y, 2 X + 2 Z < A + X, A + P < 2 X,.5 Y < Z] # Determine the differential equations for the scheme... > odes(); The differential equations are: d X(t) = k1 A Y - k2 X Y + k3 X A - 2 k4 X dt d P(t) = k1 A Y + 2 k2 X Y + k4 X dt d A(t) = - k1 A Y - k3 X A + k4 X dt 9 9/29/04
10 d ---- Y(t) = - k1 A Y - k2 X Y +.5 k5 Z dt d ---- Z(t) = 2 k3 X A - k5 Z dt # Determine the conservation laws... > laws(); The conservation laws are: 1. constant = 6.*X+4.*P+8.*A+2.*Y+Z # Check to see what species can have a zero steady state... > steady(); No species were found to have a zero steady state. # Compute, and then plot, the solution... > ksolve(); In response to the command "ksolve" the plot data window shown in Fig. 1 appears. When the window first comes up there are no entries in the Initial Value and Rate Value columns, and there are no selections in the X, Y, Z columns. The labels X, Y, Z are used here only to indicate the conventional cartesian coordinates (and not specific chemical species). In any case, it is now necessary to specify the initial values for each species. For this example, we take X(0) = 1, P(0) = 0, A(0) = 0.1, Y(0) = 1, and Z(0) = 1. The rate constants to be used in this example are k 1 = 100, k 2 = 1, k 3 = 100, k 4 = 1, and k 5 = Also, the solution is to be calculated at 200 points over the time interval 0 < t 40 (time in this example is measured in seconds). This information, as it would appear in the data window, is shown in Fig. 1. The first plot to be constructed gives Y(t) and Z(t) as functions of time t. To accomplish this, T is selected in the X column and both Y and Z are selected in the Y column (see Fig. 1). By default the window assigns a color to each species and these are 10 9/29/04
11 (in order): white, red, green, yellow, orange, blue, and violet. In the Color column it is possible to reassign the color for each curve. For this example, the curve for Y is to be drawn in green and the one for Z in orange. After entering this information the "Plot Graph" button is pushed and the result is shown in Fig. 2. This plot is done by Maple and there are several options on the plot window, which is not shown, for manipulating the plot (e.g., line styles, type of axes, printing, etc.). It is just as easy to produce 3-D plots using the data window. To illustrate, suppose one wants to plot the curve consisting of the points (X(t), Y(t), A(t)) for 0 < t 40. In regard to Fig. 1, out of the X column select Y, out of the Y column select Z and out of the Z column select A. The curve is also to drawn in violet and this is indicated in the Color column. The other aspects of the data window are not changed. Again, the "Plot Graph" button is pushed and the result is shown in Fig. 3. It is worth mentioning that the symbolic procedures used here are relatively fast. In particular, for the Oregonator model the "odes" command takes about 0.2 seconds, the "laws" command about a second, and the "steady" command less than 0.05 seconds. Moreover, these procedures will work on any computer system which has Maple. The "ksolve" procedure, because of its graphical nature, requires X-windows and the Motif X Windows libraries. Acknowledgments: We would like to thank the URP at Rensselaer for their support of this work. 11 9/29/04
12 References 12 9/29/04
13 Figure Captions Fig. 1 The X-window launched from within Maple. This window is the graphical interface used for entering data for plotting the solution of the kinetic equations. Fig. 2 Solution curves for Y(t) and Z(t) from the Oregonator model using the data values shown in Fig. 1. As indicated in Fig. 1, Y is drawn in green and Z in orange. Fig. 3 Solution curve (Y, Z, A), for 0 < t 40, from the Oregonator model using the data values shown in Fig /29/04
Athena Visual Software, Inc. 1
Athena Visual Studio Visual Kinetics Tutorial VisualKinetics is an integrated tool within the Athena Visual Studio software environment, which allows scientists and engineers to simulate the dynamic behavior
More informationComputational Study of Chemical Kinetics (GIDES)
Computational Study of Chemical Kinetics (GIDES) Software Introduction Berkeley Madonna (http://www.berkeleymadonna.com) is a dynamic modeling program in which relational diagrams are created using a graphical
More informationLAB 1: MATLAB - Introduction to Programming. Objective:
LAB 1: MATLAB - Introduction to Programming Objective: The objective of this laboratory is to review how to use MATLAB as a programming tool and to review a classic analytical solution to a steady-state
More informationMath 314 Lecture Notes Section 006 Fall 2006
Math 314 Lecture Notes Section 006 Fall 2006 CHAPTER 1 Linear Systems of Equations First Day: (1) Welcome (2) Pass out information sheets (3) Take roll (4) Open up home page and have students do same
More informationMatrices A matrix is a rectangular array of numbers. For example, the following rectangular arrays of numbers are matrices: 2 1 2
Matrices A matrix is a rectangular array of numbers For example, the following rectangular arrays of numbers are matrices: 7 A = B = C = 3 6 5 8 0 6 D = [ 3 5 7 9 E = 8 7653 0 Matrices vary in size An
More informationBCMB/CHEM 8190 Lab Exercise Using Maple for NMR Data Processing and Pulse Sequence Design March 2012
BCMB/CHEM 8190 Lab Exercise Using Maple for NMR Data Processing and Pulse Sequence Design March 2012 Introduction Maple is a powerful collection of routines to aid in the solution of mathematical problems
More information4.4 Computing π, ln 2 and e
252 4.4 Computing π, ln 2 and e The approximations π 3.1415927, ln 2 0.69314718, e 2.7182818 can be obtained by numerical methods applied to the following initial value problems: (1) y = 4, 1 + x2 y(0)
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK
ES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK What is SIMULINK? SIMULINK is a software package for modeling, simulating, and analyzing
More informationAMS 27L LAB #8 Winter 2009
AMS 27L LAB #8 Winter 29 Solving ODE s in Matlab Objectives:. To use Matlab s ODE Solvers 2. To practice using functions and in-line functions Matlab s ODE Suite Matlab offers a suite of ODE solvers including:
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
More informationLAB 2 - ONE DIMENSIONAL MOTION
Name Date Partners L02-1 LAB 2 - ONE DIMENSIONAL MOTION OBJECTIVES Slow and steady wins the race. Aesop s fable: The Hare and the Tortoise To learn how to use a motion detector and gain more familiarity
More informationUsing Spreadsheets to Teach Engineering Problem Solving: Differential and Integral Equations
Session 50 Using Spreadsheets to Teach Engineering Problem Solving: Differential and Integral Equations James P. Blanchard University of Wisconsin - Madison ABSTRACT Spreadsheets offer significant advantages
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationSymbolic Solution of higher order equations
Math 216 - Assignment 4 - Higher Order Equations and Systems of Equations Due: Monday, April 16. Nothing accepted after Tuesday, April 17. This is worth 15 points. 10% points off for being late. You may
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationChem 1 Kinetics. Objectives. Concepts
Chem 1 Kinetics Objectives 1. Learn some basic ideas in chemical kinetics. 2. Understand how the computer visualizations can be used to benefit the learning process. 3. Understand how the computer models
More informationVCell Tutorial. Building a Rule-Based Model
VCell Tutorial Building a Rule-Based Model We will demonstrate how to create a rule-based model of EGFR receptor interaction with two adapter proteins Grb2 and Shc. A Receptor-monomer reversibly binds
More information1-D Convection-Diffusion Lab
Computational Fluid Dynamics -D Convection-Diffusion Lab The lab. uses scientificworkplace symbolic calculus and maths editor software (SWP) This file Concevtion-Diffusion-Lab is available from Blackboard
More informationMATLAB BASICS. Instructor: Prof. Shahrouk Ahmadi. TA: Kartik Bulusu
MATLAB BASICS Instructor: Prof. Shahrouk Ahmadi 1. What are M-files TA: Kartik Bulusu M-files are files that contain a collection of MATLAB commands or are used to define new MATLAB functions. For the
More informationMatlab Section. November 8, 2005
Matlab Section November 8, 2005 1 1 General commands Clear all variables from memory : clear all Close all figure windows : close all Save a variable in.mat format : save filename name of variable Load
More informationMAT300/500 Programming Project Spring 2019
MAT300/500 Programming Project Spring 2019 Please submit all project parts on the Moodle page for MAT300 or MAT500. Due dates are listed on the syllabus and the Moodle site. You should include all neccessary
More informationLab 1 Uniform Motion - Graphing and Analyzing Motion
Lab 1 Uniform Motion - Graphing and Analyzing Motion Objectives: < To observe the distance-time relation for motion at constant velocity. < To make a straight line fit to the distance-time data. < To interpret
More informationVectors and Vector Arithmetic
Vectors and Vector Arithmetic Introduction and Goals: The purpose of this lab is to become familiar with the syntax of Maple commands for manipulating and graphing vectors. It will introduce you to basic
More informationA First Course on Kinetics and Reaction Engineering Example 1.4
Example 1.4 Problem Purpose This example illustrates the process of identifying reactions that are linear combinations of other reactions in a set and eliminating them until a mathematically independent
More informationSTAT 3008 Applied Regression Analysis Tutorial 1: Introduction LAI Chun Hei
STAT 3008 Applied Regression Analysis Tutorial 1: Introduction LAI Chun Hei Department of Statistics, The Chinese University of Hong Kong 1 Mathematical Background In this courses, some theorems from elementary
More informationAn IDL Based Image Deconvolution Software Package
An IDL Based Image Deconvolution Software Package F. Városi and W. B. Landsman Hughes STX Co., Code 685, NASA/GSFC, Greenbelt, MD 20771 Abstract. Using the Interactive Data Language (IDL), we have implemented
More informationA Brief Introduction To. GRTensor. On MAPLE Platform. A write-up for the presentation delivered on the same topic as a part of the course PHYS 601
A Brief Introduction To GRTensor On MAPLE Platform A write-up for the presentation delivered on the same topic as a part of the course PHYS 601 March 2012 BY: ARSHDEEP SINGH BHATIA arshdeepsb@gmail.com
More informationTOPIC 2 Computer application for manipulating matrix using MATLAB
YOGYAKARTA STATE UNIVERSITY MATHEMATICS AND NATURAL SCIENCES FACULTY MATHEMATICS EDUCATION STUDY PROGRAM TOPIC 2 Computer application for manipulating matrix using MATLAB Definition of Matrices in MATLAB
More informationMULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS. Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University
MULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University HU-1117 Budapest, Pázmány Péter sétány 1/c. karolyik@cs.elte.hu
More informationWILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
More informationAppendix B Microsoft Office Specialist exam objectives maps
B 1 Appendix B Microsoft Office Specialist exam objectives maps This appendix covers these additional topics: A Excel 2003 Specialist exam objectives with references to corresponding material in Course
More informationGreen s Functions with Reflection
Green s Functions with Reflection User s manual Alberto Cabada Fernández (USC) José Ángel Cid Araújo (UVIGO) Fernando Adrián Fernández Tojo (USC) Beatriz Máquez Villamarín (USC) Universidade de Santiago
More informationHOW TO USE MIKANA. 1. Decompress the zip file MATLAB.zip. This will create the directory MIKANA.
HOW TO USE MIKANA MIKANA (Method to Infer Kinetics And Network Architecture) is a novel computational method to infer reaction mechanisms and estimate the kinetic parameters of biochemical pathways from
More informationFirst Order Rate Laws and the Kinetics Mechanism Simulation Program
First Order Rate Laws and the Kinetics Mechanism Simulation Program The Kinetics Mechanism Simulation program (kine.html) is useful for building models in a wide variety of areas, including population
More informationKinetics of Crystal Violet Bleaching
Kinetics of Crystal Violet Bleaching Authors: V. C. Dew and J. M. McCormick* From Update March 12, 2013 with revisions Nov. 29, 2016 Introduction Chemists are always interested in whether a chemical reaction
More informationSolving Systems of Equations Row Reduction
Solving Systems of Equations Row Reduction November 19, 2008 Though it has not been a primary topic of interest for us, the task of solving a system of linear equations has come up several times. For example,
More informationClassroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft
Introduction Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft Maple has the ability to provide various series expansions and
More informationMATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction
MATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction Differential equations are a convenient way to express mathematically a change of a dependent variable (e.g. concentration
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationLinear Algebra. The analysis of many models in the social sciences reduces to the study of systems of equations.
POLI 7 - Mathematical and Statistical Foundations Prof S Saiegh Fall Lecture Notes - Class 4 October 4, Linear Algebra The analysis of many models in the social sciences reduces to the study of systems
More informationChemical Kinetics I: The Dry Lab. Up until this point in our study of physical chemistry we have been interested in
Chemical Kinetics I: The Dry Lab Up until this point in our study of physical chemistry we have been interested in equilibrium properties; now we will begin to investigate non-equilibrium properties and
More informationMichælmas 2012 Operations Research III/IV 1
Michælmas 2012 Operations Research III/IV 1 An inventory example A company makes windsurfing boards and has orders for 40, 60, 75 and 25 boards respectively over the next four production quarters. All
More informationA SAS/AF Application For Sample Size And Power Determination
A SAS/AF Application For Sample Size And Power Determination Fiona Portwood, Software Product Services Ltd. Abstract When planning a study, such as a clinical trial or toxicology experiment, the choice
More informationTP 1: Euler s Algorithm-Air Resistance-Introduction to Fortran
TP 1: Euler s Algorithm-Air Resistance-Introduction to Fortran December 10, 2009 1 References N.J.Giordano, Computational Physics. R.H.Landau, M.J.Paez, C.C.Bordeianu, Computational Physics. H.Gould, J.Tobochnick,
More informationUsing web-based Java pplane applet to graph solutions of systems of differential equations
Using web-based Java pplane applet to graph solutions of systems of differential equations Our class project for MA 341 involves using computer tools to analyse solutions of differential equations. This
More informationModeling Basics: 4. Numerical ODE Solving In Excel 5. Solving ODEs in Mathematica
Modeling Basics: 4. Numerical ODE Solving In Excel 5. Solving ODEs in Mathematica By Peter Woolf University of Michigan Michigan Chemical Process Dynamics and Controls Open Textbook version 1.0 Creative
More informationMatrix Assembly in FEA
Matrix Assembly in FEA 1 In Chapter 2, we spoke about how the global matrix equations are assembled in the finite element method. We now want to revisit that discussion and add some details. For example,
More informationMaple in Differential Equations
Maple in Differential Equations and Boundary Value Problems by H. Pleym Maple Worksheets Supplementing Edwards and Penney Differential Equations and Boundary Value Problems - Computing and Modeling Preface
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationBloomsburg University Weather Viewer Quick Start Guide. Software Version 1.2 Date 4/7/2014
Bloomsburg University Weather Viewer Quick Start Guide Software Version 1.2 Date 4/7/2014 Program Background / Objectives: The Bloomsburg Weather Viewer is a weather visualization program that is designed
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationAlgebra & Trig. I. For example, the system. x y 2 z. may be represented by the augmented matrix
Algebra & Trig. I 8.1 Matrix Solutions to Linear Systems A matrix is a rectangular array of elements. o An array is a systematic arrangement of numbers or symbols in rows and columns. Matrices (the plural
More informationAMS 27L LAB #6 Winter 2009
AMS 27L LAB #6 Winter 2009 Symbolically Solving Differential Equations Objectives: 1. To learn about the MATLAB Symbolic Solver 2. To expand knowledge of solutions to Diff-EQs 1 Symbolically Solving Differential
More informationA Symbolic Numeric Environment for Analyzing Measurement Data in Multi-Model Settings (Extended Abstract)
A Symbolic Numeric Environment for Analyzing Measurement Data in Multi-Model Settings (Extended Abstract) Christoph Richard 1 and Andreas Weber 2? 1 Institut für Theoretische Physik, Universität Tübingen,
More informationEEL2216 Control Theory CT1: PID Controller Design
EEL6 Control Theory CT: PID Controller Design. Objectives (i) To design proportional-integral-derivative (PID) controller for closed loop control. (ii) To evaluate the performance of different controllers
More informationLinear Systems of Differential Equations
Chapter 5 Linear Systems of Differential Equations Project 5. Automatic Solution of Linear Systems Calculations with numerical matrices of order greater than 3 are most frequently carried out with the
More informationJointly Learning Python Programming and Analytic Geometry
Jointly Learning Python Programming and Analytic Geometry Cristina-Maria Păcurar Abstract The paper presents an original Python-based application that outlines the advantages of combining some elementary
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationLaboratory 1. Solving differential equations with nonzero initial conditions
Laboratory 1 Solving differential equations with nonzero initial conditions 1. Purpose of the exercise: - learning symbolic and numerical methods of differential equations solving with MATLAB - using Simulink
More informationMath 308 Week 8 Solutions
Math 38 Week 8 Solutions There is a solution manual to Chapter 4 online: www.pearsoncustom.com/tamu math/. This online solutions manual contains solutions to some of the suggested problems. Here are solutions
More informationNumerical mathematics with GeoGebra in high school
Herceg 2009/2/18 23:36 page 363 #1 6/2 (2008), 363 378 tmcs@inf.unideb.hu http://tmcs.math.klte.hu Numerical mathematics with GeoGebra in high school Dragoslav Herceg and Ðorđe Herceg Abstract. We have
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationStoichiometric Reactor Simulation Robert P. Hesketh and Concetta LaMarca Chemical Engineering, Rowan University (Revised 4/8/09)
Stoichiometric Reactor Simulation Robert P. Hesketh and Concetta LaMarca Chemical Engineering, Rowan University (Revised 4/8/09) In this session you will learn how to create a stoichiometric reactor model
More informationSoPlex. and. Thorsten Koch, ZIB
SoPlex and ZIMPL Thorsten Koch, ZIB Linear program min c T x subject to Ax b with c R n x R n A R k n b R k and k n and A has full rank Linear program max c T x subject to Ax b x 0 with c R n, x R n, A
More informationLinear and nonlinear optimization with CasADi
Numerical optimization course computer exercise inear and nonlinear optimization with CasADi Prof. Dr. Moritz Diehl 1, Prof. Dr. Angelika Altmann-Dieses 2, Adrian Bürger 1,2 1 Systems Control and Optimization
More informationSolving of logic functions systems using genetic algorithm
Solving of logic functions systems using genetic algorithm V G Kurbanov,2 and M V Burakov Chair of control system of Saint-Petersburg State University of Aerospace Instrumentation, Bolshaya Morskaya, 67,
More informationTask 1: Open ArcMap and activate the Spatial Analyst extension.
Exercise 10 Spatial Analyst The following steps describe the general process that you will follow to complete the exercise. Specific steps will be provided later in the step-by-step instructions component
More informationVirtual Cell Membrane Potential Tutorial IV
Virtual Cell Membrane Potential Tutorial IV Creating the BioModel Creating the Application!" Application I -Studying voltage changes in a compartmental model!" Application II - Studying voltage, sodium,
More informationIntroduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Spring 2009
Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationThe Stoichiometry of Reactions Introduction Copyright c 2018 by Nob Hill Publishing, LLC
1 / 70 The Stoichiometry of Reactions Introduction Copyright c 2018 by Nob Hill Publishing, LLC Stoichiometry: the determination of the proportions in which chemical elements combine or are produced and
More informationAppendix F: Projecting Future Sea Level Rise with the SLRRP Model
Appendix F: Projecting Future Sea Level Rise with the SLRRP Model The Sea Level Rise Rectification Program (SLRRP) is a software package designed with a user-friendly interface to generate a suite of future
More informationThe Stoichiometry of Reactions Introduction
The Stoichiometry of Reactions Introduction Copyright c 2015 by Nob ill Publishing, LLC Stoichiometry: the determination of the proportions in which chemical elements combine or are produced and the weight
More informationFin System, Inc. Company Report. Temperature Profile Calculators. Team 1 J. C. Stewards, Lead A. B. Williams, Documentation M. D.
Fin System, Inc. Company Report Temperature Profile Calculators Team 1 J. C. Stewards, Lead A. B. Williams, Documentation M. D. Daily, Programmer Submitted in Fulfillment of Management Requirements August
More informationa11 a A = : a 21 a 22
Matrices The study of linear systems is facilitated by introducing matrices. Matrix theory provides a convenient language and notation to express many of the ideas concisely, and complicated formulas are
More informationPurpose: Materials: WARNING! Section: Partner 2: Partner 1:
Partner 1: Partner 2: Section: PLEASE NOTE: You will need this particular lab report later in the semester again for the homework of the Rolling Motion Experiment. When you get back this graded report,
More informationDouble Inverted Pendulum (DBIP)
Linear Motion Servo Plant: IP01_2 Linear Experiment #15: LQR Control Double Inverted Pendulum (DBIP) All of Quanser s systems have an inherent open architecture design. It should be noted that the following
More informationCS 246 Review of Linear Algebra 01/17/19
1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector
More informationTRANSPORTATION PROBLEMS
Chapter 6 TRANSPORTATION PROBLEMS 61 Transportation Model Transportation models deal with the determination of a minimum-cost plan for transporting a commodity from a number of sources to a number of destinations
More information9.1 - Systems of Linear Equations: Two Variables
9.1 - Systems of Linear Equations: Two Variables Recall that a system of equations consists of two or more equations each with two or more variables. A solution to a system in two variables is an ordered
More information1. Open polymath: 2. Go to Help, Contents F1 or Press F1
Polymath Tutorial Process Fluid Transport 1. Open polymath: 2. Go to Help, Contents F1 or Press F1 1 3. Read the section titled Introduction to Polymath both getting started and Variables and expressions
More informationPreparing Spatial Data
13 CHAPTER 2 Preparing Spatial Data Assessing Your Spatial Data Needs 13 Assessing Your Attribute Data 13 Determining Your Spatial Data Requirements 14 Locating a Source of Spatial Data 14 Performing Common
More informationGaussian Elimination and Back Substitution
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 4 Notes These notes correspond to Sections 31 and 32 in the text Gaussian Elimination and Back Substitution The basic idea behind methods for solving
More informationUsing SkyTools to log Texas 45 list objects
Houston Astronomical Society Using SkyTools to log Texas 45 list objects You can use SkyTools to keep track of objects observed in Columbus and copy the output into the Texas 45 observation log. Preliminary
More informationLinear Motion with Constant Acceleration
Linear Motion 1 Linear Motion with Constant Acceleration Overview: First you will attempt to walk backward with a constant acceleration, monitoring your motion with the ultrasonic motion detector. Then
More information1 One-Dimensional, Steady-State Conduction
1 One-Dimensional, Steady-State Conduction 1.1 Conduction Heat Transfer 1.1.1 Introduction Thermodynamics defines heat as a transfer of energy across the boundary of a system as a result of a temperature
More informationNONLINEAR EQUATIONS AND TAYLOR S THEOREM
APPENDIX C NONLINEAR EQUATIONS AND TAYLOR S THEOREM C.1 INTRODUCTION In adjustment computations it is frequently necessary to deal with nonlinear equations. For example, some observation equations relate
More informationSpan and Linear Independence
Span and Linear Independence It is common to confuse span and linear independence, because although they are different concepts, they are related. To see their relationship, let s revisit the previous
More informationIntroduction. How to use this book. Linear algebra. Mathematica. Mathematica cells
Introduction How to use this book This guide is meant as a standard reference to definitions, examples, and Mathematica techniques for linear algebra. Complementary material can be found in the Help sections
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationIntroduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Fall 2009
Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an
More informationFIT100 Spring 01. Project 2. Astrological Toys
FIT100 Spring 01 Project 2 Astrological Toys In this project you will write a series of Windows applications that look up and display astrological signs and dates. The applications that will make up the
More informationMath 409/509 (Spring 2011)
Math 409/509 (Spring 2011) Instructor: Emre Mengi Study Guide for Homework 2 This homework concerns the root-finding problem and line-search algorithms for unconstrained optimization. Please don t hesitate
More informationMATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.
MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If
More informationMATLAB TOOL FOR IDENTIFICATION OF NONLINEAR SYSTEMS
MATLAB TOOL FOR IDENTIFICATION OF NONLINEAR SYSTEMS M. Kalúz, Ľ. Čirka, M. Fikar Institute of Information Engineering, Automation, and Mathematics, FCFT STU in Bratislava Abstract This contribution describes
More informationUsing DERIVE to Interpret an Algorithmic Method for Finding Hamiltonian Circuits (and Rooted Paths) in Network Graphs
Liverpool John Moores University, July 12 15, 2000 Using DERIVE to Interpret an Algorithmic Method for Finding Hamiltonian Circuits (and Rooted Paths) in Network Graphs Introduction Peter Schofield Trinity
More informationO P E R A T I N G M A N U A L
OPERATING MANUAL WeatherJack OPERATING MANUAL 1-800-645-1061 The baud rate is 2400 ( 8 bits, 1 stop bit, no parity. Flow control = none) To make sure the unit is on line, send an X. the machine will respond
More informationSolving Systems of Linear Equations with the Help. of Free Technology
Solving Systems of Linear Equations with the Help of Free Technology Calin Galeriu, Ph.D. 1. Introduction The use of computer technology when teaching new math concepts, or when solving difficult math
More information