function yp = rowr(t,y)

yp = zeros(3,1);
yp(1) = -3*y(2)*y(3)-(y(1)^2+y(2)^2+y(3)^2)*y(1);
yp(2) = -3*y(1)*y(3)-(y(1)^2+y(2)^2+y(3)^2)*y(2);
yp(3) = -3*y(1)*y(2)-(y(1)^2+y(2)^2+y(3)^2)*y(3);
