diff --git a/N803_model_2.m b/N803_model_2.m
index 6410fd4..1632fee 100644
--- a/N803_model_2.m
+++ b/N803_model_2.m
@@ -1,7 +1,7 @@
 %% N803_model_2.m - solves ODE model for N-803 treatment of SIV
 %
 %     /--------------------------------------------------------------\
-%     | Date:         10/31/2021                                     |
+%     | Date:         11/05/2021                                     |
 %     | Author:       Jonathan Cody                                  |
 %     | Affiliation:  Purdue University                              |
 %     |               Weldon School of Biomedical Engineering        |
@@ -10,6 +10,8 @@
 %
 % Version 2D+ (also allowing regulation calibration)
 %
+% 12/03/2021 - added some extra terms to the full output matrix
+%
 % Nomenclature: V  = SIV virions [#/µL]
 %               E0 = resting SIV-specific CD8+ T cells [#/µL]
 %               Ea = active  SIV-specific CD8+ T cells [#/µL]
@@ -252,6 +254,18 @@
 % Ba killing proportion
 Y_OUT(:,37) = gB0*Ba./(1+gB2*R) ./ Kill ;
 
+% Average E proliferation rate
+Y_OUT(:,38) = ( p0*F(:,15).*E0 + pE*sum( Y(:,03:9) ,2) ) ./ (E0+Ea) ;
+
+% Average B proliferation rate
+Y_OUT(:,39) = ( p0*F(:,15).*B0 + pB*Y(:,12) ) ./ (B0+Ba) ;
+
+% Average E killing rate
+Y_OUT(:,40) =  gE0*Ea./(1+gE2*R) ./ (E0+Ea) ;
+
+% Average B killing rate
+Y_OUT(:,41) =  gB0*Ba./(1+gB2*R) ./ (B0+Ba) ;
+
 %% ========================================================================
 %	NESTED FUNCTIONS
 %  ========================================================================
@@ -305,7 +319,7 @@
     dY(4:9)= 2*pE*Ea(1:6) - dA*Ea(2:7) - pE*Ea(2:7) ;% dE2-7/dt
     dY(10) = 2*pE*Ea(7)   - dA*Ea(8)   - mE*Ea(8)   ;% dE8/dt
     dY(11) = p*B0 - aB*B0 - d*B0       + mB*Ba(2)   ;% dB0/dt
-    dY(12) =      + aB*B0 - dA*Ba(2)   - pB*Ba(1)   ;% dB1/dt               !!!!! dA*Ba(1) !!!!!
+    dY(12) =      + aB*B0 - dA*Ba(1)   - pB*Ba(1)   ;% dB1/dt
     dY(13) = 2*pB*Ba(1)   - dA*Ba(2)   - mB*Ba(2)   ;% dB2/dt
     dY(14) =         - ka*X ;% dX/dt
     dY(15) = ka*X/vd - ke*C ;% dC/dt