A two-compartment pharmacokinetic model
restart:with(DEtools):
The non-homogeneous linear system for the model. d(t) is the drug dosage versus time
system1:=diff(x(t),t)=-a*x(t)+d(t),diff(y(t),t)=a*x(t)-b*y(t);
NiQvLUklZGlmZkclKnByb3RlY3RlZEc2JC1JInhHNiI2I0kidEdGKkYsLCYqJkkiYUdGKiIiIkYoRjAhIiItSSJkR0YqRitGMC8tRiU2JC1JInlHRipGK0YsLCZGLkYwKiZJImJHRipGMEY3RjBGMQ==
The homogeneous version of the same system
system2:=diff(x(t),t)=-a*x(t),diff(y(t),t)=a*x(t)-b*y(t);
NiQvLUklZGlmZkclKnByb3RlY3RlZEc2JC1JInhHNiI2I0kidEdGKkYsLCQqJkkiYUdGKiIiIkYoRjAhIiIvLUYlNiQtSSJ5R0YqRitGLCwmRi5GMComSSJiR0YqRjBGNUYwRjE=
The solution of the homogeneous system. Note that the eigenvalues are -a, -b , so as time increases, the solution of the homogeneous system will converge to zero
dsolve([system2],{x(t),y(t)});
PCQvLUkieEc2IjYjSSJ0R0YmKiZJJF9DMkdGJiIiIi1JJGV4cEc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJjYjLCQqJkkiYUdGJkYrRihGKyEiIkYrLy1JInlHRiZGJywmKiZJJF9DMUdGJkYrLUYtNiMsJComSSJiR0YmRitGKEYrRjVGK0YrKiosJkY0RitGQEY1RjVGNEYrRixGK0YqRitGNQ==
Let us first consider the example of a single dosage of the drug administered from t=0 to t=2:
a:=2;b:=1.5;d:=t->1-Heaviside(t-2);
IiIj
JCIjOiEiIg==
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYiIiJGKi1JKkhlYXZpc2lkZUdGJTYjLCY5JEYqISIjRiohIiJGJUYlRiU=
plot(d(t),t=0..5,0..2,scaling=constrained);
NictJSdDVVJWRVNHNiQ3aG43JCQiIiEhIiIkIiM1ISIiNyQkIjJMTEwkM3gmKSozIiEjPCQiIzUhIiI3JCQiMm5tO0gyUCJRPyEjPCQiIzUhIiI3JCQiMkxMTGVSd1g1JCEjPCQiIzUhIiI3JCQiMU1MJDN4JTN5VCEjOyQiIzUhIiI3JCQiMW5tInolNFxZXyEjOyQiIzUhIiI3JCQiMU1MZVItL1BpISM7JCIjNSEiIjckJCItRGNtcGlzISM3JCIjNSEiIjckJCIxTUxlKik+VkIkKSEjOyQiIzUhIiI3JCQiLURKYnchUSohIzckIiM1ISIiNyQkIjJvbW1USU9vLyIhIzskIiM1ISIiNyQkIjJOTCQzXz5qVTYhIzskIiM1ISIiNyQkIjIsKytEO3YvRCIhIzskIiM1ISIiNyQkIix2PWgoZTghIzUkIiM1ISIiNyQkIix2JFs2ajkhIzUkIiM1ISIiNyQkIjJNTGUqW3ooeWIiISM7JCIjNSEiIjckJCIyb21tVFhnMG4iISM7JCIjNSEiIjckJCIyb21tbUo8Z3ciISM7JCIjNSEiIjckJCIuRDFNY3EoPSEjNyQiIzUhIiI3JCQiMkwkZTk7MD9FPiEjOyQiIzUhIiI3JCQiMm1tbTtwV2AoPiEjOyQiIzUhIiI3JCQiMkNjXi5jJTMjKT4hIzskIiM1ISIiNyQkIjImZWsuSFcjKSkpPiEjOyQiIzUhIiI3JCQiMmwhKnlMTyU+Iyo+ISM7JCIjNSEiIjckJCIyV05AeEhrYio+ISM7JCIjNSEiIjckJCIyQSFRMUtVJCopKj4hIzskIiM1ISIiNyQkIjBEMWs7L0IrIyEjOSQiIiEhIiI3JCQiMVVneC5SeTo/ISM6JCIiISEiIjckJCIyTyRlOVRPRUg/ISM7JCIiISEiIjckJCIxPGEpZTZCaTAjISM6JCIiISEiIjckJCIyLitEMWYjPSQzIyEjOyQiIiEhIiI3JCQiLUQiPUVYOCMhIzYkIiIhISIiNyQkIi52PXhwZT0jISM3JCIiISEiIjckJCIxbm0iSDI4SUgjISM6JCIiISEiIjckJCIyb20ienBTUyJSIyEjOyQiIiEhIiI3JCQiMkxMJDNfP2AoXCMhIzskIiIhISIiNyQkIjJPTGUqKT5weGcjISM7JCIiISEiIjckJCIudiRmNHQuRiEjNyQiIiEhIiI3JCQiMlBMJGUqR3N0IUchIzskIiIhISIiNyQkIitEUlc5SCEiKiQiIiEhIiI3JCQiLURKRT4+SSEjNiQiIiEhIiI3JCQiMi4rRDFSVTA3JCEjOyQiIiEhIiI3JCQiLXY9UzJMSyEjNiQiIiEhIiI3JCQiMmptbTtwKT1NTCEjOyQiIiEhIiI3JCQiMi0rK3Y9XUBXJCEjOyQiIiEhIiI3JCQiMUwkZSpbJHoqUk4hIzokIiIhISIiNyQkIjIxKytEWUtwayQhIzskIiIhISIiNyQkIjJsbSJIMnFjWlAhIzskIiIhISIiNyQkIjIvK0RKNWZGJlEhIzskIiIhISIiNyQkIjJsbW1UZy5jJlIhIzskIiIhISIiNyQkIi5EY0VzSzElISM3JCIiISEiIjckJCIxTUxMJCkqcHA7JSEjOiQiIiEhIiI3JCQiMU1MM3hlLHRVISM6JCIiISEiIjckJCIxbjtIZE89eVYhIzokIiIhISIiNyQkIjIwKytdIz4jW1olISM7JCIiISEiIjckJCIxbm1UJkchZSZlJSEjOiQiIiEhIiI3JCQiMU1MTCQpUWslbyUhIzokIiIhISIiNyQkIi5EMU1tLXolISM3JCIiISEiIjckJCIyMSt2JDQwTyIqWyEjOyQiIiEhIiI3JCQiI10hIiIkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQiIiEhIiIkIiNdISIiOyQiIiEhIiIkIiM/ISIiLSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRInQ2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lKFNDQUxJTkdHNiMlLENPTlNUUkFJTkVERy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEklISIiLSUpQk9VTkRTX1lHNiMkIiUrOCEiIi0lLUJPVU5EU19XSURUSEc2IyQiJSFcJCEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiUrOSEiIi0lKUNISUxEUkVORzYi
ics:=x(0)=0,y(0)=0;dsolve([system1,ics]);
NiQvLUkieEc2IjYjIiIhRigvLUkieUdGJkYnRig=
PCQvLUkieEc2IjYjSSJ0R0YmLCotSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiY2IywkRighIiMjISIiIiIjIyIiIkY0RjYtSSpIZWF2aXNpZGVHRiw2IywmRihGNkYxRjZGMiomRjdGNi1GKzYjLCZGKEYxIiIlRjZGNkY1Ly1JInlHRiZGJywuLUYrNiMsJEYoIyEiJEY0IyEiKSIiJEYqRjQjRjRGS0Y2RjcjRjFGS0Y7RjEqJkY3RjYtRis2IywmRihGR0ZLRjZGNiMiIilGSw==
plot(-(8/3)*exp(-(3/2)*t)+2*exp(-2*t)+2/3-(2/3)*Heaviside(t-2)-2*Heaviside(t-2)*exp(-2*t+4)+(8/3)*Heaviside(t-2)*exp(-(3/2)*t+3),t=0..10);
NiYtJSdDVVJWRVNHNiQ3YG83JCQiIiEhIiIkIjJsOkRZLUItNiIhI0s3JCQiMkxMTCQzeCYpKjMiISM8JCIyY3YkM0AqKTRaNSEjPTckJCIybW1tbVQ6KHpAISM8JCIxaiYzaDJtLnEkISM8NyQkIi1EInklKno3JCEjNyQiMShvSEZsaVMmbyEjPDckJCIyTkxMZTl1aTIlISM8JCIxRzVJTXBxWzUhIzs3JCQiLXZvTXJVXiEjNyQiMik0Q0hqUmMoWyIhIzw3JCQiMW5tbSJ6XyI0aSEjOyQiMiN6cVgxKm9uJD4hIzw3JCQiMW5tbW02bSNHKCEjOyQiMkc/UmFYUUtRIyEjPDckJCIxb21tVCZwaE4pISM7JCIyRVp6bUhqSCJHISM8NyQkIjEsK3Y9ZGRDJSohIzskIjI7J2YrJHp0cEAkISM8NyQkIjJMTCRlKj0pSFw1ISM7JCIxOT5PRSc+S2YkISM7NyQkIi12PUpOWzYhIzYkIjFhVEBBZlU6UiEjOzckJCIybm07ei8zdUMiISM7JCIyL1c7OShHXTZVISM8NyQkIjJNTCRlKm90Klw4ISM7JCIyYV0/eiJIKDNcJSEjPDckJCItREokUkRYIiEjNiQiMTpPQCU9Qk51JSEjOzckJCIyb207elInb2s7ISM7JCIxJT1NSyNIYCg9JiEjOzckJCItRDFKOnc9ISM2JCIxUj4wbnlDUGIhIzs3JCQiMk0kM3giKUhgST4hIzskIjBKNG1gOlRoJiEjOjckJCIxbjtIZEciXCk+ISM6JCIxIzMlb0J3OidvJiEjOzckJCIyTDNfXXotQCwjISM7JCIxPjEhKTRaKz5kISM7NyQkIi9EIkd0I0hSPyEjOCQiMV88WlYleilRZCEjOzckJCIxPEhkcUVbbT8hIzokIjE+c1shKjQiW3UmISM7NyQkIjJQTEwkM0VuJDQjISM7JCIxaF9kQkQwUWQhIzs3JCQiMjArXWlEbyUqPSMhIzskIjE7OyxiXFxIYyEjOzckJCIxbm07L1JFJkcjISM6JCIxNnYuNmtbP2EhIzs3JCQiMk9MJGU5cjUkUiMhIzskIjFXV1dPIkdrNSYhIzs3JCQiMi4rK11LXTRdIyEjOyQiMkZWaEV0T0p1JSEjPDckJCIqTk8jNEUhIikkIjIoKT1CeG81XU4lISM8NyQkIit2Ql88RiEiKiQiMiczdUFlNT1pUiEjPDckJCIrRGcoPSNHISIqJCIyYW5EJno+SiJmJCEjPDckJCIrdidIaSNIISIqJCIxeTk2ZiEpKmVCJCEjOzckJCIyUExla3kjKjQtJCEjOyQiMW9ccy1ndklIISM7NyQkIjFubSJ6KmV2OkohIzokIjJsVG9bRy1aayMhIzw3JCQiLkRKU1ElR0shIzckIjEqPmQ2cVg0TCMhIzs3JCQiMk5MTCQzNDdUTCEjOyQiMmpzIjQ7UUlZPyEjPDckJCIxTUwkM3h4bFYkISM6JCIyayZwPnJAY0Y9ISM8NyQkIjJQTExMak0/YCQhIzskIjJQOSVRJylRZEc7ISM8NyQkIi1EIm83VHYkISM2JCIyMWQyTSZRQk83ISM8NyQkIjFMTEwkUSpvXVIhIzokIjFEemhhL0g6JyohIzw3JCQiMSwrRCI9bGo7JSEjOiQiMG1wTTo+N0QoISM7NyQkIi12ViZSPFAlISM2JCIyOTUsP2AwT14mISM9NyQkIjFNTCRlOUVnZSUhIzokIjI7MVNEaFVSNyUhIz03JCQiMk5MJGVSIjNHeSUhIzskIjI5JFwjb0tudjkkISM9NyQkIjFubTsvVDEmKlwhIzokIjJ3O21zJilIVk0jISM9NyQkIjFubSJ6UlFiQCYhIzokIjE9JWZhZF41cyIhIzw3JCQiLXY9PlkyYSEjNiQiMDRqKmVXQTc4ISM7NyQkIjJ2bW0ielh1OWMhIzskIjFeKHo5ciw1eCohIz03JCQiKiZ5KSlHZSEiKSQiMXM0bSc0KlIiPighIz03JCQiLERFJlFRZyEjNSQiMWFHOlkvXz5gISM9NyQkIjIwK103eSUzVGkhIzskIjEjUkREW1wjb1IhIz03JCQiLHYuW2hZJyEjNSQiMi5VKD5JWztpRyEjPjckJCIxTExMJFF4JG9tISM6JCIyIzQhR1lOdjs4IyEjPjckJCIyMCsrXVArVilvISM7JCIyLThzUmw/WmIiISM+NyQkIjFtbSJ6cGUqenEhIzokIjFkIj5GKHk0bjYhIz03JCQiMSwrK0RcJ1FIKCEjOiQiMSZvb1s8JFxCJikhIz43JCQiMUxMZTlTOCZcKCEjOiQiMVpXKlwtMnVMJyEjPjckJCIxLCtEMSM9YnEoISM6JCIxXjVoTyE+a2slISM+NyQkIjFMTEwzcz82eiEjOiQiMSYzYDZYNCZHTSEjPjckJCItREpYYUUiKSEjNiQiMnUzZUNZJHkjXCMhIz83JCQiMW9tbW0qUlJMKSEjOiQiMUNcMGdnM0w9ISM+NyQkIjFvbTthPC5ZJikhIzokIjIxIj4tMUk1UTghIz83JCQiMU5MZTl0T2MoKSEjOiQiMXVxaTpfKSkqeSohIz83JCQiMSwrK11Ra1wqKSEjOiQiMVFPSm8kZVZNKCEjPzckJCIxTkwkM2RnNjwqISM6JCIxejRdX01oIkcmISM/NyQkIjFvbW1teEdwJCohIzokIjFccyMpeWMlPSRSISM/NyQkIi1EIm9LMGUqISM2JCIyWW1DKD4lZihwRyEjQDckJCIxLCt2PTVzI3kqISM6JCIyL01JQ3hLRTcjISNANyQkIiQrIiEiIiQiMXcrP3J5eU06ISM/LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQiIiEhIiIkIiQrIiEiIiUoREVGQVVMVEctJStBWEVTTEFCRUxTRzYnLUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHNiI2NVEidDYiLyUnZmFtaWx5R1EhNiIvJSVzaXplR1EjMTA2Ii8lJWJvbGRHUSZmYWxzZTYiLyUnaXRhbGljR1EldHJ1ZTYiLyUqdW5kZXJsaW5lR1EmZmFsc2U2Ii8lKnN1YnNjcmlwdEdRJmZhbHNlNiIvJSxzdXBlcnNjcmlwdEdRJmZhbHNlNiIvJStmb3JlZ3JvdW5kR1EoWzAsMCwwXTYiLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV02Ii8lJ29wYXF1ZUdRJmZhbHNlNiIvJStleGVjdXRhYmxlR1EmZmFsc2U2Ii8lKXJlYWRvbmx5R1EmZmFsc2U2Ii8lKWNvbXBvc2VkR1EmZmFsc2U2Ii8lKmNvbnZlcnRlZEdRJmZhbHNlNiIvJStpbXNlbGVjdGVkR1EmZmFsc2U2Ii8lLHBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lNnNlbGVjdGlvbi1wbGFjZWhvbGRlckdRJmZhbHNlNiIvJSxtYXRodmFyaWFudEdRJ2l0YWxpYzYiUSE2Ii0lJUZPTlRHNiUlKERFRkFVTFRHJShERUZBVUxURyIjNSUrSE9SSVpPTlRBTEclK0hPUklaT05UQUxHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkSSUhIiItJSlCT1VORFNfWUc2IyQiI10hIiItJS1CT1VORFNfV0lEVEhHNiMkIiVnTSEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiU/TiEiIi0lKUNISUxEUkVORzYi
In the next example, the drug is administered periodically:
d:=t->0.3*(cos(t)+1);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiYkIiIkISIiIiIiLCYtSSRjb3NHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IzkkRi1GLUYtRi1GJUYlRiU=
dsolve([system1,ics]);
PCQvLUkieEc2IjYjSSJ0R0YmLCotSSRjb3NHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiZGJyMiIiQiI0QtSSRzaW5HRixGJyNGMCIjXSNGMCIjPyIiIi1JJGV4cEdGLDYjLCRGKCEiIyMhI0YiJCsiLy1JInlHRiZGJywsRiojIiNDIiREJEYyIyIjVUZHI0Y4IiImRjhGOSMiI0ZGMS1GOjYjLCRGKCMhIiQiIiMjISMpKSIjbA==
plot((16/65)*cos(t)+(28/65)*sin(t)+2/3+(18/5)*exp(-2*t)-(176/39)*exp(-(3/2)*t),t=0..40);
NiYtJSdDVVJWRVNHNiQ3YXk3JCQiIiEhIiIkIiIhISIiNyQkIjJMTEwkM3gmKSozIiEjPCQiMkIoKSlwQCIzSjQjISM9NyQkIjJtbW1tVDooekAhIzwkIjEjKT5bLjljJVEoISM8NyQkIitESmRwSyEjNSQiMj0keTxcbCJvWSIhIzw3JCQiMUxMTEwzVmZWISM7JCIybFBcJ2ZIJ1FJIyEjPDckJCIxbW1tVCYpR1xhISM7JCIyYF5jL1FkQj0kISM8NyQkIipEWSJSbCEiKiQiMkNxbjE4W0kwJSEjPDckJCIxTUxMZVIrSHchIzskIjFcOSlHO1QwKVshIzs3JCQiMW1tbW07Jyk9KCkhIzskIjFbb3NJXTRTYyEjOzckJCIyTEwkZVI/YWg1ISM7JCIxSjFGQDQ2YW4hIzs3JCQiLEQiej5eNyEjNSQiMVZ6InpVI1FwdiEjOzckJCIybG07YXlgM1ciISM7JCIxYiplJUcubHchKSEjOzckJCIyTUxMJGUnNDBqIiEjOyQiMG15ZVlfJipHKSEjOjckJCIybjt6V2lKUW8iISM7JCIxKTRVNnVhKipIKSEjOzckJCIyLStEMWZgcnQiISM7JCIxVl8jRzspKioqRykhIzs3JCQiMk4kM3hjYloheiIhIzskIjFUSEthLlxnIykhIzs3JCQiMm9tO0hfKHpWPSEjOyQiMWtIP3p6RjcjKSEjOzckJCIyTEwzX1hULyY+ISM7JCIxKEcjR08pKUdqISkhIzs3JCQiLHZRJjNkPyEjNSQiMVFZaUchKlFdeSEjOzckJCIxTEwzX0tQcUEhIzokIjEiXEwhKnkwUUUoISM7NyQkIjJubW1tNmhPWyMhIzskIjEkZSpwX1YrO2whIzs3JCQiMW5tbSJ6aSQpcCMhIzokIjEpUSEpPXRNYG0mISM7NyQkIjFubW1tVzE4SCEjOiQiMTM7Q3E8LiF5JSEjOzckJCIxbm1tVGh3RkohIzokIjEqb1VqUlIsI1IhIzs3JCQiMnRtbW0ieVlVTCEjOyQiMlolSGthXTxTSiEjPDckJCIyMCsrdkdJKXBQISM7JCIyUSFcciU+Y1grIyEjPDckJCIxTExMZUY+KD4lISM6JCIyOipcL2p6Qkk7ISM8NyQkIit2Q1QkZiUhIiokIjI7bSUpZm4xXTAjISM8NyQkIjFubW0iPksnKilcISM6JCIyJipvUT8rQk48JCEjPDckJCIrdk13JT4mISIqJCIxYzgmUm1RWyhSISM7NyQkIjFNTExlWiopKlImISM6JCIyRG93UkItbilbISM8NyQkIjFubW1UZy0wYyEjOiQiMUEqPkBxdzkoZSEjOzckJCIrRHQ6NWUhIiokIjFHOCJSKFtKKSlvISM7NyQkIjFubW0iUi9CLSchIzokIjFHRGhhI1wnR3ohIzs3JCQiMU1MTGU5WE1pISM6JCIxTzNSJ29iOyIqKSEjOzckJCIxLCsrRCYpZllrISM6JCIxI1IiWyJRZE16KiEjOzckJCIxbm1tImZYKGVtISM6JCIyakhCM0JtTTAiISM7NyQkIilCQHFvISIoJCIyLkIlZlxJMTU2ISM7NyQkIjFMTEwzIXo7MyghIzokIjFhKW9lZyUpbzkiISM6NyQkIixETzd1PSghIzUkIjFFI1JvO0t0OiIhIzo3JCQiMmttbW1yWEpIKCEjOyQiMW91UkJoSGk2ISM6NyQkIi12JFI3Z00oISM2JCIyMTBmQzstRjsiISM7NyQkIjFMTCQzMnopKVIoISM6JCIyJDRIKDReQDw7IiEjOzckJCIxbW0ienVYPFgoISM6JCIyKT1IQmVwTmY2ISM7NyQkIitEQ2gvdiEiKiQiMmFVJHljXmhiNiEjOzckJCIxTUwzRj44QXghIzokIjEneVFxKnAsRTYhIzo3JCQiMW9tO0g5bFJ6ISM6JCIyY3FHci1rWjIiISM7NyQkIi1ESjQ8ZCIpISM2JCIyMTtiN1p0VSsiISM7NyQkIjFOTExML3B1JCkhIzokIjEvO1IpMyhveSIqISM7NyQkIjFvbTtIPEdtJikhIzokIjFydSFwaG8hPSQpISM7NyQkIjEtKytESSh5dikhIzokIjElcD1dKy5xUighIzs3JCQiMU9MJDNLayVcKikhIzokIjFhc3MqPSsjXGshIzs3JCQiMW9tbTtjMFQiKiEjOiQiMSR5b2xwWiQ0YiEjOzckJCIxLStdUD91YyQqISM6JCIyMSplKiozKlFPXSUhIzw3JCQiMU1MTGUlR0NkKiEjOiQiMXdvcDExOylmJCEjOzckJCIxbm07els2KXkqISM6JCIyWS9TVm11WyRHISM8NyQkIjIsKysrOCFRKzUhIzokIjFZdmlwSzpcQSEjOzckJCIyLCsrXUxQPy0iISM6JCIyS1M5QjtEcSc9ISM8NyQkIiphJXBWNSEiKCQiMShwWER4TyI0PCEjOzckJCIyLStdNyUpMyJcNSEjOiQiMlAhSC5dUSFlcSIhIzw3JCQiLEQ5QlgwIiEiKiQiMUdtUjcqNHFyIiEjOzckJCItdlZ1JCpmNSEjNSQiMU9NTUFAc1U8ISM7NyQkIitYPE5sNSEiKSQiMSVmNyRHXidHeSIhIzs3JCQiLHZNIT13NSEiKiQiMktvJCkpKio0JGYhPiEjPDckJCIqJiozcTMiISIoJCIyMlRxQyJmdyUzIyEjPDckJCIqb3p5NSIhIigkIjJDJjNYZSYqXHlEISM8NyQkIipUXShHNiEiKCQiMUVSOWNGbVxLISM7NyQkIio5QCdcNiEiKCQiMTI2TmFsN3BTISM7NyQkIiooPVxxNiEiKCQiMVZTVixjSyxdISM7NyQkIjJvbSJII1xXJSo9IiEjOiQiMWN3Q3BBaDZmISM7NyQkIjJOTCRlOXJSMzchIzokIjFIJVJBWFIqW28hIzs3JCQiLnZvdFx0QSIhIzYkIjFhdEptKVEoenghIzs3JCQiMm5tbSJmQklZNyEjOiQiMV8oKkgqKmZucScpISM7NyQkIjJNTDMtJylRKW83ISM6JCIxbiZcNTQ7VWoqISM7NyQkIi1EaGBQIkgiISM1JCIyUEE9ZkhuWi8iISM7NyQkIjJubSJIaT0iUkoiISM6JCIyYDprMWkhKnA1IiEjOzckJCIyTUxMTE9ba0wiISM6JCIxVUFyX1klcDkiISM6NyQkIjJMTCRlXFMqZk0iISM6JCIyeXhLNFI+bToiISM7NyQkIjJNTExldFJiTiIhIzokIjJLXTlEMUw9OyIhIzs3JCQiMk5MJDNBYTNsOCEjOiQiMT4lNGw+UUQ7IiEjOjckJCIyTkxMJDM2anU4ISM6JCIyL1FKRVBHKGU2ISM7NyQkIjJOTEwzW0FQUiIhIzokIjJldWwpKlJWeDgiISM7NyQkIjJOTExMJlEiR1QiISM6JCIxQiEqUSZvUycqNCIhIzo3JCQiLkQiZTstTjkhIzYkIjJIUSVlQiRcYy4iISM7NyQkIjJubTtIWUhzWCIhIzokIjFVRHkjKipbYGAqISM7NyQkIjJNTDN4RVAlejkhIzokIjFpXTEicDBMZCkhIzs3JCQiLEQyWDtdIiEiKiQiMGxkcEw6d14oISM6NyQkIjJMTDNGdS04XyIhIzokIjEjZSplLD1jWWwhIzs3JCQiMm5tO0hUZzRhIiEjOiQiMUBcZHpXOCFlJiEjOzckJCIuREozPTFjIiEjNiQiMVp6Q1EkZWJsJSEjOzckJCIyTExMTHZ2LWUiISM6JCIxOD0lZkpZJTNRISM7NyQkIi5ESkxWPWciISM2JCIxdWE8RHN5MUkhIzs3JCQiMW5tIkgiNFRCOyEjOSQiMSZlIjNyQnJ1QiEjOzckJCIyT0wzRlx5XGsiISM6JCIyZGVxXWU0OiU+ISM8NyQkIjIuKytEMllsbSIhIzokIjItK2QqZTxERjwhIzw3JCQiLkQxejlvbiIhIzYkIjIoMzdSbXpYMDwhIzw3JCQiLXYzTjMobyIhIzUkIjJoMkBgWUtmdCIhIzw3JCQiMi4rdm9BX3RwIiEjOiQiMilwbihbPmEkPT0hIzw3JCQiK1g0aTI8ISIpJCIyJVFMaDRbJj0mPiEjPDckJCIyLitdN1FlIkc8ISM6JCIxPEppXE0lZk8jISM7NyQkIix2ImVwWzwhIiokIjFCKCkqeipSemdIISM7NyQkIjJOTDN4WkMseCIhIzokIjJ2aGgxW1FxdSQhIzw3JCQiMW5tIno4YDp6IiEjOSQiMUFGbzJrJG9tJSEjOzckJCIyLitEIil6IilIIj0hIzokIjFISFIiPTgieWMhIzs3JCQiMlBMTCRlL1RNPSEjOiQiMVAyR2YnNFl0JyEjOzckJCIyTUwzeGwpM2E9ISM6JCIwPC9xdklLcSghIzo3JCQiMk5MJDNkb3d0PSEjOiQiMSh5Tk49VT1qKSEjOzckJCIyUExlazBYTSo9ISM6JCIxTWE9ajVnJVsqISM7NyQkIjJNTExlREJKIj4hIzokIjIiKlIoUkkjZkctIiEjOzckJCItdltXZGI+ISM1JCIyQjVLP0wkXEw2ISM7NyQkIjJtbW07a0QhKSo+ISM6JCIyayY9NFJpRGg2ISM7NyQkIjJrbVQ1MnQrLSMhIzokIjIwdG4yJmZyUzYhIzs3JCQiMm1tOy9dP0AvIyEjOiQiMic9cHFKZUEoNCIhIzs3JCQiMm5tInpIejtrPyEjOiQiMkcuUSpSOCpHLiIhIzs3JCQiMW5tO2ZgQCczIyEjOSQiL2piImVxIzMmKiEjOTckJCIyLitdN3IyYTUjISM6JCIxPCU0cDhBLm8pISM7NyQkIjJOTExMMStZNyMhIzokIjEmSEYnPSlHJXl4ISM7NyQkIjJtbTthVCN6VkAhIzokIjEqW0hBcDRkJG8hIzs3JCQiLHZ3JSlIOyMhIiokIjE3L3Y7SXknKWUhIzs3JCQiMm5tVE4uOFA9IyEjOiQiMSoqUiIqKmY0XipbISM7NyQkIjJMTCRlKkhUVz8jISM6JCIyJikqeUlPNEd6UiEjPDckJCIuRGNjcF5BIyEjNiQiMSlHd0ojZl15SiEjOzckJCIxbm1tSnkqZUMjISM5JCIyPFJTYnNucV8jISM8NyQkIjIuK10oZUBKbkEhIzokIjFtI3lKeFsvLyMhIzs3JCQiMk5MTGVbRSgpRyMhIzokIjJsSW5NU21edyIhIzw3JCQiMW5UZ24rMyVIIyEjOSQiMndHTno1ajV0IiEjPDckJCIudiRcT1YqSCMhIzYkIjJbIik9dUctNnIiISM8NyQkIjJOJGU5SnN5L0IhIzokIjFjIlJ3OFRgcSIhIzs3JCQiMm9tO0giMzk1QiEjOiQiMiMpKjNmbWh6ODwhIzw3JCQiMk1MZWsoeiUzSyMhIzokIjIza1ZwPTxLeCIhIzw3JCQiKjliOkwjISIoJCIyZCxfbHglbykpPSEjPDckJCIyLitdNylbXV9CISM6JCIyLSZweWtiNXJBISM8NyQkIjIuKytEaWFNUCMhIzokIjJgOmNBSlJkJUchIzw3JCQiMi0rXVBPL1dSIyEjOiQiMXgqKlFIelgoZSQhIzs3JCQiKzBUTjpDISIpJCIxIypHZSMqbyNRWSUhIzs3JCQiMi4rdm8wQ2NWIyEjOiQiMSZcaHFReVFTJiEjOzckJCItdjNTKmVYIyEjNSQiMVA3aHluaiZSJyEjOzckJCIuRDEnUjt3QyEjNiQiMDlIWyNIXClSKCEjOjckJCIyLSsrRCJSVidcIyEjOiQiMVYnKSp5RSRRciQpISM7NyQkIi5EIlEtJSo9RCEjNiQiMSVwSXYvLiNvJCohIzs3JCQiLXZqbFdURCEjNSQiMUtALUhkKEctIiEjOjckJCIudiQqKUcmUmMjISM2JCIxRURtaVInNDQiISM6NyQkIis6I2ZrZSMhIikkIjJ1ajhUal13OCIhIzs3JCQiMmxtIkgjb3FsZiMhIzokIjEuR1wiUiopNDoiISM6NyQkIjFMTGVcQG8xRSEjOSQiMjw3U15jIVFmNiEjOzckJCIyJyoqXChvaCR6O0UhIzokIjJfPSNHTiVRRjsiISM7NyQkIjJsbW1UMzBwaSMhIzokIjJgKjRxKnBHNTsiISM7NyQkIi12PSFHcmsjISM1JCIyTEc+Uid5X1U2ISM7NyQkIjFMTExgNE5uRSEjOSQiMTNhKGY1S1k1IiEjOjckJCIyKCoqKipcX0slKilvIyEjOiQiMkZOXVNnPVgvIiEjOzckJCIybW1tO2JOMHIjISM6JCIxKkdFJHolcCZvJyohIzs3JCQiMk9MTDMmeTdLRiEjOiQiMUZFMjQqUUR2KSEjOzckJCIyLSsrKzo/UHYjISM6JCIxMjFBNmxqUnghIzs3JCQiMm1tIkgjKWZHdEYhIzokIjEneU1DNE51eCchIzs3JCQiMkxMJGU5PSZHeiMhIzokIjEnUSdbRm4rNmUhIzs3JCQiLnZvazxDIkchIzYkIjJaKCplUF9JcyhbISM8NyQkIjJsbW0iek0pPiRHISM6JCIyWWl0YGRSPCwlISM8NyQkIi52PTV1TCZHISM2JCIxXkxCSCxgIz0kISM7NyQkIjJOTCRlQ1p3dUchIzokIjFyVzZjTTg3RCEjOzckJCIxbjtIWmA6JypHISM5JCIyXjgkM2V0NUo/ISM8NyQkIjIwKysrKGZhPEghIzokIjJBIilRJjR5UGg8ISM8NyQkIjJzO3pXVTR3I0ghIzokIjJrbjpYV105ciIhIzw3JCQiMk9MZSp5R25QSCEjOiQiMVk9JjR3ajtyIiEjOzckJCIvdlZManRaSCEjNyQiMi8oeU0rYyxpPCEjPDckJCIybW07enkqemRIISM6JCIxaTs5a2sqPic9ISM7NyQkIi52b3BFeihIISM2JCIxc14kW2Uxaj8jISM7NyQkIjJLTExlZ2AhKSpIISM6JCIxc29OK01wSUYhIzs3JCQiMikqKioqKlwtIzQ+SSEjOiQiMT4rZTAlcCNbTSEjOzckJCIybm1tVFdJLC8kISM6JCIydGhUbig9eDJWISM8NyQkIjJQTExMJylvNjEkISM6JCIxMl4oNDApSHJfISM7NyQkIjIuKytERzJBMyQhIzokIjAsdXAoeU4nSCchIzo3JCQiMk5MM0Y9d0Y1JCEjOiQiMVpdSmopZllKKCEjOzckJCIybm07SDNYTDckISM6JCIxKilvaSc0VmNJKSEjOzckJCIuREopUiJSOSQhIzYkIjFiJ29GcUd2QSohIzs3JCQiMUxMTCQpR1trSiEjOSQiMiRlK2AqM1hUKyIhIzs3JCQiMicqKlxpOm0sJz0kISM6JCIybF8pPm0wMXU1ISM7NyQkIjFtbSJ6TV12PyQhIzkkIjIkW1NuMSFlXjciISM7NyQkIi9EMTlzSj1LISM3JCIyJFx4JCpbXihHOSIhIzs3JCQiMUwkMy0zJTNISyEjOSQiMmtwbWY6eF06IiEjOzckJCIwRCJHOHZZTUshIzgkIjJKKyNweEsxZjYhIzs3JCQiMms7YWolNCYpUkshIzokIjJjdCU0N0ZpaDYhIzs3JCQiMSQzRiV6VkJYSyEjOSQiMmpJekwvW0Y7IiEjOzckJCIsRCJ5aF1LISIqJCIybjdyNCxPQzsiISM7NyQkIjFtOy9jdE5ySyEjOSQiMUNZL3MkcHk5IiEjOjckJCIyTUwkZSoqbzQjSCQhIzokIjJ3eU4oMyp5RTYiISM7NyQkIjIxK0RKV09HSiQhIzokIjFISyxZR1BlNSEjOjckJCIxbm1tJylmZExMISM5JCIxSyF5aiNveXMpKiEjOzckJCIxbm1UbF55YUwhIzkkIjE6O2kwcyFSKyohIzs3JCQiMW5tO1dWKmZQJCEjOSQiMXkkUmI9JkdJISkhIzs3JCQiMW5tIkhfLnNSJCEjOSQiMSxWajBLYiYqcCEjOzckJCIxbm1tLEZUPU0hIzkkIjFKaj1AQjNZZiEjOzckJCIxTSQzeERZJVJNISM5JCIyUCxcKFsib3EkXCEjPDckJCIyMCtdUCIpei9ZJCEjOiQiMjhGamBgIkgvUyEjPDckJCIxbjt6cExeIlskISM5JCIxQmp5eXYnKSk9JCEjOzckJCIxTUwkZSNwYS1OISM5JCIxY20mXG9Rbl8jISM7NyQkIi52JHpYKD1fJCEjNiQiMicqXGRYN0sneT8hIzw3JCQiMnZtO0hCLTdhJCEjOiQiMlYxOkdbI1EsPSEjPDckJCIxJTMtajlNZ2EkISM5JCIxKlwlPjxGN2c8ISM7NyQkIi92b2ZnJzNiJCEjNyQiMilRQig0QTsudCIhIzw3JCQiMTxIMnR6cGJOISM5JCIyPTshKjRkSz9yIiEjPDckJCIxTSRlaykpSDBjJCEjOSQiMjhRKVtZV0owPCEjPDckJCIxXlAlKSp6aGBjJCEjOSQiMjN4IVtIdjw1PCEjPDckJCIydjtIS3IkPnFOISM6JCIyYDtbIm8vaEU8ISM8NyQkIjElZTltaURdZCQhIzkkIjI4Jm9MLlxkYTwhIzw3JCQiMi8rKythZCl6TiEjOiQiMU8oKjRlYisleiIhIzs3JCQiMU1MMzcjND9nJCEjOSQiMURHNzUqKkg9QCEjOzckJCIydm1tVCkzO0NPISM6JCIyJT1lUDYqcFttIyEjPDckJCIxLCtEY0RKWU8hIzkkIjFhel1dXSsyTSEjOzckJCIxTUxMR1VZb08hIzkkIjEvJlIkeXdWM1YhIzs3JCQiMW5tInolcEYpbyQhIzkkIjFqZW4sM3g4XyEjOzckJCIyMCsrdm0qMzNQISM6JCIxbmQhKlEsJmY8JyEjOzckJCIxTUwzKFEtenMkISM5JCIxOEktQSRIdDooISM7NyQkIjFubW0xXnJaUCEjOSQiMSJRc3F3NSY+IikhIzs3JCQiMi8rRCIpZlIpb1AhIzokIjFDIW8mW1MhPjMqISM7NyQkIjJPTCRlKjNrKip5JCEjOiQiMVlUIXooKj1wJCoqISM7NyQkIjFuOy8iZSkzNlEhIzkkIjIqZSslKm9VbGs1ISM7NyQkIixEMjhBJFEhIiokIjIlUio+Q2lBejYiISM7NyQkIi1EMT5WX1EhIzUkIjJrY2kiXChSLDoiISM7NyQkIip1XUUoUSEiKCQiMSpmMWVuZkU7IiEjOjckJCIvdlZbYXF4USEjNyQiMiZ6K15FVGppNiEjOzckJCIyMCt2bzpnRilRISM6JCIyJls5SVk7TWg2ISM7NyQkIjEsREpsWyJ5KVEhIzkkIjEmZXVuYCZ5ZTYhIzo3JCQiMSwrdnQmcEcqUSEjOSQiMkZwPXVLc1w6IiEjOzckJCIuRDEqKnlIIVIhIzYkIjI5RXUkZW5oVjYhIzs3JCQiMi8rK3ZTKTM4UiEjOiQiMmw2SFkqM1JGNiEjOzckJCIuRGNJO1skUiEjNiQiMnYuRF9ddW8yIiEjOzckJCIyMStdUD9XbCZSISM6JCIyW11HJDQjb3ErIiEjOzckJCIudj01cyN5UiEjNiQiMD5UMkdZREAqISM6NyQkIiQrJSEiIiQiMTRIPzE3cE0jKSEjOy0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIkKyUhIiIlKERFRkFVTFRHLSUrQVhFU0xBQkVMU0c2Jy1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRInQ2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVGT05URzYlJShERUZBVUxURyUoREVGQVVMVEciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJEklISIiLSUpQk9VTkRTX1lHNiMkIiNdISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlZ00hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlP04hIiItJSlDSElMRFJFTkc2Ig==
JSFH