function f = hund(t, y, v, c)
% DGL zum schwimmenden Hund
% v ... geschwindigkeit des hundes (in Richtung Ursrung)
% c ... geschwindigkeit des flusses (in x-Richtung)

f=zeros(2,1);

help = norm(y,2);

f(1) = - v * y(1)/help + c;
f(2) = - v * y(2)/help;