Question: float [8] error in plotting graph of ODE system of equations

restart; _local(gamma); _local(GAMMA); _local(Pi); _local(pi); _local(I); _local(D);
                               I
Warning, The imaginary unit, I, has been renamed _I
b := .45; mu[t] := 1.112; delta := .181; rho := 0.2e-1; beta[k] := .123; sigma := 0.9e-1; alpha := 0.312e-1; gamma := 0.14e-2; phi := .24; xi := .134; A[h] := .576; k[1] := 0.1e-2; beta[c] := 0.1e-1; mu[c] := 0.19e-2; eta := .557; z := .636; phi[c] := 0.57e-1;
                              0.45
                             1.112
                             0.181
                              0.02
                             0.123
                              0.09
                             0.0312
                             0.0014
                              0.24
                             0.134
                             0.576
                             0.001
                              0.01
                             0.0019
                             0.557
                             0.636
                             0.057
ODE1 := diff(B(T), T) = rho*b-mu[t]*B(T)-delta*B(T);
                  d                            
                 --- B(T) = 0.0090 - 1.293 B(T)
                  dT                           
ODE2 := diff(C(T), T) = delta*B(T)-mu[t]*C(T)+1-rho*b-beta[k]*C(T)*H(T)+sigma*C(T);
  d                                                           
 --- C(T) = 0.181 B(T) - 1.022 C(T) + 0.9910 - 0.123 C(T) H(T)
  dT                                                          
ODE3 := diff(D(T), T) = beta[k]*C(T)*H(T)-(mu[t]+alpha+gamma)*C(T)-phi*C(T);
       d                                                 
      --- D(T) = 0.123 C(T) H(T) - 1.1446 C(T) - phi C(T)
       dT                                                
ODE4 := diff(E(T), T) = alpha*D(T)-xi*E(T)-mu[t]*E(T);
               d                                 
              --- E(T) = 0.0312 D(T) - 1.246 E(T)
               dT                                
ODE5 := diff(F(T), T) = phi*D(T)+xi*E(T)-mu[t]*F(T)-sigma*C(T);
    d                                                       
   --- F(T) = phi D(T) + 0.134 E(T) - 1.112 F(T) - 0.09 C(T)
    dT                                                      
ODE6 := diff(G(T), T) = (1-k[1])*A[h]-beta[c]*G(T)*H(T)-mu[c]*H(T)-eta*J(T);
 d                                                             
--- G(T) = 0.575424 - 0.01 G(T) H(T) - 0.0019 H(T) - 0.557 J(T)
 dT                                                            
ODE7 := diff(H(T), T) = k[1]*A[h]+beta[c]*G(T)*H(T)-mu[c]*H(T)-z*H(T);
        d                                                
       --- H(T) = 0.000576 + 0.01 G(T) H(T) - 0.6379 H(T)
        dT                                               
ODE8 := diff(J(T), T) = z*H(T)-(phi[c]+eta)*J(T);
                d                                
               --- J(T) = 0.636 H(T) - 0.614 J(T)
                dT                               
ans := dsolve({ODE1, ODE2, ODE3, ODE4, ODE5, ODE6, ODE7, ODE8, B(0) = B0, C(0) = C0, D(0) = D0, E(0) = E0, F(0) = F0, G(0) = G0, H(0) = H0, J(0) = J0}, numeric, output = listprocedure);
Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)
    [T = proc(T)  ...  end;, B(T) = proc(T)  ...  end;,

      C(T) = proc(T)  ...  end;, D(T) = proc(T)  ...  end;,

      E(T) = proc(T)  ...  end;, F(T) = proc(T)  ...  end;,

      G(T) = proc(T)  ...  end;, H(T) = proc(T)  ...  end;,

      J(T) = proc(T)  ...  end;]
B0 := 100; C0 := 90; D0 := 45; E0 := 38; F0 := 100; G0 := 45; H0 := 50; J0 := 20;
                              100
                               90
                               45
                               38
                              100
                               45
                               50
                               20
ans := dsolve([ODE1, ODE2, ODE3, ODE4, ODE5, ODE6, ODE7, ODE8, B(0) = B0, C(0) = C0, D(0) = D0, E(0) = E0, F(0) = F0, G(0) = G0, H(0) = H0, J(0) = J0], numeric, output = listprocedure);
Error, (in f) unable to store 'HFloat(450.486)-HFloat(90.0)*phi' when datatype=float[8]