Resolution of first order differential equations. Method of Euler

 

Let's consider a differential equation defined by dy/dx = f(x,y) with the initial value of y(0) = y₀. Let's say there isn't any method that can't solve this equation. Not knowing how to find the solution, let's then approximate it.

Let's consider two approximate points solutions of the differential equation and draw a line through these points ( see figure above). Let's calculate the slope of that line. The slope is the tangent of the angle made by by the line and the parallel to the x-axis. Let's call this angle θ. We have:

tanθ. = yⁱᵢ₊₁ - yᵢ /xᵢ₊₁-xᵢ . It is rise/run, the formula to calculate the slope.

tanθ is also the slope of the point (xᵢ, yᵢ). Then we can write f(xᵢ, yᵢ) = yⁱᵢ₊₁ - yᵢ /xᵢ₊₁-xᵢ . We can derive the formula of Euler from this but it will not be fair to do that. What I want to show here is is a reminder of the formula of the slope which is rise/run. 

Let's consider the following figure:

We still have the same differential equation dy/dx = f(x,y) with the initial condition y(0) = y₀

Let's draw a tangent line through the point (0, y₀). Knowing x₁, let y₁ be the approximate solution of the differential equation.. Its value can be calculated using the slope of the tangent line. This slope is the slope at the point (x₀, y₀) which is f(x₀, y₀). We have two points here (x₀, y₀) and (x₁, y₁). Using the formula of the slope rise/run we have:

f(x₀, y₀) = y₁ - y₀/x₁ - x₀

(fx₀, y₀)(/x₁ - x₀) = y₁ - y₀

y₁ = y₀ + (fx₀, y₀)(/x₁ - x₀). This is the approximate value of y₁ 

Now let's draw a line through the point (x₁, y₁) and let's y₂ be the second approximation of the solution knowing x₂

The slope of the line that passes through the points (x₁, y₁) and (x₂, y₂) is :

f(x₁, y₁) = y₂ -y₁/x₂ - x₁

   f(x₁, y₁)(x₂ - x₁) = y₂ -y₁

y₂ = y₁ +  f(x₁, y₁)(x₂ - x₁). This is the approximate value of y₂.

A pattern appears here based on the values of  y₁ and y₂. If we were to calculate a third approximation of the solution which is y₃, it would be: y₃ = y₂ + f(x₂, y₂)(x₃ - x₂)'

Any approximation yₙ can be calculated by yₙ = yₙ₋₁ + f(xₙ₋₁, yₙ₋₁) (xₙ - xₙ₋₁)

Let's write  xₙ - xₙ₋₁ = h, the above equation can be written as : yₙ = yₙ₋₁ + f(xₙ₋₁, yₙ₋₁)h where h is called the step size. This is the Euler formula that allows to calculate the approximate values of a function.

Example

Let's consider a differential equation with a given solution y' + 2y = 2 - e^-4t  y(0) = 1 Use the Euler method with a step size of h = 0.1 to find approximate values of the solutions at t = 0.1 t = 0.2  t = 0.3

Solution 

Let's write the equation on its explicit form: 

y' = 2 - e^-4t - 2y

Let's apply the Euler formula:   yₙ = yₙ₋₁ + f(xₙ₋₁, yₙ₋₁)h. Since the independent variable is t, we can write the formula as yₙ = yₙ₋₁ + f(tₙ₋₁, yₙ₋₁)h

We have h = 0.1 we have to find y for t = 0.1. Since this is the first value of y, we have to call it y₁Let's substitute n by 1 in the formula:

y₁ = y₀ + f(t₀, y₀).0.1

We already have y₀. Let's calculate f(t₀, y₀) or f(0,1). In order to do that, we substitute t by 0 and y by 1 in the explicit form of the differential equation

f(0,1) = 2 - e^0 - 2 = 2 - 1 - 2 = -1

Let's substitute y₀ and f(0,1):

y₁ = 1 + (-1)(0.1) = 1 - 0.1 = 0.90

Let's find y for the value of t = 0.2: In this case we have to substitute n by 2 in the Euler formula in order to find y₂. 

y₂ = y₁ + f(t₁, y₁).0.1

Let's find f(t₁, y₁):

f(0.1, 0.9) = 2 - e^-0.4 - 2(0.9) = -0.47

y₂ = 0.9 + (-0.47)(0.1) = 0.85.

Let's find another approximate value of y for t = 0.3

y₃ = y₂ + f(t₂,y₂).0.1

Let's calculate f(t₂,y₂):

f(0.2, 0.85) = 2 - e^-0.8 -2.0.85 = 2 - 1/e^0.8 - 1.70 = 2 - 0.449 - 1.70 = -0.149

y₃ = 0.85 + (-0.149)(0.1)

y₃ = 0.8351

Where the calculations of the approximate values contain decimals it's better to find the approximations with more than 2 or 3 decimals.

Interested in reviewing Calculus I check this site: Center for Integral Development 

Interval of validity of a non-linear differential equation

 In the following problem we are going to determine the interval of validity of a non-linear differential equation.

Problem. Determine the interval of validity of the initial value problem and give its dependence on the value of y₀.

y' = y²  y(0) = y₀

Let's see if f(t, y) = y² and δf/δy are continuous. It's obvious that f(t, y) =  y² is continuous. 

Let's calculate δf/δy:

δf/δy = 2y. It's continuous. Then there is a unique solution to the IVP in an interval that contains t₀ = 0

