A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Similar documents
Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

Charging of capacitor through inductor and resistor

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Midterm exam 2, April 7, 2009 (solutions)

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

CSE 245: Computer Aided Circuit Simulation and Verification

On the Speed of Heat Wave. Mihály Makai

Lecture 2: Current in RC circuit D.K.Pandey

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Elementary Differential Equations and Boundary Value Problems

A Condition for Stability in an SIR Age Structured Disease Model with Decreasing Survival Rate

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

Wave Equation (2 Week)

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

H is equal to the surface current J S

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

Smoking Tobacco Experiencing with Induced Death

EXERCISE - 01 CHECK YOUR GRASP

Impulsive Differential Equations. by using the Euler Method

Microscopic Flow Characteristics Time Headway - Distribution

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

4.3 Design of Sections for Flexure (Part II)

Transfer function and the Laplace transformation

C From Faraday's Law, the induced voltage is, C The effect of electromagnetic induction in the coil itself is called selfinduction.

Institute of Actuaries of India

Reliability Analysis of a Bridge and Parallel Series Networks with Critical and Non- Critical Human Errors: A Block Diagram Approach.

Logistic equation of Human population growth (generalization to the case of reactive environment).

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

Lecture 4: Laplace Transforms

Control System Engineering (EE301T) Assignment: 2

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

On Ψ-Conditional Asymptotic Stability of First Order Non-Linear Matrix Lyapunov Systems

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

t 2 B F x,t n dsdt t u x,t dxdt

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

EE 350 Signals and Systems Spring 2005 Sample Exam #2 - Solutions

Chapter 6 Differential Equations and Mathematical Modeling

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

