Study guide
2.4 Differential Equations and Numerical Methods
PRC Civil Engineer Licensure Examination CELE
Track your progress as you read
Finished sections get a check mark.
1 Chapter 1: Examination Overview, RA 544 & Professional Practice
1.1 CELE Structure, Format, Scoring, and 2026 Subject Order 1.2 Republic Act No. 544: The Civil Engineering Law and Board Regulations 1.3 Civil Engineering Ethics, Contracts, and Professional Practice
2 Chapter 2: Engineering Mathematics, Numerical Methods & Physics
2.1 Advanced Algebra, Progressions & Word Problems 2.2 Trigonometry and Analytical Geometry 2.3 Differential and Integral Calculus 2.4 Differential Equations and Numerical Methods 2.5 Series, Hyperbolic Functions, Multivariable Calculus, and Laplace Transforms 2.6 Numerical Methods: Linear Systems, Interpolation, Curve Fitting, and ODE Solvers 2.7 Physics for Engineers: Kinematics, Newton's Laws, Work-Energy, and Momentum
3 Chapter 3: Engineering Data Analysis & Engineering Economy
3.1 Engineering Data Analysis and Statistical Applications 3.2 Parameter Estimation, Linear Regression, Correlation, and Queuing Models 3.3 Engineering Economy: Time Value of Money and Interest Formulas 3.4 Project Evaluation, Depreciation, and Benefit-Cost Analysis
4 Chapter 4: Surveying Fundamentals, Construction Layout & Route Surveying
4.1 Errors, Corrections, and Pacing/Taping 4.2 Leveling, Earthworks, and Mass Haul Diagrams 4.3 Traverse Computations, Balancing, and Area Determinations 4.4 Route Surveying: Horizontal and Vertical Curves
5 Chapter 5: Highway & Transportation Engineering
5.1 Highway Geometric Design and Sight Distances 5.2 Traffic Engineering, Flow Fundamentals, and Pavement Design 5.3 Highway Policy, Route Location, Drainage, and Earthwork Slopes 5.4 Transportation Planning, Network Analysis, and Intermodal Systems
6 Chapter 6: Construction Materials, Methods, Management & Safety
6.1 Materials for Construction: Concrete, Formwork, Masonry, Steel, and Wood Products 6.2 Construction Equipment, Formwork, and Construction Methods 6.3 Construction Occupational Safety and Health: RA 11058 and DOLE DO 13 6.4 Construction Project Management, CPM/PERT, and Resource Scheduling 6.5 Construction Principles: Design Requirements, Project Organization, Cost Control, and Value Engineering 6.6 Quantity Surveying, Unit Cost Analysis, and Project Cost Estimates 6.7 Construction Trades and Specialties: Plumbing, Plastics, Doors and Windows, Finishes, and Electrical
7 Chapter 7: Fluid Properties, Hydrostatics & Buoyancy
7.1 Fluid Properties, Units, and Pressure Measurement 7.2 Hydrostatic Pressure on Plane and Curved Surfaces 7.3 Buoyancy, Flotation, and Accelerated Fluid Masses
8 Chapter 8: Hydrodynamics, Pipe Flow, Open Channels & Water Supply
8.1 Fluid Dynamics: Continuity and Energy Equations 8.2 Flow Measurement: Orifices, Weirs, and Venturi Meters 8.3 Fluid Flow in Pipes, Friction Losses, and Pipe Networks 8.4 Open Channel Flow, Manning's Equation, and Hydraulic Jump 8.5 Water Supply: Demand, Storage, Pumping, Distribution, and Treatment
9 Chapter 9: Soil Properties, Seepage & Stresses in Soil Mass
9.1 Soil Composition, Weight-Volume Phase Relationships 9.2 Soil Classification Systems (USCS & AASHTO) and Atterberg Limits 9.3 Soil Permeability, Seepage, Flow Nets, and Effective Stress 9.4 Stresses in Soil Mass, Stress Distribution, and Immediate and Secondary Settlement
10 Chapter 10: Consolidation, Soil Strength, Earth Pressures & Foundations
10.1 Soil Compaction and Soil Compressibility/Consolidation 10.2 Shear Strength of Soils and Mohr-Coulomb Failure Criteria 10.3 Lateral Earth Pressures (Rankine & Coulomb) and Retaining Walls 10.4 Shallow and Deep Foundation Bearing Capacity and Slope Stability
11 Chapter 11: Engineering Mechanics & Statics
11.1 Force Systems, Resultants, and Equilibrium of Rigid Bodies 11.2 Centroids, Moments of Inertia, and Friction 11.3 Cables, Trusses, and Internal Forces in Members
12 Chapter 12: Strength of Materials
12.1 Axial Stress, Strain, and Deformation 12.2 Torsion, Flanged Bolt Couplings, and Thin-Walled Pressure Vessels 12.3 Beam Stresses: Flexural and Transverse Shear 12.4 Beam Deflections, Combined Stresses, and Mohr's Circle
13 Chapter 13: Structural Analysis & NSCP Load Criteria
13.1 Indeterminate Structures: Approximate and Classical Methods 13.2 Moment Distribution Method and Slope Deflection 13.3 Influence Lines for Determinate and Indeterminate Structures 13.4 NSCP Minimum Design Loads: Dead, Live, Wind, and Earthquake
14 Chapter 14: Reinforced Concrete Design (NSCP 2015 / USD)
14.1 Ultimate Strength Design (USD) Principles and Singly Reinforced Beams 14.2 Doubly Reinforced Beams and T-Beams 14.3 Shear and Diagonal Tension, Stirrup Design, and Development Length 14.4 Reinforced Concrete Columns and Isolated Footing Design
15 Chapter 15: Steel, Timber & Prestressed Design and Materials Testing
15.1 Structural Steel Tension Members and Connections 15.2 Structural Steel Compression Members and Beam Design 15.3 Timber Design Principles, Working Stress, and Sawn Lumber Beams 15.4 Prestressed Concrete Fundamentals and Prestress Losses 15.5 Construction Materials Testing and Concrete Mix Proportioning
2.4 Differential Equations and Numerical Methods
Key Takeaways Explain This Topic AI
First-order ODEs are resolved via separation of variables, integrating factor I(x) = exp(∫ P dx) for linear forms dy/dx + Py = Q, or exactness tests (∂M/∂y = ∂N/∂x).
Newton's Law of Cooling models thermal dissipation in curing mass concrete: T(t) = T_m + (T₀ - T_m)e^(-kt), establishing thermal cracking control.
Second-order linear ODEs with constant coefficients yield three structural vibration regimes: overdamped (Δ > 0), critically damped (Δ = 0), and underdamped harmonic oscillation (Δ < 0).
Root-finding methods balance reliability and speed: Bisection guarantees linear convergence, while Newton-Raphson provides quadratic convergence x_{n+1} = x_n - f(x_n)/f'(x_n).
Simpson's 1/3 Rule requires an even number of intervals n and integrates cubic polynomials exactly with global error O(h⁴), outperforming the Trapezoidal Rule O(h²).
Last updated: October 2026
2.4 Differential Equations and Numerical Methods
Differential equations and numerical methods bridge pure mathematics and applied engineering mechanics in civil engineering. While first-order and second-order ordinary differential equations (ODEs) model hydraulic seepage, thermal dissipation in curing mass concrete, and structural dynamics, practical field problems often resist closed-form analytical solutions. Numerical methods—including root-finding iterations (Bisection, Newton-Raphson) and numerical integration (Trapezoidal, Simpson's 1/3 and 3/8 rules)—provide accurate numerical approximations required on the board examination.
First-Order Ordinary Differential Equations (ODEs)
A first-order differential equation involves independent variable x x x , dependent variable y y y , and the first derivative d y d x \frac{dy}{dx} d x d y .
1. Variable Separable ODEs
The equation can be manipulated algebraically into isolated single-variable differentials: f ( x ) d x + g ( y ) d y = 0 ⟹ ∫ f ( x ) d x + ∫ g ( y ) d y = C f(x)\,dx + g(y)\,dy = 0 \implies \int f(x)\,dx + \int g(y)\,dy = C f ( x ) d x + g ( y ) d y = 0 ⟹ ∫ f ( x ) d x + ∫ g ( y ) d y = C
2. Homogeneous Differential Equations
An equation in the differential form M ( x , y ) d x + N ( x , y ) d y = 0 M(x, y)\,dx + N(x, y)\,dy = 0 M ( x , y ) d x + N ( x , y ) d y = 0 is homogeneous of degree n n n if M ( k x , k y ) = k n M ( x , y ) M(kx, ky) = k^n M(x, y) M ( k x , k y ) = k n M ( x , y ) and N ( k x , k y ) = k n N ( x , y ) N(kx, ky) = k^n N(x, y) N ( k x , k y ) = k n N ( x , y ) .
Standard Substitution: Let y = v x ⟹ d y = v d x + x d v y = vx \implies dy = v\,dx + x\,dv y = vx ⟹ d y = v d x + x d v (or x = v y ⟹ d x = v d y + y d v x = vy \implies dx = v\,dy + y\,dv x = v y ⟹ d x = v d y + y d v ).
Substituting transforms the equation into a separable ODE in terms of v v v and x x x .
3. Exact Differential Equations
The equation M ( x , y ) d x + N ( x , y ) d y = 0 M(x, y)\,dx + N(x, y)\,dy = 0 M ( x , y ) d x + N ( x , y ) d y = 0 is exact if and only if: ∂ M ∂ y = ∂ N ∂ x \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} ∂ y ∂ M = ∂ x ∂ N The general solution is F ( x , y ) = C F(x, y) = C F ( x , y ) = C , found by: F ( x , y ) = ∫ M ( x , y ) ∂ x + g ( y ) , where g ′ ( y ) = N ( x , y ) − ∂ ∂ y [ ∫ M ( x , y ) ∂ x ] F(x, y) = \int M(x, y)\,\partial x + g(y), \quad \text{where } g'(y) = N(x, y) - \frac{\partial}{\partial y}\left[\int M(x, y)\,\partial x\right] F ( x , y ) = ∫ M ( x , y ) ∂ x + g ( y ) , where g ′ ( y ) = N ( x , y ) − ∂ y ∂ [ ∫ M ( x , y ) ∂ x ]
4. First-Order Linear Differential Equations
The standard canonical form is: d y d x + P ( x ) y = Q ( x ) \frac{dy}{dx} + P(x)y = Q(x) d x d y + P ( x ) y = Q ( x ) Where P ( x ) P(x) P ( x ) and Q ( x ) Q(x) Q ( x ) are continuous functions of x x x alone.
Integrating Factor ( I ( x ) I(x) I ( x ) ): I ( x ) = e ∫ P ( x ) d x I(x) = e^{\int P(x)\,dx} I ( x ) = e ∫ P ( x ) d x
General Closed-Form Solution: y ⋅ I ( x ) = ∫ Q ( x ) I ( x ) d x + C y \cdot I(x) = \int Q(x) I(x)\,dx + C y ⋅ I ( x ) = ∫ Q ( x ) I ( x ) d x + C
- Engineering Applications of First-Order Models
- 1. Newton's Law of Cooling (Thermal Curing of Concrete)
Mass concrete structures (such as gravity dams and bridge piers) generate substantial internal hydration heat that must dissipate to the ambient environment. Newton's Law states that the rate of temperature change of a body is proportional to the difference between its temperature T T T and the surrounding ambient temperature T m T_m T m : d T d t = − k ( T − T m ) \frac{dT}{dt} = -k(T - T_m) d t d T = − k ( T − T m ) Separating variables and integrating yields: T ( t ) = T m + ( T 0 − T m ) e − k t T(t) = T_m + (T_0 - T_m)e^{-kt} T ( t ) = T m + ( T 0 − T m ) e − k t Where T 0 T_0 T 0 is the initial temperature at t = 0 t = 0 t = 0 , and k k k is the cooling rate constant ( time − 1 \text{time}^{-1} time − 1 ).
2. Exponential Growth and Decay
Governed by d N d t = k N ⟹ N ( t ) = N 0 e k t \frac{dN}{dt} = kN \implies N(t) = N_0 e^{kt} d t d N = k N ⟹ N ( t ) = N 0 e k t . If k < 0 k < 0 k < 0 , it models radioisotope decay or contaminant breakdown; half-life is t 1 / 2 = ln 2 ∣ k ∣ t_{1/2} = \frac{\ln 2}{|k|} t 1/2 = ∣ k ∣ l n 2 .
Second-Order Linear Homogeneous ODEs with Constant Coefficients
Second-order equations describe dynamic mechanical vibrations, structural frame oscillations, and beam deflections: a d 2 y d x 2 + b d y d x + c y = 0 ( a ≠ 0 ) a \frac{d^2 y}{dx^2} + b \frac{dy}{dx} + c y = 0 \quad (a \ne 0) a d x 2 d 2 y + b d x d y + cy = 0 ( a = 0 )
The Characteristic (Auxiliary) Equation
Assuming trial solution y = e r x y = e^{rx} y = e r x yields the characteristic quadratic equation: a r 2 + b r + c = 0 a r^2 + b r + c = 0 a r 2 + b r + c = 0
Roots are obtained via r = − b ± b 2 − 4 a c 2 a r = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a} r = 2 a − b ± b 2 − 4 a c . The physical and mathematical nature of the response is classified by the discriminant Δ = b 2 − 4 a c \Delta = b^2 - 4ac Δ = b 2 − 4 a c :
Discriminant ( Δ \Delta Δ ) Nature of Characteristic Roots General Solution y ( x ) y(x) y ( x ) Physical Structural Behavior
Δ > 0 \Delta > 0 Δ > 0 Two distinct real roots r 1 ≠ r 2 r_1 \ne r_2 r 1 = r 2 y = C 1 e r 1 x + C 2 e r 2 x y = C_1 e^{r_1 x} + C_2 e^{r_2 x} y = C 1 e r 1 x + C 2 e r 2 x Overdamped: Non-oscillatory exponential decay to equilibrium
Δ = 0 \Delta = 0 Δ = 0 Repeated real root r 1 = r 2 = r = − b 2 a r_1 = r_2 = r = -\frac{b}{2a} r 1 = r 2 = r = − 2 a b y = ( C 1 + C 2 x ) e r x y = (C_1 + C_2 x) e^{r x} y = ( C 1 + C 2 x ) e r x Critically Damped: Fastest non-oscillatory return to rest without overshoot
Δ < 0 \Delta < 0 Δ < 0 Complex conjugate roots r = α ± i β r = \alpha \pm i\beta r = α ± i β y = e α x ( C 1 cos β x + C 2 sin β x ) y = e^{\alpha x}(C_1 \cos\beta x + C_2 \sin\beta x) y = e αx ( C 1 cos β x + C 2 sin β x ) Underdamped: Decaying harmonic oscillations ( α = − b 2 a , β = 4 a c − b 2 2 a \alpha = -\frac{b}{2a}, \beta = \frac{\sqrt{4ac-b^2}}{2a} α = − 2 a b , β = 2 a 4 a c − b 2 )
Numerical Root-Finding Algorithms
When non-linear transcendental equations f ( x ) = 0 f(x) = 0 f ( x ) = 0 (such as pipe friction factor formulas or open channel critical depths) cannot be solved analytically, numerical iterative algorithms are employed.
1. The Bisection Method
Based on the Intermediate Value Theorem : If a continuous function f ( x ) f(x) f ( x ) satisfies f ( a ) ⋅ f ( b ) < 0 f(a) \cdot f(b) < 0 f ( a ) ⋅ f ( b ) < 0 , at least one real root exists within the bracket [ a , b ] [a, b] [ a , b ] .
Iteration Formula: Midpoint c = a + b 2 c = \frac{a + b}{2} c = 2 a + b . If f ( a ) f ( c ) < 0 f(a)f(c) < 0 f ( a ) f ( c ) < 0 , root lies in [ a , c ] [a, c] [ a , c ] ; otherwise root lies in [ c , b ] [c, b] [ c , b ] .
Convergence Rate: Linear (error halves each iteration, convergence factor = 0.5 = 0.5 = 0.5 ).
Number of Iterations for Tolerance ϵ \epsilon ϵ : n ≥ ln ( b − a ϵ ) ln 2 n \ge \frac{\ln\left(\frac{b - a}{\epsilon}\right)}{\ln 2} n ≥ l n 2 l n ( ϵ b − a )
Strengths & Weaknesses: Absolutely guaranteed convergence, but slow compared to gradient methods.
2. The Newton-Raphson Method
Uses the local tangent line approximation at trial point x n x_n x n to project to the next root estimate: x n + 1 = x n − f ( x n ) f ′ ( x n ) x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)} x n + 1 = x n − f ′ ( x n ) f ( x n )
Convergence Rate: Quadratic near a simple root ( ∣ e n + 1 ∣ ≈ M ∣ e n ∣ 2 |e_{n+1}| \approx M |e_n|^2 ∣ e n + 1 ∣ ≈ M ∣ e n ∣ 2 ), doubling the number of significant correct decimal places with each iteration.
Common Failure Modes:
Zero Derivative ( f ′ ( x n ) = 0 f'(x_n) = 0 f ′ ( x n ) = 0 ): Horizontal tangent causes division by zero.
Oscillatory Cycles: Diverges or enters an infinite periodic loop around points of inflection.
Poor Starting Value: Converges to an unintended distant root if x 0 x_0 x 0 is poorly chosen.
Multiple Roots: If a root has multiplicity m > 1 m > 1 m > 1 , convergence slows from quadratic to linear.
Numerical Integration Methods
Numerical quadrature approximates the definite integral I = ∫ a b f ( x ) d x I = \int_a^b f(x)\,dx I = ∫ a b f ( x ) d x using n n n equal panels of step size h = b − a n h = \frac{b - a}{n} h = n b − a across discrete ordinates y i = f ( x i ) y_i = f(x_i) y i = f ( x i ) .
1. The Trapezoidal Rule
Approximates the region under each subinterval with a linear trapezoid: ∫ a b f ( x ) d x ≈ h 2 [ f ( x 0 ) + 2 ∑ i = 1 n − 1 f ( x i ) + f ( x n ) ] \int_a^b f(x)\,dx \approx \frac{h}{2}\left[f(x_0) + 2\sum_{i=1}^{n-1} f(x_i) + f(x_n)\right] ∫ a b f ( x ) d x ≈ 2 h [ f ( x 0 ) + 2 ∑ i = 1 n − 1 f ( x i ) + f ( x n ) ]
Polynomial Exactness: Degree ≤ 1 \le 1 ≤ 1 (exact for linear profiles).
Global Truncation Error: E T = − ( b − a ) h 2 12 f ′ ′ ( ξ ) = O ( h 2 ) E_T = -\frac{(b - a)h^2}{12} f''(\xi) = O(h^2) E T = − 12 ( b − a ) h 2 f ′′ ( ξ ) = O ( h 2 )
2. Simpson's 1/3 Rule
Approximates consecutive pairs of panels with quadratic parabolas.
Strict Constraint: Requires an even number of intervals n n n (an odd number of ordinates n + 1 n + 1 n + 1 ). ∫ a b f ( x ) d x ≈ h 3 [ f ( x 0 ) + 4 ∑ odd i f ( x i ) + 2 ∑ even i f ( x i ) + f ( x n ) ] \int_a^b f(x)\,dx \approx \frac{h}{3}\left[f(x_0) + 4\sum_{\text{odd } i} f(x_i) + 2\sum_{\text{even } i} f(x_i) + f(x_n)\right] ∫ a b f ( x ) d x ≈ 3 h [ f ( x 0 ) + 4 ∑ odd i f ( x i ) + 2 ∑ even i f ( x i ) + f ( x n ) ]
Polynomial Exactness: Degree ≤ 3 \le 3 ≤ 3 (exact for cubics as well as quadratics due to error cancellation!).
Global Truncation Error: E S = − ( b − a ) h 4 180 f ( 4 ) ( ξ ) = O ( h 4 ) E_S = -\frac{(b - a)h^4}{180} f^{(4)}(\xi) = O(h^4) E S = − 180 ( b − a ) h 4 f ( 4 ) ( ξ ) = O ( h 4 )
3. Simpson's 3/8 Rule
Approximates sets of three panels with cubic polynomials.
Strict Constraint: Requires the number of intervals n n n to be a multiple of 3 . ∫ a b f ( x ) d x ≈ 3 h 8 [ f ( x 0 ) + 3 f ( x 1 ) + 3 f ( x 2 ) + 2 f ( x 3 ) + 3 f ( x 4 ) + 3 f ( x 5 ) + 2 f ( x 6 ) + ⋯ + f ( x n ) ] \int_a^b f(x)\,dx \approx \frac{3h}{8}\left[f(x_0) + 3f(x_1) + 3f(x_2) + 2f(x_3) + 3f(x_4) + 3f(x_5) + 2f(x_6) + \dots + f(x_n)\right] ∫ a b f ( x ) d x ≈ 8 3 h [ f ( x 0 ) + 3 f ( x 1 ) + 3 f ( x 2 ) + 2 f ( x 3 ) + 3 f ( x 4 ) + 3 f ( x 5 ) + 2 f ( x 6 ) + ⋯ + f ( x n ) ]
Polynomial Exactness: Degree ≤ 3 \le 3 ≤ 3 .
Global Truncation Error: E 3 / 8 = − ( b − a ) h 4 80 f ( 4 ) ( ξ ) = O ( h 4 ) E_{3/8} = -\frac{(b - a)h^4}{80} f^{(4)}(\xi) = O(h^4) E 3/8 = − 80 ( b − a ) h 4 f ( 4 ) ( ξ ) = O ( h 4 )
- Step-by-Step Worked Problems
- Worked Example: Thermal Dissipation in Curing Mass Concrete
Problem: A freshly cast mass concrete bridge pier foundation is placed at an initial core temperature of T 0 = 65 ∘ C T_0 = 65^\circ\text{C} T 0 = 6 5 ∘ C . The ambient surrounding air temperature is maintained at T m = 25 ∘ C T_m = 25^\circ\text{C} T m = 2 5 ∘ C . After 4 hours 4\text{ hours} 4 hours of cooling, the core temperature drops to 55 ∘ C 55^\circ\text{C} 5 5 ∘ C . Assuming Newton's Law of Cooling governs heat transfer:
Determine the cooling rate constant k k k (in hr − 1 \text{hr}^{-1} hr − 1 ).
What will be the concrete core temperature after a total elapsed time of 12 hours 12\text{ hours} 12 hours ?
Solution:
Apply Newton's Law of Cooling formula: T ( t ) = T m + ( T 0 − T m ) e − k t T(t) = T_m + (T_0 - T_m)e^{-kt} T ( t ) = T m + ( T 0 − T m ) e − k t T ( t ) = 25 + ( 65 − 25 ) e − k t = 25 + 40 e − k t T(t) = 25 + (65 - 25)e^{-kt} = 25 + 40e^{-kt} T ( t ) = 25 + ( 65 − 25 ) e − k t = 25 + 40 e − k t
Use the condition at t = 4 hours t = 4\text{ hours} t = 4 hours ( T = 55 ∘ C T = 55^\circ\text{C} T = 5 5 ∘ C ): 55 = 25 + 40 e − 4 k 55 = 25 + 40e^{-4k} 55 = 25 + 40 e − 4 k 30 = 40 e − 4 k ⟹ e − 4 k = 30 40 = 0.75 30 = 40e^{-4k} \implies e^{-4k} = \frac{30}{40} = 0.75 30 = 40 e − 4 k ⟹ e − 4 k = 40 30 = 0.75 − 4 k = ln ( 0.75 ) ≈ − 0.28768 -4k = \ln(0.75) \approx -0.28768 − 4 k = ln ( 0.75 ) ≈ − 0.28768 k = 0.28768 4 ≈ 0.07192 hr − 1 k = \frac{0.28768}{4} \approx 0.07192\text{ hr}^{-1} k = 4 0.28768 ≈ 0.07192 hr − 1
Compute the core temperature at t = 12 hours t = 12\text{ hours} t = 12 hours : T ( 12 ) = 25 + 40 e − 0.07192 ( 12 ) = 25 + 40 e − 0.86304 T(12) = 25 + 40e^{-0.07192(12)} = 25 + 40e^{-0.86304} T ( 12 ) = 25 + 40 e − 0.07192 ( 12 ) = 25 + 40 e − 0.86304 e − 0.86304 = ( e − 4 k ) 3 = ( 0.75 ) 3 = 0.421875 e^{-0.86304} = (e^{-4k})^3 = (0.75)^3 = 0.421875 e − 0.86304 = ( e − 4 k ) 3 = ( 0.75 ) 3 = 0.421875 T ( 12 ) = 25 + 40 ( 0.421875 ) = 25 + 16.875 = 41.88 ∘ C T(12) = 25 + 40(0.421875) = 25 + 16.875 = 41.88^\circ\text{C} T ( 12 ) = 25 + 40 ( 0.421875 ) = 25 + 16.875 = 41.8 8 ∘ C
The core temperature after 12 hours 12\text{ hours} 12 hours is 41.88 ∘ C 41.88^\circ\text{C} 41.8 8 ∘ C .
CELE Board Exam Traps & Strategic Checklists
Warning
Simpson's Rule Parity Trap: Simpson's 1/3 Rule strictly requires an even number of subintervals n n n (meaning an odd count of ordinates). If given an odd number of panels (e.g., n = 5 n = 5 n = 5 ), you cannot apply Simpson's 1/3 Rule directly across the entire domain; you must combine Simpson's 1/3 for the first 4 panels and the Trapezoidal Rule for the final panel (or use Simpson's 3/8 Rule for 3 panels + 1/3 Rule for 2 panels).
Linear ODE Standard Form: Before calculating the integrating factor I ( x ) = e ∫ P ( x ) d x I(x) = e^{\int P(x)dx} I ( x ) = e ∫ P ( x ) d x , you must ensure the leading coefficient of d y d x \frac{dy}{dx} d x d y is normalized to exactly 1 1 1 . For example, if given x d y d x + 2 y = 4 x 2 x\frac{dy}{dx} + 2y = 4x^2 x d x d y + 2 y = 4 x 2 , divide by x x x first to obtain d y d x + 2 x y = 4 x \frac{dy}{dx} + \frac{2}{x}y = 4x d x d y + x 2 y = 4 x so that P ( x ) = 2 / x P(x) = 2/x P ( x ) = 2/ x .
Newton-Raphson Sign: Remember the minus sign in x n + 1 = x n − f ( x n ) f ′ ( x n ) x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)} x n + 1 = x n − f ′ ( x n ) f ( x n ) . When f ( x n ) f(x_n) f ( x n ) is negative, subtracting a negative produces an addition!