function dydt=HyChen(t,y)
global a b c d r
a=35; b=3; c=12; d=7; r=0.58;%input('r=');

dydt=[a*y(2,:)-a*y(1,:)+y(4,:)
           d*y(1,:)+c*y(2,:)-y(1,:)*y(3,:)
           y(1,:)*y(2,:)-b*y(3,:)
           y(2,:)*y(3,:)+r*y(4,:)
            ];