Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- In[122]:= Clear[dosage];
- dosage[t_, i_, d_] := d /; Mod[t, i] == 0;
- dosage[t_, i_, d_] := 0 /; Mod[t, i] != 0;
- hurt = 100;
- effektiv = 90;
- halftime = 5;
- interval = 5;
- Clear[fun];
- fun[t_, t0_, dose_] := 0 /; t < t0;
- fun[t_, t0_, dose_] := dose* (1/2)^((t - t0) /halftime) /; t >= t0;
- Clear[test];
- test[t_, i_, d_] := dosage[t, i, d] + antifun[Mod[t, i], d];
- (*DiscretePlot[dosage[t, 4, 200], {t, 0, 20}]
- DiscretePlot[dosage[t, 4, 200] + antifun[Mod[t, interval], 200], {t, \
- 0, 20}]
- Plot[test[t, 4, 200], {t, 0, 20}]*)
- Clear[antifun];
- antifun[t_, h_] := ( h) *(1/2)^(t/halftime) ;
- Plot[antifun[Mod[t, interval], hurt], {t, 0, 10}]
- NMaximize[{antifun[Mod[t, interval], hurt], t >= 0}, t]
- NMinimize[{antifun[Mod[t, interval], hurt], t >= 0}, t]
- (*Plot[NMinimize[{antifun[Mod[t, interval], hurt], t \[GreaterEqual] \
- 0, t \[LessEqual] 24}, t], {interval, 0, 24}]*)
- Out[136]= \!\(\*
- GraphicsBox[{{{}, {},
- {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
- Opacity[1.], LineBox[CompressedData["
- 1:eJwBsQNO/CFib1JlAgAAADoAAAACAAAAnX5iBjFkiz4WMyL0//9YQJ5a7DpX
- IGk/r2qOf0f9WEDZTYpyICB5P2ijxFaP+lhAd0dZDgUgiT+8G27oH/VYQEbE
- QFz3H5k/sHLOl0TqWECtgjSD8B+pP20Akh+c1FhA4WGuFu0fuT+Nub2Zg6lY
- QHtRa2DrH8k/dzZbcDJUWEAJ+ZvJbC7aP0eAVzcmn1dAccJ/jOxy4z/hhu2H
- 2vpWQC6kk9G0ruk/1/OcaB1eVkDlTmgc3zjwP9xwo0rPuFVAlxqGHaFg8z+O
- jxIX2CJVQDlyuK8DzPY/hhRmzuKEVEDD15tQbyf6P6oirJNL7lNAMYz+PhhJ
- /T/upCA+uGVTQEfmOt8wVwBAuU1fRcDVUkDnLbZF9OwBQKRqX0EeU1JAe/yJ
- M7x6A0CDnYm+fdZRQAcR6GlUKgVAyfdGLQNTUUAFzQVHC70GQMRSFnK/21BA
- +86tbJJxCEA/uInJFl5QQOVXrhkeHgpA4ObNP5PMT0BBiG5tyK0LQM2UzxRM
- 805Alf64CUNfDUBx9W4AMA5OQFscw0zc8w5AoKVR01I+TUAMwCvsIlUQQPKQ
- qg1sY0xAZTWi9VksEUALrmT8u5JLQHd+eFIg9RFABKmyhnTVSkCF6pPTzs4S
- QIAsxK73DUpAhNflyvvREkBKfhq/GQtKQITEN8Io1RJAoDUzIDwISkCCntuw
- gtsSQDZKiNSBAkpAflIjjjboEkCodsgEEfdJQHa6skieARNAim3fez7gSUBo
- itG9bTQTQB52j4bVsklAZ3cjtZo3E0CnsA2eAbBJQGZkdazHOhNA/tUzBS6t
- SUBkPhmbIUETQLTSVMKHp0lAYfJgeNVNE0BZZfT2PpxJQFpa8DI9ZxNAF6ED
- QryFSUBMKg+oDJoTQH5nYiLyWElAWchvqCmdE0DYYrUiNlZJQGZm0KhGoBNA
- nqKLbnpTSUCAopGpgKYTQNBaoegDTklAtBoUq/SyE0C4zo1lGkNJQBwLGa7c
- yxNA/qxOe1UtSUApqXmu+c4TQOJVBjCeKklANkfarhbSE0B9Vr8v5ydJQFCD
- m69Q2BNATQIVEHoiSUCE+x2xxOQTQFitbFOjF0lAkZl+seHnE0Detkxf7hRJ
- QJ4337H+6hNArHXttTkSSUC4c6CyOPETQIDSUEPRDElAxREBs1X0E0BcUgN6
- HQpJQNKvYbNy9xNAMUtW+2kHSUDfTcKzj/oTQBiwQce2BElAp+DlFLz8E0CF
- Wp1t1AJJQON3mO8=
- "]], LineBox[CompressedData["
- 1:eJwVj30w1HkAxjcreeuu42ex3pbvXrenvGXUSS9PbS++kbzs8lMu1luuCzPH
- 1NWYG66LVJvKSS5FhblTuCtvu2rTC4rdUCIdd+XEkWqFUm117o9nnvnMM/OZ
- eRyjk4Pj9Dgcjngm/3e4A/+pmMtgx5LwwQXvWRiIBquHZth1Yu3oaR0L8+Y1
- Ndw5DMpPrH409oaFILa0RmDMoHpnlcDnNQsXPYPa5XMZ/LLvp/PZkyyWFsfX
- bp7HIEulSH8wzkLS90XdzzwGyflC411jLGR79tf9Yc2Acd2gah5hkWQ1UnfH
- lsG8vqrLvGEW+yXl9YaEwWCqzYm6xywuqxcq09wYOBzIvaXrZnF7u1xZsIhB
- YVOvJqCLRbfhc2WtFwO5rcDuTCcL7ZqqBq0PAw83ecE6NQvhFY8rsb4MOk8Y
- GeVdYyGvWNzoL2OQ3ZvbsPICi/lRIllFLIOAHEVOxG8sVGZ8vU8SGDw0iFLu
- Lpvx7Xovbk+a2VUG5jXFM39W32wJSmPw9fziDrc8FrY9Qe2h+Qzqd7ukev3A
- 4gIn8a8oDQOe0VpeRiCL1pBzulRvC1za03Pc8lUYGjXJhxZyeVhaLt33IDAM
- F7dnQN3PQ9KtyKtNZ0Nh77fyq6RqSwzg3HSHTgrxx4i9RWlWGDd3erIpXorY
- 7qsv66XW6HIITCxpkWBFxo81MgEfeQdL+BVLJDiWPN5gOsyHgK9+FlIZgoH0
- 0qEilQ383Z6v/8YuBE/jFLHXs2yx7Iwk81VBMLSJZnM0MjtwxDv3pnwajEra
- pot3scdz98KzjUeDMJTgvpg7ZY9JpZzTbR0E/vmHuVkaB3TUto9EngzE3qjv
- zmr8BOj/3nlRqSgQr3XWcj+FANIFBpdOR26CqCG/oFXoiENdUn3v3QHwimna
- 9Fm2I3o/cN3Fv25EY6GFlXrKEdroa5Wz1P4QNVocCNjihHeSUPXn+v7YOHvM
- 7q7KCacsfZL2rPWD9s/BG7edCUzWFywPS9uA5EK5eXYOQSDxxkAThXqW8MWG
- owQ7/k5vU12ncN6mbDXNJYi+W1R/8irFoMdwxpHjBH01nvWhCgq2Bdr8UwSz
- PeP9Oi9QrH45oS47T1DVfKTuQS4Fz3dz5s0Wgi3EXuIeTZFSoZVl3iaIJCFc
- 60iKDrOs5b5tBGUDg2Z6ERQH+y9Ntt0hCPj9RfB9KQUnZW7MvfsEoqDN0nRK
- MXqqceXAPwSKonXh/7pTrOeG2ZQ8IfDOaTXscaEoSXj2Km6Y4JhIa9nsTBHp
- ya8cGSUQN41LSoUUXbdSbMfHCabG9vsmWFF4uBpPX5wgqN7WNRFuQXE4t/he
- 6hTBu+zYD35mFHSr5sD0NMH9rUIPD1OKshsx8cq3BJ2Pl4ULjSj0v3y7Kk1H
- YFvx7SNLAwrZ4SN2Kz4QJC8sbTbhUqgm5r/5+JHghbrfhMOh+A+axtZA
- "]]}, {}}, {{}, {}, {}}},
- AspectRatio->0.6180339887498948,
- Axes->{True, True},
- AxesLabel->{None, None},
- AxesOrigin->{0, 48.},
- DisplayFunction->Identity,
- Frame->{{False, False}, {False, False}},
- FrameLabel->{{None, None}, {None, None}},
- FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
- GridLines->{None, None},
- GridLinesStyle->Directive[
- GrayLevel[0.5, 0.4]],
- Method->{"DefaultBoundaryStyle" -> Automatic,
- "ScalingFunctions" -> None},
- PlotRange->{{0, 10}, {50.0000014145861, 99.99999717082787}},
- PlotRangeClipping->True,
- PlotRangePadding->{{
- Scaled[0.02],
- Scaled[0.02]}, {
- Scaled[0.05],
- Scaled[0.05]}},
- Ticks->{Automatic, Automatic}]\)
- Out[137]= {100., {t -> 2.28154*10^-9}}
- Out[138]= {50., {t -> 5.}}
- In[94]:=
- In[96]:=
- In[139]:=
- Plot[antifun[t, hurt], {t, 0, 24}]
- Out[139]= \!\(\*
- GraphicsBox[{{}, {},
- {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
- Opacity[1.], LineBox[CompressedData["
- 1:eJwVknk01AkAx5WixsyUZISxmBkTkdy/ssvvO6YDreeoaWaSDkde7mVyTJsc
- jWrDilRYR/VIK9tBUiKFDUm02xM5shmd7pVKk7V/fN73ff78vvcx9A33Clio
- pKTkNs//+2lrWL5HeJFj1BbZq7k5EcQqSRpPuPHk7tKyodjPIhRlkDvSuOmk
- WFM3d2pShAZBQUwIt4B0SaUZjMpFEFS7LdnF/YNMqHl6cuaRCAUHiju2cevI
- lKrHR7zPieDk1Fnpzm0nOx6rPfS3EuG0+esSN24/uYJYV+AfLARl9+jEm6Ax
- kneJy8npFOAYc//BmztmyISBE/YuIduh3OTaXNz5jewRmxwy52yDRhODFRmj
- DK/kgWVJXZ6o+lC04z5UUbr+stpPNz1whfV7IeGthmnrqzqV+u7QtxvNkvku
- Q9vAijMSFTd4Hgu4x8tUR2j0Ts/FVltRK9Ug/87TwLDkh+l9Qhdc76DdqUvS
- hGKWTxy/vQX+2cUUzb1aeCNMmTzosBkY/Fn+LUAbnyi3j47VbETxr0U1cNEF
- 2zH0ckswHzZNykyGQA+bUktygrSckHhgkZdQqA/2kaSlkwOAaXUpLzbPAKXd
- zoraIhKmc2GF6YWGsFBwR0ylDgiPefHFI5OFdnkvXy/3ezhXBOr65LDRQ4bE
- iis3IKHkAbcwj4PyKXXe6AyBr60N69JSjaDz3p8RLrEDrWJI0XiKi+CEe19f
- v7eBb1tZS6lsNQSDteKSbGsw8nkJDyTGMHtq07jT2wp7LGdmjWJNwNTRrDnP
- tMTIR6+DtkFrEMgzaLmjYoFlCe0V6RGmqF4wXdnRsxan6lf16UWaIaI1bLWi
- yxRCq6WLxL5rQYvwfa+jtQY3ArXH+/abI+SR834bD2N4tj/ppu9ch+ie9vCM
- Bi6Mo9d02XMsUFAlubRFaISJ4a/srGILiPyyM2/K2dgVcSGc950lEuObGHEX
- Wch80ZW/8rwl0tqWTrN8DBG052Wrg64VLM41Mc/BABOOHiMLsqxQuORk1MXo
- 78BX8ph8t9IaKrKjZn3+TBA9rI8m6dY4XWRyZvyMDqazZb1p6jbw/NOIF9Ww
- Cn7l7tWpmTa4ZiiPuKOlhai+D5mDFFvUqZ0NbL2gCbcxbQ/vNFtUpeolTRmv
- RKu2DtVU1Q72bW+PJwyswMfCxNrseDuUlAp6p66pY0Q1ztNIicBz+XiaVvxy
- pBznD1HiCFBzuMrqcjoa8k5Vcw8R+LSdcKb00qF0dSDV6TCBw7f8LFT+oiPu
- 2SFbaSIBzYovw4vv0xHCqkx5e4JAwNhvXZx8OjzvckyacwmcVXbVfSmggzm+
- OEx2l0C37Frv5hYaRMrbnc7XETB+sU8tq56GbMYFRm09gdzk6cbBWzTQHRzv
- /dtIoN/n0bOUEhoW/BK93P/xfEeukeKxZBqG2a9vOPUTkNCD+RtIGtiE7bHd
- LwnEqJx9WG5Hw17XZG/pPwTKc3tbOeY0dEfoL7oxTMDsQ8tWHT0a2mqFAsMx
- Aj6h7Sz2LBVLOotNHCbm/y0MPFk2ScWmoSmFaGre/VxC7d5RUUfJKMmY77jq
- kmrutm4qZvX6pVc+E/jRMPKVvIOK9ZZm7s2zBEKz+BelzVRINkrZQwoCjI6G
- 5+r1VFwXNs/MzRGgD9gnl92i4j+M9hoO
- "]]}},
- AspectRatio->0.6180339887498948,
- Axes->{True, True},
- AxesLabel->{None, None},
- AxesOrigin->{0, 0},
- DisplayFunction->Identity,
- Frame->{{False, False}, {False, False}},
- FrameLabel->{{None, None}, {None, None}},
- FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
- GridLines->{None, None},
- GridLinesStyle->Directive[
- GrayLevel[0.5, 0.4]],
- Method->{"DefaultBoundaryStyle" -> Automatic,
- "ScalingFunctions" -> None},
- PlotRange->{{0, 24}, {0., 99.99999320998704}},
- PlotRangeClipping->True,
- PlotRangePadding->{{
- Scaled[0.02],
- Scaled[0.02]}, {
- Scaled[0.05],
- Scaled[0.05]}},
- Ticks->{Automatic, Automatic}]\)
- In[98]:=
- In[99]:=
- In[100]:=
- In[101]:=
- In[102]:=
- General::stop: Further output of Set::setrpt will be suppressed during this calculation. >>
- In[103]:=
- In[148]:= Plot[halftime * (Log[x/hurt] / Log[0.5]), {x, 0, hurt},
- AxesLabel -> Automatic]
- Out[148]= \!\(\*
- GraphicsBox[{{}, {},
- {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
- Opacity[1.], LineBox[CompressedData["
- 1:eJwVznk81PkDx/GZZpRcITQYNcOQn/wyhxnXT5+3Sa5SQnxLjlxpEcWSJc1P
- fh0qZotSkg5ilaPEr34p5VgprbIdsqpt20koO46la7L7++P1eP774kYm+cfM
- otFo2//u/zqIdxqmM5jYHO8R3HfCGUc1L1fKnZgQdLmNxu91xmJNweeqnUxw
- zRo/+yQ5w0fD7uwzNTX4ffTSmyLOkKsvUTobzcaeNNOxrt+dwFazOKiSqKPw
- /XLle2cn1DLPvjTMVEeik4DctnICYXJFS1vUoZf237BufSdEMhb1h3nORVOI
- WcB/RhxRRTe1vhWsAeYFzYnIMkeIvuq35+zQgmpAOTg6zxEVmfSvjFBdnC5h
- JF7QcUDOZOfBn7J0kZPo69pIc0BEYr7p8ZO6+P1z3stn4xKww02dlg7oojhE
- i7H9qQSHpeJUKkQP6awG7s9nJNg1d8vwxfX6CFBcFGi4SLD+2IMn64IMEOM/
- uJq/UwyJ7rEYTpoBWvyHY8a3iWGwP3RyuMgAqZav3nTEivHgu2E92WMD9LRx
- 1p72E8MjjOlbHWiI6203eWaWYogsHdu++htBc7v3Ea2H9tBuKKurXMPC+bg5
- V4tF9kgz+ljcu5UFjf77NuE29niR4Z+jOsTClNnxDQKuPeqhti6gm4Uwfud3
- H3Ts4f/TN5++ehnDPTqkWGdEhOIh8YogqQkkjXop0goRzDnd/UwxGyEmt96y
- uSIc2G3ZbhfIhn7xPx5bG4sw+Sa7ZkMKG13Hf8ty1ROho0Ygq7/Exthn3mUZ
- XYQ4l2OWG5eaoaVAdjHhNyFqgiKTGxYvREXRq6h7FUKc9PgyS6DNQR7P+6ZQ
- KIRWePYJ1kIOKlxD/yi1ESIrnSak2XFQUPlynp6FEBurmJt6/DgIGxE/1zEQ
- gq2h3ZJYyEGsTuKb2CkBTt43y6w25cL9mZtLZrMApYFkwsLGHI/lTJteXwFO
- Rea8NvLkoWqcVn8ih4/i6b01z4N4UHbxubSdfBw5eCi9PJYHhlbS0JZ0PvY1
- FWsK9/CQQwuO9kngY7tGnci3g4cvstGVy4L48Lj8y+7dyy2R2rUh4J0NH6N0
- MU9JrFBZdkeSVGsH1zOD0V2O1gh72znWlboUd/LPlzB32SKCDOkcZtuiP3bT
- 4QqKj+qLxUkCrg3KFoQvrvMVwnM8fih3lzX6y+YuUCXYg6+oPxs8bIWoaxbB
- gjgJ4m6E/i+h0BKlmQmIX++IPL8+r82bePA1jMwdSnVG9nEThydLLbASnazV
- Cf/C7NDXu3SszNEifTZn+tAyKOzVNL6YcDE+pfWg1x+w9HPTD7LiwNYwxr10
- lRs2Ldlv+ahsIbhl89/sD5Ei3yJkxHmIjRtjH+u2UMtRs9Pin6xYUxjUSyrj
- v3GHu8N7sV+fMdpaNVKS01fgXkbEa3kMCyanjDbHbfdAeMbYlyumC2Cr2LEv
- J8cTiz45XZ/uMIRiUr+3QOaFI3GVAd8XGmCKcleWZXhjeUFo2ljqfLRad7v0
- 7vWBscL2WwNvfYhvH6ib3r0SGe/KVQ9d9BDVccaPFK2Ci+Xc69Xmuhg2W/Pj
- /jO+2DGxf+VGgQ5YxD66s2Q1njiGqV6s1QK9O8lu8YU1CHxas0KxQQN35Nv0
- XAP9ML8/+klwoTqaFWl/bn7gB/aFaOnI69nIT16stAhei9obUBax1dDwciL7
- 2aO1aKk1oWvJGGjKjedxQv3BZzTPyrhPh7bOTGjEgD9+uc94NCilQe72gm6+
- PgAuH7iLDAtVRFkd15z6awCO7t1iVJf+iZR/0NG8HRYIQ4dZfZ7npok4KyL5
- 1GAg1h+wZj2smCQTP6iIKnEd8hTNFWZdY4TZ3HRu3cg6HNs4rS8NGSUTDNOF
- JCUIVZnSeeqP3xIxv8TWfDwIblqMOUV3X5NzszdX7NkajLuG0tKCt8/JH/ov
- 7LP+DIbH7NzL9HM/kwLGNio+hQJnT4SN8t/txH3W++rSbymcv3ZUQrdrJx9o
- W1Q96RQSmkwajQfaSMTXiLP2WRROyRsYmyRtRPjR750ql0LogYe3vUZukyej
- fJn8KIWkhEj+9MWbZGG/srLxGoXWqejBu9FXSG9fwqfB6xTuatNT2m81kL1P
- 364yuUkh0kezoJvdQJSPfh3LbqUgD5+/Vf3xJdLa88DF6x6FvkXbZjYerCWx
- HfU9/QMU3vCKoqKCyolJu6259ksKJrXai6x0z5Ge1qpU8opCZpTPD2p3zxDH
- W6eNyxUUGCUpSyRZp4jGdXlU4ujfP/HNMtc1heTmNa2mMiWF7iNXO/Pbvycp
- V/ep945TOJK3tkmYmk8GGmU1kmkKDxWm2VdW5hL5lS8zcR8pZNWlv3dLziYr
- Gnb4l3ym4F4cnjTjnUo+XZosv6+i0MSRNUi9o0ltffL0zAyFH9M4Zs6tu5f9
- BadUDRY=
- "]]}},
- AspectRatio->0.6180339887498948,
- Axes->{True, True},
- AxesLabel->{
- FormBox[
- TagBox["x", HoldForm], TraditionalForm], None},
- AxesOrigin->{0, 0},
- DisplayFunction->Identity,
- Frame->{{False, False}, {False, False}},
- FrameLabel->{{None, None}, {None, None}},
- FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
- GridLines->{None, None},
- GridLinesStyle->Directive[
- GrayLevel[0.5, 0.4]],
- Method->{"DefaultBoundaryStyle" -> Automatic,
- "ScalingFunctions" -> None},
- PlotRange->{{0, 100}, {0., 26.604863600437422`}},
- PlotRangeClipping->True,
- PlotRangePadding->{{
- Scaled[0.02],
- Scaled[0.02]}, {
- Scaled[0.05],
- Scaled[0.05]}},
- Ticks->{Automatic, Automatic}]\)
- y'[x] == y[x]*(Log[0.5]/h - a) + a*d[x],
- y[0] == 0,
- d'[x] == y[x]*
- Log[0.5]/
- h + (2*dos/interval)*(Sin[(x - (interval/4))*2*Pi/interval] + 1),
- d[0] == 0
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement