Implementation of B-Spline Based Finite Element Analysis. for the Boundary-Value Problem. Andrew J. Marker. Master of Science
|
|
- Beverly Tucker
- 6 years ago
- Views:
Transcription
1 Implementaton of B-Splne Based Fnte Element Analyss for the Boundary-Value Problem Andrew J. Marker A project submtted to the faculty of Brgham Young Unversty n partal fulfllment of the requrements for the degree of Master of Scence Mchael A. Scott, Char Rchard J. Ballng Norman L. Jones Department of Cvl and Envronmental Engneerng Brgham Young Unversty December 2013 Copyrght 2013 Andrew Marker All Rghts Reserved
2
3 ABSTRACT Implementaton of B-Splne Based Fnte Element Analyss for the Boundary-Value Problem Andrew J. Marker Department of Cvl and Envronmental Engneerng Master of Scence Ths project wll eamne the use of B-Splnes as the bass shaped functons n fnte element (FE) analyss. The analyss wll be appled to a boundary-value problem. A program wll be wrtten to verfy that B-Splne fnte element (BSFE) analyss accurately represents the boundary-value problem. Keywords: b-splnes, fnte element
4
5 TABLE OF CONTENTS LIST OF TABLES... v LIST OF FIGURES... 1 Introducton Background... 2 Fnte Element Method Basc Concept... 3 B-Splnes Basc Concept Bass Functons Methods... 7 The Boundary Value Problem The Strong Form The Weak Form The Galerkn Form Matr Form Local Element pont of vew... 9 MATLAB Program Gaussan Quadrature Bass Values and Dervatves v
6 3.2.3 Mappng from Local to Global Error Results and Conclulsons Dscusson of Results Sources of Error Conclusons REFERENCES Append A. Plots of Solutons v
7 LIST OF TABLES Table 4-1: Error for f() = Table 4-2: Error for f() = Table 4-3: Error for f() = v
8 v
9 LIST OF FIGURES Fgure 2-1: Bézer and B-Splne Curves... 4 Fgure 2-2: Lnear B-Splne Bass over Fve Elements... 5 Fgure 2-3: Quadratc B-Splne Bass over Four Elements... 5 Fgure 2-4: Cubc B-Splne Bass over S Elements... 6 Fgure 4-1: f() = 1, n = Fgure 4-2: Enlarged pcture of f() = 1, n= Fgure 4-3: f() = 1, n = Fgure 4-4: Error for f() =
10
11 1 INTRODUCTION The fnte element method (FEM) has been wdely used for engneerng purposes snce ts ntroducton n the 1940 s and 1950 s. The age of computers accelerated the use and understandng of FEM. FEM requres the use of bass equatons to appromate the functon. Often pece-wse contnuous polynomals, LaGrange or Hermte polynomals, have been used as the bass functons n FEM. B-splnes offer an alternatve bass that s contnuous over the whole model. In nterestng applcaton of B-Splne based FEM s so-geometrc analyss (IGA). B-Splnes are used etensvely n computer modelng whch allows for B-Splnes to be a lnk between modelng and analyss. Ths project wll eamne the use of B-Splnes as the bass shaped functons n FEM analyss. The analyss wll be appled to the boundary-value problem. A program was developed to run FEM on the boundary-value problem wth frst, second, and thrd-order bass functons. The error between the actual soluton and the FEM results are calculated and compared wth FEM usng Lagrange polynomals of the same order. The rate of convergence was also calculated and used to determne the effcency of B-Splne based FEM. Ths paper provdes a bref overvew of the fnte element method and B-Splnes. It dscusses the method of the project. The results wll be presented and dscussed. 1
12 2 BACKGROUND Ths secton provdes a bref ntroducton to the fnte element method and b-splnes. It s meant to provde enough nformaton for the reader to understand the project, and s certanly not ehaustve. For each topc a bref hstory s gven as well as an ntroducton to the basc dea. Fnte Element Method The fnte method s used today for many engneerng applcatons. Much research has and contnues to be devoted to ths topc because of ts mportance n engneerng analyss, although the fnte element method s relatvely new. Its orgns are n the 1940 s and 1950 s when mathematcans and engneers were lookng to solve hgh order problems wth new numercal methods. Frst publcatons on the subject nclude R.L. Courant s Varatonal Methods for the Soluton of Problems of Equlbrum and Vbraton n whch he used pece-wse lnear appromatons on trangles. The method s more clearly outlned by Argyrs and Turner wth others n the 1950 s. The technque was frst appled to the dscplne of structures but found other uses ncludng heat conducton and fluds n the 1960 s. The development of computers has ncreased the amount of research and applcaton of the FEM. (Pepper, Henrch, 2006, ). Along wth the development of hgh-speed dgtal computers, the applcaton of the fnte element method also progressed at a very mpressve rate (Rao, 2001). FEM contnues to be a hghly regarded and used numercal method n engneerng and other dscplnes. Research contnues n order fnd new technques and applcatons. 2
13 2.1.1 Basc Concept The fnte element method (FEM) s a numercal method. Lke all numercal methods, the dea s to fnd an appromaton to a comple problem by replacng t wth a smpler one. The dea behnd FEM s to break up a larger geometry nto smaller sub-geometres, or elements. Analyss of the ndvdual elements can be completed based on the governng equatons of the problem. Appromatons are made by lnear combnatons of the bass functons. The results are than assembled nto a global soluton whch appromates the real soluton. The appromaton mproves as the mesh s refned, that s as the number of elements ncreases. Dfferent technques or greater computatonal effort can also be appled to FEM to create closer appromatons. (Huebner 2001, Reddy, 1984)). B-Splnes Smlar to the fnte element method, the hstory of B-splnes (bass splnes) s relatvely new.. B-splnes were frst ntroduced for use n statstcal data smoothng by Schoenberg n 1946 who coned the phrase. B-splnes were developed for the automotve ndustry n the 1960 s as an easy way to model the comple curves as found on automobles. Two French mathematcans, Bézer wth Renault and De Casteljau s wth Ctroën, smultaneously developed B-splnes for geometrc desgn. There are two mportant mathematcal developments n the hstory of B-splnes. Frst s the development of the recursve relaton for defnng and evaluatng B-splnes bass functons. Ths s often credted to de Boor et al. The second s the creaton of the polar form by Ramshaw whch allows for quck understandng and manpulatng of B-splnes. Today B-splnes 3
14 are used n the felds of statstcs, computer aded desgn (CAD), computer aded geometrc desgn, anmaton, graphc desgn, and fnte element analyss. Even the fonts and curves on a word processor are B-splnes Basc Concept B-splnes are a sequence of Bezer curves constructed n a way to create C n-1 contnuty, where C s the order. B-splnes are defned by the control ponts, P, assocated wth weghts, and the knot vector whch controls the parameterzaton of the splne. The polar form of B-splne mentoned above makes t easy to fnd ntermedate values between control ponts Fgure 1 shows a sngle Bézer curve and a B-splne consstng of three Bézer curves (Sederberg, 2012). P(1, 2, 3) P(2, 3, 3) P(1, 1, 2) a) Bézer Curve P(1, 1, 1) Knot Vector ( ) b) B-Splne Fgure 2-1: Bézer and B-Splne Curves P(3, 3, 3) Bass Functons An mportant aspect of FEM s choosng a bass functon, or a lnear combnaton by whch functons can be appromated. B-splne bass functons do ths and also guarantee C n-1 contnuty whch becomes convenent for hgher order problems. Usng the knot vector and recursve relaton B-splne bass functons are nepensve to evaluate as are the dervatves. AS wth any lnear combnaton, the sum of the values at any pont s one, and the sum of the dervatves at any pont 4
15 s zero. Another nnovatve reason for usng B-splne bass functons wth FEM s for the use n so-geometrc analyss, whch tres to close the gap between the graphcal model and the model used for analyss. There s currently a lot of research n ths area. Fgure 2-2, Y and Z show the lnear, quadratc and cubc bass functons. Notce there are C1 splne over each element, where C agan s the order of the splnes. Fgure 2-2: Lnear B-Splne Bass over Fve Elements Fgure 2-3: Quadratc B-Splne Bass over Four Elements 5
16 6 Fgure 2-4: Cubc B-Splne Bass over S Elements
17 3 METHODS The purpose of ths project was to mplement a B-splnes bass fnte element program to solve the boundary value problem. The man porton of the project conssted of defnng the problem, and solvng the problem wth a program n MATLAB. Also an Ecel spreadsheet was created and programmed to generate values for the bass functon for a gven order and number of elements. The Boundary Value Problem The boundary value problem (BVP) has been wdely studed as an ntroducton to fnte element method. It s a smple one dmensonal frst-order problem. It s ntroduced here brefly as n Hughes (87) and can be studed n greater detal by referrng to that book. The goal of the followng dervaton s to get to a matr form of the equaton whch s famlar and easy to solve The Strong Form The strong or classcal form of the BVP n terms of u on Ω = [0,1] s (SS) Gven ff Ω R and constants gg and h, fnd uu Ω R, ssssssh tthaaaa uu, ff = 0 oooo Ω uu(1) = gg uu, (0) = h 3.1 where a comma represent dfferentaton and f s a functon defned on Ω. The second and thrd equatons represent the boundary condtons. Both g and h are 0 for ths project. 7
18 3.1.2 The Weak Form To get the weak form of the BVP equaton we frst defne the tral soluton space, δ, and weghtng functon space, υ, as the followng δδ = {uu uu HH 1, uu(1) = gg} 3.2 Υ = {ww ww HH 1, ww(1) = 0} 3.3 where H 1 are functons who satsfy the followng equaton. 1 (uu, ) 2 dddd < Then the weak form can be wrtten as Gven ff, gg, and h, as before. Fnd uu δδ such that for all ww Υ (WW) = 1 1 ww, uu, dddd = wwww dddd ww(0)h Let 1 aa(ww, uu) = ww, uu, dddd (ww, ff) = wwww dddd Then the weak form can be wrtten as aa(ww, uu) = (ww, ff) ww(0)h The Galerkn Form Dscretzng usng Galerkn s appromaton we create subsets of the tral and watng spaces such that δδ h δδ aaaaaa Υ h Υ. Then u h s constructed such that u h = υ h g h generatng the Galerkn form: Gven ff, gg, aand h, as before, fnd uu h = υυ h gg h wwheeeeee υυ h Υ h, (GG) such that for all ww h Υ h aa(ww h, υυ h ) = (ww h, ff) ww h (0)h aa(ww h, gg h ) 3.9 8
19 3.1.4 Matr Form The Galerkn form results n a system of lnear equatons whch can be nterpolated wth the bass functons and wrtten n a matr form famlar to engneers, Kd = F, where KK = [KK AAAA ], where KK AAAA = aa(nn AA, NN BB ) 3.10 FF = {FF AA }, where FF AA = (NN AA, ff) NN AA (0)h aa(nn AA, NN nn1 )gg 3.11 Note that NA represents the bass functons, or n ths case the B-splne bass functons. Recall that g and h are both zero for ths project smplfyng the equaton for FA to FA = (NA, f). Once F and K are assembled d can be solved for whch s our soluton u Local Element pont of vew In FEM each element s analyzed and then element components are assembled nto the global problem. For each element the k e matr and f e vector s calculated wth the followng equatons. However the local lmts of ntegraton are now -1 to 1 whch requres some adjustment when assemblng K. ee kk aaaa = aa(nnnn, NNNN) 3.12 ff aa ee = (ff, NNNN) 3.13 MATLAB Program Key portons of the program most pertnent to the results wll be hghlghted here. Append A contans an outlne vew of the entre program as well as the code. The code was modfed from prevous code whch used Lagrange polynomals as the bass. The code the Lagrange bass was kept for comparson. 9
20 3.2.1 Gaussan Quadrature Gaussan quadrature was used to evaluate the ntegrals n the analyss and n fndng the error. Gaussan quadrature works by takng a weghted sum of values at f at key ntegratons ponts. 1 nn ff() ff(ξξ ll ) WW ll 1 ll= It s an effectve way to evaluate complcated ntegrals. Hughes reports that accuracy of order 2nnt s ganed by nnt ponts of ntegraton (Hughes, 1987). A nnt of 3 was used for ths project requrng that ξξ 1 = 3/5, ξξ 2 = 0, ξξ 3 = 3/5 W1 = W3 = 5/9, W2 = 8/ Bass Values and Dervatves A large porton of the new code was to evaluate the new bass. Because Gaussan quadrature was used to appromate the ntegrals, so the bass had to be evaluated at three ponts n each element. Wth Lagrange polynomals the functons are the same over each element, so only three ponts had to be evaluated. Wth the B-splnes, three ponts for each element needed to be calculated. Usng a unform knot nterval k:[ n-1 n n n] for n elements, and the recursve relaton, the equatons for the bass functons, and ther dervatves can easly be evaluated at any pont. For eample the quadratc formulas are 10
21 = elsewhere B 0, ), [, ) )( ( ) ( ), [, ) )( ( ) )( ( 1 ), [, ) )( ( ) ( ) ( Where represent the knots. Functons were created that returned values and dervatves for each of the bass equatons Mappng from Local to Global The element space s dfferent than the global space due to the way the local element s defned on [-1 1] and the global element s defned on [a b] where b a = 1. To handle ths apply a change of varables to the local equatons whch result n a dervatve term n the ntegral. In one dmenson ths s related to the frst dervatve. Pror to ntegratng, the dervatve term at each quadrature pont s determned as cc aa NN aa (ξξ) CC1 aa= Where ca are the control ponts of the B-Splnes. For convenence the control pont were chosen at the grevlle abscssae so the geometrc space and parametrc space are equvalent. 11
22 3.2.4 Error Error was calculated usng a least squared approach. Agan quadrature was appled for the followng ntegral. The rate of convergence, change n error per change n the number of elements s a good ndcator of how well a FEM s workng. 1 1/2 EE = (uu aaaatttttttt uu aaaaaaaaaaaa dddd)
23 4 RESULTS AND CONCLULSIONS Dscusson of Results In the results are ncluded numbers for both the Lagrange solutons (uh solutons) and B- Splne solutons (uhb solutons) where P s the bass order. Fgure 4-1 shows the solutons for f() =1 and n = 5. Fgure 4-2 shows s an enlarged vew of the same soluton. Fgure 4-3 s for f() = 1 wth n = Smlar plots other forcng functon and n values are n Append B. There are several thngs to pont out from these plots. Frst, the lnear functons for both Lagrange and B- splnes are the same. On Fgure 4-2, t s notceable that they are ndeed dentcal. Also, t s clear at low numbers of elements that the hgher order B-splnes are not a good appromaton. Also note that the cubc soluton s erroneously less eact than the quadratcs, and smlarly wth the quadratcs and lnears. Fnally t can be seen that the B-Splne solutons do eventually converge. Fgure 4-1: f() = 1, n = 5 13
24 Fgure 4-2: Enlarged pcture of f() = 1, n=5 Fgure 4-3: f() = 1, n =
25 The net few tables show the error and the rate of convergence for dfferent forcng functons. They confrm that the frst-order bass have the same result and that the hgh order converge at about a 1/10 of the speed of the Lagrange bass. That s, the Lagrange bass decreases error by 100 tmes wth an ncrease of n by 10. The hgher order B-Splne only ncrease at rato of 10 to 10 and the cubc B-splne soluton dverge between 500 and 5000 elements. Fgure 4-4 s the graphcal representaton of table 4-3. Table 4-1: Error for f() = 1 F() = 1 N L1 L2 L3 B1 B2 B E E E E E E E E E E E E E E E E E E E E-04 Convergence 2.00E E E E E E-01 Table 4-2: Error for f() = F() = N L1 L2 L3 B1 B2 B E E E E E E E E E E E E E E E E E E E E-04 Convergence 2.00E E E E E E-01 15
26 Table 4-3: Error for f() = 2 F() = 2 N L1 L2 L3 B1 B2 B E E E E E E E E E E E E E E E E E E E E-04 Convergence 2.00E E E E E E Error vs 1/N Error E-08 1E-10 L1 L2 L3 B1 B2 B3 1E-12 1E-14 1/n Fgure 4-4: Error for f() = 2 16
27 Sources of Error Whle there are small errors that come from the appromatons of the FEM, t s clear that the results acqured here has more error than t should. It was brefly mentoned that there s a dfference between the local and global spaces. It s most lkely that ths large error s due to the fact that the code does not properly account for mappng between the two spaces. The dervatve factor used to get the local nto global get beng ncorrect or ths mappng could be otherwse mssed n a dfferent area. Conclusons The program wrtten to mplement a B-splne bass FEM was successful n the fact that t created solutons that converged to the answer wth mesh refnement (as the number of elements ncreased). It was not successful as there s a large source of error, probably due to an ssue mappng between the two spaces, and because t converges too slowly. The program dd not successfully mplement the B-splne bass. In order to make the program a success, the source of error would need to be dentfed. Once the program worked for the boundary value problem, t could than smlarly be appled to hgher order problems. 17
28 REFERENCES Huebner, Kenneth H The Fnte Element Method for Engneers. 4th ed. New York: Wley. Hughes, Thomas J.R The Fnte Element Method: Lnear Statc and Dynamc Fnte Element Analyss. Mneola, NY: Dover Publcatons. Pepper, Darrell W., and Juan C. Henrch The Fnte Element Method: Basc Concepts and Applcatons. 2nd ed. New York: Taylor & Francs. Rao, Sngresu S., The Fnte Element Method n Engneerng. 5th ed. Burlngton, MA: Butterworth-Henemann. Reddy, J. N An Introducton to the Fnte Element Method. New York: McGraw-Hll. Sederberg, Thomas W Computer Aded Geometrc Desgn. 18
29 APPENDIX A. PLOTS OF SOLUTIONS Result plots for n = 5, 50,500, 5000 for each functon f() =1,, ^2. 19
30 20
31 21
32 22
33 23
34 24
35 25
Lecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationCubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two
World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationSolution for singularly perturbed problems via cubic spline in tension
ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationDEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND
ANZIAM J. 45(003), 195 05 DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND YOUNG JOON AHN 1 (Receved 3 August, 001; revsed 7 June, 00) Abstract In ths paper
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More information1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions
Lakshm Sreesha Ch. Non Standard Fnte Dfference Method for Sngularly Perturbed Sngular wo Pont Boundary Value Problem usng Non Polynomal Splne LAKSHMI SIREESHA CH Department of Mathematcs Unversty College
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More information36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to
ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationNUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST
Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationUnit 5: Quadratic Equations & Functions
Date Perod Unt 5: Quadratc Equatons & Functons DAY TOPIC 1 Modelng Data wth Quadratc Functons Factorng Quadratc Epressons 3 Solvng Quadratc Equatons 4 Comple Numbers Smplfcaton, Addton/Subtracton & Multplcaton
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationSystem in Weibull Distribution
Internatonal Matheatcal Foru 4 9 no. 9 94-95 Relablty Equvalence Factors of a Seres-Parallel Syste n Webull Dstrbuton M. A. El-Dacese Matheatcs Departent Faculty of Scence Tanta Unversty Tanta Egypt eldacese@yahoo.co
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationModified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls [Verson: 2000/03/22 v.0] Modfed Mass Matrces and Postvty Preservaton for Hyperbolc and
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationChapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D
Chapter Twelve Integraton 12.1 Introducton We now turn our attenton to the dea of an ntegral n dmensons hgher than one. Consder a real-valued functon f : R, where the doman s a nce closed subset of Eucldean
More informationConstruction of Serendipity Shape Functions by Geometrical Probability
J. Basc. Appl. Sc. Res., ()56-56, 0 0, TextRoad Publcaton ISS 00-0 Journal of Basc and Appled Scentfc Research www.textroad.com Constructon of Serendpty Shape Functons by Geometrcal Probablty Kamal Al-Dawoud
More informationMODELING TRAFFIC LIGHTS IN INTERSECTION USING PETRI NETS
The 3 rd Internatonal Conference on Mathematcs and Statstcs (ICoMS-3) Insttut Pertanan Bogor, Indonesa, 5-6 August 28 MODELING TRAFFIC LIGHTS IN INTERSECTION USING PETRI NETS 1 Deky Adzkya and 2 Subono
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC
COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationIsogeometric Analysis with Geometrically Continuous Functions on Multi-Patch Geometries. Mario Kapl, Florian Buchegger, Michel Bercovier, Bert Jüttler
Isogeometrc Analyss wth Geometrcally Contnuous Functons on Mult-Patch Geometres Maro Kapl Floran Buchegger Mchel Bercover Bert Jüttler G+S Report No 35 August 05 Isogeometrc Analyss wth Geometrcally Contnuous
More informationCubic Trigonometric Rational Wat Bezier Curves
Cubc Trgonometrc Ratonal Wat Bezer Curves Urvash Mshra Department of Mathematcs Mata Gujr Mahla Mahavdyalaya Jabalpur Madhya Pradesh Inda Abstract- A new knd of Ratonal cubc Bézer bass functon by the blendng
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationFuzzy Boundaries of Sample Selection Model
Proceedngs of the 9th WSES Internatonal Conference on ppled Mathematcs, Istanbul, Turkey, May 7-9, 006 (pp309-34) Fuzzy Boundares of Sample Selecton Model L. MUHMD SFIIH, NTON BDULBSH KMIL, M. T. BU OSMN
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationInstituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245
Insttuto Tecnológco de Aeronáutca FIITE ELEMETS I Class notes AE-5 Insttuto Tecnológco de Aeronáutca 5. Isoparametrc Elements AE-5 Insttuto Tecnológco de Aeronáutca ISOPARAMETRIC ELEMETS Introducton What
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v10157-01-00-y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More information: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:
764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationA LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) ,
A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS Dr. Derald E. Wentzen, Wesley College, (302) 736-2574, wentzde@wesley.edu ABSTRACT A lnear programmng model s developed and used to compare
More informationCHAPTER 4. Vector Spaces
man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationFinite Difference Method
7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use
More information( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.
3 Interpolaton {( y } Gven:,,,,,, [ ] Fnd: y for some Mn, Ma Polynomal Appromaton Theorem (Weerstrass Appromaton Theorem --- estence ε [ ab] f( P( , then there ests a polynomal
More informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationComputation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models
Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,
More informationCALCULUS CLASSROOM CAPSULES
CALCULUS CLASSROOM CAPSULES SESSION S86 Dr. Sham Alfred Rartan Valley Communty College salfred@rartanval.edu 38th AMATYC Annual Conference Jacksonvlle, Florda November 8-, 202 2 Calculus Classroom Capsules
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More information