//Created by libAntimony v2.4 // Warnings from automatic translation: // In module 'Guyton_Model__muscle_O2_delivery', the variables delivery_of_O2_to_M_tissues.PVO and PVO were unable to be set as equivalent: Loop detected: 'delivery_of_O2_to_M_tissues.PVO' may not be set to be equal to 'PVO' because PVO's definition already includes delivery_of_O2_to_M_tissues.PVO either directly or by proxy. // In module 'Guyton_Model__non_muscle_O2_delivery', the variables delivery_of_O2_to_NM_tissues.POV and POV were unable to be set as equivalent: Loop detected: 'delivery_of_O2_to_NM_tissues.POV' may not be set to be equal to 'POV' because POV's definition already includes delivery_of_O2_to_NM_tissues.POV either directly or by proxy. // In module 'Guyton_Model__kidney', the variables glomerular_filtration_rate.GLPC and GLPC were unable to be set as equivalent: Loop detected: 'glomerular_filtration_rate.GLPC' may not be set to be equal to 'GLPC' because GLPC's definition already includes glomerular_filtration_rate.GLPC either directly or by proxy. // In module 'Guyton_Model__kidney', the variables proximal_tubular_and_macula_densa_flow.GFN and GFN were unable to be set as equivalent: Loop detected: 'proximal_tubular_and_macula_densa_flow.GFN' may not be set to be equal to 'GFN' because GFN's definition already includes proximal_tubular_and_macula_densa_flow.GFN either directly or by proxy. // In module 'Guyton_Model____main', the variables circulatory_dynamics.venous_resistance.PC and PC were unable to be set as equivalent: Loop detected: 'circulatory_dynamics.venous_resistance.PC' may not be set to be equal to 'PC' because PC's definition already includes circulatory_dynamics.venous_resistance.PC either directly or by proxy. model Guyton_Model__environment(time_) // Variable initializations: time_ = ; end model Guyton_Model__temp_myogrs_and_pamk(MYOGRS, PAMK, PAMKRN) // Variable initializations: MYOGRS = 1.0; PAMK = 1.0; PAMKRN = 1.0; end model Guyton_Model__aldosterone_parameter_values(ANMALD, AMKMUL, ALDINF, ALDKNS, AMT, AM1UL, AM1LL, AMCSNS, ALDMM, AMKM, AMNAM, AMNAUL, AMNALL) // Variable initializations: ANMALD = 2.5; AMKMUL = 12; ALDINF = 0; ALDKNS = 0; AMT = 60; AM1UL = 5; AM1LL = 0; AMCSNS = 0.65; ALDMM = 2.5; AMKM = 0.5; AMNAM = 0.8; AMNAUL = 15; AMNALL = 0.04; end model Guyton_Model__angiotensin_control_of_aldosterone_secretion(ANM, ANMAL, ANMALD) // Assignment Rules: ANMAL := (ANM - 1) * ANMALD + 1; // Variable initializations: ANM = ; ANMALD = ; end model Guyton_Model__osmotic_control_of_aldosterone_secretion(CKE, OSMAL) // Assignment Rules: OSMAL := (CKE - 3.3) / 1.0; // Variable initializations: CKE = ; end model Guyton_Model__aldosterone_secretion(ANMAL, OSMAL, AMR1, AMKMUL, ALDINF, ALDKNS) // Assignment Rules: AMR1 := piecewise( ALDKNS , ALDKNS > 0 , AMR + ALDINF ); AMRBSC := ANMAL * 0.909 * OSMAL; AMRT := (AMRBSC - 1) * AMKMUL + 1; AMR := piecewise( 0 , AMRT < 0 , AMRT ); // Variable initializations: ANMAL = ; OSMAL = ; AMKMUL = ; ALDINF = ; ALDKNS = ; end model Guyton_Model__aldosterone_concentration(AMR1, time_, AMC, AMT) // Rate Rules: AMC' = (AMR1 - AMC) / AMT; // Variable initializations: AMR1 = ; time_ = ; AMC = 1.0; AMT = ; end model Guyton_Model__general_aldosterone_multiplier(AMC, AM, AM1UL, AM1LL, AMCSNS, ALDMM) // Assignment Rules: AM := (AM1 - 1) * ALDMM + 1; AM1 := AM1UL - (AM1UL - 1) / (((AM1LL - 1) / (AM1LL - AM1UL)) * (AMC - 1) * AMCSNS + 1); // Variable initializations: AMC = ; AM1UL = ; AM1LL = ; AMCSNS = ; ALDMM = ; end model Guyton_Model__aldosterone_effect_on_cell_membrane_K_transport(AM, AMK, AMKM) // Assignment Rules: AMK := piecewise( 0.2 , AMKT < 0.2 , AMKT ); AMKT := (AM - 1) * AMKM + 1; // Variable initializations: AM = ; AMKM = ; end model Guyton_Model__aldosterone_effect_on_cell_membrane_Na_transport(AM, AMNA, AMNAM, AMNAUL, AMNALL) // Assignment Rules: AMNA := piecewise( AMNALL , AMNAT < AMNALL , AMNAUL , AMNAT > AMNAUL , AMNAT ); AMNAT := (AM - 1) * AMNAM + 1; // Variable initializations: AM = ; AMNAM = ; AMNAUL = ; AMNALL = ; end model Guyton_Model__aldosterone(time_, ANM, CKE, AMK, AMNA) // Sub-modules, and any changes to those submodules: aldosterone_parameter_values: Guyton_Model__aldosterone_parameter_values(ANMALD, AMKMUL, ALDINF, ALDKNS, AMT, AM1UL, AM1LL, AMCSNS, ALDMM, AMKM, AMNAM, AMNAUL, AMNALL); angiotensin_control_of_aldosterone_secretion: Guyton_Model__angiotensin_control_of_aldosterone_secretion(ANM, ANMAL, ANMALD); osmotic_control_of_aldosterone_secretion: Guyton_Model__osmotic_control_of_aldosterone_secretion(CKE, OSMAL); aldosterone_secretion: Guyton_Model__aldosterone_secretion(ANMAL, OSMAL, AMR1, AMKMUL, ALDINF, ALDKNS); aldosterone_concentration: Guyton_Model__aldosterone_concentration(AMR1, time_, AMC, AMT); general_aldosterone_multiplier: Guyton_Model__general_aldosterone_multiplier(AMC, AM, AM1UL, AM1LL, AMCSNS, ALDMM); aldosterone_effect_on_cell_membrane_K_transport: Guyton_Model__aldosterone_effect_on_cell_membrane_K_transport(AM, AMK, AMKM); aldosterone_effect_on_cell_membrane_Na_transport: Guyton_Model__aldosterone_effect_on_cell_membrane_Na_transport(AM, AMNA, AMNAM, AMNAUL, AMNALL); end model Guyton_Model__angiotensin_parameter_values(ANXM, ANV, REK, ANGINF, ANGKNS, ANT, Z12, ANMUL, ANMLL, ANCSNS, ANULL, ANUM, ANUVM) // Variable initializations: ANXM = 0; ANV = 5000; REK = 1; ANGINF = 0; ANGKNS = 0; ANT = 12; Z12 = 5; ANMUL = 1.8; ANMLL = 0.7; ANCSNS = 0.4; ANULL = 0.8; ANUM = 6; ANUVM = 0; end model Guyton_Model__instantaneous_angiotensin_formation(MDFLW, ANGSCR) // Assignment Rules: ANGSCR := piecewise( 1 / (1 + (MDFLW3 - 1) * 72) , MDFLW3 > 1 , 10 - 9 / (1 + (1 - MDFLW3) * 8) ); MDFLW3 := MDFLW; // Variable initializations: MDFLW = ; end model Guyton_Model__time_delayed_angiotensin_formation(ANGSCR, time_, ANX1, ANXM, ANV) // Assignment Rules: ANX := (ANGSCR - 1) * ANXM; // Rate Rules: ANX1' = (ANX - ANX1) / ANV; // Variable initializations: ANGSCR = ; time_ = ; ANX1 = 0.0; ANXM = ; ANV = ; end model Guyton_Model__total_angiotensin_formation(ANGSCR, ANX1, ANPR, REK) // Assignment Rules: ANPR := piecewise( 0.00001 , ANPRT < 0.00001 , ANPRT ); ANPRT := (ANGSCR + ANX1) * REK; // Variable initializations: ANGSCR = ; ANX1 = ; REK = ; end model Guyton_Model__artificial_angiotensin_formation(ANPR, ANPR1, ANGKNS, ANGINF) // Assignment Rules: ANPR1 := piecewise( ANGKNS , ANGKNS > 0 , ANPR + ANGINF ); // Variable initializations: ANPR = ; ANGKNS = ; ANGINF = ; end model Guyton_Model__angiotensin_concentration(ANPR1, time_, ANC, ANT) // Rate Rules: ANC' = (ANPR1 - ANC) / ANT; // Variable initializations: ANPR1 = ; time_ = ; ANC = 0.859476; ANT = ; end model Guyton_Model__general_angiotensin_multiplier(ANC, ANM, ANMUL, ANMLL, ANCSNS) // Assignment Rules: ANM := ANMUL - (ANMUL - 1) / (((ANMLL - 1) / (ANMLL - ANMUL)) * (ANC - 1) * ANCSNS + 1); // Variable initializations: ANC = ; ANMUL = ; ANMLL = ; ANCSNS = ; end model Guyton_Model__angiotensin_effect_on_circulation(ANM, ANU, ANUM, ANULL) // Assignment Rules: ANU := piecewise( ANULL , ANU1 < ANULL , ANU1 ); ANU1 := (ANM - 1) * ANUM + 1; // Variable initializations: ANM = ; ANUM = ; ANULL = ; end model Guyton_Model__angiotensin_effect_on_venous_constriction(ANU, ANUVN, ANUVM) // Assignment Rules: ANUVN := (ANU - 1) * ANUVM + 1; // Variable initializations: ANU = ; ANUVM = ; end model Guyton_Model__angiotensin(time_, MDFLW, ANUVN, ANU, ANM) // Sub-modules, and any changes to those submodules: angiotensin_parameter_values: Guyton_Model__angiotensin_parameter_values(ANXM, ANV, REK, ANGINF, ANGKNS, ANT, Z12, ANMUL, ANMLL, ANCSNS, ANULL, ANUM, ANUVM); instantaneous_angiotensin_formation: Guyton_Model__instantaneous_angiotensin_formation(MDFLW, ANGSCR); time_delayed_angiotensin_formation: Guyton_Model__time_delayed_angiotensin_formation(ANGSCR, time_, ANX1, ANXM, ANV); total_angiotensin_formation: Guyton_Model__total_angiotensin_formation(ANGSCR, ANX1, ANPR, REK); artificial_angiotensin_formation: Guyton_Model__artificial_angiotensin_formation(ANPR, ANPR1, ANGKNS, ANGINF); angiotensin_concentration: Guyton_Model__angiotensin_concentration(ANPR1, time_, ANC, ANT); general_angiotensin_multiplier: Guyton_Model__general_angiotensin_multiplier(ANC, ANM, ANMUL, ANMLL, ANCSNS); angiotensin_effect_on_circulation: Guyton_Model__angiotensin_effect_on_circulation(ANM, ANU, ANUM, ANULL); angiotensin_effect_on_venous_constriction: Guyton_Model__angiotensin_effect_on_venous_constriction(ANU, ANUVN, ANUVM); end model Guyton_Model__antidiuretic_hormone_parameter_values(CNR, ADHPUL, ADHPAM, ADHINF, ADHTC, ADHVUL, ADHVLL, ADHKLL, ADHKUL) // Variable initializations: CNR = 139; ADHPUL = 85; ADHPAM = 0.3; ADHINF = 0; ADHTC = 15; ADHVUL = 2.5; ADHVLL = 0.93617; ADHKLL = 0.2; ADHKUL = 5; end model Guyton_Model__osmotic_control_of_ADH_secretion(CNA, ADHNA, CNR) // Assignment Rules: ADHNA := piecewise( 0 , ADHNA1 < 0 , ADHNA1 ); ADHNA1 := (CNA - CNR) / (142 - CNR); // Variable initializations: CNA = ; CNR = ; end model Guyton_Model__pressure_control_of_ADH_secretion(PA1, ADHPR, ADHPUL, ADHPAM) // Assignment Rules: ADHPR := power(ADHPUL - ADHPA, 2) * ADHPAM; ADHPA := piecewise( ADHPUL , PA1 > ADHPUL , PA1 ); // Variable initializations: PA1 = ; ADHPUL = ; ADHPAM = ; end model Guyton_Model__total_ADH_secretion(ADHNA, ADHPR, ADH, ADHINF) // Assignment Rules: ADH := piecewise( 0 , ADH1 < 0 , ADH1 ); ADH1 := ADHNA + ADHPR + ADHINF; // Variable initializations: ADHNA = ; ADHPR = ; ADHINF = ; end model Guyton_Model__ADH_in_blood(ADH, time_, ADHC, ADHTC) // Rate Rules: ADHC' = (ADH - ADHC) / ADHTC; // Variable initializations: ADH = ; time_ = ; ADHC = 1.0; ADHTC = ; end model Guyton_Model__ADH_effect_on_nonrenal_vascular_resistance(ADHC, ADHMV, ADHVUL, ADHVLL) // Assignment Rules: ADHMV := piecewise( ADHVLL , ADHMV1 < ADHVLL , ADHMV1 ); ADHMV1 := ADHVUL - (ADHVUL - 1) / (((ADHVLL - 1) / (ADHVLL - ADHVUL)) * (ADHC - 1) + 1); // Variable initializations: ADHC = ; ADHVUL = ; ADHVLL = ; end model Guyton_Model__ADH_effect_on_kidney(ADHC, ADHMK, ADHKLL, ADHKUL) // Assignment Rules: ADHMK := piecewise( ADHKLL , ADHMK1 < ADHKLL , ADHMK1 ); ADHMK1 := ADHKUL - (ADHKUL - 1) / (((ADHKLL - 1) / (ADHKLL - ADHKUL)) * (ADHC - 1) + 1); // Variable initializations: ADHC = ; ADHKLL = ; ADHKUL = ; end model Guyton_Model__antidiuretic_hormone(time_, CNA, PA1, ADHMK, ADHMV, ADHC) // Sub-modules, and any changes to those submodules: antidiuretic_hormone_parameter_values: Guyton_Model__antidiuretic_hormone_parameter_values(CNR, ADHPUL, ADHPAM, ADHINF, ADHTC, ADHVUL, ADHVLL, ADHKLL, ADHKUL); osmotic_control_of_ADH_secretion: Guyton_Model__osmotic_control_of_ADH_secretion(CNA, ADHNA, CNR); pressure_control_of_ADH_secretion: Guyton_Model__pressure_control_of_ADH_secretion(PA1, ADHPR, ADHPUL, ADHPAM); total_ADH_secretion: Guyton_Model__total_ADH_secretion(ADHNA, ADHPR, ADH, ADHINF); ADH_in_blood: Guyton_Model__ADH_in_blood(ADH, time_, ADHC, ADHTC); ADH_effect_on_nonrenal_vascular_resistance: Guyton_Model__ADH_effect_on_nonrenal_vascular_resistance(ADHC, ADHMV, ADHVUL, ADHVLL); ADH_effect_on_kidney: Guyton_Model__ADH_effect_on_kidney(ADHC, ADHMK, ADHKLL, ADHKUL); end model Guyton_Model__atrial_natriuretic_peptide_parameter_values(ANPKNS, ANPINF, ANPTC, ANPXUL) // Variable initializations: ANPKNS = 0; ANPINF = 0; ANPTC = 4; ANPXUL = 10; end model Guyton_Model__total_ANP_secreted(PLA, PRA, ANP) // Assignment Rules: ANP := (ANPL + ANPR2) / 3; ANPL := piecewise( 0 , (PLA - 1) * 1 < 0 , (PLA - 1) * 1 ); ANPR2 := piecewise( 0 , (PRA + 1) * 2 < 0 , (PRA + 1) * 2 ); // Variable initializations: PLA = ; PRA = ; end model Guyton_Model__ANP_into_circulation(ANP, ANP1, ANPKNS, ANPINF) // Assignment Rules: ANP1 := piecewise( ANPKNS , ANPKNS > 0 , ANP + ANPINF ); // Variable initializations: ANP = ; ANPKNS = ; ANPINF = ; end model Guyton_Model__ANP_in_plasma(ANP1, time_, ANPC, ANPTC) // Rate Rules: ANPC' = (ANP1 - ANPC) / ANPTC; // Variable initializations: ANP1 = ; time_ = ; ANPC = 1.0; ANPTC = ; end model Guyton_Model__ANP_effect_on_renal_afferent_arteriolar_resistance(ANPC, ANPX, ANPXUL) // Assignment Rules: ANPX := piecewise( -(1) , ANPX1 < -(1) , ANPX1 ); ANPX1 := ANPXUL - ANPXUL / (0.5555556 * (1 + ANPC)); // Variable initializations: ANPC = ; ANPXUL = ; end model Guyton_Model__atrial_natriuretic_peptide(time_, PLA, PRA, ANPX) // Sub-modules, and any changes to those submodules: atrial_natriuretic_peptide_parameter_values: Guyton_Model__atrial_natriuretic_peptide_parameter_values(ANPKNS, ANPINF, ANPTC, ANPXUL); total_ANP_secreted: Guyton_Model__total_ANP_secreted(PLA, PRA, ANP); ANP_into_circulation: Guyton_Model__ANP_into_circulation(ANP, ANP1, ANPKNS, ANPINF); ANP_in_plasma: Guyton_Model__ANP_in_plasma(ANP1, time_, ANPC, ANPTC); ANP_effect_on_renal_afferent_arteriolar_resistance: Guyton_Model__ANP_effect_on_renal_afferent_arteriolar_resistance(ANPC, ANPX, ANPXUL); end model Guyton_Model__autonomics_parameter_values(VV9, AUL, AULPM, AUY, AUV, AUS, O2A, AUM1, AUM2, AUDMP, AUMAX, AUMIN, AUSLP, EXC, EXCML, EXCXP, AUN1, MDMP, BAROTC, AUK, AUX, O2CHMO, AUC1, CRRFLX, EXE) // Variable initializations: VV9 = 2.51; AUL = 1.5; AULPM = 0; AUY = 0; AUV = 0.55; AUS = 1; O2A = 0.1; AUM1 = 3; AUM2 = 1; AUDMP = 0.3; AUMAX = 5.0; AUMIN = 0.4; AUSLP = 1.5; EXC = 1; EXCML = 0.01; EXCXP = 1.0; AUN1 = 0.5; MDMP = 0; BAROTC = 0.16; AUK = 0.004; AUX = 1; O2CHMO = 0.01; AUC1 = 0.3; CRRFLX = 0; EXE = 0; end model Guyton_Model__pressure_driving_autonomic_receptors(PA, PA1, CRRFLX, EXE) // Assignment Rules: PA1 := piecewise( CRRFLX , CRRFLX > 0.0000001 , PA - EXE ); // Variable initializations: PA = ; CRRFLX = ; EXE = ; end model Guyton_Model__chemoreceptors_effect_of_PA(PA1, AUC, AUC1) // Assignment Rules: AUC := piecewise( 0.005 * (80 - PA1) * AUC1 , ( PA1 < 80) && (geq(PA1, 40) ), 0.2 * AUC1 , PA1 < 40 , 0 ); // Variable initializations: PA1 = ; AUC1 = ; end model Guyton_Model__chemoreceptors_effect_of_art_PO2(PO2ART, AUC, AUC3, O2CHMO) // Assignment Rules: AUC3 := AUC + AUC2; AUC2 := piecewise( O2CHMO * (80 - PO2ART) , ( PO2ART < 80) && (geq(PO2ART, 40) ), O2CHMO * 40 , PO2ART < 40 , 0 ); // Variable initializations: PO2ART = ; AUC = ; O2CHMO = ; end model Guyton_Model__arterial_baroreceptor_reflex(PA1, time_, AU6C, AUX, AUK, BAROTC) // Assignment Rules: AU6C := AU6; AUB := piecewise( 0.016667 * (160 - PA1) , ( PA1 < 160) && (geq(PA1, 80) ), 1.3336 , PA1 < 80 , 0 ); A1B := (AUB - 1) * AUX + 1; AU6A := A1B - AU4; // Rate Rules: AU6' = (AU6A - AU6) / BAROTC; // Variable initializations: PA1 = ; time_ = ; AUX = ; AUK = ; BAROTC = ; AU4 = -0.060024; AU6 = 1.00132; end model Guyton_Model__CNS_ischemic_reflex(PA1, AUN, AUN1) // Assignment Rules: AUN := piecewise( 0.04 * (40 - PA1) * AUN1 , PA1 < 40 , 0 ); // Variable initializations: PA1 = ; AUN1 = ; end model Guyton_Model__autonomic_response_to_vasculature_pressure(PLA, PRA, PPA, AULP, AULPM) // Assignment Rules: AULP := (15 / (PLA + PRA + PPA) - 1) * AULPM + 1; // Variable initializations: PLA = ; PRA = ; PPA = ; AULPM = ; end model Guyton_Model__autonomic_response_to_exercise(AUEX, EXC, EXCXP) // Assignment Rules: AUEX := power(EXC, EXCXP); // Variable initializations: EXC = ; EXCXP = ; end model Guyton_Model__total_autonomic_stimulation(AUC3, AU6C, AUN, AULP, AUEX, AUTTL, EXCML) // Assignment Rules: AUTTL := piecewise( 0 , AUTTL1 < 0 , AUTTL1 ); AUTTL1 := (AUEX * AULP * (AUC3 + AU6C + AUN) - 1) * EXCML + 1; // Variable initializations: AUC3 = ; AU6C = ; AUN = ; AULP = ; AUEX = ; EXCML = ; end model Guyton_Model__actual_autonomic_stimulation(AUTTL, time_, AU, AUDMP, AUMAX, AUMIN, AUSLP) // Assignment Rules: AU := piecewise( AUMIN , AUT < AUMIN , AUT ); DAU := (AUTTL - AU1) / AUDMP; AUT := AUMAX - (AUMAX - 1) / exp(AUSLP * (AU1 - 1)); // Rate Rules: AU1' = DAU; // Variable initializations: AUTTL = ; time_ = ; AUDMP = ; AUMAX = ; AUMIN = ; AUSLP = ; AU1 = 1.00007; end model Guyton_Model__autonomic_drive_on_target_organs_and_tissues(AU, VVR, AUH, AUR, AOM, AUM, AVE, VV9, AUL, AUV, AUS, O2A, AUM1, AUM2, AUY) // Assignment Rules: VVR := (VV9 - AU * AUL) + AUL; AUH := AUO * AUV + 1; AUR := AUO * AUS + 1; AOM := AUO * O2A + 1; AUM := power(AUO * AUM1 + 1, AUM2); AVE := AUO * AUY + 1; AUO := AU - 1; // Variable initializations: AU = ; VV9 = ; AUL = ; AUV = ; AUS = ; O2A = ; AUM1 = ; AUM2 = ; AUY = ; end model Guyton_Model__autonomics(time_, PA, PO2ART, PLA, PRA, PPA, PA1, VVR, AUH, AUR, AUM, AOM, AVE, AU) // Sub-modules, and any changes to those submodules: autonomics_parameter_values: Guyton_Model__autonomics_parameter_values(VV9, AUL, AULPM, AUY, AUV, AUS, O2A, AUM1, AUM2, AUDMP, AUMAX, AUMIN, AUSLP, EXC, EXCML, EXCXP, AUN1, MDMP, BAROTC, AUK, AUX, O2CHMO, AUC1, CRRFLX, EXE); pressure_driving_autonomic_receptors: Guyton_Model__pressure_driving_autonomic_receptors(PA, PA1, CRRFLX, EXE); chemoreceptors_effect_of_PA: Guyton_Model__chemoreceptors_effect_of_PA(PA1, AUC, AUC1); chemoreceptors_effect_of_art_PO2: Guyton_Model__chemoreceptors_effect_of_art_PO2(PO2ART, AUC, AUC3, O2CHMO); arterial_baroreceptor_reflex: Guyton_Model__arterial_baroreceptor_reflex(PA1, time_, AU6C, AUX, AUK, BAROTC); CNS_ischemic_reflex: Guyton_Model__CNS_ischemic_reflex(PA1, AUN, AUN1); autonomic_response_to_vasculature_pressure: Guyton_Model__autonomic_response_to_vasculature_pressure(PLA, PRA, PPA, AULP, AULPM); autonomic_response_to_exercise: Guyton_Model__autonomic_response_to_exercise(AUEX, EXC, EXCXP); total_autonomic_stimulation: Guyton_Model__total_autonomic_stimulation(AUC3, AU6C, AUN, AULP, AUEX, AUTTL, EXCML); actual_autonomic_stimulation: Guyton_Model__actual_autonomic_stimulation(AUTTL, time_, AU, AUDMP, AUMAX, AUMIN, AUSLP); autonomic_drive_on_target_organs_and_tissues: Guyton_Model__autonomic_drive_on_target_organs_and_tissues(AU, VVR, AUH, AUR, AOM, AUM, AVE, VV9, AUL, AUV, AUS, O2A, AUM1, AUM2, AUY); end model Guyton_Model__capillary_dynamics_parameter_values(CFC, TRPL, CPR, LPDE, PCE, PCR, LPK, LPPR, CPK, TSSLML, TSSLTC, HYL, CMPTSS, PGHF, GCOPF, VTSF, PLDF) // Variable initializations: CFC = 0.01167; TRPL = 0; CPR = 40; LPDE = 8; PCE = 1; PCR = 15; LPK = 2.728e-14; LPPR = 0.03; CPK = 0.000253; TSSLML = 0.15; TSSLTC = 0.005; HYL = 60; CMPTSS = 2; PGHF = -2; GCOPF = 0.8092; VTSF = 6; PLDF = 4.2; end model Guyton_Model__capillary_membrane_dynamics() end model Guyton_Model__capillary_pressure(RVS, BFN, PVS, time_, PC) // Assignment Rules: PC := RVS * 1.7 * BFN + PVS; // Variable initializations: RVS = ; BFN = ; PVS = ; time_ = ; end model Guyton_Model__rate_of_fluid_out_of_capillaries(PC, PGH, PTC, PPC, VTCPL, VTC, CFC) // Assignment Rules: VTC := (((PC - PPC) - PGH) + PTC) * CFC + VTCPL; // Variable initializations: PC = ; PGH = ; PTC = ; PPC = ; VTCPL = ; CFC = ; end model Guyton_Model__plasma_volume_and_protein() end model Guyton_Model__plasma_volume(DFP, VTC, VTL, time_, VP, TRPL) // Assignment Rules: VPD := ((VTL - VTC) - DFP) + TRPL; // Rate Rules: VP' = VPD; // Variable initializations: DFP = ; VTC = ; VTL = ; time_ = ; VP = 3.00449; TRPL = ; end model Guyton_Model__plasma_protein_concentration(PRP, VP, CPP) // Assignment Rules: CPP := PRP / VP; // Variable initializations: PRP = ; VP = ; end model Guyton_Model__protein_destruction_and_formation(CPP, DLP, CPR, LPPR, LPDE, LPK) // Assignment Rules: DLP := LPPR - power(CPPD, LPDE) * LPK; CPPD := piecewise( 0 , CPP - CPR < 0 , CPP - CPR ); // Variable initializations: CPP = ; CPR = ; LPPR = ; LPDE = ; LPK = ; end model Guyton_Model__plasma_leakage(PC, VTCPL, PCR, CPK, PCE) // Assignment Rules: VTCPL := power(PRCD * CPK, PCE); PRCD := piecewise( 0 , PC - PCR < 0 , PC - PCR ); // Variable initializations: PC = ; PCR = ; CPK = ; PCE = ; end model Guyton_Model__protein_influx_into_interstitium(VTCPL, CPP, CPI, DPC) // Assignment Rules: DPC := VTCPL * CPP + (CPP - CPI) * 0.00104; // Variable initializations: VTCPL = ; CPP = ; CPI = ; end model Guyton_Model__total_plasma_protein(DPC, PPD, DPL, DLP, time_, PRP, TRPL) // Assignment Rules: DPP := (((DLP + DPL) - DPC) - PPD) + TRPL * 72; // Rate Rules: PRP' = DPP; // Variable initializations: DPC = ; PPD = ; DPL = ; DLP = ; time_ = ; PRP = 216.243; TRPL = ; end model Guyton_Model__plasma_colloid_osmotic_pressure(CPP, PPC) // Assignment Rules: PPC := 0.28 * CPP + 0.0019 * power(CPP, 2); // Variable initializations: CPP = ; end model Guyton_Model__systemic_tissue_fluid_volume_and_protein() end model Guyton_Model__total_systemic_fluid_volume(VEC, VP, VPF, VTS) // Assignment Rules: VTS := (VEC - VP) - VPF; // Variable initializations: VEC = ; VP = ; VPF = ; end model Guyton_Model__interstitial_fluid_volume(VTS, time_, VTS1, TSSLML, TSSLTC) // Assignment Rules: VTS1 := VTS - VTS2; // Rate Rules: VTS2' = ((VTS - 12) * TSSLML - VTS2) * TSSLTC; // Variable initializations: VTS = ; time_ = ; TSSLML = ; TSSLTC = ; VTS2 = 0.0; end model Guyton_Model__total_interstitial_protein(DPC, DPL, time_, TSP) // Assignment Rules: DPI := DPC - DPL; // Rate Rules: TSP' = DPI; // Variable initializations: DPC = ; DPL = ; time_ = ; TSP = 279.945; end model Guyton_Model__interstitial_protein_concentration(TSP, VTS, CPI) // Assignment Rules: CPI := TSP / VTS; // Variable initializations: TSP = ; VTS = ; end model Guyton_Model__interstitial_colloid_osmotic_pressure(CPI, PTCPR) // Assignment Rules: PTCPR := 0.28 * CPI + 0.0019 * power(CPI, 2); // Variable initializations: CPI = ; end model Guyton_Model__lymph_protein_flow(CPI, VTL, DPL) // Assignment Rules: DPL := CPI * VTL; // Variable initializations: CPI = ; VTL = ; end model Guyton_Model__tissue_gel_and_fluid_and_lymph_flow() end model Guyton_Model__hydrostatic_pressure_of_tissue_gel(VTS, PTT, CHY, PGH, HYL, CMPTSS, PGHF) // Assignment Rules: CHY := power((HYL / VTS) / 5, CMPTSS); PGH := CHY * PGHF + PTT; // Variable initializations: VTS = ; PTT = ; HYL = ; CMPTSS = ; PGHF = ; end model Guyton_Model__total_osmotic_pressure_of_tissue_gel(CHY, PTCPR, POSHYL, PTC, GCOPF) // Assignment Rules: POSHYL := CHY * 2; PTC := POSHYL * PTCPR * GCOPF; // Variable initializations: CHY = ; PTCPR = ; GCOPF = ; end model Guyton_Model__total_tissue_pressure(VTS1, PTT, VTSF) // Assignment Rules: PTT := power((VTS1 - VTSF) / VTSF, 2) * 1; // Variable initializations: VTS1 = ; VTSF = ; end model Guyton_Model__interstial_free_fluid_pressure(PGH, POSHYL, PIF) // Assignment Rules: PIF := PGH - POSHYL; // Variable initializations: PGH = ; POSHYL = ; end model Guyton_Model__interstitial_solid_tissue_pressure(PIF, PTT, PTS) // Assignment Rules: PTS := PTT - PIF; // Variable initializations: PIF = ; PTT = ; end model Guyton_Model__lymph_flow(PIF, PTT, VTL, PLDF) // Assignment Rules: VTL := piecewise( 0 , PLD < 0 , PLD * 0.02 ); PLD1 := (PIF + PLDF) - PTT; PLD := piecewise( 7 , PLD1 > 7 , PLD1 ); // Variable initializations: PIF = ; PTT = ; PLDF = ; end model Guyton_Model__interstitial_gel_volume(VTS, VG) // Assignment Rules: VG := piecewise( 0 , VTS <= 0 , 0 + ((11.4 - 0) * (VTS - 0)) / (12 - 0) , ( VTS > 0) && (VTS <= 12 ), 11.4 + ((14 - 11.4) * (VTS - 12)) / (15 - 12) , ( VTS > 12) && (VTS <= 15 ), 14 + ((16 - 14) * (VTS - 15)) / (18 - 15) , ( VTS > 15) && (VTS <= 18 ), 16 + ((17.3 - 16) * (VTS - 18)) / (21 - 18) , ( VTS > 18) && (VTS <= 21 ), 17.3 + ((18 - 17.3) * (VTS - 21)) / (24 - 21) , ( VTS > 21) && (VTS <= 24 ), 18 ); // Variable initializations: VTS = ; end model Guyton_Model__interstitial_free_fluid_volume(VTS, VG, VIF) // Assignment Rules: VIF := VTS - VG; // Variable initializations: VTS = ; VG = ; end model Guyton_Model__capillary_dynamics(time_, VEC, PPD, RVS, DFP, VPF, BFN, PVS, PC, CPP, PPC, VP) // Sub-modules, and any changes to those submodules: capillary_dynamics_parameter_values: Guyton_Model__capillary_dynamics_parameter_values(CFC, TRPL, CPR, LPDE, PCE, PCR, LPK, LPPR, CPK, TSSLML, TSSLTC, HYL, CMPTSS, PGHF, GCOPF, VTSF, PLDF); capillary_membrane_dynamics: Guyton_Model__capillary_membrane_dynamics(); capillary_pressure: Guyton_Model__capillary_pressure(RVS, BFN, PVS, time_, PC); rate_of_fluid_out_of_capillaries: Guyton_Model__rate_of_fluid_out_of_capillaries(PC, PGH, PTC, PPC, VTCPL, VTC, CFC); plasma_volume_and_protein: Guyton_Model__plasma_volume_and_protein(); plasma_volume: Guyton_Model__plasma_volume(DFP, VTC, VTL, time_, VP, TRPL); plasma_protein_concentration: Guyton_Model__plasma_protein_concentration(PRP, VP, CPP); protein_destruction_and_formation: Guyton_Model__protein_destruction_and_formation(CPP, DLP, CPR, LPPR, LPDE, LPK); plasma_leakage: Guyton_Model__plasma_leakage(PC, VTCPL, PCR, CPK, PCE); protein_influx_into_interstitium: Guyton_Model__protein_influx_into_interstitium(VTCPL, CPP, CPI, DPC); total_plasma_protein: Guyton_Model__total_plasma_protein(DPC, PPD, DPL, DLP, time_, PRP, TRPL); plasma_colloid_osmotic_pressure: Guyton_Model__plasma_colloid_osmotic_pressure(CPP, PPC); systemic_tissue_fluid_volume_and_protein: Guyton_Model__systemic_tissue_fluid_volume_and_protein(); total_systemic_fluid_volume: Guyton_Model__total_systemic_fluid_volume(VEC, VP, VPF, VTS); interstitial_fluid_volume: Guyton_Model__interstitial_fluid_volume(VTS, time_, VTS1, TSSLML, TSSLTC); total_interstitial_protein: Guyton_Model__total_interstitial_protein(DPC, DPL, time_, TSP); interstitial_protein_concentration: Guyton_Model__interstitial_protein_concentration(TSP, VTS, CPI); interstitial_colloid_osmotic_pressure: Guyton_Model__interstitial_colloid_osmotic_pressure(CPI, PTCPR); lymph_protein_flow: Guyton_Model__lymph_protein_flow(CPI, VTL, DPL); tissue_gel_and_fluid_and_lymph_flow: Guyton_Model__tissue_gel_and_fluid_and_lymph_flow(); hydrostatic_pressure_of_tissue_gel: Guyton_Model__hydrostatic_pressure_of_tissue_gel(VTS, PTT, CHY, PGH, HYL, CMPTSS, PGHF); total_osmotic_pressure_of_tissue_gel: Guyton_Model__total_osmotic_pressure_of_tissue_gel(CHY, PTCPR, POSHYL, PTC, GCOPF); total_tissue_pressure: Guyton_Model__total_tissue_pressure(VTS1, PTT, VTSF); interstial_free_fluid_pressure: Guyton_Model__interstial_free_fluid_pressure(PGH, POSHYL, PIF); interstitial_solid_tissue_pressure: Guyton_Model__interstitial_solid_tissue_pressure(PIF, PTT, PTS); lymph_flow: Guyton_Model__lymph_flow(PIF, PTT, VTL, PLDF); interstitial_gel_volume: Guyton_Model__interstitial_gel_volume(VTS, VG); interstitial_free_fluid_volume: Guyton_Model__interstitial_free_fluid_volume(VTS, VG, VIF); end model Guyton_Model__circulatory_dynamics_parameter_values(HTAUML, QRF, HSR, HSL, ANY, CV, PR1LL, RVSM, CN2, CN7, PAEX, RAR, RMULT1, RAM, FIS) // Variable initializations: HTAUML = 0.4; QRF = 0.15; HSR = 1; HSL = 1; ANY = -0.2; CV = 0.1; PR1LL = 0; RVSM = 1; CN2 = 0.0212; CN7 = 0.2; PAEX = 2; RAR = 30.52; RMULT1 = 1; RAM = 96.3; FIS = 0; end model Guyton_Model__total_blood_volume_change(VP, VRC, VVS1, VAS1, VLA1, VPA1, VRA1, VBD) // Variable initializations: VP = ; VRC = ; VVS1 = ; VAS1 = ; VLA1 = ; VPA1 = ; VRA1 = ; end model Guyton_Model__right_atrium() end model Guyton_Model__right_atrial_blood_volume(QVO, QRO, VBD, time_, VRA, VRA1) // Variable initializations: QVO = ; QRO = ; VBD = ; time_ = ; VRA1 = 0.100043; end model Guyton_Model__right_atrial_pressure(VRA, PRA) // Variable initializations: VRA = ; end model Guyton_Model__autonomic_stimulation_effect_on_right_atrial_pressure(PRA, AU, PRA1, HTAUML) // Variable initializations: PRA = ; AU = ; HTAUML = ; end model Guyton_Model__right_ventricle() end model Guyton_Model__pressure_effect_on_right_ventricular_pumping(PPA, AUH, OSA, RVM) // Variable initializations: PPA = ; AUH = ; OSA = ; end model Guyton_Model__pumping_effectiveness_of_right_ventricle(AUH, RVM, HMD, HPR, QLO, QLN, HPEF, QRF, HSR) // Variable initializations: AUH = ; RVM = ; HMD = ; HPR = ; QLO = ; QLN = ; QRF = ; HSR = ; end model Guyton_Model__right_ventricular_output(PRA1, HPEF, QRO) // Variable initializations: PRA1 = ; HPEF = ; end model Guyton_Model__pulmonary_vasculature() end model Guyton_Model__pulmonary_vasculature_blood_volume(QRO, QPO, VBD, time_, VPA, VPA1) // Variable initializations: QRO = ; QPO = ; VBD = ; time_ = ; VPA1 = 0.38131; end model Guyton_Model__pulmonary_vasculature_pressure(VPA, PPA) // Variable initializations: VPA = ; end model Guyton_Model__pulmonary_arterial_resistance(PPA, RPA) // Variable initializations: PPA = ; end model Guyton_Model__pulmonary_venous_resistance(PLA, RPV) // Variable initializations: PLA = ; end model Guyton_Model__total_pulmonary_vascular_resistance(RPV, RPA, RPT) // Variable initializations: RPV = ; RPA = ; end model Guyton_Model__pressure_gradient_through_the_lungs(PPA, PLA, PGL) // Variable initializations: PPA = ; PLA = ; end model Guyton_Model__rate_of_blood_flow_from_pulmonary_veins_to_left_atrium(PGL, RPT, QLO, QPO) // Variable initializations: PGL = ; RPT = ; QLO = ; end model Guyton_Model__left_atrium() end model Guyton_Model__left_atrial_blood_volume(QPO, QLO, VBD, time_, VLA, VLA1) // Variable initializations: QPO = ; QLO = ; VBD = ; time_ = ; VLA1 = 0.379883; end model Guyton_Model__left_atrial_pressure(VLA, PLA) // Variable initializations: VLA = ; end model Guyton_Model__autonomic_stimulation_effect_on_left_atrial_pressure(PLA, AU, PLA1, HTAUML) // Variable initializations: PLA = ; AU = ; HTAUML = ; end model Guyton_Model__left_ventricle() end model Guyton_Model__pumping_effectiveness_of_left_ventricle(PA, AUH, OSA, LVM) // Variable initializations: PA = ; AUH = ; OSA = ; end model Guyton_Model__left_ventricular_output(PLA1, LVM, AUH, HMD, HPL, PLA, PA, QLO, QLOT, QLN, HSL) // Variable initializations: PLA1 = ; LVM = ; AUH = ; HMD = ; HPL = ; PLA = ; PA = ; HSL = ; end model Guyton_Model__systemic_venous_system() end model Guyton_Model__venous_blood_volume(QAO, QVO, VBD, time_, VVS, VVS1) // Variable initializations: QAO = ; QVO = ; VBD = ; time_ = ; VVS1 = 3.28246; end model Guyton_Model__angiotensin_induced_venous_constriction(ANU, VVA, ANY) // Variable initializations: ANU = ; ANY = ; end model Guyton_Model__venous_excess_volume(VVS, VVR, VVA, VV6, VV7, ATRVFB, VVE) // Variable initializations: VVS = ; VVR = ; VVA = ; VV6 = ; VV7 = ; ATRVFB = ; end model Guyton_Model__venous_average_pressure(VVE, PVS, CV) // Variable initializations: VVE = ; CV = ; end model Guyton_Model__venous_outflow_pressure_into_heart(PRA, PR1, PR1LL) // Variable initializations: PRA = ; PR1LL = ; end model Guyton_Model__resistance_from_veins_to_right_atrium(PVS, VIM, RVG) // Variable initializations: PVS = ; VIM = ; end model Guyton_Model__rate_of_blood_flow_from_veins_to_right_atrium(PVS, PR1, RVG, QRO, QVO) // Variable initializations: PVS = ; PR1 = ; RVG = ; QRO = ; end model Guyton_Model__venous_resistance(PC, VIM, AVE, ANUVN, RVS, CN7, CN2, RVSM) // Variable initializations: PC = ; VIM = ; AVE = ; ANUVN = ; CN7 = ; CN2 = ; RVSM = ; end model Guyton_Model__NM_NR_venous_resistance(RVS, NNRVR) // Variable initializations: RVS = ; end model Guyton_Model__systemic_arterial_system() end model Guyton_Model__arterial_blood_volume(QLO, QAO, VBD, time_, VAS, VAS1) // Variable initializations: QLO = ; QAO = ; VBD = ; time_ = ; VAS1 = 0.862514; end model Guyton_Model__arterial_pressure_and_pressure_gradient(VAS, PRA, PA, PAG) // Variable initializations: VAS = ; PRA = ; end model Guyton_Model__pressure_effect_on_arterial_distention(PA, PAM, PAEX) // Variable initializations: PA = ; PAEX = ; end model Guyton_Model__non_renal_systemic_arterial_resistance_multiplier(ANU, ADHMV, AUM, VIM, PAMK, PAM, ATRRFB, R1) // Variable initializations: ANU = ; ADHMV = ; AUM = ; VIM = ; PAMK = ; PAM = ; ATRRFB = ; end model Guyton_Model__NM_NR_arterial_resistance(R1, ARM, MYOGRS, NNRAR, RAR, RMULT1) // Variable initializations: R1 = ; ARM = ; MYOGRS = ; RAR = ; RMULT1 = ; end model Guyton_Model__pressure_gradient_from_arteries_to_veins(PA, PVS, PGS) // Variable initializations: PA = ; PVS = ; end model Guyton_Model__M_systemic_resistance(R1, AMM, MYOGRS, RSM, RAM, RMULT1) // Variable initializations: R1 = ; AMM = ; MYOGRS = ; RAM = ; RMULT1 = ; end model Guyton_Model__total_NM_NR_systemic_resistance(NNRAR, NNRVR, RSN) // Variable initializations: NNRAR = ; NNRVR = ; end model Guyton_Model__blood_flow_through_M_tissues(PGS, RSM, BFM) // Variable initializations: PGS = ; RSM = ; end model Guyton_Model__blood_flow_through_NM_NR_tissues(PGS, RSN, BFN) // Variable initializations: PGS = ; RSN = ; end model Guyton_Model__blood_flow_through_AV_fistulas(PAG, FISFLO, FIS) // Variable initializations: PAG = ; FIS = ; end model Guyton_Model__systemic_blood_flow(BFM, BFN, RBF, FISFLO, SYSFLO, QAO) // Variable initializations: BFM = ; BFN = ; RBF = ; FISFLO = ; end model Guyton_Model__total_peripheral_resistance(PAG, QAO, RTP) // Variable initializations: PAG = ; QAO = ; end model Guyton_Model__circulatory_dynamics(time_, ADHMV, AMM, ANU, ANUVN, ARM, ATRRFB, ATRVFB, AU, AUH, AUM, AVE, HMD, HPL, HPR, MYOGRS, OSA, PAMK, PC, RBF, VIM, VP, VRC, VV6, VV7, VVR, BFM, BFN, PA, PLA, PPA, PRA, PVS, QAO, QRO, QLO, RPA, RPV, RVS, VVE) // Sub-modules, and any changes to those submodules: circulatory_dynamics_parameter_values: Guyton_Model__circulatory_dynamics_parameter_values(HTAUML, QRF, HSR, HSL, ANY, CV, PR1LL, RVSM, CN2, CN7, PAEX, RAR, RMULT1, RAM, FIS); total_blood_volume_change: Guyton_Model__total_blood_volume_change(VP, VRC, VVS1, VAS1, VLA1, VPA1, VRA1, VBD); right_atrium: Guyton_Model__right_atrium(); right_atrial_blood_volume: Guyton_Model__right_atrial_blood_volume(QVO, QRO, VBD, time_, VRA, VRA1); right_atrial_pressure: Guyton_Model__right_atrial_pressure(VRA, PRA); autonomic_stimulation_effect_on_right_atrial_pressure: Guyton_Model__autonomic_stimulation_effect_on_right_atrial_pressure(PRA, AU, PRA1, HTAUML); right_ventricle: Guyton_Model__right_ventricle(); pressure_effect_on_right_ventricular_pumping: Guyton_Model__pressure_effect_on_right_ventricular_pumping(PPA, AUH, OSA, RVM); pumping_effectiveness_of_right_ventricle: Guyton_Model__pumping_effectiveness_of_right_ventricle(AUH, RVM, HMD, HPR, QLO0, QLN, HPEF, QRF, HSR); right_ventricular_output: Guyton_Model__right_ventricular_output(PRA1, HPEF, QRO); pulmonary_vasculature: Guyton_Model__pulmonary_vasculature(); pulmonary_vasculature_blood_volume: Guyton_Model__pulmonary_vasculature_blood_volume(QRO, QPO, VBD, time_, VPA, VPA1); pulmonary_vasculature_pressure: Guyton_Model__pulmonary_vasculature_pressure(VPA, PPA); pulmonary_arterial_resistance: Guyton_Model__pulmonary_arterial_resistance(PPA, RPA); pulmonary_venous_resistance: Guyton_Model__pulmonary_venous_resistance(PLA, RPV); total_pulmonary_vascular_resistance: Guyton_Model__total_pulmonary_vascular_resistance(RPV, RPA, RPT); pressure_gradient_through_the_lungs: Guyton_Model__pressure_gradient_through_the_lungs(PPA, PLA, PGL); rate_of_blood_flow_from_pulmonary_veins_to_left_atrium: Guyton_Model__rate_of_blood_flow_from_pulmonary_veins_to_left_atrium(PGL, RPT, QLO0, QPO); left_atrium: Guyton_Model__left_atrium(); left_atrial_blood_volume: Guyton_Model__left_atrial_blood_volume(QPO, QLO0, VBD, time_, VLA, VLA1); left_atrial_pressure: Guyton_Model__left_atrial_pressure(VLA, PLA); autonomic_stimulation_effect_on_left_atrial_pressure: Guyton_Model__autonomic_stimulation_effect_on_left_atrial_pressure(PLA, AU, PLA1, HTAUML); left_ventricle: Guyton_Model__left_ventricle(); pumping_effectiveness_of_left_ventricle: Guyton_Model__pumping_effectiveness_of_left_ventricle(PA, AUH, OSA, LVM); left_ventricular_output: Guyton_Model__left_ventricular_output(PLA1, LVM, AUH, HMD, HPL, PLA, PA, QLO0, QLOT, QLN, HSL); systemic_venous_system: Guyton_Model__systemic_venous_system(); venous_blood_volume: Guyton_Model__venous_blood_volume(QAO0, QVO, VBD, time_, VVS, VVS1); angiotensin_induced_venous_constriction: Guyton_Model__angiotensin_induced_venous_constriction(ANU, VVA, ANY); venous_excess_volume: Guyton_Model__venous_excess_volume(VVS, VVR, VVA, VV6, VV7, ATRVFB, VVE); venous_average_pressure: Guyton_Model__venous_average_pressure(VVE, PVS, CV); venous_outflow_pressure_into_heart: Guyton_Model__venous_outflow_pressure_into_heart(PRA, PR1, PR1LL); resistance_from_veins_to_right_atrium: Guyton_Model__resistance_from_veins_to_right_atrium(PVS, VIM, RVG); rate_of_blood_flow_from_veins_to_right_atrium: Guyton_Model__rate_of_blood_flow_from_veins_to_right_atrium(PVS, PR1, RVG, QRO, QVO); venous_resistance: Guyton_Model__venous_resistance(PC, VIM, AVE, ANUVN, RVS, CN7, CN2, RVSM); NM_NR_venous_resistance: Guyton_Model__NM_NR_venous_resistance(RVS, NNRVR); systemic_arterial_system: Guyton_Model__systemic_arterial_system(); arterial_blood_volume: Guyton_Model__arterial_blood_volume(QLO, QAO0, VBD, time_, VAS, VAS1); arterial_pressure_and_pressure_gradient: Guyton_Model__arterial_pressure_and_pressure_gradient(VAS, PRA, PA, PAG); pressure_effect_on_arterial_distention: Guyton_Model__pressure_effect_on_arterial_distention(PA, PAM, PAEX); non_renal_systemic_arterial_resistance_multiplier: Guyton_Model__non_renal_systemic_arterial_resistance_multiplier(ANU, ADHMV, AUM, VIM, PAMK, PAM, ATRRFB, R1); NM_NR_arterial_resistance: Guyton_Model__NM_NR_arterial_resistance(R1, ARM, MYOGRS, NNRAR, RAR, RMULT1); pressure_gradient_from_arteries_to_veins: Guyton_Model__pressure_gradient_from_arteries_to_veins(PA, PVS, PGS); M_systemic_resistance: Guyton_Model__M_systemic_resistance(R1, AMM, MYOGRS, RSM, RAM, RMULT1); total_NM_NR_systemic_resistance: Guyton_Model__total_NM_NR_systemic_resistance(NNRAR, NNRVR, RSN); blood_flow_through_M_tissues: Guyton_Model__blood_flow_through_M_tissues(PGS, RSM, BFM); blood_flow_through_NM_NR_tissues: Guyton_Model__blood_flow_through_NM_NR_tissues(PGS, RSN, BFN); blood_flow_through_AV_fistulas: Guyton_Model__blood_flow_through_AV_fistulas(PAG, FISFLO, FIS); systemic_blood_flow: Guyton_Model__systemic_blood_flow(BFM, BFN, RBF, FISFLO, SYSFLO, QAO0); total_peripheral_resistance: Guyton_Model__total_peripheral_resistance(PAG, QAO, RTP); QLO0 is QLO; QAO0 is QAO; // Assignment Rules: BFM := PGS / RSM; BFN := PGS / RSN; PA := arterial_pressure_and_pressure_gradient.VAE / 0.00355; PLA := left_atrial_pressure.VLE / 0.01; PPA := pulmonary_vasculature_pressure.VPE / 0.0048; PRA := right_atrial_pressure.VRE / 0.005; PVS := piecewise( 0.0001 , venous_average_pressure.PVS1 < 0.0001 , venous_average_pressure.PVS1 ); QAO := SYSFLO + FISFLO; QRO := right_ventricular_output.QRN * HPEF; QLO := piecewise( QLOT + left_ventricular_output.QLO1 , left_ventricular_output.QLO1 > 0 , QLOT ); RPA := 1 / pulmonary_arterial_resistance.CPA; RPV := 1 / (pulmonary_venous_resistance.PL1 * 0.0357); RVS := AVE * venous_resistance.RV1 * VIM * ANUVN; VVE := piecewise( 0.0001 , venous_excess_volume.VVE1 < 0.0001 , venous_excess_volume.VVE1 ); VBD := ((((((VP + VRC) - VVS1) - VAS1) - VLA1) - VPA1) - VRA1) / 2; QVO := rate_of_blood_flow_from_veins_to_right_atrium.PGV / RVG; VRA := VRA1 + VBD * 0.0574; right_atrial_blood_volume.DRA := QVO - QRO; right_atrial_pressure.VRE := VRA - 0.1; PRA1 := (PRA + 8) * (HTAUML * (AU - 1) + 1) - 8; RVM := piecewise( 1.06 , pressure_effect_on_right_ventricular_pumping.PP2 <= 0 , 1.06 + ((0.97 - 1.06) * (pressure_effect_on_right_ventricular_pumping.PP2 - 0)) / (32 - 0) , ( pressure_effect_on_right_ventricular_pumping.PP2 > 0) && (pressure_effect_on_right_ventricular_pumping.PP2 <= 32 ), 0.97 + ((0.93 - 0.97) * (pressure_effect_on_right_ventricular_pumping.PP2 - 32)) / (38.4 - 32) , ( pressure_effect_on_right_ventricular_pumping.PP2 > 32) && (pressure_effect_on_right_ventricular_pumping.PP2 <= 38.4 ), 0.93 + ((0.8 - 0.93) * (pressure_effect_on_right_ventricular_pumping.PP2 - 38.4)) / (48 - 38.4) , ( pressure_effect_on_right_ventricular_pumping.PP2 > 38.4) && (pressure_effect_on_right_ventricular_pumping.PP2 <= 48 ), 0.8 + ((0.46 - 0.8) * (pressure_effect_on_right_ventricular_pumping.PP2 - 48)) / (60.8 - 48) , ( pressure_effect_on_right_ventricular_pumping.PP2 > 48) && (pressure_effect_on_right_ventricular_pumping.PP2 <= 60.8 ), 0.46 + ((0 - 0.46) * (pressure_effect_on_right_ventricular_pumping.PP2 - 60.8)) / (72 - 60.8) , ( pressure_effect_on_right_ventricular_pumping.PP2 > 60.8) && (pressure_effect_on_right_ventricular_pumping.PP2 <= 72 ), 0 ); pressure_effect_on_right_ventricular_pumping.PP2 := (PPA / AUH) / OSA; QLN := piecewise( 0.01 , PLA1 <= -(2) , 0.01 + ((3.6 - 0.01) * (PLA1 - -(2))) / (1 - -(2)) , ( PLA1 > -(2)) && (PLA1 <= 1 ), 3.6 + ((9.4 - 3.6) * (PLA1 - 1)) / (5 - 1) , ( PLA1 > 1) && (PLA1 <= 5 ), 9.4 + ((11.6 - 9.4) * (PLA1 - 5)) / (8 - 5) , ( PLA1 > 5) && (PLA1 <= 8 ), 11.6 + ((13.5 - 11.6) * (PLA1 - 8)) / (12 - 8) , ( PLA1 > 8) && (PLA1 <= 12 ), 13.5 ); HPEF := (1 - QRF) * AUH * RVM * HSR * HMD * HPR + (QRF * QLO) / QLN; right_ventricular_output.QRN := piecewise( 0 , PRA1 <= -(8) , 0 + ((0.75 - 0) * (PRA1 - -(8))) / (-(6) - -(8)) , ( PRA1 > -(8)) && (PRA1 <= -(6) ), 0.75 + ((2.6 - 0.75) * (PRA1 - -(6))) / (-(2) - -(6)) , ( PRA1 > -(6)) && (PRA1 <= -(2) ), 2.6 + ((9.8 - 2.6) * (PRA1 - -(2))) / (4 - -(2)) , ( PRA1 > -(2)) && (PRA1 <= 4 ), 9.8 + ((13.5 - 9.8) * (PRA1 - 4)) / (12 - 4) , ( PRA1 > 4) && (PRA1 <= 12 ), 13.5 ); QPO := PGL / RPT; VPA := VPA1 + VBD * 0.155; pulmonary_vasculature_blood_volume.DPA := QRO - QPO; pulmonary_vasculature_pressure.VPE := VPA - 0.30625; pulmonary_arterial_resistance.PP1T := 0.026 * PPA; pulmonary_arterial_resistance.PP1 := piecewise( 0.00001 , pulmonary_arterial_resistance.PP1T < 0.00001 , pulmonary_arterial_resistance.PP1T ); pulmonary_arterial_resistance.CPA := power(pulmonary_arterial_resistance.PP1, 0.5); pulmonary_venous_resistance.PL1 := PLA + 18; RPT := RPV + RPA; PGL := PPA - PLA; VLA := VLA1 + VBD * 0.128; left_atrial_blood_volume.DLA := QPO - QLO; left_atrial_pressure.VLE := VLA - 0.38; PLA1 := (PLA + 4) * (HTAUML * (AU - 1) + 1) - 4; LVM := piecewise( 1.04 , pumping_effectiveness_of_left_ventricle.PA2 <= 0 , 1.04 + ((1.025 - 1.04) * (pumping_effectiveness_of_left_ventricle.PA2 - 0)) / (60 - 0) , ( pumping_effectiveness_of_left_ventricle.PA2 > 0) && (pumping_effectiveness_of_left_ventricle.PA2 <= 60 ), 1.025 + ((0.97 - 1.025) * (pumping_effectiveness_of_left_ventricle.PA2 - 60)) / (125 - 60) , ( pumping_effectiveness_of_left_ventricle.PA2 > 60) && (pumping_effectiveness_of_left_ventricle.PA2 <= 125 ), 0.97 + ((0.88 - 0.97) * (pumping_effectiveness_of_left_ventricle.PA2 - 125)) / (160 - 125) , ( pumping_effectiveness_of_left_ventricle.PA2 > 125) && (pumping_effectiveness_of_left_ventricle.PA2 <= 160 ), 0.88 + ((0.59 - 0.88) * (pumping_effectiveness_of_left_ventricle.PA2 - 160)) / (200 - 160) , ( pumping_effectiveness_of_left_ventricle.PA2 > 160) && (pumping_effectiveness_of_left_ventricle.PA2 <= 200 ), 0.59 + ((0 - 0.59) * (pumping_effectiveness_of_left_ventricle.PA2 - 200)) / (240 - 200) , ( pumping_effectiveness_of_left_ventricle.PA2 > 200) && (pumping_effectiveness_of_left_ventricle.PA2 <= 240 ), 0 ); pumping_effectiveness_of_left_ventricle.PA2 := PA / (AUH * OSA); QLOT := LVM * QLN * AUH * HSL * HMD * HPL; left_ventricular_output.QLO1 := (PLA - PA) / 3; VVS := VVS1 + VBD * 0.3986; venous_blood_volume.DVS := QAO - QVO; VVA := (ANU - 1) * ANY; venous_excess_volume.VVE1 := ((((VVS - VVR) - VVA) - VV7) - VV6) - ATRVFB; venous_average_pressure.PVS1 := 3.7 + (VVE - 0.74) / CV; PR1 := piecewise( PR1LL , PRA < PR1LL , PRA ); RVG := 0.74 / power(PVS / (VIM * 3.7), 0.5); rate_of_blood_flow_from_veins_to_right_atrium.PGV := PVS - PR1; venous_resistance.CN3 := ((PC - 17) * CN7 + 17) * CN2; venous_resistance.RV1 := RVSM / venous_resistance.CN3; NNRVR := RVS * 1.79; VAS := VAS1 + VBD * 0.261; arterial_blood_volume.DAS := QLO - QAO; PAG := PA - PRA; arterial_pressure_and_pressure_gradient.VAE := VAS - 0.495; PAM := power(PA / 100, PAEX); R1 := ((ANU * ADHMV * AUM * VIM * PAMK) / PAM) / ATRRFB; NNRAR := RAR * ARM * R1 * MYOGRS * RMULT1; PGS := PA - PVS; RSM := RAM * AMM * R1 * MYOGRS * RMULT1; RSN := NNRAR + NNRVR; FISFLO := PAG * FIS; SYSFLO := BFM + BFN + RBF; RTP := PAG / QAO; // Rate Rules: VVS1' = venous_blood_volume.DVS; VAS1' = arterial_blood_volume.DAS; VLA1' = left_atrial_blood_volume.DLA; VPA1' = pulmonary_vasculature_blood_volume.DPA; VRA1' = right_atrial_blood_volume.DRA; // Variable initializations: time_ = ; ADHMV = ; AMM = ; ANU = ; ANUVN = ; ARM = ; ATRRFB = ; ATRVFB = ; AU = ; AUH = ; AUM = ; AVE = ; HMD = ; HPL = ; HPR = ; MYOGRS = ; OSA = ; PAMK = ; PC = ; RBF = ; VIM = ; VP = ; VRC = ; VV6 = ; VV7 = ; VVR = ; VRA1 = 0.100043; HTAUML = 0.4; QRF = 0.15; HSR = 1; VPA1 = 0.38131; VLA1 = 0.379883; HSL = 1; VVS1 = 3.28246; ANY = -0.2; CV = 0.1; PR1LL = 0; CN7 = 0.2; CN2 = 0.0212; RVSM = 1; VAS1 = 0.862514; PAEX = 2; RAR = 30.52; RMULT1 = 1; RAM = 96.3; FIS = 0; end model Guyton_Model__electrolytes_parameter_values(TRPL, NID, ALCLK, KID, VIDML) // Variable initializations: TRPL = 0; NID = 0.1; ALCLK = 0.3; KID = 0.08; VIDML = 0.01; end model Guyton_Model__extracellular_Na_concentration(STH, NOD, VEC, time_, CNA, NID, TRPL) // Variable initializations: STH = ; NOD = ; VEC = ; time_ = ; NID = ; TRPL = ; NAE = 2109.91; end model Guyton_Model__aldosterone_effect_on_cellular_K_distribution(AMK, AMK1, ALCLK) // Variable initializations: AMK = ; ALCLK = ; end model Guyton_Model__extracellular_K_concentration(KOD, AMK1, VEC, time_, CKE, KE, KTOT, KID) // Variable initializations: KOD = ; AMK1 = ; VEC = ; time_ = ; KTOT = 3622.54; KID = ; end model Guyton_Model__intracellular_K_concentration(KTOT, KE, VIC, CKI) // Variable initializations: KTOT = ; KE = ; VIC = ; end model Guyton_Model__intracellular_fluid_volume(CNA, CKI, time_, VIC, VID, VIDML) // Variable initializations: CNA = ; CKI = ; time_ = ; VIC = 25.0404; VIDML = ; end model Guyton_Model__total_body_water(TVD, VUD, time_, VTW) // Variable initializations: TVD = ; VUD = ; time_ = ; VTW = 39.8952; end model Guyton_Model__extracellular_fluid_volume(VIC, VTW, VEC) // Variable initializations: VIC = ; VTW = ; end model Guyton_Model__electrolytes(time_, AMK, TVD, NOD, STH, KOD, VUD, CKE, VEC, CNA, VTW) // Sub-modules, and any changes to those submodules: electrolytes_parameter_values: Guyton_Model__electrolytes_parameter_values(TRPL, NID, ALCLK, KID, VIDML); extracellular_Na_concentration: Guyton_Model__extracellular_Na_concentration(STH, NOD, VEC, time_, CNA, NID, TRPL); aldosterone_effect_on_cellular_K_distribution: Guyton_Model__aldosterone_effect_on_cellular_K_distribution(AMK, AMK1, ALCLK); extracellular_K_concentration: Guyton_Model__extracellular_K_concentration(KOD, AMK1, VEC, time_, CKE, KE, KTOT, KID); intracellular_K_concentration: Guyton_Model__intracellular_K_concentration(KTOT, KE, VIC, CKI); intracellular_fluid_volume: Guyton_Model__intracellular_fluid_volume(CNA, CKI, time_, VIC, VID, VIDML); total_body_water: Guyton_Model__total_body_water(TVD, VUD, time_, VTW); extracellular_fluid_volume: Guyton_Model__extracellular_fluid_volume(VIC, VTW, VEC); // Variable initializations: time_ = ; AMK = ; TVD = ; NOD = ; STH = ; KOD = ; VUD = ; VTW = 39.8952; TRPL = 0; NID = 0.1; ALCLK = 0.3; KID = 0.08; KTOT = 3622.54; VIC = 25.0404; VIDML = 0.01; end model Guyton_Model__heart_hypertrophy_or_deterioration_parameter_values(HSL, Z13, HSR, DHDTR) // Variable initializations: HSL = 1; Z13 = 0.625; HSR = 1; DHDTR = 0.05; end model Guyton_Model__left_ventricular_hypertrophy(QAO, PA, time_, HPL, HSL, Z13) // Variable initializations: QAO = ; PA = ; time_ = ; HPL = 1.00163; HSL = ; Z13 = ; end model Guyton_Model__right_ventricular_hypertrophy(QAO, PPA, time_, HPR, HSR, Z13) // Variable initializations: QAO = ; PPA = ; time_ = ; HPR = 1.00237; HSR = ; Z13 = ; end model Guyton_Model__heart_deterioration(POT, time_, HMD, DHDTR) // Variable initializations: POT = ; time_ = ; DHDTR = ; HMD1 = 1.0; end model Guyton_Model__heart_hypertrophy_or_deterioration(time_, QAO, PA, POT, PPA, HPL, HPR, HMD) // Sub-modules, and any changes to those submodules: heart_hypertrophy_or_deterioration_parameter_values: Guyton_Model__heart_hypertrophy_or_deterioration_parameter_values(HSL, Z13, HSR, DHDTR); left_ventricular_hypertrophy: Guyton_Model__left_ventricular_hypertrophy(QAO, PA, time_, HPL, HSL, Z13); right_ventricular_hypertrophy: Guyton_Model__right_ventricular_hypertrophy(QAO, PPA, time_, HPR, HSR, Z13); heart_deterioration: Guyton_Model__heart_deterioration(POT, time_, HMD, DHDTR); // Variable initializations: time_ = ; QAO = ; PA = ; POT = ; PPA = ; HPL = 1.00163; HPR = 1.00237; HSL = 1; Z13 = 0.625; HSR = 1; DHDTR = 0.05; end model Guyton_Model__HR_and_SV_parameter_values(PR1LL) // Variable initializations: PR1LL = 0; end model Guyton_Model__effect_of_autonomic_stimulation_on_HR(AUR, AUHR) // Variable initializations: AUR = ; end model Guyton_Model__effect_of_PRA_on_HR(PRHR, PR1LL) // Variable initializations: PR1LL = ; end model Guyton_Model__effect_of_heart_deterioration_on_HR(HMD, HDHR) // Variable initializations: HMD = ; end model Guyton_Model__heart_rate(AUHR, PRHR, HDHR, HR) // Variable initializations: AUHR = ; PRHR = ; HDHR = ; end model Guyton_Model__stroke_volume_output(QLO, HR, SVO) // Variable initializations: QLO = ; HR = ; end model Guyton_Model__heart_rate_and_stroke_volume(QLO, AUR, PRA, HMD) // Sub-modules, and any changes to those submodules: HR_and_SV_parameter_values: Guyton_Model__HR_and_SV_parameter_values(PR1LL); effect_of_autonomic_stimulation_on_HR: Guyton_Model__effect_of_autonomic_stimulation_on_HR(AUR, AUHR); effect_of_PRA_on_HR: Guyton_Model__effect_of_PRA_on_HR(PRHR, PR1LL); effect_of_heart_deterioration_on_HR: Guyton_Model__effect_of_heart_deterioration_on_HR(HMD, HDHR); heart_rate: Guyton_Model__heart_rate(AUHR, PRHR, HDHR, HR); stroke_volume_output: Guyton_Model__stroke_volume_output(QLO, HR, SVO); // Variable initializations: QLO = ; AUR = ; PRA = ; HMD = ; PR1LL = 0; end model Guyton_Model__M_autoregulatory_local_blood_flow_parameter_values(POM, A4K, AMM4, POM2, A4K2) // Variable initializations: POM = 0.04; A4K = 0.1; AMM4 = 0.005; POM2 = 2; A4K2 = 40000; end model Guyton_Model__M_autoregulatory_driving_force(PMO, PDO) // Variable initializations: PMO = ; end model Guyton_Model__M_short_term_autoregulation() end model Guyton_Model__M_ST_sensitivity_control(PDO, POE, POM) // Variable initializations: PDO = ; POM = ; end model Guyton_Model__M_ST_time_delay_and_limit(POE, time_, AMM1, A4K, AMM4) // Variable initializations: POE = ; time_ = ; A4K = ; AMM4 = ; AMM1T = 1.00269; end model Guyton_Model__M_long_term_autoregulation() end model Guyton_Model__M_LT_sensitivity_control(PDO, POF, POM2) // Variable initializations: PDO = ; POM2 = ; end model Guyton_Model__M_LT_time_delay(POF, time_, AMM2, A4K2) // Variable initializations: POF = ; time_ = ; AMM2 = 1.09071; A4K2 = ; end model Guyton_Model__global_M_blood_flow_autoregulation_output(AMM1, AMM2, AMM) // Variable initializations: AMM1 = ; AMM2 = ; end model Guyton_Model__muscle_autoregulatory_local_blood_flow_control(time_, PMO, AMM) // Sub-modules, and any changes to those submodules: M_autoregulatory_local_blood_flow_parameter_values: Guyton_Model__M_autoregulatory_local_blood_flow_parameter_values(POM, A4K, AMM4, POM2, A4K2); M_autoregulatory_driving_force: Guyton_Model__M_autoregulatory_driving_force(PMO, PDO); M_short_term_autoregulation: Guyton_Model__M_short_term_autoregulation(); M_ST_sensitivity_control: Guyton_Model__M_ST_sensitivity_control(PDO, POE, POM); M_ST_time_delay_and_limit: Guyton_Model__M_ST_time_delay_and_limit(POE, time_, AMM1, A4K, AMM4); M_long_term_autoregulation: Guyton_Model__M_long_term_autoregulation(); M_LT_sensitivity_control: Guyton_Model__M_LT_sensitivity_control(PDO, POF, POM2); M_LT_time_delay: Guyton_Model__M_LT_time_delay(POF, time_, AMM2, A4K2); global_M_blood_flow_autoregulation_output: Guyton_Model__global_M_blood_flow_autoregulation_output(AMM1, AMM2, AMM); // Variable initializations: time_ = ; PMO = ; POM = 0.04; A4K = 0.1; AMM4 = 0.005; POM2 = 2; A4K2 = 40000; AMM2 = 1.09071; end model Guyton_Model__M_O2_delivery_parameter_values(EXC, EXCXP2, OMM, PM5, PK2) // Variable initializations: EXC = 1; EXCXP2 = 0.17; OMM = 57.1; PM5 = 30; PK2 = 0.79167; end model Guyton_Model__M_O2_blood_supply(OVA, BFM, O2ARTM) // Variable initializations: OVA = ; BFM = ; end model Guyton_Model__M_venous_O2_content(O2ARTM, RMO, BFM, HM, time_, PVO, EXC, EXCXP2) // Variable initializations: O2ARTM = ; RMO = ; BFM = ; HM = ; time_ = ; EXC = ; EXCXP2 = ; end model Guyton_Model__metabolic_O2_consumption_by_M_tissue(PMO, AOM, MMO, P2O, OMM, EXC) // Variable initializations: PMO = ; AOM = ; OMM = ; EXC = ; end model Guyton_Model__delivery_of_O2_to_M_tissues(PMO, PVO, BFM, RMO, PM5) // Variable initializations: PMO = ; PVO = ; BFM = ; PM5 = ; end model Guyton_Model__volume_of_O2_in_M_tissue(RMO, MMO, time_, QOM) // Variable initializations: RMO = ; MMO = ; time_ = ; QOM1 = 48.0839; end model Guyton_Model__pressure_of_O2_in_M_tissue_cells(QOM, PMO, PK2) // Variable initializations: QOM = ; PK2 = ; end model Guyton_Model__muscle_O2_delivery(time_, BFM, OVA, HM, AOM, PMO, RMO) // Sub-modules, and any changes to those submodules: M_O2_delivery_parameter_values: Guyton_Model__M_O2_delivery_parameter_values(EXC, EXCXP2, OMM, PM5, PK2); M_O2_blood_supply: Guyton_Model__M_O2_blood_supply(OVA, BFM, O2ARTM); M_venous_O2_content: Guyton_Model__M_venous_O2_content(O2ARTM, RMO0, BFM, HM, time_, PVO, EXC, EXCXP2); metabolic_O2_consumption_by_M_tissue: Guyton_Model__metabolic_O2_consumption_by_M_tissue(PMO, AOM, MMO, P2O, OMM, EXC); delivery_of_O2_to_M_tissues: Guyton_Model__delivery_of_O2_to_M_tissues(PMO, PVO0, BFM, RMO0, PM5); volume_of_O2_in_M_tissue: Guyton_Model__volume_of_O2_in_M_tissue(RMO, MMO, time_, QOM); pressure_of_O2_in_M_tissue_cells: Guyton_Model__pressure_of_O2_in_M_tissue_cells(QOM, PMO, PK2); RMO0 is RMO; // Variable initializations: time_ = ; BFM = ; OVA = ; HM = ; AOM = ; EXC = 1; EXCXP2 = 0.17; OMM = 57.1; PM5 = 30; PK2 = 0.79167; PVO0 = ; end model Guyton_Model__NM_autoregulatory_local_blood_flow_parameter_values(POR, POK, A1K, PON, A2K, A3K, POZ, AUTOSN) // Variable initializations: POR = 35; POK = 0.1; A1K = 0.5; PON = 0.1; A2K = 60; A3K = 40000; POZ = 2; AUTOSN = 0.9; end model Guyton_Model__NM_autoregulatory_driving_force(POT, POD, POR) // Variable initializations: POT = ; POR = ; end model Guyton_Model__NM_short_term_autoregulation() end model Guyton_Model__NM_ST_sensitivity_control(POD, POB, POK) // Variable initializations: POD = ; POK = ; end model Guyton_Model__NM_ST_time_delay_and_damping(POB, time_, AR1, A1K) // Variable initializations: POB = ; time_ = ; A1K = ; AR1T = 1.02127; end model Guyton_Model__NM_intermediate_autoregulation() end model Guyton_Model__NM_I_sensitivity_control(POD, POA, PON) // Variable initializations: POD = ; PON = ; end model Guyton_Model__NM_I_time_delay_and_limit(POA, time_, AR2, A2K) // Variable initializations: POA = ; time_ = ; A2K = ; AR2T = 1.01179; end model Guyton_Model__NM_long_term_autoregulation() end model Guyton_Model__NM_LT_sensitivity_control(POD, POC, POZ) // Variable initializations: POD = ; POZ = ; end model Guyton_Model__NM_LT_time_delay_and_limit(POC, time_, AR3, A3K) // Variable initializations: POC = ; time_ = ; A3K = ; AR3T = 1.1448; end model Guyton_Model__total_NM_autoregulation(AR1, AR2, AR3, ARM1) // Variable initializations: AR1 = ; AR2 = ; AR3 = ; end model Guyton_Model__global_NM_blood_flow_autoregulation_output(ARM1, ARM, AUTOSN) // Variable initializations: ARM1 = ; AUTOSN = ; end model Guyton_Model__non_muscle_autoregulatory_local_blood_flow_control(time_, POT, ARM) // Sub-modules, and any changes to those submodules: NM_autoregulatory_local_blood_flow_parameter_values: Guyton_Model__NM_autoregulatory_local_blood_flow_parameter_values(POR, POK, A1K, PON, A2K, A3K, POZ, AUTOSN); NM_autoregulatory_driving_force: Guyton_Model__NM_autoregulatory_driving_force(POT, POD, POR); NM_short_term_autoregulation: Guyton_Model__NM_short_term_autoregulation(); NM_ST_sensitivity_control: Guyton_Model__NM_ST_sensitivity_control(POD, POB, POK); NM_ST_time_delay_and_damping: Guyton_Model__NM_ST_time_delay_and_damping(POB, time_, AR1, A1K); NM_intermediate_autoregulation: Guyton_Model__NM_intermediate_autoregulation(); NM_I_sensitivity_control: Guyton_Model__NM_I_sensitivity_control(POD, POA, PON); NM_I_time_delay_and_limit: Guyton_Model__NM_I_time_delay_and_limit(POA, time_, AR2, A2K); NM_long_term_autoregulation: Guyton_Model__NM_long_term_autoregulation(); NM_LT_sensitivity_control: Guyton_Model__NM_LT_sensitivity_control(POD, POC, POZ); NM_LT_time_delay_and_limit: Guyton_Model__NM_LT_time_delay_and_limit(POC, time_, AR3, A3K); total_NM_autoregulation: Guyton_Model__total_NM_autoregulation(AR1, AR2, AR3, ARM1); global_NM_blood_flow_autoregulation_output: Guyton_Model__global_NM_blood_flow_autoregulation_output(ARM1, ARM, AUTOSN); // Variable initializations: time_ = ; POT = ; POR = 35; POK = 0.1; A1K = 0.5; PON = 0.1; A2K = 60; POZ = 2; A3K = 40000; AUTOSN = 0.9; end model Guyton_Model__NM_O2_delivery_parameter_values(O2M) // Variable initializations: O2M = 164; end model Guyton_Model__NM_O2_blood_supply(OVA, BFN, O2ARTN) // Variable initializations: OVA = ; BFN = ; end model Guyton_Model__NM_venous_O2_content(O2ARTN, DOB, BFN, HM, time_, POV) // Variable initializations: O2ARTN = ; DOB = ; BFN = ; HM = ; time_ = ; end model Guyton_Model__O2_consumption_by_NM_tissue(POT, AOM, MO2, O2M) // Variable initializations: POT = ; AOM = ; O2M = ; end model Guyton_Model__delivery_of_O2_to_NM_tissues(POT, POV, BFN, DOB) // Variable initializations: POT = ; POV = ; BFN = ; end model Guyton_Model__volume_of_O2_in_NM_tissue(DOB, MO2, time_, QO2) // Variable initializations: DOB = ; MO2 = ; time_ = ; QO2T = 72.2362; end model Guyton_Model__pressure_of_O2_in_NM_tissue_cells(QO2, POT) // Variable initializations: QO2 = ; end model Guyton_Model__non_muscle_O2_delivery(time_, BFN, OVA, HM, AOM, POT, DOB) // Sub-modules, and any changes to those submodules: NM_O2_delivery_parameter_values: Guyton_Model__NM_O2_delivery_parameter_values(O2M); NM_O2_blood_supply: Guyton_Model__NM_O2_blood_supply(OVA, BFN, O2ARTN); NM_venous_O2_content: Guyton_Model__NM_venous_O2_content(O2ARTN, DOB0, BFN, HM, time_, POV); O2_consumption_by_NM_tissue: Guyton_Model__O2_consumption_by_NM_tissue(POT, AOM, MO2, O2M); delivery_of_O2_to_NM_tissues: Guyton_Model__delivery_of_O2_to_NM_tissues(POT, POV0, BFN, DOB0); volume_of_O2_in_NM_tissue: Guyton_Model__volume_of_O2_in_NM_tissue(DOB, MO2, time_, QO2); pressure_of_O2_in_NM_tissue_cells: Guyton_Model__pressure_of_O2_in_NM_tissue_cells(QO2, POT); DOB0 is DOB; // Variable initializations: time_ = ; BFN = ; OVA = ; HM = ; AOM = ; O2M = 164; POV0 = ; end model Guyton_Model__pulmonary_fluid_dynamics_parameter_values(CPF) // Variable initializations: CPF = 0.0003; end model Guyton_Model__pulmonary_capillary_pressure(PPA, PLA, RPV, RPA, PCP) // Variable initializations: PPA = ; PLA = ; RPV = ; RPA = ; end model Guyton_Model__fluid_filtration_into_pulmonary_interstitium(PCP, PPC, POS, PPI, PFI, CPF) // Variable initializations: PCP = ; PPC = ; POS = ; PPI = ; CPF = ; end model Guyton_Model__pulmonary_interstitial_free_fluid_volume(PFI, PLF, time_, DFP, VPF) // Variable initializations: PFI = ; PLF = ; time_ = ; VPF1 = 0.0123238; end model Guyton_Model__pulmonary_interstitial_fluid_pressure(VPF, PPI) // Variable initializations: VPF = ; end model Guyton_Model__concentration_of_protein_in_pulmonary_interstitium(PPO, PPN, VPF, time_, PPD, CPN) // Variable initializations: PPO = ; PPN = ; VPF = ; time_ = ; PPR1 = 0.419998; end model Guyton_Model__colloid_osmotic_pressure_of_pulmonary_interstitium(CPN, POS) // Variable initializations: CPN = ; end model Guyton_Model__protein_leakage_into_pulmonary_interstitium(CPP, CPN, PPN) // Variable initializations: CPP = ; CPN = ; end model Guyton_Model__lung_lymphatic_protein_flow(PPI, CPN, PLF, PPO) // Variable initializations: PPI = ; CPN = ; end model Guyton_Model__pulmonary_fluid_dynamics(time_, PPC, PPA, PLA, CPP, RPV, RPA, DFP, VPF, PPD) // Sub-modules, and any changes to those submodules: pulmonary_fluid_dynamics_parameter_values: Guyton_Model__pulmonary_fluid_dynamics_parameter_values(CPF); pulmonary_capillary_pressure: Guyton_Model__pulmonary_capillary_pressure(PPA, PLA, RPV, RPA, PCP); fluid_filtration_into_pulmonary_interstitium: Guyton_Model__fluid_filtration_into_pulmonary_interstitium(PCP, PPC, POS, PPI, PFI, CPF); pulmonary_interstitial_free_fluid_volume: Guyton_Model__pulmonary_interstitial_free_fluid_volume(PFI, PLF, time_, DFP, VPF); pulmonary_interstitial_fluid_pressure: Guyton_Model__pulmonary_interstitial_fluid_pressure(VPF, PPI); concentration_of_protein_in_pulmonary_interstitium: Guyton_Model__concentration_of_protein_in_pulmonary_interstitium(PPO, PPN, VPF, time_, PPD, CPN); colloid_osmotic_pressure_of_pulmonary_interstitium: Guyton_Model__colloid_osmotic_pressure_of_pulmonary_interstitium(CPN, POS); protein_leakage_into_pulmonary_interstitium: Guyton_Model__protein_leakage_into_pulmonary_interstitium(CPP, CPN, PPN); lung_lymphatic_protein_flow: Guyton_Model__lung_lymphatic_protein_flow(PPI, CPN, PLF, PPO); // Variable initializations: time_ = ; PPC = ; PPA = ; PLA = ; CPP = ; RPV = ; RPA = ; CPF = 0.0003; end model Guyton_Model__pulmonary_O2_uptake_parameter_values(PO2AMB, PL2, VPTISS, VNTSTM) // Variable initializations: PO2AMB = 150; PL2 = 1.8; VPTISS = 0.0175; VNTSTM = 1; end model Guyton_Model__total_O2_utilization(DOB, RMO, O2UTIL) // Variable initializations: DOB = ; RMO = ; end model Guyton_Model__alveolar_ventilation(O2UTIL, O2VAD2, O2VTS2, ALVENT, VNTSTM) // Variable initializations: O2UTIL = ; O2VAD2 = ; O2VTS2 = ; VNTSTM = ; end model Guyton_Model__alveolar_PO2(O2UTIL, ALVENT, PO2ALV, PO2AMB) // Variable initializations: O2UTIL = ; ALVENT = ; PO2AMB = ; end model Guyton_Model__respiratory_O2_diffusion_into_capillaries(PO2ALV, PO2ART, VPF, O2DFS, PL2, VPTISS) // Variable initializations: PO2ALV = ; PO2ART = ; VPF = ; PL2 = ; VPTISS = ; end model Guyton_Model__O2_volume_of_arterial_blood(O2DFS, O2UTIL, QRO, time_, OVA) // Variable initializations: O2DFS = ; O2UTIL = ; QRO = ; time_ = ; OVA = 204.497; end model Guyton_Model__arterial_PO2(OVA, HM, PO2ART, OSA) // Variable initializations: OVA = ; HM = ; end model Guyton_Model__chemoreceptor_adaptation_of_alveolar_ventilation() end model Guyton_Model__acute_chemoreceptor_adaptation_of_alveolar_ventilation(PO2ART, O2VTS2) // Variable initializations: PO2ART = ; end model Guyton_Model__progressive_chemoreceptor_adaptation_of_alveolar_ventilation(O2VTS2, time_, O2VAD2) // Variable initializations: O2VTS2 = ; time_ = ; O2VAD1 = 2.368e-07; end model Guyton_Model__pulmonary_O2_uptake(time_, VPF, DOB, QRO, RMO, HM, OVA, PO2ART, OSA) // Sub-modules, and any changes to those submodules: pulmonary_O2_uptake_parameter_values: Guyton_Model__pulmonary_O2_uptake_parameter_values(PO2AMB, PL2, VPTISS, VNTSTM); total_O2_utilization: Guyton_Model__total_O2_utilization(DOB, RMO, O2UTIL); alveolar_ventilation: Guyton_Model__alveolar_ventilation(O2UTIL, O2VAD2, O2VTS2, ALVENT, VNTSTM); alveolar_PO2: Guyton_Model__alveolar_PO2(O2UTIL, ALVENT, PO2ALV, PO2AMB); respiratory_O2_diffusion_into_capillaries: Guyton_Model__respiratory_O2_diffusion_into_capillaries(PO2ALV, PO2ART0, VPF, O2DFS, PL2, VPTISS); O2_volume_of_arterial_blood: Guyton_Model__O2_volume_of_arterial_blood(O2DFS, O2UTIL, QRO, time_, OVA); arterial_PO2: Guyton_Model__arterial_PO2(OVA, HM, PO2ART0, OSA); chemoreceptor_adaptation_of_alveolar_ventilation: Guyton_Model__chemoreceptor_adaptation_of_alveolar_ventilation(); acute_chemoreceptor_adaptation_of_alveolar_ventilation: Guyton_Model__acute_chemoreceptor_adaptation_of_alveolar_ventilation(PO2ART, O2VTS2); progressive_chemoreceptor_adaptation_of_alveolar_ventilation: Guyton_Model__progressive_chemoreceptor_adaptation_of_alveolar_ventilation(O2VTS2, time_, O2VAD2); PO2ART0 is PO2ART; // Variable initializations: time_ = ; VPF = ; DOB = ; QRO = ; RMO = ; HM = ; OVA = 204.497; VNTSTM = 1; PO2AMB = 150; PL2 = 1.8; VPTISS = 0.0175; end model Guyton_Model__red_cells_and_viscosity_parameter_values(HMK, HKM, PO2AMB, HM6, HM8, REK, RKC, TRRBC) // Variable initializations: HMK = 90; HKM = 0.53333; PO2AMB = 150; HM6 = 1850; HM8 = 4.714e-08; REK = 1; RKC = 5.8e-06; TRRBC = 0; end model Guyton_Model__blood_viscosity_calculations() end model Guyton_Model__hematocrit_fraction(VP, VRC, HM, HM1) // Variable initializations: VP = ; VRC = ; end model Guyton_Model__viscosity_due_to_RBCs(HM, VIE, HMK, HKM) // Variable initializations: HM = ; HMK = ; HKM = ; end model Guyton_Model__blood_viscosity(VIE, VIM) // Variable initializations: VIE = ; end model Guyton_Model__RBC_formation_and_destruction() end model Guyton_Model__oxygen_stimulation(HM, HM7, PO2AMB, HM6) // Variable initializations: HM = ; PO2AMB = ; HM6 = ; end model Guyton_Model__RBC_production(HM7, RC1, HM8, REK) // Variable initializations: HM7 = ; HM8 = ; REK = ; end model Guyton_Model__RBC_destruction(VRC, VIM, RC2, RKC) // Variable initializations: VRC = ; VIM = ; RKC = ; end model Guyton_Model__RBC_volume(RC1, RC2, time_, VRC, TRRBC) // Variable initializations: RC1 = ; RC2 = ; time_ = ; VRC = 2.00439; TRRBC = ; end model Guyton_Model__red_cells_and_viscosity(time_, VP, VIM, HM, HM1, VRC) // Sub-modules, and any changes to those submodules: red_cells_and_viscosity_parameter_values: Guyton_Model__red_cells_and_viscosity_parameter_values(HMK, HKM, PO2AMB, HM6, HM8, REK, RKC, TRRBC); blood_viscosity_calculations: Guyton_Model__blood_viscosity_calculations(); hematocrit_fraction: Guyton_Model__hematocrit_fraction(VP, VRC0, HM, HM1); viscosity_due_to_RBCs: Guyton_Model__viscosity_due_to_RBCs(HM, VIE, HMK, HKM); blood_viscosity: Guyton_Model__blood_viscosity(VIE, VIM); RBC_formation_and_destruction: Guyton_Model__RBC_formation_and_destruction(); oxygen_stimulation: Guyton_Model__oxygen_stimulation(HM, HM7, PO2AMB, HM6); RBC_production: Guyton_Model__RBC_production(HM7, RC1, HM8, REK); RBC_destruction: Guyton_Model__RBC_destruction(VRC0, VIM, RC2, RKC); RBC_volume: Guyton_Model__RBC_volume(RC1, RC2, time_, VRC0, TRRBC); VRC0 is VRC; // Variable initializations: time_ = ; VP = ; VRC = 2.00439; HMK = 90; HKM = 0.53333; PO2AMB = 150; HM6 = 1850; HM8 = 4.714e-08; REK = 1; RKC = 5.8e-06; TRRBC = 0; end model Guyton_Model__stress_relaxation_parameter_values(SR, SR2, SRK, SRK2) // Variable initializations: SR = 1; SR2 = 1; SRK = 5; SRK2 = 10000; end model Guyton_Model__short_term_stress_relaxation(VVE, time_, VV7, SR, SRK) // Variable initializations: VVE = ; time_ = ; VV7 = 0.00366525; SR = ; SRK = ; end model Guyton_Model__long_term_stress_relaxation(VVE, time_, VV6, SR2, SRK2) // Variable initializations: VVE = ; time_ = ; VV6 = 0.0101913; SR2 = ; SRK2 = ; end model Guyton_Model__stress_relaxation(time_, VVE, VV6, VV7) // Sub-modules, and any changes to those submodules: stress_relaxation_parameter_values: Guyton_Model__stress_relaxation_parameter_values(SR, SR2, SRK, SRK2); short_term_stress_relaxation: Guyton_Model__short_term_stress_relaxation(VVE, time_, VV7, SR, SRK); long_term_stress_relaxation: Guyton_Model__long_term_stress_relaxation(VVE, time_, VV6, SR2, SRK2); // Variable initializations: time_ = ; VVE = ; VV6 = 0.0101913; VV7 = 0.00366525; SR = 1; SRK = 5; SR2 = 1; SRK2 = 10000; end model Guyton_Model__thirst_drinking_and_salt_appetite_parameter_values(Z10, Z11, ANMSLT, AHTHM, ANMTM, DR, TVDDL) // Variable initializations: Z10 = 45; Z11 = 0.01; ANMSLT = 2; AHTHM = 2; ANMTM = 1.5; DR = 0; TVDDL = 30; end model Guyton_Model__effect_of_salt_appetite_stimulation_on_thirst(ANM, POT, STH, ANMSLT, Z10, Z11) // Variable initializations: ANM = ; POT = ; ANMSLT = ; Z10 = ; Z11 = ; end model Guyton_Model__effect_of_antidiuretic_hormone_on_thirst(ADHC, AHCM, AHTHM) // Variable initializations: ADHC = ; AHTHM = ; end model Guyton_Model__effect_of_angiotensin_on_thirst(ANM, ANMTH, ANMTM) // Variable initializations: ANM = ; ANMTM = ; end model Guyton_Model__rate_of_fluid_intake(STH, AHCM, ANMTH, time_, TVD, DR, TVDDL) // Variable initializations: STH = ; AHCM = ; ANMTH = ; time_ = ; TVD = 0.000980838; DR = ; TVDDL = ; end model Guyton_Model__thirst_drinking_and_salt_appetite(time_, ADHC, ANM, POT, TVD, STH) // Sub-modules, and any changes to those submodules: thirst_drinking_and_salt_appetite_parameter_values: Guyton_Model__thirst_drinking_and_salt_appetite_parameter_values(Z10, Z11, ANMSLT, AHTHM, ANMTM, DR, TVDDL); effect_of_salt_appetite_stimulation_on_thirst: Guyton_Model__effect_of_salt_appetite_stimulation_on_thirst(ANM, POT, STH, ANMSLT, Z10, Z11); effect_of_antidiuretic_hormone_on_thirst: Guyton_Model__effect_of_antidiuretic_hormone_on_thirst(ADHC, AHCM, AHTHM); effect_of_angiotensin_on_thirst: Guyton_Model__effect_of_angiotensin_on_thirst(ANM, ANMTH, ANMTM); rate_of_fluid_intake: Guyton_Model__rate_of_fluid_intake(STH, AHCM, ANMTH, time_, TVD, DR, TVDDL); // Variable initializations: time_ = ; ADHC = ; ANM = ; POT = ; TVD = 0.000980838; ANMSLT = 2; Z10 = 45; Z11 = 0.01; AHTHM = 2; ANMTM = 1.5; DR = 0; TVDDL = 30; end model Guyton_Model__volume_receptors_parameter_values(AH10, AH11, AH9, ATRFBM, ATRVM) // Variable initializations: AH10 = 0.333; AH11 = 1000; AH9 = 1; ATRFBM = 0; ATRVM = 0; end model Guyton_Model__effect_of_pressure_on_volume_receptors(PRA, AHZ, AH10, AH9) // Variable initializations: PRA = ; AH10 = ; AH9 = ; end model Guyton_Model__time_dependent_volume_receptor_adaptation(AHZ, time_, AHY, AH11) // Variable initializations: AHZ = ; time_ = ; AHY = 0.301963; AH11 = ; end model Guyton_Model__total_volume_nervous_feedback(AHZ, AHY, AH7) // Variable initializations: AHZ = ; AHY = ; end model Guyton_Model__volume_effect_on_arteries(AH7, ATRRFB, ATRFBM) // Variable initializations: AH7 = ; ATRFBM = ; end model Guyton_Model__volume_effect_on_unstressed_venous_volume(AH7, ATRVFB, ATRVM) // Variable initializations: AH7 = ; ATRVM = ; end model Guyton_Model__volume_receptors(time_, PRA, ATRRFB, ATRVFB) // Sub-modules, and any changes to those submodules: volume_receptors_parameter_values: Guyton_Model__volume_receptors_parameter_values(AH10, AH11, AH9, ATRFBM, ATRVM); effect_of_pressure_on_volume_receptors: Guyton_Model__effect_of_pressure_on_volume_receptors(PRA, AHZ, AH10, AH9); time_dependent_volume_receptor_adaptation: Guyton_Model__time_dependent_volume_receptor_adaptation(AHZ, time_, AHY, AH11); total_volume_nervous_feedback: Guyton_Model__total_volume_nervous_feedback(AHZ, AHY, AH7); volume_effect_on_arteries: Guyton_Model__volume_effect_on_arteries(AH7, ATRRFB, ATRFBM); volume_effect_on_unstressed_venous_volume: Guyton_Model__volume_effect_on_unstressed_venous_volume(AH7, ATRVFB, ATRVM); // Variable initializations: time_ = ; PRA = ; AH9 = 1; AH10 = 0.333; AH11 = 1000; AHY = 0.301963; ATRFBM = 0; ATRVM = 0; end model Guyton_Model__kidney_parameter_values(GBL, RAPRSP, RFCDFT, RCDFDP, RCDFPC, RNAUGN, RNAULL, RNAUUL, RNAGTC, RNAUAD, ARF, ANMAM, ANMARL, AARK, AARLL, ANPXAF, AUMK1, ANMEM, EFAFR, EARK, EARLL, REK, GPPD, GLPCA, PXTP, GFLC, GFNLL, GFNDMP, URFORM, MDFL1, RTPPR, RTPPRS, RTSPRS, MDFLKM, ANMKEM, RVRS, RFABX, RABSC, RFABDP, RFABDM, DTNAR, DTNARL, DIURET, AHMNAR, ANMNAM, RFABKM, ANMKEL, CKEEX) // Variable initializations: GBL = 0; RAPRSP = 0; RFCDFT = 0; RCDFDP = 2000; RCDFPC = 0; RNAUGN = 0.6; RNAULL = 0.3; RNAUUL = 10; RNAGTC = 15; RNAUAD = 0; ARF = 0.5; ANMAM = 1.4; ANMARL = 0.86; AARK = 1; AARLL = 4; ANPXAF = 1.5; AUMK1 = 0.3; ANMEM = 1.6; EFAFR = 0; EARK = 1; EARLL = 24; REK = 1; GPPD = 1.0; GLPCA = 1.0; PXTP = 8; GFLC = 0.0208333; GFNLL = 0.001; GFNDMP = 3; URFORM = 0.24; MDFL1 = 10; RTPPR = 0.9; RTPPRS = 15.2; RTSPRS = 6; MDFLKM = 0.667; ANMKEM = 2; RVRS = 19.167; RFABX = 0.8; RABSC = 0.5; RFABDP = 1; RFABDM = 0.3; DTNAR = 0.675; DTNARL = 1e-06; DIURET = 1; AHMNAR = 0.3; ANMNAM = 1; RFABKM = 0.03; ANMKEL = 0.3; CKEEX = 4; end model Guyton_Model__perfusion_pressure(PA, time_, PAR, GBL, RAPRSP, RFCDFT, RCDFPC, RCDFDP) // Variable initializations: PA = ; time_ = ; GBL = ; RAPRSP = ; RFCDFT = ; RCDFPC = ; RCDFDP = ; PAR1 = 103.525; end model Guyton_Model__renal_autoregulatory_feedback_factor(MDFLW, time_, RNAUG2, RNAUGN, RNAULL, RNAUUL, RNAUAD) // Variable initializations: MDFLW = ; time_ = ; RNAUGN = ; RNAULL = ; RNAUUL = ; RNAUAD = ; RNAUG3 = 0.0; end model Guyton_Model__afferent_arterial_resistance() end model Guyton_Model__autonomic_effect_on_AAR(AUM, AUMK, ARF) // Variable initializations: AUM = ; ARF = ; end model Guyton_Model__angiotensin_effect_on_AAR(ANM, ANMAR, ANMAM, ANMARL) // Variable initializations: ANM = ; ANMAM = ; ANMARL = ; end model Guyton_Model__AAR_calculation(PAMKRN, AUMK, RNAUG2, ANMAR, MYOGRS, AAR1, AARK) // Variable initializations: PAMKRN = ; AUMK = ; RNAUG2 = ; ANMAR = ; MYOGRS = ; AARK = ; end model Guyton_Model__atrial_natriuretic_peptide_effect_on_AAR(AAR1, ANPX, AAR, ANPXAF, AARLL) // Variable initializations: AAR1 = ; ANPX = ; ANPXAF = ; AARLL = ; end model Guyton_Model__efferent_arterial_resistance() end model Guyton_Model__autonomic_effect_on_EAR(AUMK, AUMK2, AUMK1) // Variable initializations: AUMK = ; AUMK1 = ; end model Guyton_Model__angiotensin_effect_on_EAR(ANM, ANMER, ANMEM) // Variable initializations: ANM = ; ANMEM = ; end model Guyton_Model__effect_of_renal_autoregulatory_feedback_on_EAR(RNAUG2, RNAUG4, EFAFR) // Variable initializations: RNAUG2 = ; EFAFR = ; end model Guyton_Model__EAR_calculation(ANMER, AUMK2, MYOGRS, RNAUG4, EAR, EARK, EARLL) // Variable initializations: ANMER = ; AUMK2 = ; MYOGRS = ; RNAUG4 = ; EARK = ; EARLL = ; end model Guyton_Model__total_renal_resistance(AAR, EAR, RR) // Variable initializations: AAR = ; EAR = ; end model Guyton_Model__normal_renal_blood_flow(RR, PAR, RFN) // Variable initializations: RR = ; PAR = ; end model Guyton_Model__actual_renal_blood_flow(RFN, RBF, REK) // Variable initializations: RFN = ; REK = ; end model Guyton_Model__glomerular_capillaries() end model Guyton_Model__glomerular_colloid_osmotic_pressure(HM1, RFN, GFN, PPC, GLPC, GPPD, GLPCA) // Variable initializations: HM1 = ; RFN = ; GFN = ; PPC = ; GPPD = ; GLPCA = ; end model Guyton_Model__glomerular_pressure(AAR, PAR, RFN, GLP) // Variable initializations: AAR = ; PAR = ; RFN = ; end model Guyton_Model__glomerular_filtration_rate(GLP, GLPC, GFN, GFR, PXTP, GFLC, GFNLL, REK) // Variable initializations: GLP = ; GLPC = ; PXTP = ; GFLC = ; GFNLL = ; REK = ; end model Guyton_Model__glomerular_urea_concentration(UROD, time_, PLUR, URFORM) // Variable initializations: UROD = ; time_ = ; PLUR = 159.549; URFORM = ; end model Guyton_Model__proximal_tubular_and_macula_densa_flow(GFN, MDFLW, MDFL1) // Variable initializations: GFN = ; MDFL1 = ; end model Guyton_Model__renal_tissue_osmotic_pressure(GLPC, RTSPPC, RTPPR, RTPPRS) // Variable initializations: GLPC = ; RTPPR = ; RTPPRS = ; end model Guyton_Model__urea_handling() end model Guyton_Model__plasma_urea_concentration(PLUR, VTW, PLURC) // Variable initializations: PLUR = ; VTW = ; end model Guyton_Model__renal_peritubular_capillaries() end model Guyton_Model__peritubular_capillary_pressure(RFN, RCPRS, RFABX, RVRS) // Variable initializations: RFN = ; RFABX = ; RVRS = ; end model Guyton_Model__peritubular_capillary_reabsorption_factor(RCPRS, GLPC, RTSPPC, RFABD, RTSPRS, RABSC, RFABDP, RFABDM) // Variable initializations: RCPRS = ; GLPC = ; RTSPPC = ; RTSPRS = ; RABSC = ; RFABDP = ; RFABDM = ; end model Guyton_Model__sodium_and_potassium_handling() end model Guyton_Model__distal_tubular_Na_delivery(MDFLW, CNA, DTNAI) // Variable initializations: MDFLW = ; CNA = ; end model Guyton_Model__Na_reabsorption_into_distal_tubules(ADHMK, AMNA, RFABD, DTNARA, DTNAR, DIURET, AHMNAR, DTNARL) // Variable initializations: ADHMK = ; AMNA = ; RFABD = ; DTNAR = ; DIURET = ; AHMNAR = ; DTNARL = ; end model Guyton_Model__angiotensin_induced_Na_reabsorption_into_distal_tubules(ANM, DTNANG, ANMNAM) // Variable initializations: ANM = ; ANMNAM = ; end model Guyton_Model__distal_tubular_K_delivery(DTNAI, CNA, CKE, DTKI) // Variable initializations: DTNAI = ; CNA = ; CKE = ; end model Guyton_Model__effect_of_physical_forces_on_distal_K_reabsorption(RFABD, RFABK, RFABKM) // Variable initializations: RFABD = ; RFABKM = ; end model Guyton_Model__effect_of_fluid_flow_on_distal_K_reabsorption(MDFLW, MDFLK, MDFLKM) // Variable initializations: MDFLW = ; MDFLKM = ; end model Guyton_Model__K_reabsorption_into_distal_tubules(KODN, VUDN, time_, DTKA) // Variable initializations: KODN = ; VUDN = ; time_ = ; DTKA = 0.0367573; end model Guyton_Model__K_secretion_from_distal_tubules(CKE, AMK, MDFLK, ANM, DTKSC, ANMKEM, ANMKEL, CKEEX) // Variable initializations: CKE = ; AMK = ; MDFLK = ; ANM = ; ANMKEM = ; ANMKEL = ; CKEEX = ; end model Guyton_Model__urinary_excretion() end model Guyton_Model__normal_Na_excretion(DTNAI, DTNANG, DTNARA, NODN) // Variable initializations: DTNAI = ; DTNANG = ; DTNARA = ; end model Guyton_Model__normal_K_excretion(DTKI, RFABK, DTKSC, DTKA, KODN) // Variable initializations: DTKI = ; RFABK = ; DTKSC = ; DTKA = ; end model Guyton_Model__normal_urea_excretion(GFN, PLURC, DTURI) // Variable initializations: GFN = ; PLURC = ; end model Guyton_Model__normal_osmolar_and_water_excretion(DTURI, NODN, KODN, OSMOPN1, OSMOPN) // Variable initializations: DTURI = ; NODN = ; KODN = ; end model Guyton_Model__normal_urine_volume(OSMOPN, OSMOPN1, ADHMK, VUDN) // Variable initializations: OSMOPN = ; OSMOPN1 = ; ADHMK = ; end model Guyton_Model__actual_Na_excretion_rate(NODN, NOD, REK) // Variable initializations: NODN = ; REK = ; end model Guyton_Model__actual_K_excretion_rate(KODN, KOD, REK) // Variable initializations: KODN = ; REK = ; end model Guyton_Model__actual_urea_excretion_rate(DTURI, UROD, REK) // Variable initializations: DTURI = ; REK = ; end model Guyton_Model__actual_urine_volume(VUDN, VUD, REK) // Variable initializations: VUDN = ; REK = ; end model Guyton_Model__kidney(time_, ADHMK, AMK, AMNA, ANM, ANPX, AUM, CKE, CNA, HM1, MYOGRS, PA, PAMKRN, PPC, VTW, RBF, MDFLW, NOD, KOD, VUD) // Sub-modules, and any changes to those submodules: kidney_parameter_values: Guyton_Model__kidney_parameter_values(GBL, RAPRSP, RFCDFT, RCDFDP, RCDFPC, RNAUGN, RNAULL, RNAUUL, RNAGTC, RNAUAD, ARF, ANMAM, ANMARL, AARK, AARLL, ANPXAF, AUMK1, ANMEM, EFAFR, EARK, EARLL, REK, GPPD, GLPCA, PXTP, GFLC, GFNLL, GFNDMP, URFORM, MDFL1, RTPPR, RTPPRS, RTSPRS, MDFLKM, ANMKEM, RVRS, RFABX, RABSC, RFABDP, RFABDM, DTNAR, DTNARL, DIURET, AHMNAR, ANMNAM, RFABKM, ANMKEL, CKEEX); perfusion_pressure: Guyton_Model__perfusion_pressure(PA, time_, PAR, GBL, RAPRSP, RFCDFT, RCDFPC, RCDFDP); renal_autoregulatory_feedback_factor: Guyton_Model__renal_autoregulatory_feedback_factor(MDFLW, time_, RNAUG2, RNAUGN, RNAULL, RNAUUL, RNAUAD); afferent_arterial_resistance: Guyton_Model__afferent_arterial_resistance(); autonomic_effect_on_AAR: Guyton_Model__autonomic_effect_on_AAR(AUM, AUMK, ARF); angiotensin_effect_on_AAR: Guyton_Model__angiotensin_effect_on_AAR(ANM, ANMAR, ANMAM, ANMARL); AAR_calculation: Guyton_Model__AAR_calculation(PAMKRN, AUMK, RNAUG2, ANMAR, MYOGRS, AAR1, AARK); atrial_natriuretic_peptide_effect_on_AAR: Guyton_Model__atrial_natriuretic_peptide_effect_on_AAR(AAR1, ANPX, AAR, ANPXAF, AARLL); efferent_arterial_resistance: Guyton_Model__efferent_arterial_resistance(); autonomic_effect_on_EAR: Guyton_Model__autonomic_effect_on_EAR(AUMK, AUMK2, AUMK1); angiotensin_effect_on_EAR: Guyton_Model__angiotensin_effect_on_EAR(ANM, ANMER, ANMEM); effect_of_renal_autoregulatory_feedback_on_EAR: Guyton_Model__effect_of_renal_autoregulatory_feedback_on_EAR(RNAUG2, RNAUG4, EFAFR); EAR_calculation: Guyton_Model__EAR_calculation(ANMER, AUMK2, MYOGRS, RNAUG4, EAR, EARK, EARLL); total_renal_resistance: Guyton_Model__total_renal_resistance(AAR, EAR, RR); normal_renal_blood_flow: Guyton_Model__normal_renal_blood_flow(RR, PAR, RFN); actual_renal_blood_flow: Guyton_Model__actual_renal_blood_flow(RFN, RBF, REK); glomerular_capillaries: Guyton_Model__glomerular_capillaries(); glomerular_colloid_osmotic_pressure: Guyton_Model__glomerular_colloid_osmotic_pressure(HM1, RFN, GFN, PPC, GLPC, GPPD, GLPCA); glomerular_pressure: Guyton_Model__glomerular_pressure(AAR, PAR, RFN, GLP); glomerular_filtration_rate: Guyton_Model__glomerular_filtration_rate(GLP, GLPC0, GFN, GFR, PXTP, GFLC, GFNLL, REK); glomerular_urea_concentration: Guyton_Model__glomerular_urea_concentration(UROD, time_, PLUR, URFORM); proximal_tubular_and_macula_densa_flow: Guyton_Model__proximal_tubular_and_macula_densa_flow(GFN0, MDFLW, MDFL1); renal_tissue_osmotic_pressure: Guyton_Model__renal_tissue_osmotic_pressure(GLPC, RTSPPC, RTPPR, RTPPRS); urea_handling: Guyton_Model__urea_handling(); plasma_urea_concentration: Guyton_Model__plasma_urea_concentration(PLUR, VTW, PLURC); renal_peritubular_capillaries: Guyton_Model__renal_peritubular_capillaries(); peritubular_capillary_pressure: Guyton_Model__peritubular_capillary_pressure(RFN, RCPRS, RFABX, RVRS); peritubular_capillary_reabsorption_factor: Guyton_Model__peritubular_capillary_reabsorption_factor(RCPRS, GLPC, RTSPPC, RFABD, RTSPRS, RABSC, RFABDP, RFABDM); sodium_and_potassium_handling: Guyton_Model__sodium_and_potassium_handling(); distal_tubular_Na_delivery: Guyton_Model__distal_tubular_Na_delivery(MDFLW, CNA, DTNAI); Na_reabsorption_into_distal_tubules: Guyton_Model__Na_reabsorption_into_distal_tubules(ADHMK, AMNA, RFABD, DTNARA, DTNAR, DIURET, AHMNAR, DTNARL); angiotensin_induced_Na_reabsorption_into_distal_tubules: Guyton_Model__angiotensin_induced_Na_reabsorption_into_distal_tubules(ANM, DTNANG, ANMNAM); distal_tubular_K_delivery: Guyton_Model__distal_tubular_K_delivery(DTNAI, CNA, CKE, DTKI); effect_of_physical_forces_on_distal_K_reabsorption: Guyton_Model__effect_of_physical_forces_on_distal_K_reabsorption(RFABD, RFABK, RFABKM); effect_of_fluid_flow_on_distal_K_reabsorption: Guyton_Model__effect_of_fluid_flow_on_distal_K_reabsorption(MDFLW, MDFLK, MDFLKM); K_reabsorption_into_distal_tubules: Guyton_Model__K_reabsorption_into_distal_tubules(KODN, VUDN, time_, DTKA); K_secretion_from_distal_tubules: Guyton_Model__K_secretion_from_distal_tubules(CKE, AMK, MDFLK, ANM, DTKSC, ANMKEM, ANMKEL, CKEEX); urinary_excretion: Guyton_Model__urinary_excretion(); normal_Na_excretion: Guyton_Model__normal_Na_excretion(DTNAI, DTNANG, DTNARA, NODN); normal_K_excretion: Guyton_Model__normal_K_excretion(DTKI, RFABK, DTKSC, DTKA, KODN); normal_urea_excretion: Guyton_Model__normal_urea_excretion(GFN, PLURC, DTURI); normal_osmolar_and_water_excretion: Guyton_Model__normal_osmolar_and_water_excretion(DTURI, NODN, KODN, OSMOPN1, OSMOPN); normal_urine_volume: Guyton_Model__normal_urine_volume(OSMOPN, OSMOPN1, ADHMK, VUDN); actual_Na_excretion_rate: Guyton_Model__actual_Na_excretion_rate(NODN, NOD, REK); actual_K_excretion_rate: Guyton_Model__actual_K_excretion_rate(KODN, KOD, REK); actual_urea_excretion_rate: Guyton_Model__actual_urea_excretion_rate(DTURI, UROD, REK); actual_urine_volume: Guyton_Model__actual_urine_volume(VUDN, VUD, REK); // Variable initializations: time_ = ; ADHMK = ; AMK = ; AMNA = ; ANM = ; ANPX = ; AUM = ; CKE = ; CNA = ; HM1 = ; MYOGRS = ; PA = ; PAMKRN = ; PPC = ; VTW = ; GBL = 0; RAPRSP = 0; RFCDFT = 0; RCDFPC = 0; RCDFDP = 2000; RNAUGN = 0.6; RNAULL = 0.3; RNAUUL = 10; RNAUAD = 0; ARF = 0.5; ANMAM = 1.4; ANMARL = 0.86; AARK = 1; ANPXAF = 1.5; AARLL = 4; AUMK1 = 0.3; ANMEM = 1.6; EFAFR = 0; EARK = 1; EARLL = 24; REK = 1; GPPD = 1.0; GLPCA = 1.0; PXTP = 8; GFLC = 0.0208333; GFNLL = 0.001; MDFL1 = 10; RTPPR = 0.9; RTPPRS = 15.2; URFORM = 0.24; PLUR = 159.549; RFABX = 0.8; RVRS = 19.167; RTSPRS = 6; RABSC = 0.5; RFABDP = 1; RFABDM = 0.3; DTNAR = 0.675; DIURET = 1; AHMNAR = 0.3; DTNARL = 1e-06; ANMNAM = 1; RFABKM = 0.03; MDFLKM = 0.667; DTKA = 0.0367573; ANMKEM = 2; ANMKEL = 0.3; CKEEX = 4; RNAGTC = 15; GFNDMP = 3; GLPC0 = ; GFN0 = ; end model *Guyton_Model____main() // Sub-modules, and any changes to those submodules: environment: Guyton_Model__environment(time_); temp_myogrs_and_pamk: Guyton_Model__temp_myogrs_and_pamk(MYOGRS, PAMK, PAMKRN); aldosterone: Guyton_Model__aldosterone(time_, ANM, CKE, AMK, AMNA); angiotensin: Guyton_Model__angiotensin(time_, MDFLW, ANUVN, ANU, ANM); antidiuretic_hormone: Guyton_Model__antidiuretic_hormone(time_, CNA, PA1, ADHMK, ADHMV, ADHC); atrial_natriuretic_peptide: Guyton_Model__atrial_natriuretic_peptide(time_, PLA, PRA, ANPX); autonomics: Guyton_Model__autonomics(time_, PA, PO2ART, PLA, PRA, PPA, PA1, VVR, AUH, AUR, AUM, AOM, AVE, AU); capillary_dynamics: Guyton_Model__capillary_dynamics(time_, VEC, PPD, RVS, DFP, VPF, BFN, PVS, PC, CPP, PPC, VP); circulatory_dynamics: Guyton_Model__circulatory_dynamics(); electrolytes: Guyton_Model__electrolytes(time_, AMK, TVD, NOD, STH, KOD, VUD, CKE, VEC, CNA, VTW); heart_hypertrophy_or_deterioration: Guyton_Model__heart_hypertrophy_or_deterioration(time_, QAO, PA, POT, PPA, HPL, HPR, HMD); heart_rate_and_stroke_volume: Guyton_Model__heart_rate_and_stroke_volume(QLO, AUR, PRA, HMD); muscle_autoregulatory_local_blood_flow_control: Guyton_Model__muscle_autoregulatory_local_blood_flow_control(time_, PMO, AMM); muscle_O2_delivery: Guyton_Model__muscle_O2_delivery(time_, BFM, OVA, HM, AOM, PMO, RMO); non_muscle_autoregulatory_local_blood_flow_control: Guyton_Model__non_muscle_autoregulatory_local_blood_flow_control(time_, POT, ARM); non_muscle_O2_delivery: Guyton_Model__non_muscle_O2_delivery(time_, BFN, OVA, HM, AOM, POT, DOB); pulmonary_fluid_dynamics: Guyton_Model__pulmonary_fluid_dynamics(time_, PPC, PPA, PLA, CPP, RPV, RPA, DFP, VPF, PPD); pulmonary_O2_uptake: Guyton_Model__pulmonary_O2_uptake(time_, VPF, DOB, QRO, RMO, HM, OVA, PO2ART, OSA); red_cells_and_viscosity: Guyton_Model__red_cells_and_viscosity(time_, VP, VIM, HM, HM1, VRC); stress_relaxation: Guyton_Model__stress_relaxation(time_, VVE, VV6, VV7); thirst_drinking_and_salt_appetite: Guyton_Model__thirst_drinking_and_salt_appetite(time_, ADHC, ANM, POT, TVD, STH); volume_receptors: Guyton_Model__volume_receptors(time_, PRA, ATRRFB, ATRVFB); kidney: Guyton_Model__kidney(time_, ADHMK, AMK, AMNA, ANM, ANPX, AUM, CKE, CNA, HM1, MYOGRS, PA, PAMKRN, PPC, VTW, RBF, MDFLW, NOD, KOD, VUD); circulatory_dynamics.QLO is QLO; circulatory_dynamics.PRA is PRA; circulatory_dynamics.PLA is PLA; circulatory_dynamics.PA is PA; circulatory_dynamics.PPA is PPA; circulatory_dynamics.AVE is AVE; circulatory_dynamics.VVR is VVR; circulatory_dynamics.AUH is AUH; circulatory_dynamics.AUM is AUM; circulatory_dynamics.AU is AU; circulatory_dynamics.ANU is ANU; circulatory_dynamics.ANUVN is ANUVN; circulatory_dynamics.ADHMV is ADHMV; circulatory_dynamics.PAMK is PAMK; circulatory_dynamics.MYOGRS is MYOGRS; circulatory_dynamics.ARM is ARM; circulatory_dynamics.VIM is VIM; circulatory_dynamics.VRC is VRC; circulatory_dynamics.ATRVFB is ATRVFB; circulatory_dynamics.ATRRFB is ATRRFB; circulatory_dynamics.VP is VP; circulatory_dynamics.RVS is RVS; circulatory_dynamics.BFN is BFN; circulatory_dynamics.PVS is PVS; circulatory_dynamics.HPR is HPR; circulatory_dynamics.HMD is HMD; circulatory_dynamics.HPL is HPL; circulatory_dynamics.QAO is QAO; circulatory_dynamics.OSA is OSA; circulatory_dynamics.QRO is QRO; circulatory_dynamics.AMM is AMM; circulatory_dynamics.time_ is time_; circulatory_dynamics.BFM is BFM; circulatory_dynamics.RPV is RPV; circulatory_dynamics.RPA is RPA; circulatory_dynamics.VVE is VVE; circulatory_dynamics.VV6 is VV6; circulatory_dynamics.VV7 is VV7; circulatory_dynamics.RBF is RBF; // Variable initializations: VRC = 2.00439; HPR = 1.00237; HPL = 1.00163; TVD = 0.000980838; OVA = 204.497; VV6 = 0.0101913; VV7 = 0.00366525; VTW = 39.8952; end