MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL

Size: px
Start display at page:

Download "MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL"

Transcription

1 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1: Specl Issue No 1 MATLAB PROGRAM FOR THE NUMERICAL SOLUTION OF DUHAMEL CONVOLUTION INTEGRAL M BOTIŞ 1 Astrct: I te ler lyss of structures troug modl superposto fter decouplg te moto equtos, tere must e solved te moto equtos for ec degree of dymc freedom, troug umercl tegrto, y usg te Dumel covoluto tegrl To fd structurl respose to ertquke cto t s lso eeded to determe te mxmum spectrl vlues dsplcemet, speed d ccelerto for y ccelerogrm of te ste; tese vlues re determed y umercl tegrto wt Dumel tegrl I ts pper t s preseted Mtl progrm for tegrtg umercl of Dumel covoluto tegrl wt trpezods metod tt c e used ot for modl lyss d spectrl lyss of structures Key words: Dumel tegrl covoluto, dymc lyss y modl superposto, spectrl lyss ertquke egeerg, dymc respose of te structures for rdom forces 1 Itroducto Dymc lyss of cvl structures s ecessry we ctos re tme-vryg If dymc lyss s ler, modl superposto prcple c e used To use modl superposto prcple we must frst ceve modl decomposto d te mke modl superposto After decomposto modl vrto modes re oted, wc glgt modl deformed spes ssocted wt ec mode of vrto d ertl mss correspodg to ec mode of vrto Troug modl decomposto two mportt results re oted, oe s qulttve d s represeted y te form of ow modes of vrto d oter s qutttvely d dctes level of prtcpto of te ertl msses te movemet of dymc system Te lyss of te vrto mode s very mportt ecuse ec vrto mode s cotruto for te moto of te dymc system After spectrl decomposto, modl superposto c e mde, wc volves composg sgfct vrto modes moto dymc system I cse of system wt degrees of freedom (fg1 dymc equto of moto s: [ M ]{ x} [ C]{ x} [ K]{ x} { F} (1 were: -[M ] ertl mtrx of dymc system; -[C] dmpg mtrx of dymc system; -[K] stffess mtrx of dymc system; 1 Dept of Cvl Egeerg, Trslv Uversty of Brşov

2 46 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 ẋ -{ } ccelertos vector for te degrees of freedom dymc; -{ẋ} velocty vector for te degrees of freedom dymc; -{x} dsplcemet vector for te degrees of freedom dymc; -{F} force vector for te degrees of freedom dymc Istteous dymc respose o te degree of freedom k s oted y ler superposto of stteous dymc resposes correspodg to te modes of vrto fg: k x ( x ( t ( 1 k, To dscoect equtos of moto te followg ler trsformto s used: k 1 k, ( k, ( 1 x ( x (3 Reltosp etwee turl coordtes dymc system x k (, k 1 d prcpl coordtes or modl (, 1 c e expressed y te modl mtrx s follows: Fg 1 Te system wt degrees of freedom dymc x x x x 1 1 1,1,1 1,1,1 1, 1,,, 1, 1 1, 1, 1, 1 1, 1,,, 1 1 Dsplcemet vector velocty d ccelerto o te dymc degrees of freedom k 1 ecome: Fg Vrtos mode d egevlues for system wt degrees of freedom dymc { } [ ]{ ( }; { } [ ]{ ( }; (4 { } [ ]{ ( } were: -[] modl mtrx -{ ( } dsplcemets vector turl coordtes -{ ( } velocty vector turl coordtes

3 M BOTIŞ: Mtl progrm for te umercl soluto of Dumel covoluto tegrl 47 -{ ( } ccelertos vector turl coordtes Replcg te system equtos of moto (1 rel coordtes wt modl coordtes, te equtos of moto s relesed Due to ortogolty of vectors ssocted vrto modes system of equtos (1 s reduced to decoupled equtos of moto for ec degree of freedom dymc: As te system of dfferetl equtos of moto s decoupled, we c wrte te correspodg equto of moto for ec degree, s follows: F { ( } { ( } { ( } (5 M 1 were: T - F { F( } correspodg equvlet force mode ; T - M M equvlet mss mode ; - correspodg dmpg mode -[ C] [ M ] [ K] Ryleg dmpg mtrx To solve decoupled equtos of moto, Dumel covoluto tegrl c e used Dsplcemet respose oted usg Te covoluto tegrl Dumel tegrl o ec dymc degrees of freedom s: t 1 ( t ( F ( e s ( t d M 1 (6 were, - 1 dmped pulsto If you kow te swer dsplcemets o ec dymc degree of freedom te we c clculte veloctes ccelertos d ertl forces for dymc system: { } [ ]{ ( }; { } [ ]{ ( }; { } [ ]{ ( } Te lyss of te dymcs of system wt degrees of freedom dymcs c e see s solvg te system of equtos s reduced to clculto covoluto tegrl Dumel For ts reso ts pper presets Mtl progrm for umercl tegrto of covoluto tegrl Dumel If te t pply H mpulse of te system, ts s removed from equlrum (7 posto wt velocty x H / m, (fg3 d performs free dmped vrtos descred y te equto : H e m H ( t e m ( t s ( t s ( t for t (8 Fg 3 Free respose dsplcemet due to mpulse fte H If te dymc system rdom dsturg force cts,(fg4 equto

4 48 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 of moto of te dymc system s: F( m Acto of dsturg force ppled dymc system t rest c e represeted s cotuous elemetry mpulses, o te tme tervl from t t After pplyg elemetry mpulse dsplcemet respose of te dymc system dymc o degree of freedom s: dh ( t d e s ( t d m (9 F( ( t e s ( t d m Dymc dsplcemet respose of te system s ceved y superposto of effects s tegrl mout of moto geerted y elemetry pulse successo o te legt of tme tt force cts [1]: F( ( t e s ( t d m f ( ( t F e s ( t d m Fg 4 Dsplcemet respose oted y superposto of te effects of ftesml mpulses dh If dsturg force s ppled to dymc system wle t s moto due to tl dsplcemet d velocty codtos, te dsplcemet respose s expressed: e [( x x cost s st] f ( ( t F e s ( t d m were, - x tl dsplcemet for t= - x tl velocty for t= It c e see tt f te dymc system wt sgle degree of freedom dymc system to determe dymc respose Dumel covoluto tegrl must e clculted Mtemtcl spects of umercl soluto of te Dumel tegrl To tegrte umerclly covoluto tegrl Dumel, furter s preseted trpezods metod wc s te smplest metod used prctce Trpezods metod s sed o pproxmtg te re delmted y te grp of fucto f :[, ] R wt te re of trpezod If te tervl [, ] s dvded to two sutervls [, x1] d [ x 1, ] te tegrl of fucto f(x c e pproxmted y te relto: dx x1 dx dx [ f ( f ( x1] [ f ( f ( x1] [ f ( f ( x1 f ( ] [ f ( f ( x1 f ( ] 4 To reduce te error of tegrto, te tervl [,] my e dvded sutervls wt dvso ( /, x1

5 M BOTIŞ: Mtl progrm for te umercl soluto of Dumel covoluto tegrl 49 tegrl of fucto f(x c e pproxmted y te relto: dx x1 dx x11 dx [ f ( f ( x1] [ f ( f ( x 1] [ f ( f ( x1 f ( x 1 f ( ] [ f ( f ( x1 f ( x 1 f ( ] Error clcultg troduced y te metod of trpezods depeds o te tegrto step d s te followg expresso o tervl [ x x+]:: 3,, E ( f ( 1 Full error s: 1 3 1,, ET ( E ( f ( 1 I 1 I 1 3 Implemetto of Mtl progrm for clculto Dumel tegrl For umercl clculto of te covoluto tegrl Dumel strts from respose dsplcemet for system wt oe degree of freedom: F( ( t e s( t d m F( t e cos d e st m F( t e s d e cot m t t A( e st B( e cost were, F( A( e cos d m F( B( e s d m To solve te Dumel tegrl trpezods metod wll e used to clculte tegrls A( d B( wt Mtl progrm [] % Solvg covoluto tegrl Dumel %y umercl tegrto metod of %trpezods %Aswer dplcemet, velocty d %ccelerto for dymc sstem oe degree %of freedom % force rmp cler;clf; % Tme step used =; % Mss system wt 1GLD m=175; % Rgdty of te system wt 1GLD k=7; % Frcto of crtcl dmpg u=5/1; % Pulsto system 1GLD wtout %dmpg omeg=sqrt(k/m; % Pulsto system 1GLD wt %dmpg omegd=omegsqrt(1-u^; % Tme tervl tt s computed %respose system wt 1 DOF for %dsplcemet, velocty d ccelerto tmx=14; % Force rmp fucto vrto cure=[ ]; % Numer of pots of fucto curves pt=msze(cure; t=; f=terp1([cure(1,:p]', [cure(,:p]',; y=f1exp(uomeg cos(omegd(1momegd^-1; y=f1exp(uomeg s(omegd(1momegd^-1; % Numer of steps for umercl tegrto ps=tmx/; suplot(,1,1 for =1:ps t=t+; % Clculte te force t tme t fort=terp1([cure(1,:p]', [cure(,:p]', f1=fort5;

6 5 Bullet of te Trslv Uversty of Brşov Vol 5 (54-1 Seres 1 f(=f1; tmp(=t; % Vrto represetto force %fucto tme p=plot(t,f1,'--rs','lewdt',, 'MrkerEdgeColor','', 'MrkerFceColor','', 'MrkerSze',; old o; grd o; % Clculto of vlues y d y %from Dumel covoluto tegrl y(=f11exp(uomeg cos(omegd(1momegd^-1; y(=f11exp(uomeg s(omegd(1momegd^-1; ed % Fl vlues y=[y y(1,:' ]; y=[y y(1,:']; % Metod trpezods t= suplot(,1, for =:ps+1 t=t+; A(1=; B(1=; A(=A(-1+(y(+y(-1/; B(=B(-1+(y(+y(-1/; % Respose moto for rmp sgl %cse 1GLD dymc system d(=(a(exp(-uomeg s(omegd- B(exp(-uomeg cos(omegd(k/5; % Grpcl represetto of te rto of %dymc dsplcemet d sttc %dsplcemet plot(t,d(,'--rs','lewdt',, 'MrkerEdgeColor','r', 'MrkerFceColor','r', 'MrkerSze',; old o; grd o; ed 4 Results d cocluso Te progrm proposed te pper s ee tested te cse of dymcl system wt oe degree of freedom dymc Dymc lyss ws performed for two steps of tegrto = d =1 Te results fter umercl tegrto re preseted fg5 d fg6 Fg 5 Te respose dsplcemet for dymc system wt 1 degree of freedom, tegrto step, Fg 6 Te respose dsplcemet for dymc system wt 1 degree of freedom, tegrto step,1 Te result oted y umercl tegrto ws compred wt te exct soluto, for te tegrto step =, ot error 13,7 % d f tegrto step s = error ws,3% Altoug te progrm preseted s smple, t c e successfully used f te tegrto step s lower t were T/1, were T s ow perod for dymc system Refereces 1 Ifrm, M: Dmc structurlor s gere sesmc Edtur ddctc s pedgogc, 1973 Adelw Kr, Rold B Gueter,: Itroducto to Numercl Metods: A MATLAB Approc CRC Press Tylor &Frcs Group, 1

MTH 146 Class 7 Notes

MTH 146 Class 7 Notes 7.7- Approxmte Itegrto Motvto: MTH 46 Clss 7 Notes I secto 7.5 we lered tht some defte tegrls, lke x e dx, cot e wrtte terms of elemetry fuctos. So, good questo to sk would e: How c oe clculte somethg

More information

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is Mth Sprg 08 L Approxmtg Dete Itegrls I Itroducto We hve studed severl methods tht llow us to d the exct vlues o dete tegrls However, there re some cses whch t s ot possle to evlute dete tegrl exctly I

More information

The z-transform. LTI System description. Prof. Siripong Potisuk

The z-transform. LTI System description. Prof. Siripong Potisuk The -Trsform Prof. Srpog Potsuk LTI System descrpto Prevous bss fucto: ut smple or DT mpulse The put sequece s represeted s ler combto of shfted DT mpulses. The respose s gve by covoluto sum of the put

More information

Numerical differentiation (Finite difference)

Numerical differentiation (Finite difference) MA umercl Alyss Week 4 To tk d wrte portle progrms we eed strct model of ter computtol evromet. Ftful models do est ut tey revel tt evromet to e too dverse forcg portle progrmmers to lot eve te smplest

More information

Chapter Simpson s 1/3 Rule of Integration. ( x)

Chapter Simpson s 1/3 Rule of Integration. ( x) Cpter 7. Smpso s / Rule o Itegrto Ater redg ts pter, you sould e le to. derve te ormul or Smpso s / rule o tegrto,. use Smpso s / rule t to solve tegrls,. develop te ormul or multple-segmet Smpso s / rule

More information

Mathematically, integration is just finding the area under a curve from one point to another. It is b

Mathematically, integration is just finding the area under a curve from one point to another. It is b Numerl Metods or Eg [ENGR 9] [Lyes KADEM 7] CHAPTER VI Numerl Itegrto Tops - Rem sums - Trpezodl rule - Smpso s rule - Rrdso s etrpolto - Guss qudrture rule Mtemtlly, tegrto s just dg te re uder urve rom

More information

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Chapter 2 Intro to Math Techniques for Quantum Mechanics Wter 3 Chem 356: Itroductory Qutum Mechcs Chpter Itro to Mth Techques for Qutum Mechcs... Itro to dfferetl equtos... Boudry Codtos... 5 Prtl dfferetl equtos d seprto of vrbles... 5 Itroducto to Sttstcs...

More information

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ]

this is the indefinite integral Since integration is the reverse of differentiation we can check the previous by [ ] Atervtves The Itegrl Atervtves Ojectve: Use efte tegrl otto for tervtves. Use sc tegrto rules to f tervtves. Aother mportt questo clculus s gve ervtve f the fucto tht t cme from. Ths s the process kow

More information

Chapter 12-b Integral Calculus - Extra

Chapter 12-b Integral Calculus - Extra C - Itegrl Clulus Cpter - Itegrl Clulus - Etr Is Newto Toms Smpso BONUS Itroduto to Numerl Itegrto C - Itegrl Clulus Numerl Itegrto Ide s to do tegrl smll prts, lke te wy we preseted tegrto: summto. Numerl

More information

UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM

UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 7/8 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 7 MASA : 3 JAM KOD KURSUS : KKKQ33/KKKF33 TAJUK : PENGIRAAN BERANGKA ARAHAN :.

More information

CS321. Numerical Analysis

CS321. Numerical Analysis CS Numercl Alyss Lecture 4 Numercl Itegrto Proessor Ju Zg Deprtmet o Computer Scece Uversty o Ketucky Legto, KY 456 6 Octoer 6, 5 Dete Itegrl A dete tegrl s tervl or tegrto. For ed tegrto tervl, te result

More information

On Several Inequalities Deduced Using a Power Series Approach

On Several Inequalities Deduced Using a Power Series Approach It J Cotemp Mth Sceces, Vol 8, 203, o 8, 855-864 HIKARI Ltd, wwwm-hrcom http://dxdoorg/02988/jcms2033896 O Severl Iequltes Deduced Usg Power Seres Approch Lored Curdru Deprtmet of Mthemtcs Poltehc Uversty

More information

CS321. Introduction to Numerical Methods

CS321. Introduction to Numerical Methods CS Itroducto to Numercl Metods Lecture Revew Proessor Ju Zg Deprtmet o Computer Scece Uversty o Ketucky Legto, KY 6 6 Mrc 7, Number Coverso A geerl umber sould be coverted teger prt d rctol prt seprtely

More information

Integration by Parts for D K

Integration by Parts for D K Itertol OPEN ACCESS Jourl Of Moder Egeerg Reserc IJMER Itegrto y Prts for D K Itegrl T K Gr, S Ry 2 Deprtmet of Mtemtcs, Rgutpur College, Rgutpur-72333, Purul, West Begl, Id 2 Deprtmet of Mtemtcs, Ss Bv,

More information

i+1 by A and imposes Ax

i+1 by A and imposes Ax MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 09.9 NUMERICAL FLUID MECHANICS FALL 009 Mody, October 9, 009 QUIZ : SOLUTIONS Notes: ) Multple solutos

More information

under the curve in the first quadrant.

under the curve in the first quadrant. NOTES 5: INTEGRALS Nme: Dte: Perod: LESSON 5. AREAS AND DISTANCES Are uder the curve Are uder f( ), ove the -s, o the dom., Prctce Prolems:. f ( ). Fd the re uder the fucto, ove the - s, etwee,.. f ( )

More information

Module 2: Introduction to Numerical Analysis

Module 2: Introduction to Numerical Analysis CY00 Itroducto to Computtol Chemtr Autum 00-0 Module : Itroducto to umercl Al Am of the preet module. Itroducto to c umercl l. Developg mple progrm to mplemet the umercl method opc of teret. Iterpolto:

More information

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ] Are d the Defte Itegrl 1 Are uder Curve We wt to fd the re uder f (x) o [, ] y f (x) x The Prtto We eg y prttog the tervl [, ] to smller su-tervls x 0 x 1 x x - x -1 x 1 The Bsc Ide We the crete rectgles

More information

Evaluation of the Singular Stress Field. of the Crack Tip of the Infinite Elastic Plate

Evaluation of the Singular Stress Field. of the Crack Tip of the Infinite Elastic Plate ICCM05, 4-7 th July, Auckld, NZ Evluto of the Sgulr Stress Feld of the Crck Tp of the Ifte Elstc lte * J.Y. Lu,, G. L, Z.C.Su 3,. Zhg d Y.Y. Wg Fculty of Ifrstructure Egeerg, Dl Uversty of Techology, Dl

More information

Spectrum Estimation by Several Interpolation Methods

Spectrum Estimation by Several Interpolation Methods IJCSNS Itertol Jourl o Computer Scece d Network Securty VOL.6 No.A Februry 006 05 Spectrum Estmto by Severl Iterpolto Metods Mbu Isr Oym Ntol College o Tecology Oym-S Tocg 33-0806 JAPAN Summry I ts pper

More information

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

INTERPOLATION(2) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek ELM Numerl Alss Dr Murrem Merme INTEROLATION ELM Numerl Alss Some of te otets re dopted from Luree V. Fusett Appled Numerl Alss usg MATLAB. rete Hll I. 999 ELM Numerl Alss Dr Murrem Merme Tod s leture

More information

EFFECT OF CRACK PARAMETERS ON FREE VIBRATIONS OF THE BERNOULLI-EULER BEAM

EFFECT OF CRACK PARAMETERS ON FREE VIBRATIONS OF THE BERNOULLI-EULER BEAM Jourl of Appled themtcs d Computtol echcs 05 () 67-7 www.mcm.pcz.pl p-in 99-9965 DOI: 0.75/jmcm.05..7 e-in 353-0588 EFFECT OF CRACK PARAETER ON FREE VIBRATION OF THE BERNOULLI-EULER BEA Izbel Zmorsk Dwd

More information

CS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department

CS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department CS473-Algorthms I Lecture 3 Solvg Recurreces Cevdet Aykt - Blket Uversty Computer Egeerg Deprtmet Solvg Recurreces The lyss of merge sort Lecture requred us to solve recurrece. Recurreces re lke solvg

More information

Free Vibration Analysis of Rectangular Plates Using Galerkin-Based Finite Element Method

Free Vibration Analysis of Rectangular Plates Using Galerkin-Based Finite Element Method Itertol Jourl of Mechcl Egeerg ISSN : 77-709 Volume Issue http://.meourl.com/ https://stes.google.com/ste/ourlme/ Free Vbrto Alyss of Rectgulr Pltes Usg Glerk-Bsed Fte Elemet Method Nefft M. Werfll, Abobker

More information

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11: Soo Kg Lm 1.0 Nested Fctorl Desg... 1.1 Two-Fctor Nested Desg... 1.1.1 Alss of Vrce... Exmple 1... 5 1.1. Stggered Nested Desg for Equlzg Degree of Freedom... 7 1.1. Three-Fctor Nested Desg... 8 1.1..1

More information

14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y

14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y 4. Le Itegrls I ths secto we defe tegrl tht s smlr to sgle tegrl except tht sted of tegrtg over tervl [ ] we tegrte over curve. Such tegrls re clled le tegrls lthough curve tegrls would e etter termology.

More information

DESIGNING AND CONSTRUCTION OF ROADS, SUBWAYS, AIRFIELDS, BRIDGES AND TRANSPORT TUNNELS

DESIGNING AND CONSTRUCTION OF ROADS, SUBWAYS, AIRFIELDS, BRIDGES AND TRANSPORT TUNNELS Issue (), ISSN 5- DESIGNING AND CONSTRUCTION OF ROADS, SUBWAYS, AIRFIELDS, BRIDGES AND TRANSPORT TUNNELS UDC 59.. A. Kozlov THE DEFLECTED MODE OF MULTI COHERENT PRISMATIC CONSTRUCTIE ELEMENTS OF BRIDGE

More information

CURVE FITTING LEAST SQUARES METHOD

CURVE FITTING LEAST SQUARES METHOD Nuercl Alss for Egeers Ger Jord Uverst CURVE FITTING Although, the for of fucto represetg phscl sste s kow, the fucto tself ot be kow. Therefore, t s frequetl desred to ft curve to set of dt pots the ssued

More information

Differential Method of Thin Layer for Retaining Wall Active Earth Pressure and Its Distribution under Seismic Condition Li-Min XU, Yong SUN

Differential Method of Thin Layer for Retaining Wall Active Earth Pressure and Its Distribution under Seismic Condition Li-Min XU, Yong SUN Itertol Coferece o Mechcs d Cvl Egeerg (ICMCE 014) Dfferetl Method of Th Lyer for Retg Wll Actve Erth Pressure d Its Dstrbuto uder Sesmc Codto L-M XU, Yog SUN Key Lbortory of Krst Evromet d Geologcl Hzrd

More information

Numerical Solution of Second order Integro- Differential Equations(Ides) with Different Four Polynomials Bases Functions

Numerical Solution of Second order Integro- Differential Equations(Ides) with Different Four Polynomials Bases Functions Numercl Soluto of Secod order Itegro- Dfferetl Equtos(Ides) wth Dfferet Four Polyomls Bses Fuctos wo O A R M Deprtmet of Mthemtcs Uersty of Ilor Deprtmet of Mthemtcs d Sttstcs he Poly Id Astrct: - I ths

More information

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini DATA FITTING Itesve Computto 3/4 Als Mss Dt fttg Dt fttg cocers the problem of fttg dscrete dt to obt termedte estmtes. There re two geerl pproches two curve fttg: Iterpolto Dt s ver precse. The strteg

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Mtr Trsformtos usg Egevectors September 8, Mtr Trsformtos Usg Egevectors Lrry Cretto Mechcl Egeerg A Semr Egeerg Alyss September 8, Outle Revew lst lecture Trsformtos wth mtr of egevectors: = - A ermt

More information

Level-2 BLAS. Matrix-Vector operations with O(n 2 ) operations (sequentially) BLAS-Notation: S --- single precision G E general matrix M V --- vector

Level-2 BLAS. Matrix-Vector operations with O(n 2 ) operations (sequentially) BLAS-Notation: S --- single precision G E general matrix M V --- vector evel-2 BS trx-vector opertos wth 2 opertos sequetlly BS-Notto: S --- sgle precso G E geerl mtrx V --- vector defes SGEV, mtrx-vector product: r y r α x β r y ther evel-2 BS: Solvg trgulr system x wth trgulr

More information

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Chapter 2 Intro to Math Techniques for Quantum Mechanics Fll 4 Chem 356: Itroductory Qutum Mechcs Chpter Itro to Mth Techques for Qutum Mechcs... Itro to dfferetl equtos... Boudry Codtos... 5 Prtl dfferetl equtos d seprto of vrbles... 5 Itroducto to Sttstcs...

More information

Roberto s Notes on Integral Calculus Chapter 4: Definite integrals and the FTC Section 2. Riemann sums

Roberto s Notes on Integral Calculus Chapter 4: Definite integrals and the FTC Section 2. Riemann sums Roerto s Notes o Itegrl Clculus Chpter 4: Defte tegrls d the FTC Secto 2 Rem sums Wht you eed to kow lredy: The defto of re for rectgle. Rememer tht our curret prolem s how to compute the re of ple rego

More information

19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in

19 22 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y Write the given sum or difference as a single integral in SECTION. THE DEFINITE INTEGRAL. THE DEFINITE INTEGRAL A Clck here for swers. S Clck here for solutos. Use the Mdpot Rule wth the gve vlue of to pproxmte the tegrl. Roud the swer to four decml plces. 9

More information

20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2

20 23 Evaluate the integral by interpreting it in terms of areas. (1 x ) dx. y 2 2 SECTION 5. THE DEFINITE INTEGRAL 5. THE DEFINITE INTEGRAL A Clck here for swers. S Clck here for solutos. 7 Use the Mdpot Rule wth the gve vlue of to pproxmte the tegrl. Roud the swer to four decml plces.

More information

ScienceDirect. About Verification of Discrete-Continual Finite Element Method of Structural Analysis. Part 2: Three-Dimensional Problems

ScienceDirect. About Verification of Discrete-Continual Finite Element Method of Structural Analysis. Part 2: Three-Dimensional Problems Avlle ole t wwwscecedrectcom SceceDrect Proced Egeerg 9 (04 4 9 XXIII R-S-P semr heoretcl Foudto of Cvl Egeerg (RSP (FoCE 04 Aout Verfcto of Dscrete-Cotul Fte Elemet Method of Structurl Alyss Prt : hree-dmesol

More information

Analytical Approach for the Solution of Thermodynamic Identities with Relativistic General Equation of State in a Mixture of Gases

Analytical Approach for the Solution of Thermodynamic Identities with Relativistic General Equation of State in a Mixture of Gases Itertol Jourl of Advced Reserch Physcl Scece (IJARPS) Volume, Issue 5, September 204, PP 6-0 ISSN 2349-7874 (Prt) & ISSN 2349-7882 (Ole) www.rcourls.org Alytcl Approch for the Soluto of Thermodymc Idettes

More information

Hopf Bifurcation of Maglev System with Coupled Elastic Guideway

Hopf Bifurcation of Maglev System with Coupled Elastic Guideway Hopf Bfurcto of Mglev System wth Coupled Elstc udewy No Loghu She, Ho Wg, Dogsheg Zou, Zhzhou Zhg* d Wese Chg Egeerg Reserch Ceter of Mglev Techology, Chgsh Isttute of Techology, Chgsh 7, Hu, PRCh sheloghu@yhoocomc,

More information

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl Alyss for Egeers Germ Jord Uversty ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl soluto of lrge systems of ler lgerc equtos usg drect methods such s Mtr Iverse, Guss

More information

Sequences and summations

Sequences and summations Lecture 0 Sequeces d summtos Istructor: Kgl Km CSE) E-ml: kkm0@kokuk.c.kr Tel. : 0-0-9 Room : New Mleum Bldg. 0 Lb : New Egeerg Bldg. 0 All sldes re bsed o CS Dscrete Mthemtcs for Computer Scece course

More information

Itō Calculus (An Abridged Overview)

Itō Calculus (An Abridged Overview) Itō Clculus (A Abrdged Overvew) Arturo Ferdez Uversty of Clfor, Berkeley Sttstcs 157: Topcs I Stochstc Processes Semr Aprl 14, 211 1 Itroducto I my prevous set of otes, I troduced the cocept of Stochstc

More information

Preliminary Examinations: Upper V Mathematics Paper 1

Preliminary Examinations: Upper V Mathematics Paper 1 relmr Emtos: Upper V Mthemtcs per Jul 03 Emer: G Evs Tme: 3 hrs Modertor: D Grgortos Mrks: 50 INSTRUCTIONS ND INFORMTION Ths questo pper sts of 0 pges, cludg swer Sheet pge 8 d Iformto Sheet pges 9 d 0

More information

Strategies for the AP Calculus Exam

Strategies for the AP Calculus Exam Strteges for the AP Clculus Em Strteges for the AP Clculus Em Strtegy : Kow Your Stuff Ths my seem ovous ut t ees to e metoe. No mout of cochg wll help you o the em f you o t kow the mterl. Here s lst

More information

St John s College. UPPER V Mathematics: Paper 1 Learning Outcome 1 and 2. Examiner: GE Marks: 150 Moderator: BT / SLS INSTRUCTIONS AND INFORMATION

St John s College. UPPER V Mathematics: Paper 1 Learning Outcome 1 and 2. Examiner: GE Marks: 150 Moderator: BT / SLS INSTRUCTIONS AND INFORMATION St Joh s College UPPER V Mthemtcs: Pper Lerg Outcome d ugust 00 Tme: 3 hours Emer: GE Mrks: 50 Modertor: BT / SLS INSTRUCTIONS ND INFORMTION Red the followg structos crefull. Ths questo pper cossts of

More information

Non-Polynomial Spline Method for the Solution of Problems in Calculus of Variations

Non-Polynomial Spline Method for the Solution of Problems in Calculus of Variations World Acdemy of Scece Egeerg d Tecology No-Polyoml Sple etod for te Solto of Prolems Clcls of Vrtos. Zre. Hosyr. Sedgt Astrct I ts pper mercl solto sed o opolyoml cc sple fctos s sed for fdg te solto of

More information

General Method for Calculating Chemical Equilibrium Composition

General Method for Calculating Chemical Equilibrium Composition AE 6766/Setzma Sprg 004 Geeral Metod for Calculatg Cemcal Equlbrum Composto For gve tal codtos (e.g., for gve reactats, coose te speces to be cluded te products. As a example, for combusto of ydroge wt

More information

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS Pweł Pędzch Jerzy Blcerz Wrsw Uversty of Techology Fculty of Geodesy d Crtogrphy Astrct Usully to pproto of

More information

Least Squares Method For Solving Integral Equations With Multiple Time Lags

Least Squares Method For Solving Integral Equations With Multiple Time Lags Eg. & Tech. Jourl, Vol.8, No., Lest Squres Method For Solvg Itegrl Equtos Wth Multple Tme Lgs Dr.Suh N. Sheh*, Hyt Adel Al* & Hl Mohmmed Ysee* Receved o:6//9 Accepted o:// Astrct The m purpose of ths work

More information

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS PubH 745: REGRESSION ANALSIS REGRESSION IN MATRIX TERMS A mtr s dspl of umbers or umercl quttes ld out rectgulr rr of rows d colums. The rr, or two-w tble of umbers, could be rectgulr or squre could be

More information

MATRIX AND VECTOR NORMS

MATRIX AND VECTOR NORMS Numercl lyss for Egeers Germ Jord Uversty MTRIX ND VECTOR NORMS vector orm s mesure of the mgtude of vector. Smlrly, mtr orm s mesure of the mgtude of mtr. For sgle comoet etty such s ordry umers, the

More information

A Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares

A Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares Itertol Jourl of Scetfc d Reserch Publctos, Volume, Issue, My 0 ISSN 0- A Techque for Costructg Odd-order Mgc Squres Usg Bsc Lt Squres Tomb I. Deprtmet of Mthemtcs, Mpur Uversty, Imphl, Mpur (INDIA) tombrom@gml.com

More information

FORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS

FORMULAE FOR FINITE DIFFERENCE APPROXIMATIONS, QUADRATURE AND LINEAR MULTISTEP METHODS Jourl o Mtemtcl Sceces: Advces d Alctos Volume Number - Pges -9 FRMULAE FR FINITE DIFFERENCE APPRXIMATINS QUADRATURE AND LINEAR MULTISTEP METHDS RAMESH KUMAR MUTHUMALAI Dertmet o Mtemtcs D G Vsv College

More information

INTRODUCTION ( ) 1. Errors

INTRODUCTION ( ) 1. Errors INTRODUCTION Numercl lyss volves the study, developmet d lyss of lgorthms for obtg umercl solutos to vrous mthemtcl problems. Frequetly umercl lyss s clled the mthemtcs of scetfc computg. Numercl lyss

More information

A METHOD FOR THE RAPID NUMERICAL CALCULATION OF PARTIAL SUMS OF GENERALIZED HARMONICAL SERIES WITH PRESCRIBED ACCURACY

A METHOD FOR THE RAPID NUMERICAL CALCULATION OF PARTIAL SUMS OF GENERALIZED HARMONICAL SERIES WITH PRESCRIBED ACCURACY UPB c Bull, eres D, Vol 8, No, 00 A METHOD FOR THE RAPD NUMERAL ALULATON OF PARTAL UM OF GENERALZED HARMONAL ERE WTH PRERBED AURAY BERBENTE e roue o etodă ouă etru clculul rd l suelor rţle le serlor roce

More information

Numerical Differentiation and Integration

Numerical Differentiation and Integration Numerl Deretto d Itegrto Overvew Numerl Deretto Newto-Cotes Itegrto Formuls Trpezodl rule Smpso s Rules Guss Qudrture Cheyshev s ormul Numerl Deretto Forwrd te dvded deree Bkwrd te dvded deree Ceter te

More information

Chapter 4: Distributions

Chapter 4: Distributions Chpter 4: Dstrbutos Prerequste: Chpter 4. The Algebr of Expecttos d Vrces I ths secto we wll mke use of the followg symbols: s rdom vrble b s rdom vrble c s costt vector md s costt mtrx, d F m s costt

More information

4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula

4. Runge-Kutta Formula For Differential Equations. A. Euler Formula B. Runge-Kutta Formula C. An Example for Fourth-Order Runge-Kutta Formula NCTU Deprme o Elecrcl d Compuer Egeerg Seor Course By Pro. Yo-Pg Ce. Ruge-Ku Formul For Derel Equos A. Euler Formul B. Ruge-Ku Formul C. A Emple or Four-Order Ruge-Ku Formul

More information

POWERS OF COMPLEX PERSYMMETRIC ANTI-TRIDIAGONAL MATRICES WITH CONSTANT ANTI-DIAGONALS

POWERS OF COMPLEX PERSYMMETRIC ANTI-TRIDIAGONAL MATRICES WITH CONSTANT ANTI-DIAGONALS IRRS 9 y 04 wwwrppresscom/volumes/vol9issue/irrs_9 05pdf OWERS OF COLE ERSERIC I-RIIGOL RICES WIH COS I-IGOLS Wg usu * Q e Wg Hbo & ue College of Scece versty of Shgh for Scece d echology Shgh Ch 00093

More information

EXACT STOCHASTIC ANALYSIS OFCOUPLED BENDING-TORSION BEAMS WITH IN-SPAN SUPPORTS AND MASSES

EXACT STOCHASTIC ANALYSIS OFCOUPLED BENDING-TORSION BEAMS WITH IN-SPAN SUPPORTS AND MASSES Meccc de Mterl e delle Strutture Vol. VI 06), o., pp. 5-3 ISS: 035-679X Dprtmeto d Igeger Cvle, Ambetle, Aerospzle, De Mterl DICAM EXACT STOCHASTIC AALYSIS OFCOUPLED BEDIG-TORSIO BEAMS WITH I-SPA SUPPORTS

More information

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation PGE 30: Formulato ad Soluto Geosystems Egeerg Dr. Balhoff Iterpolato Numercal Methods wth MATLAB, Recktewald, Chapter 0 ad Numercal Methods for Egeers, Chapra ad Caale, 5 th Ed., Part Fve, Chapter 8 ad

More information

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III

Outline. Finite Difference Grids. Numerical Analysis. Finite Difference Grids II. Finite Difference Grids III Itrodcto to Nmercl Alyss Mrc, 9 Nmercl Metods or PDEs Lrry Cretto Meccl Egeerg 5B Semr Egeerg Alyss Mrc, 9 Otle Revew mdterm soltos Revew bsc mterl o mercl clcls Expressos or dervtves, error d error order

More information

4. Runge-Kutta Formula For Differential Equations

4. Runge-Kutta Formula For Differential Equations NCTU Deprme o Elecrcl d Compuer Egeerg 5 Sprg Course by Pro. Yo-Pg Ce. Ruge-Ku Formul For Derel Equos To solve e derel equos umerclly e mos useul ormul s clled Ruge-Ku ormul

More information

Linear Algebra Concepts

Linear Algebra Concepts Ler Algebr Cocepts Nuo Vscocelos (Ke Kreutz-Delgdo) UCSD Vector spces Defto: vector spce s set H where ddto d sclr multplcto re defed d stsf: ) +( + ) = (+ )+ 5) H 2) + = + H 6) = 3) H, + = 7) ( ) = (

More information

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n 0. Sere I th ecto, we wll troduce ere tht wll be dcug for the ret of th chpter. Wht ere? If we dd ll term of equece, we get whch clled fte ere ( or jut ere) d deoted, for hort, by the ymbol or Doe t mke

More information

Computer Programming

Computer Programming Computer Progrmmg I progrmmg, t s ot eough to be vetve d geous. Oe lso eeds to be dscpled d cotrolled order ot be become etgled oe's ow completes. Hrl D. Mlls, Forwrd to Progrmmg Proverbs b Her F. Ledgrd

More information

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions HK Km Slghtly moded //9 /8/6 Frstly wrtte t Mrch 5 Itertve Methods er Systems: Guss-Sedel Noler Systems Cse Study: Chemcl Rectos Itertve or ppromte methods or systems o equtos cosst o guessg vlue d the

More information

Rendering Equation. Linear equation Spatial homogeneous Both ray tracing and radiosity can be considered special case of this general eq.

Rendering Equation. Linear equation Spatial homogeneous Both ray tracing and radiosity can be considered special case of this general eq. Rederg quto Ler equto Sptl homogeeous oth ry trcg d rdosty c be cosdered specl cse of ths geerl eq. Relty ctul photogrph Rdosty Mus Rdosty Rederg quls the dfferece or error mge http://www.grphcs.corell.edu/ole/box/compre.html

More information

Advanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University

Advanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University Advced Algorthmc Prolem Solvg Le Arthmetc Fredrk Hetz Dept of Computer d Iformto Scece Lköpg Uversty Overvew Arthmetc Iteger multplcto Krtsu s lgorthm Multplcto of polyomls Fst Fourer Trsform Systems of

More information

An Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix

An Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix Itertol Jourl of Scetfc d Reserch Pulctos, Volue 3, Issue 6, Jue 3 ISSN 5-353 A Altertve Method to Fd the Soluto of Zero Oe Iteger Ler Frctol Progrg Prole wth the Help of -Mtr VSeeregsy *, DrKJeyr ** *

More information

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ. Численные методы. Учебно-методическое пособие

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ. Численные методы. Учебно-методическое пособие МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Нижегородский государственный университет им. Н.И. Лобачевского Численные методы К.А.Баркалов Учебно-методическое пособие Рекомендовано методической

More information

Some Unbiased Classes of Estimators of Finite Population Mean

Some Unbiased Classes of Estimators of Finite Population Mean Itertol Jourl O Mtemtcs Ad ttstcs Iveto (IJMI) E-IN: 3 4767 P-IN: 3-4759 Www.Ijms.Org Volume Issue 09 etember. 04 PP-3-37 ome Ubsed lsses o Estmtors o Fte Poulto Me Prvee Kumr Msr d s Bus. Dertmet o ttstcs,

More information

On Solution of Min-Max Composition Fuzzy Relational Equation

On Solution of Min-Max Composition Fuzzy Relational Equation U-Sl Scece Jourl Vol.4()7 O Soluto of M-Mx Coposto Fuzzy eltol Equto N.M. N* Dte of cceptce /5/7 Abstrct I ths pper, M-Mx coposto fuzzy relto equto re studed. hs study s geerlzto of the works of Ohsto

More information

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices. 4.3 - Modal Aalyss Physcal coordates are ot always the easest to work Egevectors provde a coveet trasformato to modal coordates Modal coordates are lear combato of physcal coordates Say we have physcal

More information

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Joural of Appled Matematcs ad Computatoal Mecacs 04, 3(4), 7-34 FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM Ata Cekot, Stasław Kukla Isttute of Matematcs, Czestocowa Uversty of Tecology Częstocowa,

More information

Systems of second order ordinary differential equations

Systems of second order ordinary differential equations Ffth order dgolly mplct Ruge-Kutt Nystrom geerl method solvg secod Order IVPs Fudzh Isml Astrct A dgolly mplct Ruge-Kutt-Nystróm Geerl (SDIRKNG) method of ffth order wth explct frst stge for the tegrto

More information

Quantization in Dynamic Smarandache Multi-Space

Quantization in Dynamic Smarandache Multi-Space Quatzato Dyamc Smaradache Mult-Space Fu Yuhua Cha Offshore Ol Research Ceter, Beg, 7, Cha (E-mal: fuyh@cooc.com.c ) Abstract: Dscussg the applcatos of Dyamc Smaradache Mult-Space (DSMS) Theory. Supposg

More information

Numerical Analysis Formulae Booklet

Numerical Analysis Formulae Booklet Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx

More information

Chapter Unary Matrix Operations

Chapter Unary Matrix Operations Chpter 04.04 Ury trx Opertos After redg ths chpter, you should be ble to:. kow wht ury opertos mes,. fd the trspose of squre mtrx d t s reltoshp to symmetrc mtrces,. fd the trce of mtrx, d 4. fd the ermt

More information

Method for Hydrodynamic Coupling of Concentric Cylindrical Shells and Beams

Method for Hydrodynamic Coupling of Concentric Cylindrical Shells and Beams ethod for Hydrodymc Couplg of Cocetrc Cyldrcl Shells d Bems tthew D. Syder Westghouse Electrc Compy, LLC Abstrct Dymc lyses of pressurzed wter rector vessels, core support structures, d terls requre ccurte

More information

Introduction to local (nonparametric) density estimation. methods

Introduction to local (nonparametric) density estimation. methods Itroducto to local (oparametrc) desty estmato methods A slecture by Yu Lu for ECE 66 Sprg 014 1. Itroducto Ths slecture troduces two local desty estmato methods whch are Parze desty estmato ad k-earest

More information

Vibration Analysis Of a Self-Excited Elastic Beam

Vibration Analysis Of a Self-Excited Elastic Beam Vbrto lyss Of Self-Excted Elstc Bem M.. Brró-Mez Deprtmeto de Mterles Uversdd utóom Metropolt zcpotzlco v. S Pblo 80, ol. Reyos-Tmulps.P. 000, Méxco, D.F., MÉXIO bmm@correo.zc.um.mx BSTRT The vbrto behvor

More information

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers Regresso By Jugl Klt Bsed o Chpter 7 of Chpr d Cle, Numercl Methods for Egeers Regresso Descrbes techques to ft curves (curve fttg) to dscrete dt to obt termedte estmtes. There re two geerl pproches two

More information

6.6 Moments and Centers of Mass

6.6 Moments and Centers of Mass th 8 www.tetodre.co 6.6 oets d Ceters of ss Our ojectve here s to fd the pot P o whch th plte of gve shpe lces horzotll. Ths pot s clled the ceter of ss ( or ceter of grvt ) of the plte.. We frst cosder

More information

INTEGRATION TECHNIQUES FOR TWO DIMENSIONAL DOMAINS

INTEGRATION TECHNIQUES FOR TWO DIMENSIONAL DOMAINS IJRET: Itertol Jourl of Reserch Egeerg d Techolog eissn: 9-6 pissn: -78 INTEGRATION TECHNIQUES FOR TWO DIMENSIONAL DOMAINS Logh Peruml Lecturer, Fcult of Egeerg d Techolog, Multmed Uverst, Mlcc, Mls Astrct

More information

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II CEE49b Chapter - Free Vbrato of Mult-Degree-of-Freedom Systems - II We ca obta a approxmate soluto to the fudametal atural frequecy through a approxmate formula developed usg eergy prcples by Lord Raylegh

More information

Numerical Analysis Topic 4: Least Squares Curve Fitting

Numerical Analysis Topic 4: Least Squares Curve Fitting Numerl Alss Top 4: Lest Squres Curve Fttg Red Chpter 7 of the tetook Alss_Numerk Motvto Gve set of epermetl dt: 3 5. 5.9 6.3 The reltoshp etwee d m ot e ler. Fd futo f tht est ft the dt 3 Alss_Numerk Motvto

More information

SUM PROPERTIES FOR THE K-LUCAS NUMBERS WITH ARITHMETIC INDEXES

SUM PROPERTIES FOR THE K-LUCAS NUMBERS WITH ARITHMETIC INDEXES Avlble ole t http://sc.org J. Mth. Comput. Sc. 4 (04) No. 05-7 ISSN: 97-507 SUM PROPERTIES OR THE K-UCAS NUMBERS WITH ARITHMETIC INDEXES BIJENDRA SINGH POOJA BHADOURIA AND OMPRAKASH SIKHWA * School of

More information

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem Joural of Amerca Scece ;6( Cubc Nopolyomal Sple Approach to the Soluto of a Secod Order Two-Pot Boudary Value Problem W.K. Zahra, F.A. Abd El-Salam, A.A. El-Sabbagh ad Z.A. ZAk * Departmet of Egeerg athematcs

More information

The definite Riemann integral

The definite Riemann integral Roberto s Notes o Itegrl Clculus Chpter 4: Defte tegrls d the FTC Secto 4 The defte Rem tegrl Wht you eed to kow lredy: How to ppromte the re uder curve by usg Rem sums. Wht you c ler here: How to use

More information

Asymptotic Formulas Composite Numbers II

Asymptotic Formulas Composite Numbers II Iteratoal Matematcal Forum, Vol. 8, 203, o. 34, 65-662 HIKARI Ltd, www.m-kar.com ttp://d.do.org/0.2988/mf.203.3854 Asymptotc Formulas Composte Numbers II Rafael Jakmczuk Dvsó Matemátca, Uversdad Nacoal

More information

Lecture 3-4 Solutions of System of Linear Equations

Lecture 3-4 Solutions of System of Linear Equations Lecture - Solutos of System of Ler Equtos Numerc Ler Alger Revew of vectorsd mtrces System of Ler Equtos Guss Elmto (drect solver) LU Decomposto Guss-Sedel method (tertve solver) VECTORS,,, colum vector

More information

DISCRETE TIME MODELS OF FORWARD CONTRACTS INSURANCE

DISCRETE TIME MODELS OF FORWARD CONTRACTS INSURANCE G Tstsshvl DSCRETE TME MODELS OF FORWARD CONTRACTS NSURANCE (Vol) 008 September DSCRETE TME MODELS OF FORWARD CONTRACTS NSURANCE GSh Tstsshvl e-ml: gurm@mdvoru 69004 Vldvosto Rdo str 7 sttute for Appled

More information

Mathematics HL and further mathematics HL formula booklet

Mathematics HL and further mathematics HL formula booklet Dplom Progrmme Mthemtcs HL d further mthemtcs HL formul boolet For use durg the course d the emtos Frst emtos 04 Publshed Jue 0 Itertol Bcclurete Orgzto 0 5048 Mthemtcs HL d further mthemtcs formul boolet

More information

Proceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition IMECE2011 November 11-17, 2011, Denver, Colorado, USA

Proceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition IMECE2011 November 11-17, 2011, Denver, Colorado, USA Proceedgs of the ASME Itertol Mechcl Egeerg ogress & Exosto IMEE ovember -7 Dever olordo USA IMEE-6869 ALORITHM OF LAUHED VEHILES MASS ALULATIO AT THE EARLY STAE OF DESII IMEE-6869 A. AYUTDIOVA Dertmet

More information

Mathematics HL and further mathematics HL formula booklet

Mathematics HL and further mathematics HL formula booklet Dplom Progrmme Mthemtcs HL d further mthemtcs HL formul boolet For use durg the course d the emtos Frst emtos 04 Edted 05 (verso ) Itertol Bcclurete Orgzto 0 5048 Cotets Pror lerg Core 3 Topc : Algebr

More information

An Analysis of Queueing Model with a Multiple System Governed by a Quasi-Birth Death Process and PH Service and PH Repair

An Analysis of Queueing Model with a Multiple System Governed by a Quasi-Birth Death Process and PH Service and PH Repair Itertol Jourl of Scece d Reserch (IJSR) ISSN (Ole): 39-764 Idex Copercus Vlue (3): 64 Impct Fctor (4): 56 A Alyss of Queueg Model wth Multple System Govered y Qus-Brth Deth Process d PH Servce d PH Repr

More information

Chapter 1. Introduction. Fundamental Concepts. Introduction. Historical background. Historical background. Fundamental Concepts

Chapter 1. Introduction. Fundamental Concepts. Introduction. Historical background. Historical background. Fundamental Concepts Chpter udmetl Cocepts Lecture Notes Dr Mohd Afed Uverst Mlys Perls N67 te lemet Alyss Itroducto A or sometmes referred s M, hs ecome powerful tool for umercl soluto of wde rge of egeerg prolems A s computtol

More information

Solutions Manual for Polymer Science and Technology Third Edition

Solutions Manual for Polymer Science and Technology Third Edition Solutos ul for Polymer Scece d Techology Thrd Edto Joel R. Fred Uer Sddle Rver, NJ Bosto Idols S Frcsco New York Toroto otrel Lodo uch Prs drd Cetow Sydey Tokyo Sgore exco Cty Ths text s ssocted wth Fred/Polymer

More information