function dy=raindropf(r,y);

%epsilon=7.5e-4;
epsilon=.1;
delta=.0124;
g=9.8;

dy=-6*g/epsilon+6*delta/epsilon./r.^1.5.*y.^.75 +6./r.*y;
