function xdot=cfunc_reg_perm(t,x)
xdot=zeros(3,1);
xdot(1)=x(2);
xdot(2)=x(3);
xdot(3)=-x(1)-2*x(2)-x(3)+sin(t);