Fmd1[x_, rho_] := Module[{n=0}, T=Floor[x]; t=x-T; v=RecMd1[T,rho]; If[ (t == 0), Plus @@ v, E^(rho t) Sum[(-rho t)^n/n! (Plus @@ Take[v,T-n+1]), {n,0,T}] ] ]