Runge-Kutta method, dP/dt = (0.25*cos(t))P(1-P)-0.228, P(0) = 10, from 1 to 100