function res=circ_dae(x,dx,t)
    res=[dx(1)-x(3);dx(2)-x(1)-x(3);2*x(3)+x(2)+x(1)];
endfunction
