$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}{\| #1 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

# 3.2E: The Improved Euler Method and Related Methods (Exercises)

• • Contributed by William F. Trench
• Andrew G. Cowles Distinguished Professor Emeritus (Mathamatics) at Trinity University

$$\newcommand{\vecs}{\overset { \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

Most of the following numerical exercises involve initial value problems considered in the exercises in Section 3.1. You’ll find it instructive to compare the results that you obtain here with the corresponding results that you obtained in Section 3.1.

In Exercises [exer:3.2.1}– [exer:3.2.5} use the improved Euler method to find approximate values of the solution of the given initial value problem at the points $$x_i=x_0+ih$$, where $$x_0$$ is the point where the initial condition is imposed and $$i=1$$, $$2$$, $$3$$.

[exer:3.2.1} $$y'=2x^2+3y^2-2,\quad y(2)=1;\quad h=0.05$$

[exer:3.2.2} $$y'=y+\sqrt{x^2+y^2},\quad y(0)=1;\quad h=0.1$$

[exer:3.2.3} $$y'+3y=x^2-3xy+y^2,\quad y(0)=2;\quad h=0.05$$

[exer:3.2.4} $$y'= {1+x\over1-y^2},\quad y(2)=3;\quad h=0.1$$

[exer:3.2.5} $$y'+x^2y=\sin xy,\quad y(1)=\pi;\quad h=0.2$$

[exer:3.2.6} Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of the initial value problem $y'+3y=7e^{4x},\quad y(0)=2\nonumber$ at $$x=0$$, $$0.1$$, $$0.2$$, $$0.3$$, …, $$1.0$$. Compare these approximate values with the values of the exact solution $$y=e^{4x}+e^{-3x}$$, which can be obtained by the method of Section 2.1. Present your results in a table like Table [table:3.2.2}.

[exer:3.2.7} Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of the initial value problem $y'+{2\over x}y={3\over x^3}+1,\quad y(1)=1\nonumber$ at $$x=1.0$$, $$1.1$$, $$1.2$$, $$1.3$$, …, $$2.0$$. Compare these approximate values with the values of the exact solution $y={1\over3x^2}(9\ln x+x^3+2)\nonumber$ which can be obtained by the method of Section 2.1. Present your results in a table like Table [table:3.2.2}.

[exer:3.2.8} Use the improved Euler method with step sizes $$h=0.05$$, $$h=0.025$$, and $$h=0.0125$$ to find approximate values of the solution of the initial value problem $y'={y^2+xy-x^2\over x^2},\quad y(1)=2,\nonumber$ at $$x=1.0$$, $$1.05$$, $$1.10$$, $$1.15$$, …, $$1.5$$. Compare these approximate values with the values of the exact solution $y={x(1+x^2/3)\over1-x^2/3}\nonumber$ obtained in Example [example:2.4.3}. Present your results in a table like Table [table:3.2.2}.

[exer:3.2.9} In Example [example:3.2.2} it was shown that $y^5+y=x^2+x-4\nonumber$ is an implicit solution of the initial value problem $y'={2x+1\over5y^4+1},\quad y(2)=1. \eqno{\rm(A)}\nonumber$ Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of (A) at $$x=2.0$$, $$2.1$$, $$2.2$$, $$2.3$$, …, $$3.0$$. Present your results in tabular form. To check the error in these approximate values, construct another table of values of the residual $R(x,y)=y^5+y-x^2-x+4\nonumber$ for each value of $$(x,y)$$ appearing in the first table.

[exer:3.2.10} You can see from Example [example:2.5.1} that $x^4y^3+x^2y^5+2xy=4\nonumber$ is an implicit solution of the initial value problem $y'=-{4x^3y^3+2xy^5+2y\over3x^4y^2+5x^2y^4+2x},\quad y(1)=1. \eqno{\rm(A)}\nonumber$ Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of (A) at $$x=1.0$$, $$1.14$$, $$1.2$$, $$1.3$$, …, $$2.0$$. Present your results in tabular form. To check the error in these approximate values, construct another table of values of the residual $R(x,y)=x^4y^3+x^2y^5+2xy-4\nonumber$ for each value of $$(x,y)$$ appearing in the first table.

[exer:3.2.11} Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of the initial value problem $(3y^2+4y)y'+2x+\cos x=0, \quad y(0)=1 \mbox{\ (Exercise~2.2.~\hspace*{-3pt}\ref{exer:2.2.13})}\nonumber$ at $$x=0$$, $$0.1$$, $$0.2$$, $$0.3$$, …, $$1.0$$.

[exer:3.2.12} Use the improved Euler method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of the initial value problem $y'+{(y+1)(y-1)(y-2)\over x+1}=0, \quad y(1)=0 \mbox{\ (Exercise~2.2.~\hspace*{-3pt}\ref{exer:2.2.14})}\nonumber$ at $$x=1.0$$, $$1.1$$, $$1.2$$, $$1.3$$, …, $$2.0$$.

[exer:3.2.13} Use the improved Euler method and the improved Euler semilinear method with step sizes $$h=0.1$$, $$h=0.05$$, and $$h=0.025$$ to find approximate values of the solution of the initial value problem $y'+3y=e^{-3x}(1-2x),\quad y(0)=2,\nonumber$ at $$x=0$$, $$0.1$$, $$0.2$$, $$0.3$$, …, $$1.0$$. Compare these approximate values with the values of the exact solution $$y=e^{-3x}(2+x-x^2)$$, which can be obtained by the method of Section 2.1. Do you notice anything special about the results? Explain.

[exer:3.2.14} $$y'-2y= {1\over1+x^2},\quad y(2)=2$$; $$h=0.1,0.05,0.025$$ on $$[2,3}$$

[exer:3.2.15} $$y'+2xy=x^2,\quad y(0)=3$$; $$h=0.2,0.1,0.05$$ on $$[0,2}$$ (Exercise 2.1. [exer:2.1.38})

[exer:3.2.16} $${y'+{1\over x}y={\sin x\over x^2},\quad y(1)=2}$$, $$h=0.2,0.1,0.05$$ on $$[1,3}$$ (Exercise 2.1. [exer:2.1.39})

[exer:3.2.17} $${y'+y={e^{-x}\tan x\over x},\quad y(1)=0}$$; $$h=0.05,0.025,0.0125$$ on $$[1,1.5}$$ (Exercise 2.1. [exer:2.1.40}),

[exer:3.2.18} $${y'+{2x\over 1+x^2}y={e^x\over (1+x^2)^2}, \quad y(0)=1}$$; $$h=0.2,0.1,0.05$$ on $$[0,2}$$ (Exercise 2.1. [exer:2.1.41})

[exer:3.2.19} $$xy'+(x+1)y=e^{x^2},\quad y(1)=2$$; $$h=0.05,0.025,0.0125$$ on $$[1,1.5}$$ (Exercise 2.1. [exer:2.1.42})

[exer:3.2.20} $$y'+3y=xy^2(y+1),\quad y(0)=1$$; $$h=0.1,0.05,0.025$$ on $$[0,1}$$

[exer:3.2.21} $${y'-4y={x\over y^2(y+1)},\quad y(0)=1}$$; $$h=0.1,0.05,0.025$$ on $$[0,1}$$

[exer:3.2.22} $${y'+2y={x^2\over1+y^2},\quad y(2)=1}$$; $$h=0.1,0.05,0.025$$ on $$[2,3}$$

[exer:3.2.23} Do Exercise [exer:3.2.7} with “improved Euler method” replaced by “midpoint method.”

[exer:3.2.24} Do Exercise [exer:3.2.7} with “improved Euler method” replaced by “Heun’s method.”

[exer:3.2.25} Do Exercise [exer:3.2.8} with “improved Euler method” replaced by “midpoint method.”

[exer:3.2.26} Do Exercise [exer:3.2.8} with “improved Euler method” replaced by “Heun’s method.”

[exer:3.2.27} Do Exercise [exer:3.2.11} with “improved Euler method” replaced by “midpoint method.”

[exer:3.2.28} Do Exercise [exer:3.2.11} with “improved Euler method” replaced by “Heun’s method.”

[exer:3.2.29} Do Exercise [exer:3.2.12} with “improved Euler method” replaced by “midpoint method.”

[exer:3.2.30} Do Exercise [exer:3.2.12} with “improved Euler method” replaced by “Heun’s method.”

[exer:3.2.31} Show that if $$f$$, $$f_x$$, $$f_y$$, $$f_{xx}$$, $$f_{yy}$$, and $$f_{xy}$$ are continuous and bounded for all $$(x,y)$$ and $$y$$ is the solution of the initial value problem $y'=f(x,y),\quad y(x_0)=y_0,\nonumber$ then $$y''$$ and $$y'''$$ are bounded.

[exer:3.2.32} Numerical Quadrature (see Exercise 3.1. [exer:3.1.23}).

1. Derive the quadrature formula $\int_a^bf(x)\,dx \approx 0.5h(f(a)+f(b))+h\sum_{i=1}^{n-1}f(a+ih) \tag{A}\nonumber$ (where $$h=(b-a)/n)$$ by applying the improved Euler method to the initial value problem $y'=f(x),\quad y(a)=0.\nonumber$ The quadrature formula (A) is called the trapezoid rule. Draw a figure that justifies this terminology.
2. For several choices of $$a$$, $$b$$, $$A$$, and $$B$$, apply (A) to $$f(x)=A+Bx$$, with $$n = 10,20,40,80,160,320$$. Compare your results with the exact answers and explain what you find.
3. For several choices of $$a$$, $$b$$, $$A$$, $$B$$, and $$C$$, apply (A) to $$f(x)=A+Bx+Cx^2$$, with $$n=10$$, $$20$$, $$40$$, $$80$$, $$160$$, $$320$$. Compare your results with the exact answers and explain what you find.