
% Wie am Anfang was zu setzen ist, ist in dgl1.txt
% erklärt.
%______________________________________________________
%
% nicht-lineares Beispiel
%
% gestörte harmonische Schwingung

 function xp = dgl3e(t, x)
 
 eps = .1;
 xp = [x(2), eps*x(1)^3 - x(1)]';

%  alpha = .01;
%  beta  = .02;  
%  xp = diag([1 - alpha*x(2),-1 + beta*x(1)])*x;
%