46. Let y = ln r. Then dy = dr, and so. = [ sin (ln r) cos (ln r)

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

Modelling of three dimensional liquid steel flow in continuous casting process

System of Linear Differential Equations

3.4 Repeated Roots; Reduction of Order

Differential Equations

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

Midterm Examination (100 pts)

Stability Analysis of a delayed HIV/AIDS Epidemic Model with Saturated Incidence

Asymptotic Solutions of Fifth Order Critically Damped Nonlinear Systems with Pair Wise Equal Eigenvalues and another is Distinct

Double Slits in Space and Time

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

Discussion 06 Solutions

The Optimal Timing of Transition to New Environmental Technology in Economic Growth

4.6 One Dimensional Kinematics and Integration

( ) ( ) + = ( ) + ( )

10. If p and q are the lengths of the perpendiculars from the origin on the tangent and the normal to the curve

EE 434 Lecture 22. Bipolar Device Models

Math 334 Fall 2011 Homework 11 Solutions

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

DISCRETE GRONWALL LEMMA AND APPLICATIONS

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Chapter 12 Introduction To The Laplace Transform

On General Solutions of First-Order Nonlinear Matrix and Scalar Ordinary Differential Equations

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

6.2 Transforms of Derivatives and Integrals.

CHAPTER 2 Signals And Spectra

Physics 160 Lecture 3. R. Johnson April 6, 2015

Chapter 3 Boundary Value Problem

Math 333 Problem Set #2 Solution 14 February 2003

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

Reliability of Technical Systems

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

ME 391 Mechanical Engineering Analysis

Modelling of repairable items for production inventory with random deterioration

2. The Laplace Transform

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

Simulation of multi-component mass diffusion pellet model using least squares spectral element method for steam methane reforming process

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

Solutionbank Edexcel AS and A Level Modular Mathematics

EXERCISES FOR SECTION 1.5

Complex Dynamic Models of Star and Delta Connected Multi-phase Asynchronous Motors

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Transcription:

A THREE COPARTENT ATHEATICAL ODEL OF LIVER V. An N. Ch. Paabhi Ramacharyulu Faculy of ahmaics, R D collgs, Hanamonda, Warangal, India Dparmn of ahmaics, Naional Insiu of Tchnology, Warangal, India E-ail: valupadas@gmail.com ABSTRACT ahmaical modling of pharmacoinics is an imporan growing fild in drug dvlopmn. Pharmacoinics concrns wih h disribuion of drugs, chmicals or racrs by a fluid among h various comparmn of human body. In his wor w discuss h comparmn mahmaical modl of livr funcion basd on fundamnal biological pharmacological principls. Hr w prsn bhavior of hyroxin, iodin bil ovr a priod of im. Kywords: livr, mahmaical modl, diffrnial quaions, pharmacoinics, comparmn modl. INTRODUCTION Pharmacoinics concrns wih h disribuion of drugs, chmicals or racrs by a fluid among h various comparmn of human body []. Th comparmns could b ficiious spacs hrough which biomarials pass hrough various locaions (comparmns of h body. Th prsn invsigaion is on a hr comparmn modl rlad o h livr in human body [,, 4]. Whn a chmical, hyroxin is injcd ino h blood sram i is carrid o h livr. Th livr convrs hyroxin o iodin, which is absorbd ino h bil [5]. Howvr, nihr h convrsion nor h absorpion of which ino h bil, would occur insananously. Som of h hyroxin (unconvrd rnrs ino h blood sram gs rcirculad. Th ahmaical modl is composd of hr comparmns: Comparmn I, Comparmn II comparmn III which rprsn blood vssls, livr bil rspcivly as shown in h Figur-. Figur-. ahmaical modl for livr. Noaions adopd: x ( = h quaniy of hyroxin in h blood vssl a h insan im x ( = h quaniy of iodin in h livr x ( = h quaniy of iodin absorbd in o h bil = Th ra of convrsion of hyroxin ino iodin = Th ra of h quaniy of unabsorbd hyroxin sn ou for rcycling from Comparmn II o comparmn I = ra of absorpion of Iodin from comparmn II ino bil comparmn III x 0, x 0 x 0 ar iniial h valus of x, x, x rspcivly h ra consans x, x x ar all posiiv. I is assumd ha h ra x ( a which hyroxin is convrd o iodin as i ransfrrd from comparmn I o comparmn II is proporional o concnraion x ( of hyroxin in h comparmn I. ODEL BLOCK - DIAGRAS AND ODEL EQUATIONS dx x x, x(0 x0 ( 0

dx x x, x(0 x0 ( ( 0 4 (8 Evidnly is ral lss han. Hnc ar boh ral ngaiv roos. Th soluions of quaions (, (, ( can b wrin as dx, x (0 x 0 ( x Th quaions (, (, ( can b pu ino h of marix form d whr X AX (4 A 0 ( 0 x 0 X x 0 x L X X 0 b a rial soluion wih iniial condiions X 0 [x 0 x 0 x 0 ] T Th xponn λ saisfis h characrisic quaion of A: d AI 0 i, ; ( 0 (6 h roos of which ar λ=0, (7 whr = (5 x AB C (9 ( x ( A B( C( B( C( x( A E (0 ( whr A, B, C E ar arbirary consans. Using h iniial condiions, w g h valus of h consans A = 0 ( x0 ( x0 B ( x x 0 ( 0 C (4 E ( x0 x0 x0 (5 Subsiuing hs valus in (., (., (. w g afr som simplificaion x ( x ( 0 0 x x0 ( x0 x0 x0( sinh ( x0cosh + (6 0

x0( x0 x0x0( x( ( ( ( x0 x0 sinh ( x0cosh (7 x0 ( x0 x0 x0( x( ( ( ( x0 x0 x0 ( ( x0 ( x0 x0 sinh ( x0 x0cosh ( x0 x0 x0 (8 Th variaions of x (, x ( x ( Vrs ar illusrad for a slc rang of valus of,, (vid Figur- o Figur-5 for h iniial valus x0 50, x0 5 x0 65 of hyroxin, iodin bil, rspcivly. Figur-. Variaion of hyroxin, iodin bil for for h ransfr cofficins =0.85 =0.05 =0.000. Figur-. Variaion of hyroxin, iodin, bil for h ransfr cofficins =0.85 =0.056 =0.006. Figur-4. Variaion of hyroxin, iodin, bil for h ransfr cofficins =0.85, =0.05 =0.000. 04

Figur-5. Variaion of hyroxin, iodin, bil for h ransfr cofficins =0.85, =0.055 =0.0005. CONCLUSIONS Thyroxin monoonically dcrass raching zro lvl. Iodin iniially riss o rach a maximum falls asympoically o approach o zro as is h cas wih hyroxin. Th bil incras monoonically approachs asympoic lvl. Figur-6. =0.85, =0.05, =0.006. Incipin variaions of x (, x ( x ( (i..; variaions for small rang im I is now ( cosh O( sinh O(! Nglcing rm of O (, w g xx0 ( 0 x x0 ( ( x0 0 x (9 Figur-7. =0.85, =0.05, =0.007. x x ( x x x ( x 0 0 0 0 0 (0 x x0( x0 ( x0 ( x0 ( Th variaions of x (, x ( x ( Vrs small im ar illusrad for a slc rang of valus of,, (vid Figur-6 o Figur-9 for h iniial valus x0 50, x0 5 x0 65 of hyroxin, iodin bil, rspcivly. Figur-8. =0.85, =0.05, =0.008. 05

[] Jacquz John A. 97. Comparmn analysis in biology mdicin. Nw Yor:Elsvir Scinific publishing company. [4] J.. Wa Andrw Young. 96. An amp o simula h livr on a compur. Compur. Journal. 5: -7. [5] Goffry Gordon. 999. Sysm simulaion. Prnichal India. pp. 4-6. Figur-9. =0.85, =0.05, =0.000. Asympoic variaion of x (, x ( x ( (variaion for larg rang im In his cas dominas ovr (sinc >0 cosh sinh Hnc h asympoic xprssion for x (, x (, x ( ar ( K x ( x0 x0( K ( K x ( x0 x0 ( K ( x K ( x ( x x ( K ( x x 0 0 0 0 0 x0 ( ( (4 CONCLUSIONS For larg rang im, ach of x (, x ( x ( dcras xponnially wih h Characrisic im. REFERENCES [] J.N. Kapur. 985. ahmaical modls in Biology dicin. Affiliad as-ws prss pv. Ld. pp. 6-7. [] Evr C.F.F. Ral. 970. Formulaion compuaion of comparmn modls. J. pharm. Sci. 9(: 0-4. 06