Differential equations
|
|
- Roy Horn
- 5 years ago
- Views:
Transcription
1 Differential equations First order differential equation Imlicit F(, y, y ) = or Elicit y = f ( y, ) In general = hysical laws involving a rate of change (with time) of an event = function
2 Solution of DE = Function y= h Derivable y h = and satisfy DE on interval a< < b Conditions for DE to have a solution are general Many simle DE do not have a solution E. ( y ) = no solution E. ( y ) y + y = General solution: y= c c family of straight lines; Singular solution: y = arabola 4 General solution: many solutions differing by constant family of curves Singular solution: additional solution not obtained from general solution E. Radioactive decay Observed decay with time of radioactive substance to amount resent dy ky dt = dy k is a constant (half life) deending on radioactive substance; since dt < k < General solution: y kt = ce y = cke kt = ky t Verification: Particular solution: at t =, y = Y y = ce = c c= Y yt = Ye kt
3 Geometrical interretations y of y- y = f ( y, ) is the sloe of y( ) if y( ) solution assing through (, ) lane then sloe at that oint is f (, y ) f (, ) solution curves Isoclines: (, ) y = f y = k= constant E. y = y, the isoclines are hyerbolas y = direction field (or sloe field) of General solution eist: y = ce ; verification: = = y ce ( ) y E. van der Pol equation of electronics y.( ) = y Direction fields give all solutions with limited accuracy More accurate method = numeric method Allows to solve DE for which we have no solution E Euler-Cauchy method and Runge- Kutta method 3
4 Methods to find general solutions Case : Searable DE g y y = f g y dy = f d dy Solution: g( y) d= f d+ C d Integrals eist rovided g and f continuous dy d = dy g y dy = f d + C d but since E. 9yy + 4 = 9ydy+ 4d = Solution: or 9ydy+ 4d+ C = y + = C 4 9 family of ellises centered on (,) E. dy y = + y = d+ C + y Solution: arctan y = + C y = tan( + C) Note y = tan + C is not a solution E. 3 y = ky with k > dy k+ C Solution: = k d+ C ln y = k+ C y = e = Ce y C C Where C =+ e for y > and C = e y < k 4
5 E. 4 Initial value roblem initial condition (IC) y DE: y = IC: y = Solution: dy = d ln y = ln + C = ln + C C y = y C IC = C = y = E. 5 Initial value roblem DE: y y Solution: y = = IC: dy = ln = + = = y c e =+ C for y > c Where e = C for y < C = for y = IC: y y e + C d y C y e Ce = = Bell-shaed curves 5
6 Case : Reduction to searable form y y = g Solution use transformation y u = y = u y = u + u u + u= g u u = g u u du d = g u u E. y y yy = y y = = y y u + u + u = u u = u+ = u u u udu d = ln+ u = ln + C + u C y C + u = + = + y = C C C + y = 4 Family of circles, tangent to origin, with centers on -ais 6
7 Construction of mathematical models Radio carbon dating 4 Based on isotoe C 6 Rendered radioactive by cosmic rays bombardment 4 In the air the ratio of the isotoes C C is constant Princile of method: When living organism dies the absortion (breathing and eating) of The ratio of C 4 ends 6 4 C C measured comared to value in air to date of death Half life of 4 6C is 573 yrs 4 ± yrs (error of 7%) Mathematical model: kt y = ky yt = ye where y is the initial amount of C 4 6 From definition of half life time at which the level is ½ the original one: k 573 k ye = y e = k = ln. yr 573 Measurement: if remains contain 5% of original level.t ln.5 ye =.5 y t = yrs 46± 8yrs. Comarisons with other dating methods suggest t 3 yrs value given by model is too low, why? 7
8 Miing roblem A thank contains gal of water in which 4 lb of salt is dissolved. Each minute 5gal of water, with lb/gal of salt is added, and the same amount is retrieved yt lb Problem: amount of salt as a function of time: Simlification: the miing rocess is instantaneous and erfect y t = dy = dt Mathematical model: rate of change inflow outflow Inflow: lb/min Outflow: since one gal of water contains y( t) lb gal yt 5yt lb 5 gal coming out will contain = =.5yt 4 min y =.5yt Rate of change: IC: y = 4 lb Solution: y =.5( y 4) dy =.5dt ( y 4) = + = +.5t ln y 4.5t C y Ce 4 IC: y = 4 C = 36.5t yt 4 36e = Solution increases with time from 4 to 4. Equilibrium reached when amount going in = amount going out y = ; 8
9 Heating roblem (Newton s law of cooling) The time rate of change of temerature of a body is described by the following equation: dt T = Tbody TMedium dt Model assuming ambient temerature constant at T medium = 3 dt dt = k( T 3 body ) = kdt dt T 3 ( body ) ln T 3 = kt + C T = 3+ Ce body body kt F IC: T() = T = 66 F 3 F+ C = 66 F C = 34 F T t = 3 F+ 34 F e kt Observation: after t = h T = 63 F body 63 3 T = + e = e = k ( h ) = ln.9765=.4687 h k k h 3 F 34 F 63 F.9765 Rate of change of temerature decreases with time as the body temerature aroach temerature of ambient Equilibrium = temerature with the ambient t = h For T e.4687 h = 3 F+ 34 F = 53.4 F 9
10 Velocity escae from Earth MODEL = Newton s law of gravitationma = GMm r At the surface of Earth with radius R : At any radius a( r) gr = r a = g = GM R Alying definition of acceleration and chain rule dr v R Searation of variable DE vdv= gr = g + C r r v IC: on surface of the Earth r = R and v= v C = gr gr General solution: v = + v gr r dv dv dr dv R a r = = = v= g dt dr dt dr r When v = the rojectile stos and reverse trajectory (fall back to Earth) v However, if we take v = gr, then gr = and ositive reaching only at the infinite v gr = > and v remains r Escae velocity: v = gr
11 Eact DE General form: M( yd ) Solution: u(, y), + N, y dy = = C, such that u u du = d+ dy= Md+ Ndy = y Condition (necessary and sufficient): M and N have continuous first artial derivatives N u u M and = = = y y y Solution: if M( yd, ) + N(, y) dy = then u = Md+ k ( y) role of constant to determine k( y ) derive u y Note: we could have started with u = Ndy+ l 3 3 E. + 3y d+ 3y + y dy = M N Verification: = = 6y y where k( y ) lays the to get dk and integrate. dy 3 First integration: u = ( 3 + 3y ) d+ k ( y) 4 + y + k( y) Determination of k: General solution: 4 4 u dk 3 dk y = + = + = = + y y y y k C y dy dy y + y = C 4 4
12 E. (sin cosh yd ) (cossinh ydy ) Verification: y = 3 =, with IC M N = sinsinh y = y Integration: u = ( sin cosh y) d+ k ( y) = coscosh y + k ( y) u dk Determination of k: = cossinh y + = N = cossinh y y dy dk = k = C dy te u, y = coscosh y+ C and since du = u = C coscosh y = C IC: coscosh3= C.7 coscosh y =.7
13 Reduction to eact form using integrating factors If P( y) d, + Q, y dy = not EDE Multilication by integrating factor (, ) Note: may eist more than one integrating factor F y makes DE eact FPd+ FQdy = Conditio n: FP = FQ FP y + FPy = FQ + FQ y df F = F Fy = F = d Easy to solve if and df df P Q FPy = Q+ FQ = R d F d Q = y Th. if right side deends only on then F = e R d Note: similarly we could have worked with F( y ) for convenience sake E. + = with IC: y sin y d ycosy dy = π Verification: M y N = 4ycosy = ycos y not eact DE We have P = sin y P = 4y cos y and y Q= ycos y Q = ycos y Since R P Q 3ycosy 3 = Q y R = 4ycosy ycos y ycosy = = ycos y Therefore F = e R d e e[ 3ln ] 3 F = d = = 3 3
14 3 4 Multilication by integrating factor: Verification: M sin y d+ ycosy dy = N 3 = 4ycos y = y, eact Solution: u = 3 sin( y ) d+ k ( y) = 4 sin y + k( y) = + = = = u 4 cos 4 cos y y y k y y y k y k y C 4 u(, y) = sin y = C π = = = = IC: 4 4 sin 8 C u, y siny 8 4
15 Linear DE Linear DE: y + y= r Homogeneous linear DE: r( ) = Solution: searation of variables If non homogeneous: ( y r) d+ dy = Pd+ Qdy = df P Q = F e d F d Q = = y Integrating factor: d d d e ( y + y) = e y e = r d d e y = e rd+ C Putting h = h h d a general solution is: y = e erd+ C E. y y = e = and = + e y e e d C r = e h = d = { } { } y = e e e d+ C = e e + C = Ce + e More simly using integration factor = y ye = e e ye = e ye = e + C y= e + Ce ( ) e 5
16 E. y + y= e ( 3sin+ cos) = h= = + + = + = + 3 y e e e 3sin cos d C e e sin C Ce e sin Integration: + = + By art: udv = uv vdu The second integral yields: e 3sin cos d e 3sind e cosd 3 3 = = 3 and = cos = sin u e du e dv v + = + = e 3sind e cosd e 3sind e sin e 3sind e sin = + = + 3 y e e sin C Ce e sin E. 3 y tan sin y = + = with IC: tan tan ln sec h h = h = d = e = sec e = cos r = sin= sincos h er = sec sin cos = sin = + = IC: y cos sin d C Ccos cos = = C C 3 = y 3cos cos Reduction to linear form Bernoulli equation: + = y y g y Where a is any real number for a = or a = yields linear DE a Transformation: u = y a u ( a) y a = y = ( a) y a ( gy a y ) u = ( a)( g y a ) ( ) ( ) u + a u = a g 6
17 E. Verhulst equation (logistic oulation model) y Ay = By where A and B are ositive constant Bernoulli form with a = u = y u = y y = y By + Ay = B Ay u + ua= B, which is linear Using revious method: = A h= A r = B B A A A A A A u = e Be d+ C = e e + C = Ce + General solution: y = = u B + Ce A A B A Interretation: utting For B = the oulation grows eonentially The breaking term Initially small oulations y oulations y = t we find the logistic law of oulation growth At y = e (law of Malthus) C By revent the oulation to grow without bound A < < increases monotone to A B B A > decrease to same value B while large 7
18 Autonomous DE In the logistic equation, time, the indeendent variable, does not aear elicitly y = f ( t, y) = y f ( y) = = autonomous DE Autonomous DE has constant solutions = equilibrium solutions (or equilibrium oints), te determined by zeros (critical oints) of f ( y ), because f ( y) = y = y = C An equilibrium solution can be stable solutions close to it for some t remain close to it for all further t (e y = 4 in logistic equation) unstable solutions initially close to it do not remain close to it as t increases (e. y = in logistic equation) E. hase line lot y = y y Stable equilibrium solution at y = Unstable equilibrium at y = We can condense the direction field (A) to a hase line lot (C) giving y and y and direction of arrows (B) showing equilibrium 8
19 Orthogonal trajectories Imortant class of roblem in hysics = find family of curves that intersect a given family of curves at right angle = orthogonal trajectories E. Isotherms, heat flows, equiotential curves, curves of steeest descent, etc. Ste : find DE for which the given curves are solution curves y = f ( y, ) Ste : write down DE of orthogonal trajectories Ste 3: solve DE of orthogonal trajectories y = f ( y, ) E. one arameter family curves Fyc,, = y c = or y = c family of arabola Ste : y = c not good because we must eliminate arameter y y = c y = y = y 3 Ste : DE of orthogonal trajectories y = y Ste 3: orthogonal trajectories ellises + y = c = = + one arameter family of ydy d y c 9
20 Eistence and uniqueness of solutions For DE, in most cases, general solutions eists For an IVrb there eist 3 ossibilities: ) No solution ) Only one solution 3) Infinity of solution Problem of eistence: conditions for at least one solution Problem of uniqueness: conditions for at most one solution Th. eistence theorem If f ( y, ) continuous at all oints (, ) y in R: < a and y y < b and bounded in R: f (, y) K then the IVrb has at least one solution defined for all in interval < α where α smaller than two numbers a and b K Th. Uniqueness theorem If f (, ) f M y y and f y continuous for all (, ) y in R and bounded f K and then the IVrb has at most one solution y defined at least for all in interval < α. a) b a α = a; b) b < a α = b, for larger or smaller s nothing can be deduced K K K
21 Linear DE of nd order Many imortant alications in hysics Theory tyical of that of linear DE in general etension to higher order is straightforward General form: y + y + q y= r Homogeneous: y + y + q y= Solution: function h( ) in interval I: a< < b and with continuous derivatives h and h in I E. y y = with y= e or y= e as solutions General solution = linear combination: y = ce + ce Suerosition or linearity rincile Th. Fundamental theorem for the homogeneous linear DE nd For homogeneous linear DE of nd any linear combinations of solutions (or sum and constant multiles) on an oen interval are solutions
22 IVrb - General solution as basis Solution: y= cy + cy IC conditions: y( ) = k and y = k E. y y = with IC: y = 4 and Alying IC to solution: y= ce + ce y = c + c = 4 y = ce ce y = c c = c = and c = 3 y = General solution: y= e + 3e A general solution of homogeneous linear DE of nd order is a solution of the tye y= cy + cy, where y ky The two solutions form a basis y ( ) and y ky ky + = imlies k = k = are linearly indeendent E. 3 y y = basis: y = e and y= e solution: y = ce + ce y + y = basis: y = cos and y = sin solution: y= ccos+ csin
23 Reduction of order Method to find y when y is known y + y + qy = Set y = uy y= uy+ uy y = uy+ uy+ uy Substitute in DE: y + y + qy = ( ) ( ) uy u ( y y ) u( y y qy ) uy + uy + uy + uy + uy + quy = = y + y u + u = y y Or U + + U = y du y Searation of variables: = d lnu ln y d U e + = U y = y Because U = u and y = uy y = y Ud y Since u Ud y = =, the ratio is not constant and the two solutions form a basis. E. 4 y y + = or y + = y y y By insection y Then = since y = and y = = d= ln ln And U = e = y = Ud = ln General solution: y = c + c ln d A articular solution describes the unique behavior of a given hysical system. If, q and r are continuous function on I then there always eist a solution on I which is unique (no singular solution eists) 3
24 Homogeneous linear DE of nd order with constant coefficients General form: y + ay + by = As a solution we try: y= e y = λe y = λ e λ λ λ Relacing in the equation: ( λ λ ) Solution ( λ aλ b) + a + b e λ = + + = (characteristic equation) Roots: ( a a 4b) ( a a 4b) λ = + λ = So the basis are: y = e( λ ) and y = e( λ ) Deending on the sign of the discriminant ) > distinct real roots ) = real double root 3) < comle conjugate roots a 4b 3 solutions ossible: λ λ Case distinct real roots y = ce + ce E. DE: y y y + = IC: y = 4 and λ = ( + 9) = λ = 9 = y = ce + ce y = ce ce y = 5 IC: y = c + c = 4 y = c c = 5 c = and c = 3 y= e + 3e 4
25 Case real double root a λ = λ = λ = A first solution y e a = Second solution: y = uy y = uy + uy y = uy + uy + uy ( ) ( ) uy + uy + uy + auy + uy + buy =, which reduces to uy + u y + ay = a But since y = y uy = u = u = c + c Possible second solution is y y = and general solution is = ( + ) y c ce a E. y + 8y + 6y = Characteristic equation: General solution: = ( + ) λ y c ce + 8λ+ 6= double root: λ = 4 4 E. 3 IVrob: DE: y 4y 4y Characteristic equation: + = CI: y = 3 and λ 4λ + 4 = λ λ = General solution: = ( + ) = + ( + ) y c ce y ce c ce = + = = 5 IC: y = c = 3 y c c c y= 3 5 e Particular solution: y = 5
26 Case 3: a E. Basis: 4b < comle roots y + y = λ + = λ = λ =± =± i i e and e i Using Euler formula: i cos = e + e i sin = e e i i i i e = cos+ isin i e = cos isin cos and sin are also solutions Comle eonential function ( cos sin ) z s+ it s it s e = e = ee = e t+ i t Obtained from Maclaurin series of 3 4 e it = + it + it + it + it +! 3! 4! t t t t t t e it = i t + +! 4! 6! 3! 5! 7! e Since 4 6 t t t cost = + + and! 4! 6! t t t sint = t + + 3! 5! 7! it e = cost+ isin t 6
27 Case 3 (continue): comle roots To make radicand ositive: λ = a+ a 4b = a+ i 4b a = a + i b a 4 λ = a+ iω and λ = a iω where ω = b a 4 a a + iω ( ω ω ) λ e = e = e cos + isin a a iω ( ω ω ) λ e = e = e cos isin a Adding and dividing by : y = e cosω a Subtracting and dividing byi : y = e sinω General solution: y= e ( A cos ω+ B sin ω) a E. DE: y +.y + 4.y = IC: y = y = Characteristic equation: λ +.λ+ 4.= Comle roots:.± i ω = General solution:. y= e Acos+ Bsin IC: y = A= (.sin cos ) = +. y Be y = B = B = Particular solution:. y= e sin Interretation: damed vibration amlitude of oscillation decreases eonentially 7
28 Summary of cases to 3 Case Roots Basis General solution Distinct real λ, λ e, e y = ce + ce λ λ λ λ Real double a a λ = a e, e = ( + ) y c c e a Comle conjugate 3 λ = a+ iω λ = a iω e e a a cosω sinω ( cosω sinω ) a y= e A + B Note: in mechanics or electrical circuits, these 3 cases corresond to 3 different forms of motions or flows current 8
29 Differential oerator The differential can be seen as a transformation writing DE as Dy = y differential oerator D that oerates on y E. D = or D( sin) = cos D Dy = D y = y or Higher derivative: Second order DE oerator: L y = D + ad+ b y = y + ay + by Dy L= P D = D + ad+ b The oerator is linear: L[ αy+ βw] = αl( y) + βl( w) Homogeneous differential equation: L( y ) = Since λ D e = e D e = λ e λ λ λ λ = ( λ + λ+ ) = ( λ) P D e a b e P e Since λ λ λ e λ is solution P ( λ ) = Case : P ( λ ) different real roots Case : P ( λ ) double real root needs second indeendent solution 9
30 Start with = ( λ + λ+ ) = ( λ) P D e a b e P e Differentiating: Because P dp dλ For double root: λ λ λ = ( λ) + ( λ) = λ λ λ λ P D e P e P e P D e = ( λ) + ( λ) = ( λ+ ) + ( λ + λ+ ) P e P e a e a b e λ λ λ λ ( λ λ ) ( λ ) = + e + a + e + be = D + ad+ b e a e λ is second solution λ λ λ λ λ = P( λ) P ( λ) P( D)( e λ ) Since P ( λ ) is a olynomial P( D) = = = is olynomial oerator E. Factorize P( D) = D + D 6 and solve P( D)[ y ] = D D D D + 6= + 3 D y= y y D+ 3 D y= D+ 3 y y = y + y 6y Since D+ 3 y= y + 3y D D+ 3 y= D y + 3y = y + y 6y Similarly: 3 Solutions: ( D+ 3y) = y = e and D y = y = e + 3
31 Modeling: free oscillation (mass-sring system) Physical assumtions: mass of body >> mass of sring Orientation: ositive dislacements downward Hooke s law: F = ks where k > sring modulus and s stretch the smaller k the larger s Assume static equilibrium: F balance the weight W Alying Newton law: F + W = ks+ mg = = mg Putting y = at osition of equilibrium y measures dislacement of body from equilibrium osition Eeriment: ulling weight downward restoring force F Daming: dissiation of energy of mechanical system = ky 3
32 k Ignoring daming: my + ky= y + y = m yt = Acosω t+ Bsinω t where ω = Solution: Harmonic oscillation Putting B = sinδ and A = cosδ Solution transformed into yt = Ccos( ω t δ) k m B C = A + B and tanδ = A Period: π ω ω Frequency: π (Hertz) 3 ossible cases deending on initial velocity: y = or y =± V 3
33 Including daming dy daming force y = F = cy dt Newton s law: c k my + cy + ky= y + y + y = m m Characteristic equation: Roots: c λ λ, = ± c m m c k + λ + = m m 4mk c Putting α = and m β = c m λ = α+ β and λ = α β 4mk Formal solution deends on amount of daming: Case : over-daming Case : critical daming c > 4km c = 4km Case 3: under-daming c < 4km 33
34 Solution in case (Over-daming): ( α β ) t ( α β ) yt = ce + ce + t Both terms aroach zero as t no oscillation 34
35 Solution in case (Critical daming): = ( + ) Term ( c ct) yt c ct e αt + allows one zero system ass at most by y = osition Solution in case 3 (Under-daming): αt αt = ( cosω + sinω ) = cos ( ω δ) y t e A t B t Ce t Where k c m m 4m ω = 4mk c = Amlitude varies between curves y=± Ce αt. ω The smaller is c, the higher the frequency. At the limit c ω ω π 35
36 Euler-Cauchy equation General form: + + = where a and b are constant y ay by Trying the solution: = = = ( ) m m m mm + am + b = y y m y m m m m m m m m m m + am + b = Droing m + a m+ b= m factor m m Case : distinct real roots y = and y = m m General solution: y = c + c E..5. = y y y m m = 3.5. m =.5 and m = 4 c General solution: y= + c 4 36
37 Case : double root ( a) m ( a) = m + ( a ) m+ b 4 where b = ( a) First solution: y = a for second solution, substituting y = uy in y+ ay+ by = ( ) ( ) = uy uy uy a uy uy buy uy + u y + ay + u y + ay + by = y = a = a y y + ay = y ( a) Since last term is zero and ( ) ( ) u ( u + u ) y = = ln u = ln u = u= ln u Second solution: y = y ln General solution: y= ( c + c ) ( ) ln a E = y y y ( a) = m =, and solution = ( + ) y c cln 37
38 Case 3 comle conjugate roots: m = µ + iν and m ln i iν ν iνln = e = e = µ iν Using Euler formula: ( cos( νln ) sin( ν ln )) m µ iν µ iν ln µ = = e = + i ( cos( ν ln ) sin( νln )) m µ iν µ iν ln µ = = e = i Usual transformation yields real solutions: µ µ = ( ν ) and y = ( ν ) y cos ln sin ln E = y y y + + = = ± = ± m 6m 3 m, i 3 General solution: y = Acos( ln) + Bsin( ln) 38
39 Eistence and Uniqueness theory: Wronskian Homogeneous linear DE of nd : y + y + q y= General solution: y= cy + cy IVrob: y( ) = K and y = K Th. Eistence and uniqueness for IVrob If ( ) and q continuous on I and in I then eist solution y( ) on I Linear indeendent solutions Wronskian ky + ky = k = k = Linear indeendence: Wronskian determinant: (, ) y y W y y = = yy yy Th. linear deendence and indeendence Linear deendent W = for = Linear indeendent W Th. 3 general solution y y If ( ) and q continuous on oen interval I eists general solution Th. 4 Form of general solution Under conditions Th. 3, general solution Y cy cy basis and coefficients are suitable constants. No singular solution eists. = + where y and y are 39
40 Nonhomogeneous DE of nd General form: y + y + q y= r Th. Relation between nonhomogeneous DE and homogeneous DE a) The difference of solutions of an nonhomogeneous DE on some oen interval I is a solution of the equivalent homogenous DE on I b) The sum of a solution of a nonhomogeneous DE on I and of a solution of the equivalent homogenous DE on I is a solution of the nonhomogeneous DE General solution of nonhomogeneous DE on oen interval I is of the form: y = y + y h Where y = cy + cy = general solution of homogeneous DE and h any solution of nonhomogeneous DE without arbitrary constants Particular solution: c and c are values defined in y ( ) h y = A general solution includes all solutions If coefficients of nonhomogeneous DE of nd and r( ) are continuous in I Solution eists because both yh and y Once y established, IVrb solution is unique eist on I Th. Suose that coefficient and r continuous on I then every solution of nonhomogeneous DE is obtained by assigning suitable values to arbitrary constants in general solution of nonhomogeneous DE Practical conclusion: to solve a nonhomogeneous DE of nd order or IVrb: ) Solve the homogenous DE ) Find articular solution y ( ) 4
41 E. DE: y y y.4e + + = IC: y =. and y =.9 Ste solve homogenous DE λ λ λ, + + = = ± i ( cos sin ) = + yh e A B Ste determine articular solution of nonhomogeneous DE We try y = Ce + + Ce =.4e C =. y= y + y = e Acos+ Bsin +.e General solution: h Ste 3 articular solution y = A+.=. A= y = e cos Bsin sin+ Bcos +.e y = + B+.=.9 B= Particular solution: y= e cos+.e 4
42 Solution by undetermined coefficients Considering the form y + ay + by = r Choose for substituting in solution. Rules: y a form similar to r involving unknown coefficients to be determined by A) Basic If r( ) is one of the function in table below choose of table and determine coefficients by substituting y in second column y + derivatives in DE Terms in r( ) Choice for y ke γ n k ( n =,,,... ) ke ke kcosω ksinω α α cosω sinω Ce γ n K + K K+ K n n n K cosω+ Msinω ( cosω + sinω ) α e K M B) Modification rule if a term in y haens to be a solution of homogeneous DE then multily by and (double root) C) Sum rule if r( ) is a sum of a functions in table then use sum of y 4
43 E. y + 4y= 8 y = K + K+ K y = K Substituting: ( ) K + 4 K + K+ K = 8 Equating coefficients: 4K = 8, 4K =, K + 4K = K =, K = and K = y = General solution: y= y + y = A + B + h cos sin E. y 3y + y= e Characteristic equation: λ λ λ λ 3 + = = and = y = ce + ce h Choice for y cannot be Ce because this is solution of homogenous DE Trying (rule B) y = Ce y = C ( e + e ) and y = C( e + e ) Substituting: C + e 3C + e + Ce = e C= General solution: y = ce + ce e 43
44 E. 3 IVrb DE: y + y + y= ( D+ ) y = e y = and Characteristic equation: ( λ ) h y = c + ce Rule B (double root) + = λ = (double root) y = C e y = C( ) e = ( 4 + ) Substituting: and terms dro out y C e C 4+ e + C e + C e = Ce = e C = y c ce e General solution: = ( + ) + IC: y = c = and y c c c = y e = = = y = 44
45 E. 4 DE:.5 y + y + 5y=.5e + 4cos4 55sin4 IC: y =. y = 6. Characteristic equation: λ + λ + 5= λ = ± i ( cos sin ) = + yh e A B, From the table:.5 y = Ce + Kcos4+ Msin4 = +.5 y.5 Ce 4Ksin4 4M cos4 =.5 y.5 Ce 6Kcos4 6Msin4.5 Substituting: ( + + ) + ( + + ) + ( + ).5 5 Ce 6K 8M 5K cos4 6M 8K 5M sin4 Equating coefficients: 6.5C =.5 C =. K+ 8M = 4 K = and M = 5 8K M = 55.5 y= e Acos+ Bsin +.e + 5sin4 General solution: y = A+.=. A= IC: = y e Bsin Bcos.e cos4 y = B+.+. = 6. B= Particular solution:.5 y= e sin+.e + 5sin4 45
46 Solution by variation of arameters Standard form: y + y + q y= r With, q, and r arbitrary variable functions continuous in interval I yr yr y y d y d W W Particular solution: = + Where y and y solutions of homogeneous equation and W = yy yy E. y + y = sec Solution to homogeneous equation: y + y= y = cos and y = sin W = coscos sin sin = y = cos sinsecd+ sin cossecd= cos ln cos + sin General solution: = ( + ) + ( + ) y c ln cos cos c sin 46
47 Modeling: forced oscillations. Resonance Free motion of mass on sring homogenous DE: my + cy + ky = Driving force (or inut) r( t ) nonhomogeneous DE: my + cy + ky= r( t) Particular interest: r t = F cosωt Solution = outut or resonse: y = yh + y Method of underdetermined coefficients = cosω + sin sin cos y t = ωa ωt+ ωb ωt = y t ω acosωt ω bsinωt y t a t b ωt Substituting: ( ω ) ω ω ( ω ) k m a+ cb= F ca + k m b = k mω a ωcb cosωt ωca k mω b sinωt = F cosωt Linear system of algebraic equations in two unknowns Solution by elimination (Cha 6) a = F k mω ( ω ) k m + ω c and b = F ωc ( ω ) k m + ω c k Putting ω = with ω > m a = F m m( ω ω ) ( ) + ω ω ω c and b = F m ωc ( ) ω ω + ω c 47
48 Case I: undamed forced oscillation No daming c = F F = = m ( ω ω ) ω k ω y t cosωt cosωt F y= Ccos( ω t δ ) + cosωt ω k ω Interretation: outut = suerosition of harmonic oscillations, one with natural ω frequency π and the other with the frequency of the inut ω π F Maimum amlitude: a = k ρ where ρ = ω ω Ratio of amlitude of ρ a y and inut: = k F is the resonance factor 48
49 Resonance henomenon: For ω ω ρ a F In case of resonance DE: y + ω y= cosωt m From modification rule: F y t = tsinωt becomes larger and larger mω Without daming, vibrations could destroy the system 49
50 Beats: but very close, for IC: y y If ω ω Particular solution: Or ( ω ω ) = = F yt = t t m ( cosω cosω ) ω ω F ω + ω ω ω = m yt sin t sin t Since ω ω small, eriod of last sine becomes very large this roduces the beat henomenon (amlitude varies harmonically) 5
51 Case II: damed forced oscillation Daming c > h αt ( cosω sinω ) y t = e A t+ b t with c ω = m The solution aroach zero after a sufficiently long time Solution of nonhomogeneous DE = transient solution: y = yh + y Aroach steady state solution harmonic oscillation y after sufficiently long time outut becomes This is the normal behavior of real hysical systems no daming is ideal case Practical resonance: in damed case amlitude stay finite as ω ω, becomes maimum for some frequency ω ( c) some inut may destroy the system Amlitude of y cos( ω η) y t = C t where C = a + b = F m ( ) ω ω + ω c b and tanη = = a m dc dω ωc ( ω ω ) Maimum m ( ) ( ω ω ) m = c = + = ω ω c ω For large daming ω increases c > m ω = mk no solution C decreases in monotone way as 5
52 If c m ω one solution c For ω = ωma = ω increases as c decreases aroaching ω as c m Amlitude C ( ω ) = Since C ( ω) ma mf c 4m ω c always finite for c > increases as c decreases and aroaches as c Phase lag η : π ω < ω η < π ω = ω η = π ω > ω η > Note: Whereas DE of nd order have various alications, higher order DE occur much more rarely in hysics (E. binding of elastic beams) The theory is much similar to that of DE of nd order 5
53 ) Transformation v = ay+ b+ k 4y+ 5 y + y+ 3= E. 7 The terms 4y and y v= y y= v y = ( v ) = Substitution: ( v) ( v ) ( v) [ v+ 5] ( v ) + v + 3= ( v+ 5) ( v+ 5) v = v 6 4v+ 4v+ ( v+ 5) v = 4v+ dv d dv d 4v = 4v = + + dv= d+ C v ln 4v+ = + C 4v + 4 Substitution ( ) v= y y ln 4 8y+ = + C 4 4 8y ln 4 8y+ = 8+ C 4+ 8y+ ln 4 8y+ = C 53
f ax ; a 0 is a periodic function b is a periodic function of x of p b. f which is assumed to have the period 2 π, where
(a) (b) If () Year - Tutorial: Toic: Fourier series Time: Two hours π π n Find the fundamental eriod of (i) cos (ii) cos k k f a ; a is a eriodic function b is a eriodic function of of b. f is a eriodic
More informationChapter 2 Second-Order Linear ODEs
Chapter 2 Second-Order Linear ODEs Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2.1 Homogeneous Linear ODEs of Second Order 2 Homogeneous Linear ODEs
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationOrdinary Differential Equations (ODEs)
c01.tex 8/10/2010 22: 55 Page 1 PART A Ordinary Differential Equations (ODEs) Chap. 1 First-Order ODEs Sec. 1.1 Basic Concepts. Modeling To get a good start into this chapter and this section, quickly
More informationBrief Notes on Differential Equations
Brief Notes on Differential Equations (A) Searable First-Order Differential Equations Solve the following differential equations b searation of variables: (a) (b) Solution ( 1 ) (a) The ODE becomes d and
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationSecond-Order Linear ODEs
C0.tex /4/011 16: 3 Page 13 Chap. Second-Order Linear ODEs Chapter presents different types of second-order ODEs and the specific techniques on how to solve them. The methods are systematic, but it requires
More informationENGI 2422 First Order ODEs - Separable Page 3-01
ENGI 4 First Order ODEs - Separable Page 3-0 3. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with
More informationName Class. 5. Find the particular solution to given the general solution y C cos x and the. x 2 y
10 Differential Equations Test Form A 1. Find the general solution to the first order differential equation: y 1 yy 0. 1 (a) (b) ln y 1 y ln y 1 C y y C y 1 C y 1 y C. Find the general solution to the
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More information3 More applications of derivatives
3 More alications of derivatives 3.1 Eact & ineact differentials in thermodynamics So far we have been discussing total or eact differentials ( ( u u du = d + dy, (1 but we could imagine a more general
More informationAPPLIED MATHEMATICS. Part 1: Ordinary Differential Equations. Wu-ting Tsai
APPLIED MATHEMATICS Part 1: Ordinary Differential Equations Contents 1 First Order Differential Equations 3 1.1 Basic Concepts and Ideas................... 4 1.2 Separable Differential Equations................
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
More informationME scope Application Note 16
ME scoe Alication Note 16 Integration & Differentiation of FFs and Mode Shaes NOTE: The stes used in this Alication Note can be dulicated using any Package that includes the VES-36 Advanced Signal Processing
More informationF(p) y + 3y + 2y = δ(t a) y(0) = 0 and y (0) = 0.
Page 5- Chater 5: Lalace Transforms The Lalace Transform is a useful tool that is used to solve many mathematical and alied roblems. In articular, the Lalace transform is a technique that can be used to
More informationShort Solutions to Practice Material for Test #2 MATH 2421
Short Solutions to Practice Material for Test # MATH 4 Kawai (#) Describe recisely the D surfaces listed here (a) + y + z z = Shere ( ) + (y ) + (z ) = 4 = The center is located at C (; ; ) and the radius
More informationTHE FIRST LAW OF THERMODYNAMICS
THE FIRST LA OF THERMODYNAMIS 9 9 (a) IDENTIFY and SET UP: The ressure is constant and the volume increases (b) = d Figure 9 Since is constant, = d = ( ) The -diagram is sketched in Figure 9 The roblem
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More information1 Entropy 1. 3 Extensivity 4. 5 Convexity 5
Contents CONEX FUNCIONS AND HERMODYNAMIC POENIALS 1 Entroy 1 2 Energy Reresentation 2 3 Etensivity 4 4 Fundamental Equations 4 5 Conveity 5 6 Legendre transforms 6 7 Reservoirs and Legendre transforms
More informationSecond-Order Linear ODEs
Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationENGI 3424 First Order ODEs Page 1-01
ENGI 344 First Order ODEs Page 1-01 1. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with respect
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More informationMA 266 Review Topics - Exam # 2 (updated)
MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential
More informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationIntroductory Differential Equations
Introductory Differential Equations Lecture Notes June 3, 208 Contents Introduction Terminology and Examples 2 Classification of Differential Equations 4 2 First Order ODEs 5 2 Separable ODEs 5 22 First
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationBifurcation-Based Stability Analysis of Electrostatically Actuated Micromirror as a Two Degrees of Freedom System
Coyright 8 Tech Science Press CMES, vol.4, no.,.6-76, 8 Bifurcation-Based Stability Analysis of Electrostatically Actuated Micromirror as a Two Degrees of Freedom System Kuntao Ye, *, Yan Luo and Yingtao
More informationf(x) p(x) =p(b)... d. A function can have two different horizontal asymptotes...
Math Final Eam, Fall. ( ts.) Mark each statement as either true [T] or false [F]. f() a. If lim f() =and lim g() =, then lim does not eist......................!5!5!5 g() b. If is a olynomial, then lim!b
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More informationThe Second Law: The Machinery
The Second Law: The Machinery Chater 5 of Atkins: The Second Law: The Concets Sections 3.7-3.9 8th Ed, 3.3 9th Ed; 3.4 10 Ed.; 3E 11th Ed. Combining First and Second Laws Proerties of the Internal Energy
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationSecond order linear equations
Second order linear equations Samy Tindel Purdue University Differential equations - MA 266 Taken from Elementary differential equations by Boyce and DiPrima Samy T. Second order equations Differential
More informationME 375 System Modeling and Analysis. Homework 11 Solution. Out: 18 November 2011 Due: 30 November 2011 = + +
Out: 8 November Due: 3 November Problem : You are given the following system: Gs () =. s + s+ a) Using Lalace and Inverse Lalace, calculate the unit ste resonse of this system (assume zero initial conditions).
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationChapter 1. Introduction
I. Classical Physics Chater 1. Introduction Classical Mechanics (Newton): It redicts the motion of classical articles with elegance and accuracy. d F ma, mv F: force a: acceleration : momentum q: osition
More informationAPPM 2360: Midterm 3 July 12, 2013.
APPM 2360: Midterm 3 July 12, 2013. ON THE FRONT OF YOUR BLUEBOOK write: (1) your name, (2) your instructor s name, (3) your recitation section number and (4) a grading table. Text books, class notes,
More information2 nd ORDER O.D.E.s SUBSTITUTIONS
nd ORDER O.D.E.s SUBSTITUTIONS Question 1 (***+) d y y 8y + 16y = d d d, y 0, Find the general solution of the above differential equation by using the transformation equation t = y. Give the answer in
More informationVisualizing Free and Forced Harmonic Oscillations using Maple V R5*
Int. J. Engng Ed. Vol. 15, No. 6,. 437±455, 1999 0949-149X/91 $3.00+0.00 Printed in Great Britain. # 1999 TEMPUS Publications. Visualizing Free and Forced Harmonic Oscillations using Male V R5* HARALD
More informationPhase transition. Asaf Pe er Background
Phase transition Asaf Pe er 1 November 18, 2013 1. Background A hase is a region of sace, throughout which all hysical roerties (density, magnetization, etc.) of a material (or thermodynamic system) are
More informationMath 2214 Solution Test 1D Spring 2015
Math 2214 Solution Test 1D Spring 2015 Problem 1: A 600 gallon open top tank initially holds 300 gallons of fresh water. At t = 0, a brine solution containing 3 lbs of salt per gallon is poured into the
More informationUniversity of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009
University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental
More informationSetting up the Mathematical Model Review of Heat & Material Balances
Setting u the Mathematical Model Review of Heat & Material Balances Toic Summary... Introduction... Conservation Equations... 3 Use of Intrinsic Variables... 4 Well-Mixed Systems... 4 Conservation of Total
More informationPartial Differential Equations (PDEs)
c12.qxd 7/19/11 2:6 PM Page 238 CHAPTER 12 Partial Differential Equations (PDEs) SECTION 12.1 Basic Concets of PDEs, age 54 Purose. To familiarize the student with the following: Concet of solution, verification
More informationSecond In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 2011
Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 211 (1) [6] Give the interval of definition for the solution of the initial-value problem d 4 y dt 4 + 7 1 t 2 dy dt
More informationElectronic properties of Graphene and 2-D materials. Part 2
Electronic roerties of Grahene and -D materials z Part y The laws of Quantum Mechanics continued Time indeendent Schroedinger equation Eamles Particle in a well Harmonic oscillator Quantum Tunneling The
More informationMath 211. Substitute Lecture. November 20, 2000
1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +
More informationI Poles & zeros. I First-order systems. I Second-order systems. I E ect of additional poles. I E ect of zeros. I E ect of nonlinearities
EE C28 / ME C34 Lecture Chater 4 Time Resonse Alexandre Bayen Deartment of Electrical Engineering & Comuter Science University of California Berkeley Lecture abstract Toics covered in this resentation
More informationMATH 4B Differential Equations, Fall 2016 Final Exam Study Guide
MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions
More informationE p,rms = 240 V E rms = 120 V N p N s C. f = 60 Hz R = 3.8 L
Discussion Question 1A P1, Week 1 Power in AC Circuits An electronic device, consisting of a simle C circuit, is designed to be connected to an American-standard ower outlet delivering an EMF of 1 V at
More informationChapter 2 Notes, Kohler & Johnson 2e
Contents 2 First Order Differential Equations 2 2.1 First Order Equations - Existence and Uniqueness Theorems......... 2 2.2 Linear First Order Differential Equations.................... 5 2.2.1 First
More informationODE. Philippe Rukimbira. Department of Mathematics Florida International University PR (FIU) MAP / 92
ODE Philippe Rukimbira Department of Mathematics Florida International University PR (FIU) MAP 2302 1 / 92 4.4 The method of Variation of parameters 1. Second order differential equations (Normalized,
More informationpp physics, RWTH, WS 2003/04, T.Hebbeker
1. PP TH 03/04 Accelerators and Detectors 1 hysics, RWTH, WS 2003/04, T.Hebbeker 2003-12-03 1. Accelerators and Detectors In the following, we concentrate on the three machines SPS, Tevatron and LHC with
More informationQ ABS (x,t s ) = S o /4 S(x)a p (x,t s ).
Lecture 14 ICE The feedback of exanding and contracting ice sheets has often been offered as a lausible exlanation for how the contrasting climates of the glacialinterglacial times can be both (relatively)
More informationMath Spring 2014 Homework 2 solution
Math 3-00 Spring 04 Homework solution.3/5 A tank initially contains 0 lb of salt in gal of weater. A salt solution flows into the tank at 3 gal/min and well-stirred out at the same rate. Inflow salt concentration
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.
More informationSelected Topics in Physics a lecture course for 1st year students by W.B. von Schlippe Spring Semester 2007
Selected Topics in Physics a lecture course for st year students by W.B. von Schlippe Spring Semester 7 Lecture : Oscillations simple harmonic oscillations; coupled oscillations; beats; damped oscillations;
More informationForced Mechanical Vibrations
Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State
More informationCentral Force Motion Challenge Problems
Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has
More informationAnswer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.
Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.
More informationThe most up-to-date version of this collection of homework exercises can always be found at bob/math365/mmm.pdf.
Millersville University Department of Mathematics MATH 365 Ordinary Differential Equations January 23, 212 The most up-to-date version of this collection of homework exercises can always be found at http://banach.millersville.edu/
More informationODE Math 3331 (Summer 2014) June 16, 2014
Page 1 of 12 Please go to the next page... Sample Midterm 1 ODE Math 3331 (Summer 2014) June 16, 2014 50 points 1. Find the solution of the following initial-value problem 1. Solution (S.O.V) dt = ty2,
More informationSolution Week 75 (2/16/04) Hanging chain
Catenary Catenary is idealized shae of chain or cable hanging under its weight with the fixed end oints. The chain (cable) curve is catenary that minimizes the otential energy PHY 322, Sring 208 Solution
More informationChapter 6: Sound Wave Equation
Lecture notes on OPAC0- ntroduction to Acoustics Dr. Eser OLĞAR, 08 Chater 6: Sound Wave Equation. Sound Waves in a medium the wave equation Just like the eriodic motion of the simle harmonic oscillator,
More informationSolutions to Problems in Goldstein, Classical Mechanics, Second Edition. Chapter 9
Solutions to Problems in Goldstein, Classical Mechanics, Second Edition Homer Reid October 29, 2002 Chater 9 Problem 9. One of the attemts at combining the two sets of Hamilton s equations into one tries
More informationCSE 599d - Quantum Computing When Quantum Computers Fall Apart
CSE 599d - Quantum Comuting When Quantum Comuters Fall Aart Dave Bacon Deartment of Comuter Science & Engineering, University of Washington In this lecture we are going to begin discussing what haens to
More information+++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic +++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic
Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic Daniel Leitner, Johannes Krof, Günther Zauner, TU Vienna, Austria, dleitner@osiris.tuwien.ac.at Yuri Karov, Yuri Senichenkov, Yuri
More informationAP Calculus Testbank (Chapter 10) (Mr. Surowski)
AP Calculus Testbank (Chater 1) (Mr. Surowski) Part I. Multile-Choice Questions 1. The grah in the xy-lane reresented by x = 3 sin t and y = cost is (A) a circle (B) an ellise (C) a hyerbola (D) a arabola
More informationA Numerical Method for Critical Buckling Load for a Beam Supported on Elastic Foundation
A Numerical Method for Critical Buckling Load for a Beam Suorted on Elastic Foundation Guo-ing Xia Institute of Bridge Engineering, Dalian University of Technology, Dalian, Liaoning Province, P. R. China
More informationAdam Paweł Zaborski. 8 Plasticity. reloading. 1. Bauschinger s effect. 2. unchanged yielding limit. 3. isotropic hardening
8 lasticity Introduction Definitions loading/unloading/reloading words commonly used in lasticity lastic strain a ermanent strain that doesn t vanish after unloading, lastically assive rocess a rocess
More informationDrill Exercise Differential equations by abhijit kumar jha DRILL EXERCISE - 1 DRILL EXERCISE - 2. e x, where c 1
DRILL EXERCISE -. Find the order and degree (if defined) of the differential equation 5 4 3 d y d y. Find the order and degree (if defined) of the differential equation n. 3. Find the order and degree
More informationIntroduction to Landau s Fermi Liquid Theory
Introduction to Landau s Fermi Liquid Theory Erkki Thuneberg Deartment of hysical sciences University of Oulu 29 1. Introduction The rincial roblem of hysics is to determine how bodies behave when they
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationHIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations
HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations David Levermore Department of Mathematics University of Maryland 22 October 28 Because the presentation of this
More informationSimplifications to Conservation Equations
Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,
More informationApplications of First Order Differential Equation
Dr Mansoor Alshehri King Saud University MATH204-Differential Equations Center of Excellence in Learning and Teaching 1 / 39 Orthogonal Trajectories How to Find Orthogonal Trajectories Growth and Decay
More informationHIGHER-ORDER LINEAR DIFFERENTIAL EQUATIONS
CHAPTER 5 HIGHER-ORDER LINEAR DIFFERENTIAL EQUATIONS SECTION 5. INTRODUCTION: SECOND-ORDER LINEAR EQUATIONS In this section the central ideas of the theory of linear differential equations are introduced
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationAdiabatic Shear Bands in Simple and Dipolar Plastic Materials
Adiabatic Shear Bands in Simle and Diolar Plastic Materials T W \-1RIGHT us Army Ballistic Research Laboratory Aberdeen Proving Ground, MD 215 R C BATRA University of Missouri-Rolla Rolla, Missouri 6541
More informationwhether a process will be spontaneous, it is necessary to know the entropy change in both the
93 Lecture 16 he entroy is a lovely function because it is all we need to know in order to redict whether a rocess will be sontaneous. However, it is often inconvenient to use, because to redict whether
More informationUniform Law on the Unit Sphere of a Banach Space
Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a
More informationMath 266 Midterm Exam 2
Math 266 Midterm Exam 2 March 2st 26 Name: Ground Rules. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use one 4-by-6
More informationFourier Series Tutorial
Fourier Series Tutorial INTRODUCTION This document is designed to overview the theory behind the Fourier series and its alications. It introduces the Fourier series and then demonstrates its use with a
More informationWeek 8 lectures. ρ t +u ρ+ρ u = 0. where µ and λ are viscosity and second viscosity coefficients, respectively and S is the strain tensor:
Week 8 lectures. Equations for motion of fluid without incomressible assumtions Recall from week notes, the equations for conservation of mass and momentum, derived generally without any incomressibility
More informationWork sheet / Things to know. Chapter 3
MATH 251 Work sheet / Things to know 1. Second order linear differential equation Standard form: Chapter 3 What makes it homogeneous? We will, for the most part, work with equations with constant coefficients
More informationDifferential Equations
Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown
More informationMath Applied Differential Equations
Math 256 - Applied Differential Equations Notes Existence and Uniqueness The following theorem gives sufficient conditions for the existence and uniqueness of a solution to the IVP for first order nonlinear
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationChapter 4. Systems of ODEs. Phase Plane. Qualitative Methods
Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationMath K (24564) - Lectures 02
Math 39100 K (24564) - Lectures 02 Ethan Akin Office: NAC 6/287 Phone: 650-5136 Email: ethanakin@earthlink.net Spring, 2018 Contents Second Order Linear Equations, B & D Chapter 4 Second Order Linear Homogeneous
More information6. Kinematics-Dynamics of Rotational Rest-Motion of a Material Point A potential-kinetic vector of circular motion
3 6. Kinematics-Dynamics of Rotational Rest-Motion of a Material Point 6.1. A otential-inetic vector of circular motion Kinematics and dynamics of rectilinear and rotational motions are still at the stage
More informationMath 104B: Number Theory II (Winter 2012)
Math 104B: Number Theory II (Winter 01) Alina Bucur Contents 1 Review 11 Prime numbers 1 Euclidean algorithm 13 Multilicative functions 14 Linear diohantine equations 3 15 Congruences 3 Primes as sums
More informationSplit the integral into two: [0,1] and (1, )
. A continuous random variable X has the iecewise df f( ) 0, 0, 0, where 0 is a ositive real number. - (a) For any real number such that 0, rove that the eected value of h( X ) X is E X. (0 ts) Solution:
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More information