function yp=kani(t,y)
yp=[1, -.01*y(1);...
   .005*y(2), -0.5]*y;
end