Practical 5
Practical 5
Section 1
Case of a single cold pill and a course of cold pills
Question:-(1) Solve drug single cold pill model and plot its graph.
Solution:-
Needs"PlotLegends`"
eqn1 x 't k1 xt
eqn2 y 't k1 xt k2 yt
sol DSolveeqn1, eqn2, x0 1, y0 0, x, y, t
solx xt . sol .
k1 1.3860, k2 0.1386, k1 0.6931, k2 0.0231, k1 0.6931, k2 0.08
soly yt . sol . k1 1.3860, k2 0.1386,
k1 0.6931, k2 0.0231, k1 0.6931, k2 0.08
Plotsolx1, solx2, t, 0, 15,
PlotStyle Thickness0.02, Thickness0.03,
PlotLegend "Decogestant", "Antihistamine in GITract"
Plotsoly1, soly2, soly3, t, 0, 15,
PlotStyle Thickness0.02, Thickness0.03, Thick,
PlotLegend "Decogestant", "Antihistamine", "Antihistamine in blood"
x t k1 xt
0.20
0.15
0.10
0.05
Decogestant
2 4 6 8 10 12 14
Antihistamine in GI-T
0.8
0.6
0.4
Decogestant
0.2
Antihistamine
2 4 6 8 10 12 14
Antihistamine in blo
Question:-(2) Solve course of cold pills model and plot its graph.
Solution:-
Practical_5[1].nb 3
Needs"PlotLegends`"
eqn1 x 't A k1 xt
eqn2 y 't k1 xt k2 yt
A1
sol DSolveeqn1, eqn2, x0 0, y0 0, x, y, t
solx xt . sol . k1 1.3860, k2 0.1386, k1 0.6931, k2 0.0231
soly yt . sol . k1 1.3860, k2 0.1386, k1 0.6931, k2 0.0231
Plotsolx1, solx2, t, 0, 60,
PlotStyle Thickness0.02, Thickness0.03,
PlotLegend "Decogestant", "Antihistamine in GITract"
Plotsoly1, soly2, t, 0, 60,
PlotStyle Thickness0.02, Thickness0.03, Thick,
PlotLegend "Decogestant", "Antihistamine"
x t A k1 xt
k1 t 1 k1 t
x Functiont, , y
Functiont, k1 tk2 t k1 t k1 k1 tk2 t k1 k2 t k2 k1 tk2 t k2 k2 k1 k2
k1
5.78404 1.5246 t 0.1386 0.1386 t 1.386 1.386 t 1.2474 1.5246 t ,
64.612 0.7162 t 0.0231 0.0231 t 0.6931 0.6931 t 0.67 0.7162 t
1.4
1.2
1.0
0.8
0.6
0.4
0.2 Decogestant
10 20 30 40 50 60
Antihistamine in GI-T
4 Practical_5[1].nb
30
25
20
15
10
5 Decogestant
10 20 30 40 50 60
Antihistamine
Section-2
Limited growth of population(with and without
harvesting)
Solution:-
Practical_5[1].nb 5
r 1; k 1000;
eqn x 't r xt 1 xt k
sol DSolveeqn, x0 1, xt, t
Plotxt . sol, t, 0, 10, PlotStyle Red, Thickness0.03,
PlotRange All, AxesLabel t, xt
x t 1
xt
xt
1000
6ROYH
LIXQ
,QYHUVH
IXQFWLRQV
DUHEHLQJXVHGE\ 6ROYH
VRVRPH VROXWLRQV
PD\ QRWEH IRXQGXVH5HGXFHIRUFRPSOHWHVROXWLRQ
LQIRUPDWLRQ
à
xt
1000 t
999 t
xt
800
600
400
200
t
2 4 6 8 10
Solution
6 Practical_5[1].nb
r1
k 10
h 9 10
eqn x 't r k xt ^ 2 k xt k h r
sol DSolveeqn, x0 2, xt, t
Plotxt . sol, t, 0, 8, PlotStyle Magenta, Thickness0.02,
PlotRange All, AxesLabel t, xt
10
9
10
6ROYH
LIXQ
,QYHUVH
IXQFWLRQV
DUHEHLQJXVHGE\ 6ROYH
VRVRPH VROXWLRQV
PD\ QRWEH IRXQGXVH5HGXFHIRUFRPSOHWHVROXWLRQ
LQIRUPDWLRQ
à
xt
7 9 4 t5
7 4 t5
xt
9
t
2 4 6 8