Initial Value Problems - Numerical Methods in Engineering with MATLAB - page 291

The initialconditions are

0
0

y (0) =

We solved the problemwith the function bulStoer using the increment H =

0 . 5 s. The following program utilizes the plotting facilities of MATLAB:

% Example 7.11 (Bulirsch-Stoer integration)
[xSol,ySol] = bulStoer(@fex7

_

11,0,[0 0],10,0.5);

plot(xSol,ySol(:,2),'k:o')
grid on
xlabel('Time (s)')
ylabel('Current (A)')

Recall that in each interval H (the spacing of open circles) the integrationwas per-
formedbythe modified midpoint method and refinedbyRichardson'sextrapolation.

PROBLEM SET 7.2

1. Derive the analytical solution of the problem

y + y − y (0) 380 y = 0 y (0) = 1 =− 20

Wouldyou expect difficulties in solving this problem numerically?