diff --git a/scripts/Models/M1.txt b/scripts/Models/M1.txt
index 84cf631d421b3478ed67aa963f97b68ec99afabc..5236f8094be2b1cac7ea59634a028698a5f013d7 100644
--- a/scripts/Models/M1.txt
+++ b/scripts/Models/M1.txt
@@ -31,12 +31,11 @@ d/dt(GLUT4) = -v7f+v7b %%24
d/dt(S6K) = v9b1-v9f1 %%26
d/dt(S6Kp) = v9f1-v9b1 %%27
d/dt(S6) = v9b2-v9f2 %%28
-d/dt(S6p) = v9f2-v9b2 %%29
-
+d/dt(S6p) = v9f2-v9b2 %%29
d/dt(G_p) = EGP+Ra-E-U_ii-k_1*G_p+k_2*G_t
-d/dt(G_t) = (-U_id)+k_1*G_p-k_2*G_t
-d/dt(I_p) = (-m_2*I_p)-m_4*I_p+m_1*I_l
+d/dt(G_t) = (-U_id)+k_1*G_p-k_2*G_t
d/dt(I_l) = (-m_1*I_l)-m_3*I_l+m_2*I_p+S
+d/dt(I_p) = (-m_2*I_p)-m_4*I_p+m_1*I_l
d/dt(Q_sto1) = -k_gri*Q_sto1
d/dt(Q_sto2) = (-k_empt*Q_sto2)+k_gri*Q_sto1
d/dt(Q_gut) = (-k_abs*Q_gut)+k_empt*Q_sto2
@@ -74,8 +73,7 @@ GLUT4(0) = 100
S6K(0) = 100
S6Kp(0) = 0
S6(0) = 100
-S6p(0) = 0
-
+S6p(0) = 0
G_p(0) = 178
G_t(0) = 130
I_l(0) = 4.5
@@ -137,7 +135,7 @@ k9b2 = 30.9967 %%42
km9 = 5872.68 %%43
n9 = 0.985466 %%44
kbf = 1e+06 %%45
-nC = 2.1e-06 %%46
+scaleModel = 2.1e-06 %%46
V_G = 1.88 %%47
k_1 = 0.065 %%48
@@ -166,38 +164,33 @@ k_p3 = 0.009 %%70
k_p4 = 0.0618 %%71
k_i = 0.0079 %%72
-V_m0 = 2.5 %%73
-V_mX = 0.047 %%74
-K_m0 = 225.59 %%75
-V_f0 = 2.5 %%76
-V_fX = 0.047 %%77
-K_f0 = 225.59 %%78
-p_2U = 0.0331 %%79
-part = 0.2 %%80
-K = 2.3 %%81
-alpha = 0.05 %%82
-beta = 0.11 %%83
-gamma = 0.5 %%84
-k_e1 = 0.0005 %%85
-k_e2 = 339 %%86
-D = 78000 %%87
+V_m0 = 2.5 %%74
+V_mX = 0.047 %%75
+K_m0 = 225.59 %%76
+V_f0 = 2.5 %%77
+V_fX = 0.047 %%78
+K_f0 = 225.59 %%79
+p_2U = 0.0331 %%80
+part = 0.2 %%81
+K = 2.3 %%82
+alpha = 0.05 %%83
+beta = 0.11 %%84
+gamma = 0.5 %%85
+k_e1 = 0.0005 %%86
+k_e2 = 339 %%87
+D = 78000 %%88
-K_l0 = 225.59 %%88
-V_l0 = 2.5 %%89
-V_lX = 0.047 %%90
-pf = 34
-
-U_ii = 0.8
-part_m = 0.5
-part_f = 0.5
-part_l = 1.7
+K_l0 = 225.59
+V_l0 = 2.5
+V_lX = 0.047
********** MODEL VARIABLES
aa = 5/2/(1-b)/D
cc = 5/2/d/D
EGP = k_p1-k_p2*G_p-k_p3*I_d-k_p4*I_po
V_mmax = (1-part)*(V_m0+V_mX*INS)
-V_fmax = part*(V_f0+V_fX*INS)
+V_fmax = part*(V_f0+V_fX*INS)
+V_lmax = (1-part)*(V_l0+V_lX*G_t)
Ee = k_e1*(G_p - k_e2)
S = gamma*I_po
I = I_p/V_I
@@ -207,11 +200,13 @@ m_3 = HE*m_1/(1-HE)
Q_sto = Q_sto1+Q_sto2
Ra = f*k_abs*Q_gut/BW
k_empt = k_min+(k_max-k_min)/2*(tanh(aa*(Q_sto-b*D))-tanh(cc*(Q_sto-d*D))+2)
-U_idm = part_m*V_mmax*G_t/(K_m0+G_t)
-U_idf = part_f*nC*(k8*GLUT4m*G_t/pf + glut1*G_t/pf + kbf*(INS+5)*1e-3)
+U_idm = 0.5*V_mmax*G_t/(K_m0+G_t)
+U_idf = 0.5*scaleModel*(k8*GLUT4m*G_t/34 + glut1*G_t/34 + kbf*(INS+5)*1e-3)
+U = U_ii+U_id
S_po = Y+K*(EGP+Ra-E-U_ii-k_1*G_p+k_2*G_t)/V_G+S_b
-U_idl = part_l*U_idm
-U_id = U_idm+U_idf+U_idl
+U_liver = 45*U_idm/27
+U_ii = 23*U_liver/45
+U_id = U_idm+U_idf+U_liver
********** MODEL REACTIONS
v1a = IR*k1a*(INS+5)*1e-3
@@ -259,4 +254,4 @@ glucoseuptake = k8*GLUT4m*(G_t*5/170) + glut1*(G_t*5/170) + kbf*((INS+5)*1e-3)
********** MODEL EVENTS
event1 = gt(G_p,339),E,Ee
event2 = lt(G_p,339),E,0
-********** MODEL MATLAB FUNCTIONS
\ No newline at end of file
+********** MODEL MATLAB FUNCTIONS