Advertisement
Guest User

Untitled

a guest
Apr 28th, 2015
211
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 10.78 KB | None | 0 0
  1. In[122]:= Clear[dosage];
  2. dosage[t_, i_, d_] := d /; Mod[t, i] == 0;
  3. dosage[t_, i_, d_] := 0 /; Mod[t, i] != 0;
  4.  
  5. hurt = 100;
  6. effektiv = 90;
  7. halftime = 5;
  8. interval = 5;
  9.  
  10. Clear[fun];
  11. fun[t_, t0_, dose_] := 0 /; t < t0;
  12. fun[t_, t0_, dose_] := dose* (1/2)^((t - t0) /halftime) /; t >= t0;
  13.  
  14. Clear[test];
  15. test[t_, i_, d_] := dosage[t, i, d] + antifun[Mod[t, i], d];
  16.  
  17.  
  18.  
  19. (*DiscretePlot[dosage[t, 4, 200], {t, 0, 20}]
  20. DiscretePlot[dosage[t, 4, 200] + antifun[Mod[t, interval], 200], {t, \
  21. 0, 20}]
  22. Plot[test[t, 4, 200], {t, 0, 20}]*)
  23.  
  24. Clear[antifun];
  25. antifun[t_, h_] := ( h) *(1/2)^(t/halftime) ;
  26.  
  27. Plot[antifun[Mod[t, interval], hurt], {t, 0, 10}]
  28. NMaximize[{antifun[Mod[t, interval], hurt], t >= 0}, t]
  29. NMinimize[{antifun[Mod[t, interval], hurt], t >= 0}, t]
  30.  
  31. (*Plot[NMinimize[{antifun[Mod[t, interval], hurt], t \[GreaterEqual] \
  32. 0, t \[LessEqual] 24}, t], {interval, 0, 24}]*)
  33.  
  34.  
  35.  
  36. Out[136]= \!\(\*
  37. GraphicsBox[{{{}, {},
  38. {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
  39. Opacity[1.], LineBox[CompressedData["
  40. 1:eJwBsQNO/CFib1JlAgAAADoAAAACAAAAnX5iBjFkiz4WMyL0//9YQJ5a7DpX
  41. IGk/r2qOf0f9WEDZTYpyICB5P2ijxFaP+lhAd0dZDgUgiT+8G27oH/VYQEbE
  42. QFz3H5k/sHLOl0TqWECtgjSD8B+pP20Akh+c1FhA4WGuFu0fuT+Nub2Zg6lY
  43. QHtRa2DrH8k/dzZbcDJUWEAJ+ZvJbC7aP0eAVzcmn1dAccJ/jOxy4z/hhu2H
  44. 2vpWQC6kk9G0ruk/1/OcaB1eVkDlTmgc3zjwP9xwo0rPuFVAlxqGHaFg8z+O
  45. jxIX2CJVQDlyuK8DzPY/hhRmzuKEVEDD15tQbyf6P6oirJNL7lNAMYz+PhhJ
  46. /T/upCA+uGVTQEfmOt8wVwBAuU1fRcDVUkDnLbZF9OwBQKRqX0EeU1JAe/yJ
  47. M7x6A0CDnYm+fdZRQAcR6GlUKgVAyfdGLQNTUUAFzQVHC70GQMRSFnK/21BA
  48. +86tbJJxCEA/uInJFl5QQOVXrhkeHgpA4ObNP5PMT0BBiG5tyK0LQM2UzxRM
  49. 805Alf64CUNfDUBx9W4AMA5OQFscw0zc8w5AoKVR01I+TUAMwCvsIlUQQPKQ
  50. qg1sY0xAZTWi9VksEUALrmT8u5JLQHd+eFIg9RFABKmyhnTVSkCF6pPTzs4S
  51. QIAsxK73DUpAhNflyvvREkBKfhq/GQtKQITEN8Io1RJAoDUzIDwISkCCntuw
  52. gtsSQDZKiNSBAkpAflIjjjboEkCodsgEEfdJQHa6skieARNAim3fez7gSUBo
  53. itG9bTQTQB52j4bVsklAZ3cjtZo3E0CnsA2eAbBJQGZkdazHOhNA/tUzBS6t
  54. SUBkPhmbIUETQLTSVMKHp0lAYfJgeNVNE0BZZfT2PpxJQFpa8DI9ZxNAF6ED
  55. QryFSUBMKg+oDJoTQH5nYiLyWElAWchvqCmdE0DYYrUiNlZJQGZm0KhGoBNA
  56. nqKLbnpTSUCAopGpgKYTQNBaoegDTklAtBoUq/SyE0C4zo1lGkNJQBwLGa7c
  57. yxNA/qxOe1UtSUApqXmu+c4TQOJVBjCeKklANkfarhbSE0B9Vr8v5ydJQFCD
  58. m69Q2BNATQIVEHoiSUCE+x2xxOQTQFitbFOjF0lAkZl+seHnE0Detkxf7hRJ
  59. QJ4337H+6hNArHXttTkSSUC4c6CyOPETQIDSUEPRDElAxREBs1X0E0BcUgN6
  60. HQpJQNKvYbNy9xNAMUtW+2kHSUDfTcKzj/oTQBiwQce2BElAp+DlFLz8E0CF
  61. Wp1t1AJJQON3mO8=
  62. "]], LineBox[CompressedData["
  63. 1:eJwVj30w1HkAxjcreeuu42ex3pbvXrenvGXUSS9PbS++kbzs8lMu1luuCzPH
  64. 1NWYG66LVJvKSS5FhblTuCtvu2rTC4rdUCIdd+XEkWqFUm117o9nnvnMM/OZ
  65. eRyjk4Pj9Dgcjngm/3e4A/+pmMtgx5LwwQXvWRiIBquHZth1Yu3oaR0L8+Y1
  66. Ndw5DMpPrH409oaFILa0RmDMoHpnlcDnNQsXPYPa5XMZ/LLvp/PZkyyWFsfX
  67. bp7HIEulSH8wzkLS90XdzzwGyflC411jLGR79tf9Yc2Acd2gah5hkWQ1UnfH
  68. lsG8vqrLvGEW+yXl9YaEwWCqzYm6xywuqxcq09wYOBzIvaXrZnF7u1xZsIhB
  69. YVOvJqCLRbfhc2WtFwO5rcDuTCcL7ZqqBq0PAw83ecE6NQvhFY8rsb4MOk8Y
  70. GeVdYyGvWNzoL2OQ3ZvbsPICi/lRIllFLIOAHEVOxG8sVGZ8vU8SGDw0iFLu
  71. Lpvx7Xovbk+a2VUG5jXFM39W32wJSmPw9fziDrc8FrY9Qe2h+Qzqd7ukev3A
  72. 4gIn8a8oDQOe0VpeRiCL1pBzulRvC1za03Pc8lUYGjXJhxZyeVhaLt33IDAM
  73. F7dnQN3PQ9KtyKtNZ0Nh77fyq6RqSwzg3HSHTgrxx4i9RWlWGDd3erIpXorY
  74. 7qsv66XW6HIITCxpkWBFxo81MgEfeQdL+BVLJDiWPN5gOsyHgK9+FlIZgoH0
  75. 0qEilQ383Z6v/8YuBE/jFLHXs2yx7Iwk81VBMLSJZnM0MjtwxDv3pnwajEra
  76. pot3scdz98KzjUeDMJTgvpg7ZY9JpZzTbR0E/vmHuVkaB3TUto9EngzE3qjv
  77. zmr8BOj/3nlRqSgQr3XWcj+FANIFBpdOR26CqCG/oFXoiENdUn3v3QHwimna
  78. 9Fm2I3o/cN3Fv25EY6GFlXrKEdroa5Wz1P4QNVocCNjihHeSUPXn+v7YOHvM
  79. 7q7KCacsfZL2rPWD9s/BG7edCUzWFywPS9uA5EK5eXYOQSDxxkAThXqW8MWG
  80. owQ7/k5vU12ncN6mbDXNJYi+W1R/8irFoMdwxpHjBH01nvWhCgq2Bdr8UwSz
  81. PeP9Oi9QrH45oS47T1DVfKTuQS4Fz3dz5s0Wgi3EXuIeTZFSoZVl3iaIJCFc
  82. 60iKDrOs5b5tBGUDg2Z6ERQH+y9Ntt0hCPj9RfB9KQUnZW7MvfsEoqDN0nRK
  83. MXqqceXAPwSKonXh/7pTrOeG2ZQ8IfDOaTXscaEoSXj2Km6Y4JhIa9nsTBHp
  84. ya8cGSUQN41LSoUUXbdSbMfHCabG9vsmWFF4uBpPX5wgqN7WNRFuQXE4t/he
  85. 6hTBu+zYD35mFHSr5sD0NMH9rUIPD1OKshsx8cq3BJ2Pl4ULjSj0v3y7Kk1H
  86. YFvx7SNLAwrZ4SN2Kz4QJC8sbTbhUqgm5r/5+JHghbrfhMOh+A+axtZA
  87.  
  88. "]]}, {}}, {{}, {}, {}}},
  89. AspectRatio->0.6180339887498948,
  90. Axes->{True, True},
  91. AxesLabel->{None, None},
  92. AxesOrigin->{0, 48.},
  93. DisplayFunction->Identity,
  94. Frame->{{False, False}, {False, False}},
  95. FrameLabel->{{None, None}, {None, None}},
  96. FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  97. GridLines->{None, None},
  98. GridLinesStyle->Directive[
  99. GrayLevel[0.5, 0.4]],
  100. Method->{"DefaultBoundaryStyle" -> Automatic,
  101. "ScalingFunctions" -> None},
  102. PlotRange->{{0, 10}, {50.0000014145861, 99.99999717082787}},
  103. PlotRangeClipping->True,
  104. PlotRangePadding->{{
  105. Scaled[0.02],
  106. Scaled[0.02]}, {
  107. Scaled[0.05],
  108. Scaled[0.05]}},
  109. Ticks->{Automatic, Automatic}]\)
  110.  
  111. Out[137]= {100., {t -> 2.28154*10^-9}}
  112.  
  113. Out[138]= {50., {t -> 5.}}
  114.  
  115. In[94]:=
  116.  
  117.  
  118.  
  119.  
  120.  
  121. In[96]:=
  122.  
  123. In[139]:=
  124. Plot[antifun[t, hurt], {t, 0, 24}]
  125.  
  126. Out[139]= \!\(\*
  127. GraphicsBox[{{}, {},
  128. {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
  129. Opacity[1.], LineBox[CompressedData["
  130. 1:eJwVknk01AkAx5WixsyUZISxmBkTkdy/ssvvO6YDreeoaWaSDkde7mVyTJsc
  131. jWrDilRYR/VIK9tBUiKFDUm02xM5shmd7pVKk7V/fN73ff78vvcx9A33Clio
  132. pKTkNs//+2lrWL5HeJFj1BbZq7k5EcQqSRpPuPHk7tKyodjPIhRlkDvSuOmk
  133. WFM3d2pShAZBQUwIt4B0SaUZjMpFEFS7LdnF/YNMqHl6cuaRCAUHiju2cevI
  134. lKrHR7zPieDk1Fnpzm0nOx6rPfS3EuG0+esSN24/uYJYV+AfLARl9+jEm6Ax
  135. kneJy8npFOAYc//BmztmyISBE/YuIduh3OTaXNz5jewRmxwy52yDRhODFRmj
  136. DK/kgWVJXZ6o+lC04z5UUbr+stpPNz1whfV7IeGthmnrqzqV+u7QtxvNkvku
  137. Q9vAijMSFTd4Hgu4x8tUR2j0Ts/FVltRK9Ug/87TwLDkh+l9Qhdc76DdqUvS
  138. hGKWTxy/vQX+2cUUzb1aeCNMmTzosBkY/Fn+LUAbnyi3j47VbETxr0U1cNEF
  139. 2zH0ckswHzZNykyGQA+bUktygrSckHhgkZdQqA/2kaSlkwOAaXUpLzbPAKXd
  140. zoraIhKmc2GF6YWGsFBwR0ylDgiPefHFI5OFdnkvXy/3ezhXBOr65LDRQ4bE
  141. iis3IKHkAbcwj4PyKXXe6AyBr60N69JSjaDz3p8RLrEDrWJI0XiKi+CEe19f
  142. v7eBb1tZS6lsNQSDteKSbGsw8nkJDyTGMHtq07jT2wp7LGdmjWJNwNTRrDnP
  143. tMTIR6+DtkFrEMgzaLmjYoFlCe0V6RGmqF4wXdnRsxan6lf16UWaIaI1bLWi
  144. yxRCq6WLxL5rQYvwfa+jtQY3ArXH+/abI+SR834bD2N4tj/ppu9ch+ie9vCM
  145. Bi6Mo9d02XMsUFAlubRFaISJ4a/srGILiPyyM2/K2dgVcSGc950lEuObGHEX
  146. Wch80ZW/8rwl0tqWTrN8DBG052Wrg64VLM41Mc/BABOOHiMLsqxQuORk1MXo
  147. 78BX8ph8t9IaKrKjZn3+TBA9rI8m6dY4XWRyZvyMDqazZb1p6jbw/NOIF9Ww
  148. Cn7l7tWpmTa4ZiiPuKOlhai+D5mDFFvUqZ0NbL2gCbcxbQ/vNFtUpeolTRmv
  149. RKu2DtVU1Q72bW+PJwyswMfCxNrseDuUlAp6p66pY0Q1ztNIicBz+XiaVvxy
  150. pBznD1HiCFBzuMrqcjoa8k5Vcw8R+LSdcKb00qF0dSDV6TCBw7f8LFT+oiPu
  151. 2SFbaSIBzYovw4vv0xHCqkx5e4JAwNhvXZx8OjzvckyacwmcVXbVfSmggzm+
  152. OEx2l0C37Frv5hYaRMrbnc7XETB+sU8tq56GbMYFRm09gdzk6cbBWzTQHRzv
  153. /dtIoN/n0bOUEhoW/BK93P/xfEeukeKxZBqG2a9vOPUTkNCD+RtIGtiE7bHd
  154. LwnEqJx9WG5Hw17XZG/pPwTKc3tbOeY0dEfoL7oxTMDsQ8tWHT0a2mqFAsMx
  155. Aj6h7Sz2LBVLOotNHCbm/y0MPFk2ScWmoSmFaGre/VxC7d5RUUfJKMmY77jq
  156. kmrutm4qZvX6pVc+E/jRMPKVvIOK9ZZm7s2zBEKz+BelzVRINkrZQwoCjI6G
  157. 5+r1VFwXNs/MzRGgD9gnl92i4j+M9hoO
  158. "]]}},
  159. AspectRatio->0.6180339887498948,
  160. Axes->{True, True},
  161. AxesLabel->{None, None},
  162. AxesOrigin->{0, 0},
  163. DisplayFunction->Identity,
  164. Frame->{{False, False}, {False, False}},
  165. FrameLabel->{{None, None}, {None, None}},
  166. FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  167. GridLines->{None, None},
  168. GridLinesStyle->Directive[
  169. GrayLevel[0.5, 0.4]],
  170. Method->{"DefaultBoundaryStyle" -> Automatic,
  171. "ScalingFunctions" -> None},
  172. PlotRange->{{0, 24}, {0., 99.99999320998704}},
  173. PlotRangeClipping->True,
  174. PlotRangePadding->{{
  175. Scaled[0.02],
  176. Scaled[0.02]}, {
  177. Scaled[0.05],
  178. Scaled[0.05]}},
  179. Ticks->{Automatic, Automatic}]\)
  180.  
  181. In[98]:=
  182.  
  183. In[99]:=
  184.  
  185. In[100]:=
  186.  
  187. In[101]:=
  188.  
  189.  
  190. In[102]:=
  191.  
  192. General::stop: Further output of Set::setrpt will be suppressed during this calculation. >>
  193.  
  194. In[103]:=
  195.  
  196. In[148]:= Plot[halftime * (Log[x/hurt] / Log[0.5]), {x, 0, hurt},
  197. AxesLabel -> Automatic]
  198.  
  199. Out[148]= \!\(\*
  200. GraphicsBox[{{}, {},
  201. {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
  202. Opacity[1.], LineBox[CompressedData["
  203. 1:eJwVznk81PkDx/GZZpRcITQYNcOQn/wyhxnXT5+3Sa5SQnxLjlxpEcWSJc1P
  204. fh0qZotSkg5ilaPEr34p5VgprbIdsqpt20koO46la7L7++P1eP774kYm+cfM
  205. otFo2//u/zqIdxqmM5jYHO8R3HfCGUc1L1fKnZgQdLmNxu91xmJNweeqnUxw
  206. zRo/+yQ5w0fD7uwzNTX4ffTSmyLOkKsvUTobzcaeNNOxrt+dwFazOKiSqKPw
  207. /XLle2cn1DLPvjTMVEeik4DctnICYXJFS1vUoZf237BufSdEMhb1h3nORVOI
  208. WcB/RhxRRTe1vhWsAeYFzYnIMkeIvuq35+zQgmpAOTg6zxEVmfSvjFBdnC5h
  209. JF7QcUDOZOfBn7J0kZPo69pIc0BEYr7p8ZO6+P1z3stn4xKww02dlg7oojhE
  210. i7H9qQSHpeJUKkQP6awG7s9nJNg1d8vwxfX6CFBcFGi4SLD+2IMn64IMEOM/
  211. uJq/UwyJ7rEYTpoBWvyHY8a3iWGwP3RyuMgAqZav3nTEivHgu2E92WMD9LRx
  212. 1p72E8MjjOlbHWiI6203eWaWYogsHdu++htBc7v3Ea2H9tBuKKurXMPC+bg5
  213. V4tF9kgz+ljcu5UFjf77NuE29niR4Z+jOsTClNnxDQKuPeqhti6gm4Uwfud3
  214. H3Ts4f/TN5++ehnDPTqkWGdEhOIh8YogqQkkjXop0goRzDnd/UwxGyEmt96y
  215. uSIc2G3ZbhfIhn7xPx5bG4sw+Sa7ZkMKG13Hf8ty1ROho0Ygq7/Exthn3mUZ
  216. XYQ4l2OWG5eaoaVAdjHhNyFqgiKTGxYvREXRq6h7FUKc9PgyS6DNQR7P+6ZQ
  217. KIRWePYJ1kIOKlxD/yi1ESIrnSak2XFQUPlynp6FEBurmJt6/DgIGxE/1zEQ
  218. gq2h3ZJYyEGsTuKb2CkBTt43y6w25cL9mZtLZrMApYFkwsLGHI/lTJteXwFO
  219. Rea8NvLkoWqcVn8ih4/i6b01z4N4UHbxubSdfBw5eCi9PJYHhlbS0JZ0PvY1
  220. FWsK9/CQQwuO9kngY7tGnci3g4cvstGVy4L48Lj8y+7dyy2R2rUh4J0NH6N0
  221. MU9JrFBZdkeSVGsH1zOD0V2O1gh72znWlboUd/LPlzB32SKCDOkcZtuiP3bT
  222. 4QqKj+qLxUkCrg3KFoQvrvMVwnM8fih3lzX6y+YuUCXYg6+oPxs8bIWoaxbB
  223. gjgJ4m6E/i+h0BKlmQmIX++IPL8+r82bePA1jMwdSnVG9nEThydLLbASnazV
  224. Cf/C7NDXu3SszNEifTZn+tAyKOzVNL6YcDE+pfWg1x+w9HPTD7LiwNYwxr10
  225. lRs2Ldlv+ahsIbhl89/sD5Ei3yJkxHmIjRtjH+u2UMtRs9Pin6xYUxjUSyrj
  226. v3GHu8N7sV+fMdpaNVKS01fgXkbEa3kMCyanjDbHbfdAeMbYlyumC2Cr2LEv
  227. J8cTiz45XZ/uMIRiUr+3QOaFI3GVAd8XGmCKcleWZXhjeUFo2ljqfLRad7v0
  228. 7vWBscL2WwNvfYhvH6ib3r0SGe/KVQ9d9BDVccaPFK2Ci+Xc69Xmuhg2W/Pj
  229. /jO+2DGxf+VGgQ5YxD66s2Q1njiGqV6s1QK9O8lu8YU1CHxas0KxQQN35Nv0
  230. XAP9ML8/+klwoTqaFWl/bn7gB/aFaOnI69nIT16stAhei9obUBax1dDwciL7
  231. 2aO1aKk1oWvJGGjKjedxQv3BZzTPyrhPh7bOTGjEgD9+uc94NCilQe72gm6+
  232. PgAuH7iLDAtVRFkd15z6awCO7t1iVJf+iZR/0NG8HRYIQ4dZfZ7npok4KyL5
  233. 1GAg1h+wZj2smCQTP6iIKnEd8hTNFWZdY4TZ3HRu3cg6HNs4rS8NGSUTDNOF
  234. JCUIVZnSeeqP3xIxv8TWfDwIblqMOUV3X5NzszdX7NkajLuG0tKCt8/JH/ov
  235. 7LP+DIbH7NzL9HM/kwLGNio+hQJnT4SN8t/txH3W++rSbymcv3ZUQrdrJx9o
  236. W1Q96RQSmkwajQfaSMTXiLP2WRROyRsYmyRtRPjR750ql0LogYe3vUZukyej
  237. fJn8KIWkhEj+9MWbZGG/srLxGoXWqejBu9FXSG9fwqfB6xTuatNT2m81kL1P
  238. 364yuUkh0kezoJvdQJSPfh3LbqUgD5+/Vf3xJdLa88DF6x6FvkXbZjYerCWx
  239. HfU9/QMU3vCKoqKCyolJu6259ksKJrXai6x0z5Ge1qpU8opCZpTPD2p3zxDH
  240. W6eNyxUUGCUpSyRZp4jGdXlU4ujfP/HNMtc1heTmNa2mMiWF7iNXO/Pbvycp
  241. V/ep945TOJK3tkmYmk8GGmU1kmkKDxWm2VdW5hL5lS8zcR8pZNWlv3dLziYr
  242. Gnb4l3ym4F4cnjTjnUo+XZosv6+i0MSRNUi9o0ltffL0zAyFH9M4Zs6tu5f9
  243. BadUDRY=
  244. "]]}},
  245. AspectRatio->0.6180339887498948,
  246. Axes->{True, True},
  247. AxesLabel->{
  248. FormBox[
  249. TagBox["x", HoldForm], TraditionalForm], None},
  250. AxesOrigin->{0, 0},
  251. DisplayFunction->Identity,
  252. Frame->{{False, False}, {False, False}},
  253. FrameLabel->{{None, None}, {None, None}},
  254. FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  255. GridLines->{None, None},
  256. GridLinesStyle->Directive[
  257. GrayLevel[0.5, 0.4]],
  258. Method->{"DefaultBoundaryStyle" -> Automatic,
  259. "ScalingFunctions" -> None},
  260. PlotRange->{{0, 100}, {0., 26.604863600437422`}},
  261. PlotRangeClipping->True,
  262. PlotRangePadding->{{
  263. Scaled[0.02],
  264. Scaled[0.02]}, {
  265. Scaled[0.05],
  266. Scaled[0.05]}},
  267. Ticks->{Automatic, Automatic}]\)
  268.  
  269. y'[x] == y[x]*(Log[0.5]/h - a) + a*d[x],
  270. y[0] == 0,
  271. d'[x] == y[x]*
  272. Log[0.5]/
  273. h + (2*dos/interval)*(Sin[(x - (interval/4))*2*Pi/interval] + 1),
  274. d[0] == 0
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement