Milne s Spreadsheet Calculator Using VBA with Excel Programming for solving Ordinary Differential Equations

Similar documents
Fourth Order RK-Method

9.6 Predictor-Corrector Methods

Chapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods

Applied Numerical Analysis

Solution of First Order Initial Value Problem by Sixth Order Predictor Corrector Method

The Fourth Version Of Richardson s Extrapolation Spreadsheet Calculator Using Vba Programming For Numerical Differentiations

Differential Equations

Numerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.

Unit I (Testing of Hypothesis)

NUMERICAL ANALYSIS SYLLABUS MATHEMATICS PAPER IV (A)

A Modified Predictor-Corrector Formula For Solving Ordinary Differential Equation Of First Order And First Degree

ECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations

Numerical Solution of Differential Equations

Spreadsheet Implementation for Momentum Representation of Gaussian Wave Packet and Uncertainty Principle

Excel for Scientists and Engineers Numerical Method s. E. Joseph Billo

MATHEMATICAL METHODS INTERPOLATION

PARTIAL DIFFERENTIAL EQUATIONS

STATISTICS AND NUMERICAL METHODS

Physics 584 Computational Methods

16.7 Multistep, Multivalue, and Predictor-Corrector Methods

Mathematics for chemical engineers. Numerical solution of ordinary differential equations

multistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):

