MATLAB Mathematics - SERC - Index of

5 Differential Equations
dydt(i+1,:) = 3*y(i,:) - y(i+1,:).*y(i,:).^2 + ...
c*(y(i-1,:)-2*y(i+1,:)+3);
end % End nested function f
end % End function brussode
% --------------------------------------------------------------
function S = jpattern(N)
B = ones(2*N,5);
B(2:2:2*N,2) = zeros(N,1);
B(1:2:2*N-1,4) = zeros(N,1);
S = spdiags(B,-2:2,2*N,2*N);
end;
The Brusselator for N = 20
3
2.5
2
solution u
1.5
1
0.5
0
10
8
6
time
4
2
0
0
0.2
0.4
space
0.6
0.8
1
Example: Simple Event Location
ballode models the motion of a bouncing ball. This example illustrates the
event location capabilities of the ODE solvers.
5-28