Question: float in loop

Dear guys! I've written some expressions and then the final equations which I want to solve as below:

> Omega0 := 0.26: H0 := 0.8: q0 := -0.61: j0 := 1: s0 := -0.18: l0 := 2.61:

> f0 := 6*H0^2*(Omega0-2):      f10 := 1:    f20 := (6*H0^2)^(-1)*(-3*Omega0/(4*(1+q0))+1/2):     f30 := (6*H0^2)^(-2)*(-3*Omega0*(3*q0^2+6*q0+j0+2)/(8*(1+q0)^3)+3/4):     f40 := (6*H0^2)^(-3)*(-3*Omega0*(s0*(1+q0)+j0*(3*j0+6*q0^2+17*q0+5)+3*q0*(5*q0^3+20*q0^2+29*q0+16)+9)/(16*(1+q0)^5)+15/8):     f50 := (6*H0^2)^(-4)*(-3*Omega0*(l0*(1+q0)^2+s0*(10*q0^3+43*q0^2+46*q0+13)+10*j0*s0*(1+q0)+5*j0^2*(3*j0+7+6*q0^2+22*q0)+j0*(45*q0^4+225*q0^3+412*q0^2+219*q0+32)+3*q0*(35*q0^5+210*q0^4+518*q0^3+666*q0^2+448*q0+150)+60)/(32*(1+q0)^7)+105/16):

> T := -6*H^2:     T0 := -6*H0^2:     f := f0+f10*(T-T0)+f20*(T-T0)^2/2+f30*(T-T0)^3/6+f40*(T-T0)^4/24+f50*(T-T0)^5/120:     f1 := f10+f20*(T-T0)+f30*(T-T0)^2/2+f40*(T-T0)^3/6+f50*(T-T0)^4/24:     Omega := H0^2*Omega0*(1+z)^3/H^2:    

> eq := unapply(H^2=(-1/(12*f1))*(T*Omega+f),z):    

> Y := z->if not type(z,numeric) then 'procname(z)' else fsolve(eq(z), H=H0) end  if:

> l := dsolve({D(L)(z) = L(z)/(1+z)+(1+z)/Y(z),L(0)=0}, type=numeric, range=0..5, known=Y):

Until here anything is right, but when I want to inter all of these into a loop for the "Omega0, H0, q0, j0, s0, l0" I get an error like this: " Error, (in f) unable to store 'HFloat(1.0)/fsolve ... ". I don't really understand what the problem is. Can anyone explain it for me? Thanks a lot.