2.1 NUMERICAL SOLUTION OF SIMULTANEOUS FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS. differential equations with the initial values y(x 0. ; l.

Chabot College Fall Course Outline for Mathematics 25 COMPUTATIONAL METHODS FOR ENGINEERS AND SCIENTISTS

Numerical Analysis. Introduction to. Rostam K. Saeed Karwan H.F. Jwamer Faraidun K. Hamasalh

Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester

Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations

Math 128A Spring 2003 Week 11 Solutions Burden & Faires 5.6: 1b, 3b, 7, 9, 12 Burden & Faires 5.7: 1b, 3b, 5 Burden & Faires 5.

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK

ODE Runge-Kutta methods

Romberg Integration and Gaussian Quadrature

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

Applied Math for Engineers

SME 3023 Applied Numerical Methods

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

M.SC. PHYSICS - II YEAR

SKMM 3023 Applied Numerical Methods

Review. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn

Jim Lambers MAT 772 Fall Semester Lecture 21 Notes

HIGHER ORDER METHODS. There are two principal means to derive higher order methods. b j f(x n j,y n j )

SOLUTION OF EQUATION AND EIGENVALUE PROBLEMS PART A( 2 MARKS)

Solving Ordinary Differential Equations

Modified Milne Simpson Method for Solving Differential Equations

Computational Methods

NUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group

Simulation of 1D Gaussian Wave Packet and 2D Waves in Square Membrane Using Excel Spreadsheet with Visual Basic for Teaching Activity

Enabling Advanced Problem Solving in Spreadsheets with Access to Physical Property Data

Numerical Methods for Engineers

INTRODUCTION TO COMPUTER METHODS FOR O.D.E.

The integration of response surface method in microsoft excel with visual basic application

Chapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods

Virtual University of Pakistan

Mathematical Methods for Numerical Analysis and Optimization

16.7 Multistep, Multivalue, and Predictor-Corrector Methods

A Review: Utilizing VBA to Solve Burgers Equation by The Method of Lines

Integration of Ordinary Differential Equations

Numerical Methods. Scientists. Engineers

NORTH MAHARASHTRA UNIVERSITY JALGAON.

Preliminary Examination in Numerical Analysis

Numerical Differential Equations: IVP

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

BHARATHIAR UNIVERSITY, COIMBATORE. B.Sc. Mathematics CA (Revised papers with effect from onwards)

Using Spreadsheets to Teach Engineering Problem Solving: Differential and Integral Equations

VEER NARMAD SOUTH GUJARAT UNIVERSITY, SURAT. SYLLABUS FOR B.Sc. (MATHEMATICS) Semester: III, IV Effective from July 2015

Previous Year Questions & Detailed Solutions

Multistep Methods for IVPs. t 0 < t < T

USHA RAMA COLLEGE OF ENGINEERING & TECHNOLOGY

Do not turn over until you are told to do so by the Invigilator.

Fundamental Numerical Methods for Electrical Engineering

Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 1 : Introduction

Name of the Student: Unit I (Solution of Equations and Eigenvalue Problems)

The family of Runge Kutta methods with two intermediate evaluations is defined by

CS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations

Applied Linear Algebra in Geoscience Using MATLAB

+ h4. + h5. 6! f (5) i. (C.3) Since., it yields

On A General Formula of Fourth Order Runge-Kutta Method

Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point

A First Course on Kinetics and Reaction Engineering Supplemental Unit S5. Solving Initial Value Differential Equations

Linear Multistep Methods I: Adams and BDF Methods

Reactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package

Department of Mathematics Faculty of Natural Science, Jamia Millia Islamia, New Delhi-25

Index. function 292. Excep for Scientists and Engineers: Numerical Methods by E. Joseph Billo Copyright John Wiley & Sons, Inc.

Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations

Mathematics for Chemists

A CLASS OF CONTINUOUS HYBRID LINEAR MULTISTEP METHODS FOR STIFF IVPs IN ODEs

Initial value problems for ordinary differential equations

International Journal of Scientific & Engineering Research, Volume 7, Issue 12, December ISSN

Exam in TMA4215 December 7th 2012

Exact and Approximate Numbers:

UNIVERSITY OF PUNE, PUNE BOARD OF STUDIES IN MATHEMATICS S.Y. B. Sc. (MATHEMATICS) SYLLABUS. S.Y.B.Sc. MT:211 Linear Algebra MT:221

Journal of Business & Economics Research September, 2010 Volume 8, Number 9

Exam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20

Bindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.

446 CHAP. 8 NUMERICAL OPTIMIZATION. Newton's Search for a Minimum of f(x,y) Newton s Method

COSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods

Ph 22.1 Return of the ODEs: higher-order methods

Solving Delay Differential Equations (DDEs) using Nakashima s 2 Stages 4 th Order Pseudo-Runge-Kutta Method

NUMERICAL METHODS FOR ENGINEERING APPLICATION

16.1 Runge-Kutta Method

Transcription:

Milne s Spreadsheet Calculator Using VBA with Excel Programming for solving Ordinary Differential Equations Miss. S. P. Hingmire 1, Dr. M.R. Gosavi 2, Dr. N. A.Patil 3 1 Research Scholar, Asst.Professor, SVPM s College of Commerce, Science & Computer Education, Malegaon (Bk), Tal Baramati, Dist-Pune, (Maharashtra) (India) 2 Professor &Head Of Mathematics Dept.,Maharashtra Mahavidyalaya, Nilanga,Dist :Latur, (Maharashtra), India 3 Professor & Head ofapplied Sciences & Humanities Department, Shri Sant Gajanan Maharaj College of Engineering, Shegaon, (Maharashtra), India ABSTRACT In this paper we have shown a calculator for solving differential equations that is developed in Microsoft Excel. We have selected Milne s method. We have built a VBA program for Milne s method for solving ODE numerically. In this program, computer accepts input form and sufficient information to solve the ODE via Milne s method. On the other hand, they are needed to input the interval for independent variable, initial value for dependent variable, step size h in this program.the computational procedures were translated into Visual Basic for Application. Keywords: Milne s Method, Ordinary Differential Equations.,VBA Programme. I. INTRODUCTION Numerical solution of ordinary differential equation with integer order has been developed by long time and has been a standard topic in numerical and computational Mathematics. These solutions can be divided into the exact solution method and the numerical method. The exact solution of the nonlinear ODEs is difficult hence it is necessary to develop a tool in order to solve the ODEs easily. Most of the Numerical methods used for solving ordinary differential equations are based on approximation. VBA stands for Visual Basic for Applications an event driven programming language from Microsoft that is now predominantly used with Microsoft office applications such as MS-Excel, MS-Word and MS-Access. It helps techies to build customized applications and solutions to enhance the capabilities of those applications. The advantage of this facility is that we NEED NOT have visual basic installed on our PC but installing office will implicitly help us to achieve the purpose. Visual Basic for Application (VBA) is the scripting language used to create all macros in this application. The macros were written either from scratch or by editing pre-recorded macros. 204 P a g e

Tay, et, al (2009b) solved systems of nonlinear equations using spread sheet. Kek and Tay (2009c) developed Newton s divided interpolation solver in spread sheet. Lagrange interpolating solver in spread sheet was developed by Tay et, al(2010) whereas bivariate approximation spread sheet calculator was developed by Tay et, al(2014). All the interpolation solvers and spread sheet calculators were developed without applying VBA programming. The spread sheet calculator for the numerical differentiations is developed by Tay et, al (2013a,2013b,2014). In this paper we are using Milne-Simpson method with VBA programming for solving ordinary differential equations. Here we are comparing the result with output values by VBA programming with numerical solutions of a given function. Milne-Simpson Method Details of method can be found in many textbooks for ex.[5,10]. The essences of the method given below.we have discussed many methods for obtaining the solutions of Ordinary Differential Equationswith initial condition. All of them use information only from the computed point (x i, y i ) to compute the next point (x i+1, y i+1 ).These all are called single step methods. They do not make use of the information available at the earlier steps, y i-1, y i-2 etc., even when they are available. Milne-Simpsons method is multistep method since this method use information points from more than one previous points to compute the next point. It is also called Predictor-Corrector method. Milne s formula used as a predictor and Simpson s formula as a corrector. These are based fundamental theorem of calculus. y(x i+1 ) =y(x j ) + xi+1 xi f(x, y) dx Ex:Given equation y (x)= 2y/x with y(1)=2. Estimate y(2) using Milne- Simpsons Predictor Corrector Method. Mines formula at i=3 is y p 4 =y 0 +(4h/3)(2f 1 -f 2 +2f 3 ) Simpson s Fomula at i=3 is y c 4 =y 2 +h/3(f 2 +4f 3 +f p 4 ) Using RK4, where fi=f(xi,yi) y 1 =y(1.25)=3.13, y 2 = y(1.5)= 4.50,y 3 = y(1.75)= 6.13 f 1 =5.01, f 2 = 6, f 3 = 7.01 eq.(1) Y p 4=8.01 f p 4 =8.01 eq.(2) y c 4 =4.5+0.25/3(6+4*7.0+8.01)=8 f 4 =2*8/2=8 y c 4 =4.5+0.25/3(6.00+4*7.07+8)= 8 y c 4 =8 In this section, first frame shows excel spreadsheet for Milne-Simpson method in which a button for Milne- Simpson Method and relevant information are inserted. 205 P a g e

In the second frame, we have shown a program code for Milne-Simpson Method. In the third frame, we get output by providing the relevant information such as initial values and step length. 206 P a g e

The table shows comparison between actual calculated value and output of the program. Actual calculated Value Output of the Program 8 8.002667 IV. CONCLUSION In this paper, we have discussed Milne-Simpson Method for solving Ordinary Differential Equations was developed by VBA programming with Excel. Since users need only relevant information to compute the full solution of the ODEs efficiently, it is very user friendly and time consuming. In this paper we have compared the actual value with output in the third frame of the program designed, which are very close to each other s. REFERENCES [1.] Tay, K. G. and Kek, S. L. (2008), Solver for System of Linear Equations, Proceeding of the National Symposium on Application of Science Mathematics 2008 (SKASM2008),pp. 605-615 [2.] Tay, K. G. and Kek, S. L. (2009), Design of Spreadsheet Solver for Polynomials Interpolation, National Seminar on Science and Technology 2009(PKPST 2009), pp. 69-73 [3.] Tay, K. G. and Kek, S. L. (2009), Fourth Order Runge Kutta Method Using Spreadsheet Excel, Proceeding of the 4 th International Conference on Research and Education in Mathematics(ICREM4), pp. 666-672 207 P a g e

[4.] Tay, K. G. and Kek, S. L., Cheong, T.H., Rosmila. A.K.(2014), The Forth Order Runge Kutta Spreadsheet Calculator Using VBA Programming For Ordinary Differential Equations, 4 th World Congress on Technical and Vocational Education and Training(WoCTVET), pp.231-239 [5.] E Balagurusamy, Numerical Methods, TataMcGRAW HILL Publications. [6.] Dr. J.S.Chitode, Numerical Methods and C Programming, Technical Publications. [7.] G. Haribaskaran (2 nd edition 2008) Numerical Methods. [8.] C.Xavier (1 st edition 1999)) C Language and Numerical Methods New Age International. [9.] S. S. Sastry, Introductory Methods of Numerical Analysis, Fourth Edition Prentice Hall Publication. 208 P a g e