% Example of Euler Integration
% AJS

clear all
deltat=0.01;

t=[0:deltat:10];

% a scalar example
x(:,1)=0;

% % an example with three states
% x(:,1)=[0 0 0]';

for k=1:length(t)-1
    K1 = diffeq(t(k),x(:,k));
    x(:,k+1)=x(:,k)+ K1*deltat;
end