function [dydt]=rhs(t,y)
dydt=0.5*y*(1-y/100);
