I have adjusted limits. it shows the below Warning.

```
┌ Warning: First function call produced NaNs. Exiting.
└ @ OrdinaryDiffEq C:\Users\julia\.julia\packages\OrdinaryDiffEq\HO8vN\src\initdt.jl:137
┌ Warning: Automatic dt set the starting dt as NaN, causing instability.
└ @ OrdinaryDiffEq C:\Users\julia\.julia\packages\OrdinaryDiffEq\HO8vN\src\solve.jl:481
┌ Warning: NaN dt detected. Likely a NaN value in the state, parameters, or derivative value caused this outcome.
└ @ DiffEqBase C:\Users\julia\.julia\packages\DiffEqBase\3iigH\src\integrator_interface.jl:323
```

Adjusted Upper bounds

```
phenbm = @model begin
@param begin
tvka ∈ RealDomain(lower=0,init =0.225, upper = 2)
tvk23 ∈ RealDomain(lower=0,init = 5.02, upper = 10)
tvQ ∈ RealDomain(lower=0,init = 60 , upper = 150 )
tvvc ∈ RealDomain(lower=0,init = 223, upper = 500)
tvvp ∈ RealDomain(lower=0,init = 2.8, upper = 10)
tvCL ∈ RealDomain(lower=0,init = 7 ,upper =40)
tvbio1 ∈ RealDomain(lower=0,init =0.63 , upper =1)
tvbio2 ∈ RealDomain(lower=0,init = 0.66, upper =1)
Ω ∈ PDiagDomain(init=[3.37,0.09,0.04,0.72,0.04,0.23,0.74,0.04]) # CL = 0.23,Vc = 0.72,KA = 3.37,vp = 0.04,k23 = 0.09
σ ∈ RealDomain(lower=0,init = 0.04,upper = 0.1)
end
@random begin
η ~ MvNormal(Ω)
end
@pre begin
ka = tvka * exp(η[1])
k23 = tvk23 * exp(η[2])
Q = tvQ * exp(η[3])
vc = tvvc * exp(η[4])
vp = tvvp * exp(η[5])
CL = tvCL * exp(η[6])
bioav = (Depot1 = tvbio1* exp(η[7]), Depot2 = tvbio2* exp(η[8]))
end
@dynamics begin
Depot1' = -ka * Depot1
Depot2' = -k23 * Depot2
Central' = ka * Depot1 + k23 * Depot2 - (Q/vc)*(Central) + (Q/vp)*Peripheral - (CL/vc)*Central
Peripheral' = (Q/vc) * Central - (Q/vp) * Peripheral
end
@derived begin
cp = @. Central/vc
DV ~ @. Normal(cp,abs(cp)*σ)
end
end
```