function y=sprucebudworm_exact(t,a,K,P)

y = K*P./(P+(K-P)*exp(-t*a*K));