% driver program for the pendulum


clear all;


% index of the DAE
index = 1;

% Init
[tspan,y0,options] = mpend_Idx123([],[],'init',index);

Tol = 1e-6;
options = odeset(options,'RelTol',Tol,'AbsTol',Tol);

options.Debug = 1;
options.Predict = 'on';
options.Complex = 'off';

tspan = [0 100];
y0    = [1 0 0 0 0]';

% 
[tout,yout,stats] = radau5('mpend_Idx123',tspan,y0,options,index);


% Ausgabe ...

