
% Wie am Anfang was zu setzen ist, ist in dgl1.txt
% erklärt.
%______________________________________________________
%
% 3-tes Beispiel
%
% erzwungene gedämpfte/ungedämpfte Schwingung

 function xp = dgl3(t, x)
 
 delta = 0.0;
 %xp = [x(2), -delta*x(2) - 2*x(1)+1]';
 %xp = [x(2), -delta*x(2) - x(1)]';
 A = [0  1; -1 -delta];
 xp = A*x;


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