D(rx)=rx(PD/100)+DO where, parameter PD=round(100P*C(P)) and parameter DO=round[HCT0*P*(1???C(P))???(W/2)???L]; and, where P represents a value for px corresponding to the desired cardiac phase, and round is a function which rounds off its argument to a given place.