We have to determine this solution in order to find the interval of validity. 

Let's solve the differential equation

dy/dt = :y² 

dy = y² dt

dy/ y² = dt

y^-2 dy = dt

∫y^-2 dy = ∫dt

y^-2+1/-2 +1 = t + c

y^-1/-1 = t

-1/y = t + c

y = -1/t + c

Let's determine the value of c. Let's substitute t by 0 and y by  y₀ as it's given in the initial condition

y₀ = -1/c c = -1/y₀

Let's substitute c in y we have:

y = -y₀ /y₀ t -1

The interval of validity of the IVP is the interval where y is defined. 

Let's determine this interval :of validity.

The function y isn't defined if the denominator is equal to 0.

If y₀ t -1 = 0, t = 1/y₀ . Therefore y isn't defined for that value of t.

Y is defined in the following intervals: -∞<t<1/y₀  and 1//y₀<t<+∞

Let's determine the interval of validity for different values of y₀:

If y₀ = 0, then the interval of validity is -∞<t<+∞

If y₀<0  1/y₀<0 t = 0 doesn't belong to  -∞<t<1/y₀ but belongs to 1//y₀<t<+∞, which is an interval of validity

If y₀>0 then 1/y₀>0  -∞<t<1/y₀

The problem that we just solved demonstrates the difference between the interval of validity of linear differential equations and that of non-linear differential equations. The interval of validity of linear differential equations doesn't depend the value of y₀. The interval of validity of non-linear differential equations depends on the value of  y₀.

Homogeneous Differential Equations

 An homogeneous differential equation is one that can be written as y' = F(y/x). The equation is not written in this form most of the times. You have to transform it in order to write it under this form. This type of equation is solved by substitution. That's another type of equation that's solved by substitution. Remember that in the case of a Bernoulli differential equation we use the substitution v = y¹⁻ⁿ. In the case of an homogenous differential equation we are going to see what type of substitution we can use.

Let's bring this homogeneous differential equation to us in order to solve it. Here it comes as y' = F(y/x). It seems natural that we can substitute y/x by something. Let's call it v. We can write v =  y/x. Then we have: y = vx.

Let's derive both sides of the equation:

y' = v'x + v

Remember also that the differential equation is y' = F(y/x) or y' = F(v)

 Equaling both expressions of y' we have:

v'x + v = F(v)

We have a differential equation with a new dependent variable v. This equation can be solved by the method of separation of variables. Let's continue;

v'x = F(v) - v

(dv/dx)x = F(v) - v

x/dx = F(v) - v/dv

dx/x  = dv/F(v) - v

Now that the variables are separated we can integrate both sides of the equation in order to find v. Once we find v we can use the substitution to find y.

Example. Solve the differential equation xyy' + 4x² + y² = 0

Let's bring the expression y/x in the equation by dividing by x. Let's divide both sides by x²:

 xyy'/x² + 4x²/x² +  y²/x²  = 0

yy'/x + 4 + (y/x)² = 0

y/xy' + 4 + (y/x)² = 0

Let's substitute y/x by v:

vy' + 4 + v² = 0  (1)

Let's find y' in order to have an equation with the new variable v:

y/x = v 

y = vx

y' = v'x + v

Let's substitute y' in equation (1)

v(v'x + v) + 4 + v² = 0

vv'x + v² + 4 +  v² = 0

vv'x + 2v² + 4 = 0

vv'x  = -2v² - 4 

Let's divide both sides by v in order to separate the variables:

v'x  = -2v² - 4/v

(dv/dx)x = -2v² - 4/v

x/dx = -2v² - 4/vdv

dx/x = -vdv/ 2v² + 4

Now that the variable are separated let's integrate both sides in order to find v;

∫dx/x = ∫(-v/ 2v² + 4)dv = -∫(v/ 2v² + 4)dv (2)

Let's solve the integral in the second side by doing the substitution u =  2v² +4

du = 4vdv dv = du/4v

The integral at the second side becomes:

∫(v/ 2v² +4)dv = ∫(v/u)(du/4v) = ∫du/4u = 1/4∫du/u = 1/4lnu = 1/4ln (2v² +4) + k

Let's substitute the value of the integral we just calculated in (2)

∫dx/x = -1/4ln (2v² + 4) + k

lnx + K' = -1/4ln (2v² +4) + k

1/4ln (2v² + 4) = -lnx - k' + k

Let's call c the expression-k'+k 

1/4ln (2v² + 4) = -lnx + c

1/4ln (2v² + 4) = -lnx + c

ln(2v² + 4)^1/4  =  -lnx + c

(2v² + 4)^1/4 = e^(-lnx + c)

                     = (e^-lnx)(e^c)

                     = (e^lnx^-1)(d) We call e^c d

  e^lnx^-1 being equal to x^-1, we have:

(2v² + 4)^1/4 = (x^-1)(d) 

Let's raise both sides to the power of 4:

2v² + 4 = x⁻⁴d⁴

Let's call d⁴ g:

2v² + 4 = x⁻⁴g

            = g/x⁴

v²  = g-4x⁴/2x⁴

v = ∓(⎷2/2) (⺁g-4x⁴)/x²

By substituting v in the expression v = y/x, we obtain:

y = ∓⎷2/2⺁g-4x⁴/x

Practice. Solve the differential equation for the initial condition y(2) = -7 and x>0