# exact solution of B-S for price of Euro (also American) call option
CSt_formula := S*N(d[1]) - K*exp(-r*(T-t))*N(d[2]); NiM+SSxDU3RfZm9ybXVsYUc2IiwmKiZJIlNHRiUiIiItSSJOR0YlNiMmSSJkR0YlNiNGKUYpRikqKEkiS0dGJUYpLUkkZXhwRzYkSSpwcm90ZWN0ZWRHRjVJKF9zeXNsaWJHRiU2IywkKiZJInJHRiVGKSwmSSJUR0YlRilJInRHRiUhIiJGKUY+RiktRis2IyZGLjYjIiIjRilGPg== N_function := x->(1/2)*(1+erf(x/sqrt(2)));
plot(N_function(x),x=-5..5); NiM+SStOX2Z1bmN0aW9uRzYiZio2I0kieEdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYjIiIiIiIjRi4tSSRlcmZHRiU2IyomOSRGLi1JJXNxcnRHRiU2I0YvISIiRi1GJUYlRiU= LSUlUExPVEc2JS0lJ0NVUlZFU0c2JDdTNyQkISImIiIhJCIzTlMtbz1kXm1HISNDNyQkITNZTExMZSVHP3klISM8JCIzRWRlKikqPWJubilGLzckJCEzUW1tVCZlc0JmJUYzJCIzI1IjUSpSTnY2PiMhI0I3JCQhM0FMTCQzcyUzelZGMyQiM0o+aF5NKVIqZWZGOzckJCEzXkxMJGUvJFFrVEYzJCIzIlxQbjkueTRjIiEjQTckJCEzb21tVDU9cV1SRjMkIjNnXF8+QSE9aCpRRkY3JCQhM0lMTDNfPmZfUEYzJCIzVzZNUEYteV0oKUZGNyQkITNLKyt2bzFZWk5GMyQiM3N6JDRmNj5bJT4hI0A3JCQhMztMTDMtT0pOTEYzJCIzLS9FRmhFO2dVRlY3JCQhM3AqKipcUCpvJVE3JEYzJCIzJ0doJjMmcGxeIyopRlY3JCQhM0ttbW0iUkZqIUhGMyQiM1NmRCNwKkhcRz0hIz83JCQhMzNMTCRlNE9aciNGMyQiMzwnZmRgOUlrSiRGYG83JCQhM3cqKioqKlxuXCEqXCNGMyQiM0tybnYzRk1FaUZgbzckJCEzJikqKioqKlxpeENHI0YzJCIzcUR5aF15MEI2ISM+NyQkITMiKioqKioqXEtxUDIjRjMkIjNFYnIlKT5CLjA+RmBwNyQkITM6TEwzLVRDJSk9RjMkIjNHKjNWJUgjKWZ3SEZgcDckJCEzWm1tbSI0eillO0YzJCIzb2YmUi0keitkW0ZgcDckJCEzTG1tbW1gJ3pZIkYzJCIzRHIjSF5gIm8wckZgcDckJCEzIyoqKipcKD10KWVDIkYzJCIzJUg/KW9KKUdTMSIhIz03JCQhMyJvbW1taDUkXDVGMyQiM0pJYl0xZjxxOUZqcTckJCEzUyQqKipcKD1bakwpRmpxJCIzUihlOUsqUlZBP0ZqcTckJCEzKmYqKipcaVhnI0cnRmpxJCIzImZRIiopZW07XEVGanE3JCQhM25kbW1UJlEoUlRGanEkIjNHdWUqW0ZtV1IkRmpxNyQkITMkXG1tVGc9PjwjRmpxJCIzOVZycFFUSFNURmpxNyQkITN1RE1MTGUqZSRcRmBvJCIzLz81Yy8oMy4pXEZqcTckJCIzenJtO3pSUWJARmpxJCIzI1JOJWUzQEVgZUZqcTckJCIzXywrXSg9PlkyJUZqcSQiM3lmJ1suQGM7ZSdGanE3JCQiM3J1bW0ielh1OSdGanEkIjNzak1BRkhRMXRGanE3JCQiMyQ0KysrXXkpKUcpRmpxJCIzKlxWY1RrZ1QnekZqcTckJCIzSCsrXWlfUVE1RjMkIjNzdlo5Z2phLyYpRmpxNyQkIjNiKytEInklM1Q3RjMkIjNtQypRXiV5N0YqKUZqcTckJCIzLCsrXVAhW2hZIkYzJCIzKT1kZFI+Z3BHKkZqcTckJCIzaktMTCRReCRvO0YzJCIzXWYhKnBDZXpCJipGanE3JCQiM1orKyt2LkklKT1GMyQiM29oTks7I3lCcSpGanE3JCQiM0BtbSJ6cGUqej9GMyQiMy5YcUQkUmBCIikqRmpxNyQkIjM7LCsrRFwnUUgjRjMkIjM/MDliMSgzNSopKkZqcTckJCIzJ0hMJGU5UzgmXCNGMyQiMzF5VGY7Xy9QKipGanE3JCQiM3MrK0QxIz1icSNGMyQiMz9hZjNYNSFmJyoqRmpxNyQkIjMiSExMJDNzPzZIRjMkIjMiPnIicFQjKSo+KSoqRmpxNyQkIjNhKioqXDdgV2w3JEYzJCIzLSlvUFVKYzYqKipGanE3JCQiM2VubW1tKlJSTCRGMyQiM3csSD5WKD1kKioqRmpxNyQkIjMlem1tVHZKZ2EkRjMkIjMncHMwWWdXISkqKipGanE3JCQiM15NTGU5dE9jUEYzJCIzK0U4Jz4sUSIqKioqRmpxNyQkIjMxLCsrXVFrXFJGMyQiM2docj1pJzMnKioqKkZqcTckJCIzIVtMTDNkZzY8JUYzJCIzayF6IypRWlspKioqKkZqcTckJCIzJnltbW13KEdwVkYzJCIzdSw1enR3JCoqKioqRmpxNyQkIjNDKytEIm9LMGUlRjMkIjM7d0xNNW8oKioqKipGanE3JCQiMzYsK3Y9NXMjeSVGMyQiMm9jJzRgOCoqKioqKkYzNyQkIiImRiwkIjI/RyVbTHIqKioqKipGMy0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiRGLEZnW2xGaFtsLSUrQVhFU0xBQkVMU0c2JFEieDYiUSFGXVxsLSUlVklFV0c2JDskISNdRmdbbCQiI11GZ1tsOyQhMmNfTyMpPXEqKio+RmpxJCIyX08jKT1xKioqPjUhIzs= d1_function := t -> (log(S/K) + (r + (1/2)*sigma^2)*(T-t))/(sigma*sqrt(T-t));
d2_function := t -> (log(S/K) + (r - (1/2)*sigma^2)*(T-t))/(sigma*sqrt(T-t)); NiM+SSxkMV9mdW5jdGlvbkc2ImYqNiNJInRHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSooLCYtSSRsb2dHRiU2IyomSSJTR0YlIiIiSSJLR0YlISIiRjMqJiwmSSJyR0YlRjMqJEkmc2lnbWFHRiUiIiMjRjNGO0YzLCZJIlRHRiVGMzkkRjVGM0YzRjNGOkY1LUklc3FydEc2JEkqcHJvdGVjdGVkR0ZDSShfc3lzbGliR0YlNiNGPUY1RiVGJUYl NiM+SSxkMl9mdW5jdGlvbkc2ImYqNiNJInRHRiVGJTYkSSlvcGVyYXRvckdGJUkmYXJyb3dHRiVGJSooLCYtSSRsb2dHRiU2IyomSSJTR0YlIiIiSSJLR0YlISIiRjMqJiwmSSJyR0YlRjMqJEkmc2lnbWFHRiUiIiMjRjVGO0YzLCZJIlRHRiVGMzkkRjVGM0YzRjNGOkY1LUklc3FydEc2JEkqcHJvdGVjdGVkR0ZDSShfc3lzbGliR0YlNiNGPUY1RiVGJUYl CSt := subs(d[1]=d1_function(t),d[2]=d2_function(t),CSt_formula);
CSt := eval(subs(N=N_function,%)); NiM+SSRDU3RHNiIsJiomSSJTR0YlIiIiLUkiTkdGJTYjKigsJi1JI2xuRzYkSSpwcm90ZWN0ZWRHRjJJKF9zeXNsaWJHRiU2IyomRihGKUkiS0dGJSEiIkYpKiYsJkkickdGJUYpKiRJJnNpZ21hR0YlIiIjI0YpRj1GKSwmSSJUR0YlRilJInRHRiVGN0YpRilGKUY8RjdGPyNGN0Y9RilGKSooRjZGKS1JJGV4cEdGMTYjLCQqJkY6RilGP0YpRjdGKS1GKzYjKigsJkYvRikqJiwmRjpGKUY7RkJGKUY/RilGKUYpRjxGN0Y/RkJGKUY3 NiM+SSRDU3RHNiIsJiomSSJTR0YlIiIiLCYjRikiIiNGKS1JJGVyZkc2JEkqcHJvdGVjdGVkR0YwSShfc3lzbGliR0YlNiMsJCoqLCYtSSNsbkdGLzYjKiZGKEYpSSJLR0YlISIiRikqJiwmSSJyR0YlRikqJEkmc2lnbWFHRiVGLEYrRiksJkkiVEdGJUYpSSJ0R0YlRjtGKUYpRilGQEY7RkEjRjtGLEYsRitGK0YrRilGKSooRjpGKS1JJGV4cEdGLzYjLCQqJkY+RilGQUYpRjtGKSwmRitGKS1GLjYjLCQqKiwmRjZGKSomLCZGPkYpRj9GREYpRkFGKUYpRilGQEY7RkFGREYsRitGK0YrRilGOw== # Exercise L10: 6 month call, r = 0.04
# L.10 Use a binomial lattice to find the value of a value of a six-month call option with strike price # # $52, Current price $55, interest rate 4%, standard deviation 0.25, and Delta t= 1 month.
# L.11 Use a binomial lattice to find the value of a value of a six-month American put and a six-month # # European put option with the same parameters as in exercise L. 10. Fri 11/13
# Read: Sec 13.1, 13.2, 13.3, 13.5.
# Homework: 13.2, 13.4
# L.12 Use the Black-Scholes formulas to compute the values of a European call option and European put # # option with the same parameters as in Exercise L.10. (Show your work.)
T_value := evalf(6/12);
sigma_value := 0.25;
r_value := 0.04;
params := {K = 52, sigma = sigma_value, r = r_value, T = T_value}; NiM+SShUX3ZhbHVlRzYiJCIrKysrK10hIzU= NiM+SSxzaWdtYV92YWx1ZUc2IiQiI0QhIiM= NiM+SShyX3ZhbHVlRzYiJCIiJSEiIw== NiM+SSdwYXJhbXNHNiI8Ji9JIktHRiUiI18vSSZzaWdtYUdGJSQiI0QhIiMvSSJyR0YlJCIiJUYuL0kiVEdGJSQiKysrKytdISM1 CSt_params := subs(params,N=N_function,CSt); NiM+SStDU3RfcGFyYW1zRzYiLCYqJkkiU0dGJSIiIiwmI0YpIiIjRiktSSRlcmZHNiRJKnByb3RlY3RlZEdGMEkoX3N5c2xpYkdGJTYjLCQqKCwoLUkjbG5HRi82IywkRigjRikiI19GKSQiKysrXWlOISM2RilJInRHRiUkISsrKytEckY+RiksJiQiKysrKytdISM1RilGPyEiIiNGRkYsRixGKyQiKysrKys/ISIqRitGKUYpKiYtSSRleHBHRi82IywmJCErKysrKz9GPkYpRj8kIiIlISIjRiksJkYrRiktRi42IywkKigsKEY2RikkIisrKyt2ViEjN0YpRj8kISorKyt2KUY+RilGQkZHRixGK0ZIRitGKSEjXw== plot0 := plot(subs(t=0,CSt_params),S=40..60):
plot1 := plot(subs(t=1/12,CSt_params),S=40..60):
plot2 := plot(subs(t=2/12,CSt_params),S=40..60):
plot3 := plot(subs(t=3/12,CSt_params),S=40..60):
plot4 := plot(subs(t=4/12,CSt_params),S=40..60):
plot5 := plot(subs(t=5/12,CSt_params),S=40..60):
plot6 := plot(subs(t=6/12,CSt_params),S=40..60):
with(plots):
display([plot0,plot1,plot2,plot3,plot4,plot5,plot6]);
subs(S=55,t=0,CSt_params);
evalf(%); LSUlUExPVEc2JS0lJ0NVUlZFU0c2KjdTNyQkIiNTIiIhJCIzIlEheltBKDMtNiQhIz03JCQiM1xMTEwzVmZWUyEjOyQiM01ePnQpPTglcE5GLzckJCIzIXBtO0hbRDozJUYzJCIzWy5jNXZUeTRTRi83JCQiM09MTCRlMCQ9Q1RGMyQiM0RxZFY8IVtKYiVGLzckJCIzWkxMJDNSQnI7JUYzJCIzaF1PI1s9Wlc6JkYvNyQkIjNcbW0iempmKTRVRjMkIjN2S09HJFtWKjRlRi83JCQiM2hMTGU0O1tcVUYzJCIzOk8xYiVwMTVaJ0YvNyQkIjNHKytEbXldIUglRjMkIjMlKXlFQj4rKT5AKEYvNyQkIjNqTExlenMkSEwlRjMkIjM1SC1tcmY/VCEpRi83JCQiM0MrK0RAMUJ2VkYzJCIzOi4ySidwaFEkKilGLzckJCIzNG5tbUBYdD1XRjMkIjMhXDNPR0RXTCMqKkYvNyQkIjM3TEwkM3lfcVglRjMkIjM9aj1mbi1tJjMiISM8NyQkIjNLKysrbCs+K1hGMyQiM1c8PXcpSEx5PiJGYm83JCQiM0krKyt2V11WWEYzJCIzZ1hrYyJlSiM9OEZibzckJCIzUCsrK05mQyZlJUYzJCIza0FuNiNwJHlUOUZibzckJCIzWUxMZXo2OkJZRjMkIjMxKTNbMywqXGc6RmJvNyQkIjN6bW1tIj1DI29ZRjMkIjMlPXcwM1ZtKTQ8RmJvNyQkIjNRbW1tRXBTMVpGMyQiM0A5K0RHQllWPUZibzckJCIzJCoqKipcaWBBM3YlRjMkIjMvYDYoUjwvcisjRmJvNyQkIjNzbW1td3k4IXolRjMkIjM9OCM+Qk8vJWZARmJvNyQkIjM4KytET0lGTFtGMyQiM3AjZT9PXCVmTUJGYm83JCQiMyEqKioqXCgzek11W0YzJCIzRGJiSlowSjRERmJvNyQkIjNubW07SF8/PFxGMyQiMz0nW1d1OEEpKnAjRmJvNyQkIjMpb207emlobCZcRjMkIjM/Wy5pTW45IylHRmJvNyQkIjM5TEwkMyNHLCoqXEYzJCIzPUszbSkqR3AnMyRGYm83JCQiM1dMTGV6dzVWXUYzJCIzK3d6JlxYb3hJJEZibzckJCIzUisrdiRRI1wiMyZGMyQiM0AjKiopNCI+KUgyTkZibzckJCIzbkxMJGUiKltINyZGMyQiM0x2cEkrUzdJUEZibzckJCIzPysrK3F2eGxeRjMkIjMuYCdRJ1IiPiNvUkZibzckJCIzMSsrXV9xbjJfRjMkIjNYY0pFbiMqejNVRmJvNyQkIjNZKytEY3BAW19GMyQiM3UoUT1kNGQnW1dGYm83JCQiM2sqKioqXDInSEtIJkYzJCIzPm0oWzoxKilIcyVGYm83JCQiM2BtbW13YW5MYEYzJCIzRkNWInolR2B3XEZibzckJCIzWCsrK3YrJ29QJkYzJCIzJ1FgdWRPUldEJkZibzckJCIzQ0xMZVI8KmZUJkYzJCIzKDNIMWliSkReJkZibzckJCIzQisrKyYpSHhlYUYzJCIzUGwtUDMsTSxlRmJvNyQkIjNTbW0iSCFvLSpcJkYzJCIzR2prRXZ6RHpnRmJvNyQkIjNLKytEVE81VGJGMyQiM1A+OUg6ZChmUCdGYm83JCQiM3dtbW1UOUMjZSZGMyQiMykqNHhmSCpHP24nRmJvNyQkIjNiKioqXGkhKjNgaSZGMyQiM2NrJFs4aXkhKSlwRmJvNyQkIjMoUUxMTCp6eW1jRjMkIjMrIz0xbTwkPSlIKEZibzckJCIzJVJMTDNOMSM0ZEYzJCIzIVs8cigpKT0hNGkoRmJvNyQkIjNFbm0iSFl0N3YmRjMkIjNRJVsqbyE9SWolekZibzckJCIzQCsrK3EoRyoqeSZGMyQiMzgoUkxgI2YiKlwjKUZibzckJCIzS25tOzlAQk1lRjMkIjNhKlE4bjk3SWcpRmJvNyQkIjNATExMYHYmUShlRjMkIjNPYykqWy1dSUIqKUZibzckJCIzcCoqKlxpYDFoImZGMyQiM3lqMlx6KG8jcCMqRmJvNyQkIjNlKyt2LlVhY2ZGMyQiMz1BSk5CaWAvJypGYm83JCQiI2dGLCQiM2dzSUhqczhwKipGYm83UzckRiokIjMtdj5BTHZzRT9GLzckRjEkIjNhMFBwWy1qdEJGLzckRjckIjMyeiM0emQqMzdGRi83JEY8JCIzTShbKSlwWiRmT0pGLzckRkEkIjNRJypSNSM+MFVoJEYvNyRGRiQiM3oiUSZbSGA6VlRGLzckRkskIjMpcHllKj0lb1ZvJUYvNyRGUCQiMyEzJkhiRltFKkgmRi83JEZVJCIzJyo+ODZIbGonKmZGLzckRlokIjNdOSx1ZXQyZG5GLzckRmluJCIzYidmKmUzW101d0YvNyRGXm8kIjMrZjRwWlBdQyUpRi83JEZkbyQiMypbaHBxJnpLOCUqRi83JEZpbyQiMyU+U01RIVJnWzVGYm83JEZecCQiMzpsUkk6RHlmNkZibzckRmNwJCIzMnVKbyJIX3ZFIkZibzckRmhwJCIzQTIlR0lzSVZTIkZibzckRl1xJCIzYEdhPzllbkY6RmJvNyRGYnEkIjNbbE1PKSopRyp6O0ZibzckRmdxJCIzNlE7IypcUG5BPUZibzckRlxyJCIzdk9hVyxQKnopPkZibzckRmFyJCIzO1gycndiJVI6I0ZibzckRmZyJCIzYHhtRSE9MmdMI0ZibzckRltzJCIzYUY+W3o7QTZERmJvNyRGYHMkIjMwTG5eV3IiKTNGRmJvNyRGZXMkIjMnb3AtVF16TSNIRmJvNyRGanMkIjN4SWkkeWkvIj1KRmJvNyRGX3QkIjNudSNvTENZakwkRmJvNyRGZHQkIjMjb05HViE9W3FORmJvNyRGaXQkIjM2PD8sa2AleiFRRmJvNyRGXnUkIjMoKjNZZjAqKltYU0ZibzckRmN1JCIzd3YsQl9bLz1WRmJvNyRGaHUkIjNuJzNYZGcjb3FYRmJvNyRGXXYkIjNIQ2FJNiU9JFtbRmJvNyRGYnYkIjMiKSozJylHUGluNSZGYm83JEZndiQiMyUzR2YqKUd2bFImRmJvNyRGXHckIjNrWmpKVCQqKmZuJkZibzckRmF3JCIzJyo0Jz1VRFJbKGZGYm83JEZmdyQiM3kkSGQ1L3VNRidGYm83JEZbeCQiMy9hJnA5YkNGZidGYm83JEZgeCQiMydHa0U8XmJqIXBGYm83JEZleCQiMyQqUTxfXDY0THNGYm83JEZqeCQiM1BcbmxGUihHYyhGYm83JEZfeSQiMzVzWXknKkh2cXlGYm83JEZkeSQiM2BMTyVHOnUhSCMpRmJvNyRGaXkkIjM6O0k/NjZHYSYpRmJvNyRGXnokIjNSJXl1bykpKXAwKilGYm83JEZjeiQiMypIZEhuRGhqQypGYm83JEZoeiQiMykpMzUkPTg6cGgqRmJvN1M3JEYqJCIzI2ZzV2ZjYjw3IkYvNyRGMSQiM2tWWyhcIT1TYDhGLzckRjckIjMwQUV4KW9fX2UiRi83JEY8JCIzIXo5TDBgdEopPUYvNyRGQSQiM0NaY09nZW5FQUYvNyRGRiQiM2pXV2VjJ3pnaCNGLzckRkskIjNcPG9jWikqNEJJRi83JEZQJCIzbyFHIiplMFRbXCRGLzckRlUkIjNLKmZgMCJ5WFNTRi83JEZaJCIzL0Z6PHBbeFlZRi83JEZpbiQiM24lRykpb18rKFJgRi83JEZebyQiMzx6Vi9AUWA2Z0YvNyRGZG8kIjNQcTk4YmRUU29GLzckRmlvJCIzWGUxYCF6NFB2KEYvNyRGXnAkIjNhZCopKUg4TFNyKUYvNyRGY3AkIjNzTStBaSNwb2wqRi83JEZocCQiM1dmJ1whSFAnbzMiRmJvNyRGXXEkIjN5UXdwW0hXKD4iRmJvNyRGYnEkIjNzND9mKHpmYUwiRmJvNyRGZ3EkIjMnPjZqPSR6QW05RmJvNyRGXHIkIjMhPiNSJSk+KGUiPjtGYm83JEZhciQiMzpHKFwxPDNUeCJGYm83JEZmciQiM2VxY1V6bWRYPkZibzckRltzJCIzeVtPRDZgIT42I0ZibzckRmBzJCIzdEVbenInZTNJI0ZibzckRmVzJCIzLGp1dy9GaDJERmJvNyRGanMkIjM3JSk+Rm8lZWlwI0ZibzckRl90JCIzJypSJFttdikqKjNIRmJvNyRGZHQkIjNKLCJ6c08nXFFKRmJvNyRGaXQkIjNVaHYoSEpXQ1AkRmJvNyRGXnUkIjMnKkckR2xPTnZnJEZibzckRmN1JCIzYkI4VnZPVnlRRmJvNyRGaHUkIjMkR1FoJ2U/XklURmJvNyRGXXYkIjN4JD5wIyopUVszV0ZibzckRmJ2JCIzKTQmR1ROeC1vWUZibzckRmd2JCIzNiM0JmVIRicpZlxGYm83JEZcdyQiMz1YXXlFXCM+QyZGYm83JEZhdyQiM0Q0PTwhZktVYSZGYm83JEZmdyQiM045UWs4aSFwJWVGYm83JEZbeCQiM1N0JWVEWjc1PCdGYm83JEZgeCQiM10jeSdRRCZvKSpbJ0ZibzckRmV4JCIzMyEqekMvUldBb0ZibzckRmp4JCIzSjI5ZSdIYSVlckZibzckRl95JCIzInpdZi9gKVFzdUZibzckRmR5JCIzcylmJjR3Qyl6JHlGYm83JEZpeSQiM0gjKVF2L1slKnAiKUZibzckRl56JCIzPy8uXydbdChHJilGYm83JEZjeiQiM3Q2c2QpZSRvdykpRmJvNyRGaHokIjNeeF9UNFs5YiMqRmJvN1M3JEYqJCIzT3VFVjhqZXFYISM+NyRGMSQiM1QmRzozOXJ1eSZGaF1tNyRGNyQiM18jUkwncEhpZHFGaF1tNyRGPCQiM0omND13O05ndilGaF1tNyRGQSQiMzcnRz5RPmolejVGLzckRkYkIjNoWSdlSDN0Jj44Ri83JEZLJCIzJilwQWs8LFl6OkYvNyRGUCQiMyFIRzcqMyFbMSo9Ri83JEZVJCIzUXpkd3ZPR2lBRi83JEZaJCIzVTdyJHllJD0pbyNGLzckRmluJCIzYCVSeiZ6ITMmKj0kRi83JEZebyQiMzMscVUtdW0pbyRGLzckRmRvJCIzTF5MNUYheixLJUYvNyRGaW8kIjMrWmonR3AuUC4mRi83JEZecCQiMyM9byZcV21oLGVGLzckRmNwJCIzeE9xKjREZTZkJ0YvNyRGaHAkIjNLcndATVBMIWUoRi83JEZdcSQiMzZAMj5VIlEpPSYpRi83JEZicSQiMzp3aiQzWVo1cipGLzckRmdxJCIzNEolei9Ba2YzIkZibzckRlxyJCIzKTRNcHFbb0JBIkZibzckRmFyJCIzdkFxZXJYZ2k4RmJvNyRGZnIkIjNBJXBKViY9ISo+OkZibzckRltzJCIzRl9oRzBsT3U7RmJvNyRGYHMkIjNEJVxPakFPPSY9RmJvNyRGZXMkIjNbZ3YqXF91Ils/RmJvNyRGanMkIjNzLCN5QCkqXCFIQUZibzckRl90JCIzQyw8VEc0eU1DRmJvNyRGZHQkIjNUaTJCUXViZUVGYm83JEZpdCQiM0suJG9fOSpSKSlHRmJvNyRGXnUkIjM1YlFKNzshNDckRmJvNyRGY3UkIjNrV0ksMD9bIVIkRmJvNyRGaHUkIjMjR015Inpec1VPRmJvNyRGXXYkIjNnZSM9c3E9QSNSRmJvNyRGYnYkIjNnInprTzF5VT0lRmJvNyRGZ3YkIjM5Sy9lLjAtIVslRmJvNyRGXHckIjMneWcudCRHeW1aRmJvNyRGYXckIjMmUllaYHF0XDImRmJvNyRGZnckIjNGUFEqPidlRSVRJkZibzckRlt4JCIzZmJyeixwNjtkRmJvNyRGYHgkIjNMJlFJa2VASi8nRmJvNyRGZXgkIjMpKXphZDwkR1lRJ0ZibzckRmp4JCIzaSxdUD9ZKipIbkZibzckRl95JCIzJmYhSHIiMyIpRzAoRmJvNyRGZHkkIjN3NiJ6emllIUh1RmJvNyRGaXkkIjN2aj0oelUqcHF4RmJvNyRGXnokIjMocGdPYk0oKSpSIilGYm83JEZjeiQiM1czK1UoKUcoelwpRmJvNyRGaHokIjM7S2F5didlcykpKUZibzdTNyRGKiQiMylvP3ZOPVJlOCohIz83JEYxJCIzKj4mZlV0Oio0RiJGaF1tNyRGNyQiMy9pOipHRitpbiJGaF1tNyRGPCQiM2BoJ0dRRkxCRSNGaF1tNyRGQSQiMzFXZVInZjBRLSRGaF1tNyRGRiQiM3MlcGp5aiFwISpSRmhdbTckRkskIjNGNlV3ZEBPNl5GaF1tNyRGUCQiMyNRXiE0J1tDPGEnRmhdbTckRlUkIjN5K0pRL0xQZyQpRmhdbTckRlokIjN5eWx3YCNvdDAiRi83JEZpbiQiM2FnVSkpUk1KTDhGLzckRl5vJCIzMitAPUIqNEVpIkYvNyRGZG8kIjNkd3B5c0kobysjRi83JEZpbyQiM3YhUjo8LzhFWSNGLzckRl5wJCIzU3lGLmAuX3ZIRi83JEZjcCQiM0t3eihmaCwsXiRGLzckRmhwJCIzdWpcKSo0I28nUVVGLzckRl1xJCIzX2ltJDNlLDUlXEYvNyRGYnEkIjM7JTQmUSE9cU0nZUYvNyRGZ3EkIjMjPnNFTE41MHknRi83JEZcciQiM0VuUW0+bkgsekYvNyRGYXIkIjM2YCRva2RCXzMqRi83JEZmciQiM1JkMzdmIXBZLyJGYm83JEZbcyQiM295cnJlaVIiPSJGYm83JEZgcyQiMyEqby4lKipRUjxNIkZibzckRmVzJCIzdl1YRiw4bkE6RmJvNyRGanMkIjMnSGBpP25FQXAiRmJvNyRGX3QkIjNeJil5W3ckKip6KT1GYm83JEZkdCQiM0AieWlMJikzUzUjRmJvNyRGaXQkIjMmcFdbKilmUihHQkZibzckRl51JCIzRSJmP1dtNSdlREZibzckRmN1JCIzNlNbOno/J3kjR0ZibzckRmh1JCIzZE9hRyV5Uj8zJEZibzckRl12JCIzSXYjKUg0YHpsTEZibzckRmJ2JCIzJGUmWydmP1NOaiRGYm83JEZndiQiM29rWiQpUUdKUFJGYm83JEZcdyQiMyUqenglRzAhPkxVRmJvNyRGYXckIjMlSGNITlpLQmIlRmJvNyRGZnckIjMmZXUyTyNmYXRbRmJvNyRGW3gkIjN6cGpgISpHJio9X0ZibzckRmB4JCIzYip6QnB0XylmYkZibzckRmV4JCIzRyczek52P2kiZkZibzckRmp4JCIzVmloOTRqeXdpRmJvNyRGX3kkIjM1c3FYPTIqUWgnRmJvNyRGZHkkIjMhby4pNCl5KlsxcUZibzckRml5JCIzNGdHa1tYeml0RmJvNyRGXnokIjM6VW9FIjNndnUoRmJvNyRGY3okIjNSLSU+ei0pND8iKUZibzckRmh6JCIzNSxNdXAtakMmKUZibzdTNyRGKiQiMyZ6VnUhKkhAPU0iISNANyRGMSQiMyUqZVxVYnJ1Y0NGYmBuNyRGNyQiMz0peVpcUiVmcFNGYmBuNyRGPCQiMydbPzNVLDFXLChGYmBuNyRGQSQiM01kR255VmolPSJGXWdtNyRGRiQiM3ZPX1hOIWYsJj5GXWdtNyRGSyQiMzBeN3N2JmZbLiRGXWdtNyRGUCQiM3krOCVmSG9acSVGXWdtNyRGVSQiM1R3RyhRO09jRChGXWdtNyRGWiQiM1t4XDp5UGEmNCJGaF1tNyRGaW4kIjNCTD82PCQpKTNrIkZoXW03JEZebyQiM1pgN0dQKilwL0JGaF1tNyRGZG8kIjMmUkAkKXloTCo+TEZoXW03JEZpbyQiMyY+WEllRDpXcSVGaF1tNyRGXnAkIjMmW0FRQyE9VXZrRmhdbTckRmNwJCIzKSkqXD9DLm8nUiYpRmhdbTckRmhwJCIzNCEpPT5iLSJ5OyJGLzckRl1xJCIzZyQpUk4ncDVEXSJGLzckRmJxJCIzME9fRD5kPCYpPkYvNyRGZ3EkIjNJSzI1Q2NiM0RGLzckRlxyJCIzJUdYXEghXHkrS0YvNyRGYXIkIjM+Vy9BRiIzeilSRi83JEZmciQiM3YkNFNBbV1kJlxGLzckRltzJCIzKG8hKmY8PEhzKWZGLzckRmBzJCIzdyEqZjhtaVJqc0YvNyRGZXMkIjNpKj53JUg6T3koKUYvNyRGanMkIjNmXCNSai1SZy0iRmJvNyRGX3QkIjMxKClwKnlbKGYuN0ZibzckRmR0JCIzJD5zJW9sJFFqUyJGYm83JEZpdCQiM1Ftck5UYGxCO0ZibzckRl51JCIzJnAkelNGQGJePUZibzckRmN1JCIzUGpUIil6JSpRQ0BGYm83JEZodSQiMzB4RCF5LGZtUSNGYm83JEZddiQiM0U/UT9yKSlvJG8jRmJvNyRGYnYkIjNoMSFHWFlLcidIRmJvNyRGZ3YkIjNnIWYpZiYqcFsiSCRGYm83JEZcdyQiM1YuPnhsT1U0T0ZibzckRmF3JCIzX2V5KHBgLFEmUkZibzckRmZ3JCIzVEQkUUNhWjdJJUZibzckRlt4JCIzZSwrLmcpZV5uJUZibzckRmB4JCIzcTdNKDQ6dFIvJkZibzckRmV4JCIzTk85dTxHJilHYUZibzckRmp4JCIzVio+NGFqQ3MiZUZibzckRl95JCIzMSkqXEliQjp6aEZibzckRmR5JCIzY1AzYyt3LipmJ0ZibzckRml5JCIzQSd5eEsiPlt5cEZibzckRl56JCIzcFhjJ1FgT2tRKEZibzckRmN6JCIzTTUqbyF6TWp6eEZibzckRmh6JCIzTGNaRVooR1k/KUZibzdmbjckRiokRixGLDckRjFGZWluNyRGN0ZlaW43JEY8RmVpbjckRkFGZWluNyRGRkZlaW43JEZLRmVpbjckRlBGZWluNyRGVUZlaW43JEZaRmVpbjckRmluRmVpbjckRl5vRmVpbjckRmRvRmVpbjckRmlvRmVpbjckRl5wRmVpbjckRmNwRmVpbjckRmhwRmVpbjckRl1xRmVpbjckRmJxRmVpbjckRmdxRmVpbjckRlxyRmVpbjckRmFyRmVpbjckRmZyRmVpbjckRltzRmVpbjckRmBzRmVpbjckRmVzRmVpbjckRmpzRmVpbjckRl90RmVpbjckJCIzJHBtO0hDalY5JkYzRmVpbjckRmR0RmVpbjckJCIzOytdaVNDRHdeRjNGZWluNyQkIjM4KytENnRzJz0mRjNGZWluNyQkIjNbK0RjWVonPj4mRjNGZWluNyQkIjMzK10oPT0tcz4mRjNGZWluNyQkIjNjXDdgKiozIykqPiZGM0ZlaW43JCQiM3MqXCg9PCdSQz8mRjMkIjN5PygqXCg9PCdSQ0ZoXW03JCQiMyopXFAlW0xlXT8mRjMkIjMmKSopKVxQJVtMZV1GaF1tNyRGaXQkIjMhKmUrKytEMHh3RmhdbTckJCIzRStdUC9xJXpBJkYzJCIzO0UrXVAvcSV6I0YvNyRGXnUkIjNYWSsrRGNwQFtGLzckRmN1JCIzXGsqKioqXDInSEsqRi83JEZodSQiM0VsbW1tWnZPOEZibzckRl12JCIzXS8rK10yZ288RmJvNyRGYnYkIjNVS0wkZVI8KmZARmJvNyRGZ3YkIjNMLSsrXSlIeGUjRmJvNyRGXHckIjM3a207SCFvLSpIRmJvNyRGYXckIjNBLitdN2suNk1GYm83JEZmdyQiM2ZubW07V1RBUUZibzckRlt4JCIzXyYqKipcaSEqM2BVRmJvNyRGYHgkIjNzUUxMTCp6eW0lRmJvNyRGZXgkIjNXUkxMM04xIzQmRmJvNyRGangkIjNjc207SFl0N2JGYm83JEZfeSQiMzctKysreEcqKmVGYm83JEZkeSQiMzp0bW1UNktVakZibzckRml5JCIzOktMTExiZFFuRmJvNyRGXnokIjMjcCoqKlxpYDFockZibzckRmN6JCIzeDArXVA/V2x2RmJvNyRGaHokIiIpRiwtJSZDT0xPUkc2OCUkUkdCRyQiIzUhIiIkRixGZWFvRmZhb0ZjYW9GZmFvRmZhb0ZjYW9GZmFvRmZhb0ZjYW9GZmFvRmZhb0ZjYW9GZmFvRmZhb0ZjYW9GZmFvRmZhb0ZjYW9GZmFvRmZhby0lK0FYRVNMQUJFTFNHNidRIlM2IlEhRltiby0lJUZPTlRHNiQlKkhFTFZFVElDQUdGZGFvJStIT1JJWk9OVEFMR0ZhYm8tJSVWSUVXRzYkO0YqRmh6OyQhMldoZUVYRlEqPkZibyQiMk0qZSYzK19vLCIhIzo= NiMsKCMiI2IiIiMiIiItSSRlcmZHNiRJKnByb3RlY3RlZEdGK0koX3N5c2xpYkc2IjYjLCQqJiwmLUkjbG5HRio2IyNGJSIjX0YnJCIrKytdaU4hIzZGJ0YnRiYjRidGJiQiK0NyVUdHISIqRiQqJi1JJGV4cEdGKjYjJCErKysrKz9GOUYnLCZGOkYnLUYpNiMsJComLCZGMkYnJCIrKysrdlYhIzdGJ0YnRiZGOkY7RjpGJyEjXw== NiMkIitOKFxnMychIio= # Check by substitution that CSt(S,t) actually satisfies Black-Sholes equation
BS_pde := diff(f(S,t),t)+r*S*diff(f(S,t),S) + (1/2)*sigma^2*S^2*diff(diff(f(S,t),S),S) - r*f(S,t) = 0;
print(` If f(S,t) = CSt is a solution the following should give 0 = 0`);
simplify(subs(f(S,t)=CSt,BS_pde)); NiM+SSdCU19wZGVHNiIvLCotSSVkaWZmR0kqcHJvdGVjdGVkR0YqNiQtSSJmR0YlNiRJIlNHRiVJInRHRiVGMCIiIiooSSJyR0YlRjFGL0YxLUYpNiRGLEYvRjFGMSooSSZzaWdtYUdGJSIiI0YvRjgtRik2JEYsLUkiJEdGKjYkRi9GOEYxI0YxRjgqJkYzRjFGLEYxISIiIiIh NiNJaW5+SWZ+ZihTLHQpfj1+Q1N0fmlzfmF+c29sdXRpb25+dGhlfmZvbGxvd2luZ35zaG91bGR+Z2l2ZX4wfj1+MEc2Ig== NiMvIiIhRiQ= # Use put-call parity to determine price of Euro. put
PSt := CSt + K*exp(-r*(T-t)) - S; NiM+SSRQU3RHNiIsKiomSSJTR0YlIiIiLCYjRikiIiNGKS1JJGVyZkc2JEkqcHJvdGVjdGVkR0YwSShfc3lzbGliR0YlNiMsJCoqLCYtSSNsbkdGLzYjKiZGKEYpSSJLR0YlISIiRikqJiwmSSJyR0YlRikqJEkmc2lnbWFHRiVGLEYrRiksJkkiVEdGJUYpSSJ0R0YlRjtGKUYpRilGQEY7RkEjRjtGLEYsRitGK0YrRilGKSooRjpGKS1JJGV4cEdGLzYjLCQqJkY+RilGQUYpRjtGKSwmRitGKS1GLjYjLCQqKiwmRjZGKSomLCZGPkYpRj9GREYpRkFGKUYpRilGQEY7RkFGREYsRitGK0YrRilGOyomRjpGKUZGRilGKUYoRjs= print(` Uf PSt is a solution the following sould give 0 = 0`);
simplify(subs(f(S,t)=PSt,BS_pde)); NiNJVX5VZn5QU3R+aXN+YX5zb2x1dGlvbn50aGV+Zm9sbG93aW5nfnNvdWxkfmdpdmV+MH49fjBHNiI= NiMvIiIhRiQ= PSt_params := subs(params,PSt); NiM+SStQU3RfcGFyYW1zRzYiLCoqJkkiU0dGJSIiIiwmI0YpIiIjRiktSSRlcmZHNiRJKnByb3RlY3RlZEdGMEkoX3N5c2xpYkdGJTYjLCQqKCwoLUkjbG5HRi82IywkRigjRikiI19GKSQiKysrXWlOISM2RilJInRHRiUkISsrKytEckY+RiksJiQiKysrKytdISM1RilGPyEiIiNGRkYsRixGKyQiKysrKys/ISIqRitGKUYpKiYtSSRleHBHRi82IywmJCErKysrKz9GPkYpRj8kIiIlISIjRiksJkYrRiktRi42IywkKigsKEY2RikkIisrKyt2ViEjN0YpRj8kISorKyt2KUY+RilGQkZHRixGK0ZIRitGKSEjX0ZMRjtGKEZG plot0 := plot(subs(t=0,PSt_params),S=40..60):
plot1 := plot(subs(t=1/12,PSt_params),S=40..60):
plot2 := plot(subs(t=2/12,PSt_params),S=40..60):
plot3 := plot(subs(t=3/12,PSt_params),S=40..60):
plot4 := plot(subs(t=4/12,PSt_params),S=40..60):
plot5 := plot(subs(t=5/12,PSt_params),S=40..60):
plot6 := plot(subs(t=6/12,PSt_params),S=40..60):
with(plots):
display([plot0,plot1,plot2,plot3,plot4,plot5,plot6]);
subs(S=55,t=0,PSt_params);
evalf(%); Error, numeric exception: division by zero
LSUlUExPVEc2JS0lJ0NVUlZFU0c2KjdTNyQkIiNTIiIhJCIzTDErVSk9TiJHNiEjOzckJCIzXExMTDNWZlZTRi8kIjNlOCJcWiNIOCozIkYvNyQkIjMhcG07SFtEOjMlRi8kIjNEOyZIK1kwYzAiRi83JCQiM09MTCRlMCQ9Q1RGLyQiM2deaGBEOlE9NUYvNyQkIjNaTEwkM1JCcjslRi8kIjNvLS41QT1hOSkqISM8NyQkIjNcbW0iempmKTRVRi8kIjM2cEg2IikqR0ZYKkZENyQkIjNoTExlNDtbXFVGLyQiM00sSWQmZTpFNypGRDckJCIzRysrRG15XSFIJUYvJCIzWFNYUF4uWCd5KUZENyQkIjNqTExlenMkSEwlRi8kIjNFY1JHOCl5XVcpRkQ3JCQiM0MrK0RAMUJ2VkYvJCIzZ1VCM3A0VDYiKUZENyQkIjM0bm1tQFh0PVdGLyQiM1w2I28wQT9geChGRDckJCIzN0xMJDN5X3FYJUYvJCIzci8pNDdaamFbKEZENyQkIjNLKysrbCs+K1hGLyQiMyZcMzgybWppOyhGRDckJCIzSSsrK3ZXXVZYRi8kIjNVOXheVnleYG9GRDckJCIzUCsrK05mQyZlJUYvJCIzMCMqejFhYGxmbEZENyQkIjNZTExlejY6QllGLyQiM11IZydwQT0kKkgnRkQ3JCQiM3ptbW0iPUMjb1lGLyQiM09tLjRFYyZ6KmZGRDckJCIzUW1tbUVwUzFaRi8kIjNYQVlgdFNzXGRGRDckJCIzJCoqKipcaWBBM3YlRi8kIjNvQ0NVQik0I3BhRkQ3JCQiM3NtbW13eTgheiVGLyQiM0VAUWcybU5HX0ZENyQkIjM4KytET0lGTFtGLyQiM2lkPTJWXj5zXEZENyQkIjMhKioqKlwoM3pNdVtGLyQiMy5Lb3dyQztPWkZENyQkIjNubW07SF8/PFxGLyQiMz8jNEh4JjM1KVwlRkQ3JCQiMylvbTt6aWhsJlxGLyQiMypRJlxTbjknb0clRkQ3JCQiMzlMTCQzI0csKipcRi8kIjMlM3h5QXEmKm8xJUZENyQkIjNXTExlenc1Vl1GLyQiMzY3ZjJyRS1aUUZENyQkIjNSKyt2JFEjXCIzJkYvJCIzRmYtYmxgcWlPRkQ3JCQiM25MTCRlIipbSDcmRi8kIjM0NFwjUiZlJzRaJEZENyQkIjM/KysrcXZ4bF5GLyQiM0pGKiplXld4IUckRkQ3JCQiMzErK11fcW4yX0YvJCIzKmVVOVV2ZkI1JEZENyQkIjNZKytEY3BAW19GLyQiM2RkJ3BeYT1vJEhGRDckJCIzayoqKipcMidIS0gmRi8kIjNLVitdKSpSLWhGRkQ3JCQiM2BtbW13YW5MYEYvJCIzLkgqKT4kNDQsaCNGRDckJCIzWCsrK3YrJ29QJkYvJCIzUip6RHhpcGhYI0ZENyQkIjNDTExlUjwqZlQmRi8kIjMuS1VLc14lSEsjRkQ3JCQiM0IrKysmKUh4ZWFGLyQiMzFMOktxNyVSPSNGRDckJCIzU21tIkghby0qXCZGLyQiM151NTBlNEtmP0ZENyQkIjNLKytEVE81VGJGLyQiMyd6b1VaSnFfJD5GRDckJCIzd21tbVQ5QyNlJkYvJCIzPzlCKVtfWCo+PUZENyQkIjNiKioqXGkhKjNgaSZGLyQiM2hVJyp6cTBLMDxGRDckJCIzKFFMTEwqenltY0YvJCIzMjpUQWJVais7RkQ3JCQiMyVSTEwzTjEjNGRGLyQiM0kxIipRI1JwIipcIkZENyQkIjNFbm0iSFl0N3YmRi8kIjM9KDN1TWNFUlMiRkQ3JCQiM0ArKytxKEcqKnkmRi8kIjNlb1lHUCNmNEsiRkQ3JCQiM0tubTs5QEJNZUYvJCIzJyoqKXoqcCxBNUIiRkQ3JCQiM0BMTExgdiZRKGVGLyQiMy0nejI2W2ddOiJGRDckJCIzcCoqKlxpYDFoImZGLyQiM0BVPyUqR1dgeTVGRDckJCIzZSsrdi5VYWNmRi8kIjMrIVIveT5EJTQ1RkQ3JCQiI2dGLCQiMydSV1ZDdiNvJVIqISM9N1M3JEYqJCIzJFFlKD1OKT1WOCJGLzckRjEkIjNNbiozU2IkPiU0IkYvNyRGNiQiM3UqeTxGKHBrZjVGLzckRjskIjNGOj56UVdCQDVGLzckRkAkIjNSLzJWNEFxSSkqRkQ3JEZGJCIzWjZiTF9aQmMlKkZENyRGSyQiM3IlPTtZTU9UNipGRDckRlAkIjNMSiozKT1NT2woKUZENyRGVSQiM1pQOWpiayEzVClGRDckRlokIjMpPmxGPjc8UjEpRkQ3JEZpbiQiM1dEZnVpM0E5eEZENyRGXm8kIjNMLiUqKVw7UUNUKEZENyRGY28kIjNaKzFFV3glKnpxRkQ3JEZobyQiM2ROISlRX1oyYW5GRDckRl1wJCIzei53JlF3UXlXJ0ZENyRGYnAkIjNJIVsuVzRjbDwnRkQ3JEZncCQiM2wiUTpcXS9FJ2VGRDckRlxxJCIzeDJDNFlANy9jRkQ3JEZhcSQiM3AwclRNIj5BSiZGRDckRmZxJCIzKXlnMz1lNT0xJkZENyRGW3IkIjNbcyEqXFAqeWR6JUZENyRGYHIkIjMmKilRa3gzIyk0YiVGRDckRmVyJCIzKjRsYHJbcVdJJUZENyRGanIkIjN6KiopbykpKjQ3JzMlRkQ3JEZfcyQiM3NRcXRNWD9mUUZENyRGZHMkIjMpMytCb0k9SGokRkQ3JEZpcyQiMyUpbykqKSkpUSdwVk1GRDckRl50JCIzKil6JillJG9zdEMkRkQ3JEZjdCQiM1kiKj4pR3FASzAkRkQ3JEZodCQiMyxhY2NQL3ByR0ZENyRGXXUkIjNmViNbOyVmJFFxI0ZENyRGYnUkIjNvPFF5dlZFRURGRDckRmd1JCIzeWY/alBNV3VCRkQ3JEZcdiQiM0dkIWYpZktCP0FGRDckRmF2JCIzbCdSMWNkZ3QzI0ZENyRGZnYkIjM3OkheUDVPXD5GRDckRlt3JCIzWURMcTVwQ0U9RkQ3JEZgdyQiMypmQ3MtVz1WcSJGRDckRmV3JCIzY2tVJUhBdjpmIkZENyRGanckIjMicD5CdjNeLFsiRkQ3JEZfeCQiMyZwJXAlcDwiKip5OEZENyRGZHgkIjNLVT91UkthIkciRkQ3JEZpeCQiM0c5UC8oKltsIT4iRkQ3JEZeeSQiM3o3JFFgKjMpPjYiRkQ3JEZjeSQiMyd5ZkkoNCdvcy0iRkQ3JEZoeSQiM0NvS0JrPEBpJipGW1tsNyRGXXokIjNlXFVHSDdcXikpRltbbDckRmJ6JCIzJVsxS3k8W1ZAKUZbW2w3JEZneiQiM1YiW1lRSTJWZChGW1tsN1M3JEYqJCIzLXBSMihSV0I5IkYvNyRGMSQiMyVwJDQ4XmwxLDZGLzckRjYkIjM6ImUmUSYpUVhsNUYvNyRGOyQiMyUpPmwpM192ZC0iRi83JEZAJCIzZlolWz00LUYnKSpGRDckRkYkIjNzJSpwaisrR3UlKkZENyRGSyQiM0drJm9JSWkoPSIqRkQ3JEZQJCIzK1VWQmRRbmIoKUZENyRGVSQiM1tTc08qUlRmUSlGRDckRlokIjNlM0ljKW9STy0pRkQ3JEZpbiQiM0B3dDtdSyp5bChGRDckRl5vJCIzJz1LO3opUio9TShGRDckRmNvJCIzV2okZSdSJDRNKnBGRDckRmhvJCIzd3ojKT5WWWZebUZENyRGXXAkIjMhRzZXdVA3LUwnRkQ3JEZicCQiM3gpKVFgV05XWGdGRDckRmdwJCIzaDUjR2xLISplciZGRDckRlxxJCIzKWU+d2g0VVlXJkZENyRGYXEkIjM2RnNCXEddUV5GRDckRmZxJCIzemY7TXp2NndbRkQ3JEZbciQiMzdTIipbcm5wKGYlRkQ3JEZgciQiM0JXXEgoXCgqPU0lRkQ3JEZlciQiM1hAVSE+IUd6JTMlRkQ3JEZqciQiM2EqPktpV2R2JlFGRDckRl9zJCIzLTluZ3gpKSo+aSRGRDckRmRzJCIzd1YkeklLL3lRJEZENyRGaXMkIjMjKik+PFwvLkU+JEZENyRGXnQkIjNPPS1ZNyF5MipIRkQ3JEZjdCQiMyEqPVZVIjMqKT56I0ZENyRGaHQkIjMkM3lBQD9VcGcjRkQ3JEZddSQiM3ZXTjw9VWpPQ0ZENyRGYnUkIjN4WGxkOWdTZEFGRDckRmd1JCIzY0wqUmhxRF01I0ZENyRGXHYkIjNGMFdUYDo6Xj5GRDckRmF2JCIzZE1ac2AoeSQ+PUZENyRGZnYkIjN5MC50JEgsTW8iRkQ3JEZbdyQiM3MqZmo8SkRIYyJGRDckRmB3JCIzIVEtPD1ma1dXIkZENyRGZXckIjN3aUI3Ni13TjhGRDckRmp3JCIzbSZwLlUjPT5IN0ZENyRGX3gkIjNhaCcpPjFxREw2RkQ3JEZkeCQiM1VvKWYrIilbOy8iRkQ3JEZpeCQiM24qWypmOjMpKXAmKkZbW2w3JEZeeSQiMyYqKj5aZ1cjcFYpKUZbW2w3JEZjeSQiM3E6OXUmWyhIcCEpRltbbDckRmh5JCIzRHp3bGJvUEV1RltbbDckRl16JCIzLzpfbCJRbCgqeSdGW1tsNyRGYnokIjMjPj9DQWwqNERpRltbbDckRmd6JCIzI0cmXGdOQThrY0ZbW2w3UzckRiokIjMhcDU0PXNIRzoiRi83JEYxJCIzeGJKSSlIXy82IkYvNyRGNiQiM2NTJ1tiRiJ6dDVGLzckRjskIjNuZSpIWzZLRy4iRi83JEZAJCIzI1xMMGd3LiQ+KipGRDckRkYkIjMnNCVmMyVHX2ZeKkZENyRGSyQiMzRNdykzRkBkOSpGRDckRlAkIjN4bXpNJFx4bHcpRkQ3JEZVJCIzLiUpKipwRXBXeiQpRkQ3JEZaJCIzN1svQyxONSoqekZENyRGaW4kIjNwMiFbaSVwPjl3RkQ3JEZebyQiM0E0aCdvR0k0RyhGRDckRmNvJCIzdXFxbTInM0YicEZENyRGaG8kIjNGckxDdXAiNGInRkQ3JEZdcCQiM24valNwT0g1aUZENyRGYnAkIjMoSDZCVVEmPjNmRkQ3JEZncCQiM1NRPXI2SFFlYkZENyRGXHEkIjN3WSI0RCFmU3FfRkQ3JEZhcSQiM0x6LmFRMlpYXEZENyRGZnEkIjNDLyZwKDNvPG5ZRkQ3JEZbciQiM2wkM0UmeiVIQVAlRkQ3JEZgciQiM1xuUC9Sb3IsVEZENyRGZXIkIjN2cjxpPDRXSVFGRDckRmpyJCIzcUBpMiJlVDhmJEZENyRGX3MkIjNVLSpmSFAqSFdMRkQ3JEZkcyQiM21uNDclMypvKjQkRkQ3JEZpcyQiM19RXGoqXDxuKkdGRDckRl50JCIzJnA1Tl03JCl5byNGRDckRmN0JCIzVy12PSQ0dExbI0ZENyRGaHQkIjNTWl1Bdio+VUgjRkQ3JEZddSQiM1gjZnFaU0I4NyNGRDckRmJ1JCIzaiF6cFxHeDIlPkZENyRGZ3UkIjNjOiVvdXdoJil5IkZENyRGXHYkIjMpeSpcPDckNGlqIkZENyRGYXYkIjMxKT4peUE/JnBdIkZENyRGZnYkIjNZcXJgMz8pW1AiRkQ3JEZbdyQiM14kbyQ0amg1ZjdGRDckRmB3JCIzOy5VIXlrR2w5IkZENyRGZXckIjNNNVJHK0dXVzVGRDckRmp3JCIzSTIhUkQlPj5jJSpGW1tsNyRGX3gkIjNbIil5YCEzSSR5JilGW1tsNyRGZHgkIjMiKkcpKSk+a2g6dihGW1tsNyRGaXgkIjNLKW9dO1k4JikqcEZbW2w3JEZeeSQiM0shWydwbXQlPU8nRltbbDckRmN5JCIzPUc9cDclKUckcCZGW1tsNyRGaHkkIjMoZnBfZlxfcjkmRltbbDckRl16JCIzekNOJCp6TDg6WUZbW2w3JEZieiQiMztedXdbP0FeVEZbW2w3JEZneiQiM1c1PFUyLl0pcCRGW1tsN1M3JEYqJCIzJHlKRmtAaWo7IkYvNyRGMSQiM3pfIVxkSkRKNyJGLzckRjYkIjNjJj1sKVwkKmYmMyJGLzckRjskIjN6YiVcKjR6X1Y1Ri83JEZAJCIzJEctc3MvXDgrIkYvNyRGRiQiMzN2Oklfb3omZipGRDckRkskIjMjRyJSJDMmUXk1IypGRDckRlAkIjNTMjBYMlwjWyIpKUZENyRGVSQiMyg0NStCRTwoMyUpRkQ3JEZaJCIzWmtpJ3okcCJ6KylGRDckRmluJCIzIylmOFRfQ1ord0ZENyRGXm8kIjNsJlt1ITQmPmlDKEZENyRGY28kIjN0LSRvQilIRmBvRkQ3JEZobyQiMz1XNkV6R3Fsa0ZENyRGXXAkIjN3ekdSNSF6Jio0J0ZENyRGYnAkIjNIclMmM3AlKVJ4JkZENyRGZ3AkIjM3UDlVSDg2J1ImRkQ3JEZccSQiM3E2bV07c14lMyZGRDckRmFxJCIzOzwiRzEpemdLWkZENyRGZnEkIjNzNmN2JGZlNlYlRkQ3JEZbciQiMzUqKWZiTWMpPTYlRkQ3JEZgciQiMykqUWtCJmZIJj5RRkQ3JEZlciQiMyZwelZENC5yXyRGRDckRmpyJCIzLjwsay9qRXFLRkQ3JEZfcyQiM21TbXAxdjQxSUZENyRGZHMkIjMoPiMzYEkzM1lGRkQ3JEZpcyQiMyJcOF8nZiIqeUpERkQ3JEZedCQiM3JgVENWbCpISiNGRDckRmN0JCIzRiRRXyV5JT4yNSNGRDckRmh0JCIzc14hUSEqUmJrIT5GRDckRl11JCIzdiQ+NXFVRjR0IkZENyRGYnUkIjNDXldDSEIwXTpGRDckRmd1JCIzbXckM0ZNcigqUiJGRDckRlx2JCIzaCEpeVElKTNvXjdGRDckRmF2JCIzYUk2QU4iNCJHNkZENyRGZnYkIjNnb1YjUkpwUysiRkQ3JEZbdyQiM0dac3EoWyQzdSopRltbbDckRmB3JCIzKD1rIj5oUSN5JnpGW1tsNyRGZXckIjNORm9JPyR5aDAoRltbbDckRmp3JCIzJT0kKWY3YCxOPydGW1tsNyRGX3gkIjNSKHArb0ciZmthRltbbDckRmR4JCIzTC9OTi5kVid5JUZbW2w3JEZpeCQiMz9fNHBdK1EmPSVGW1tsNyRGXnkkIjM9aW5ZTk0qM3AkRltbbDckRmN5JCIzdzIoNF9yXGw9JEZbW2w3JEZoeSQiMyMzTiIqUlNgcXkjRltbbDckRl16JCIzJGVdayZRLyIpNENGW1tsNyRGYnokIjNHYis0YkpVIjQjRltbbDckRmd6JCIzJW8vSSRbZjsieiJGW1tsN1M3JEYqJCIzcih5ODwlKjNGPSJGLzckRjEkIjMyPSFSRXlEIlI2Ri83JEY2JCIzenohR2x0NTc1IkYvNyRGOyQiM2swd3M2RWVlNUYvNyRGQCQiMzxKYDUnZiFwOjVGLzckRkYkIjNbZm4oWyw0LnQqRkQ3JEZLJCIzT2MudikqUjxOJCpGRDckRlAkIjNqOFIhM00iZUUqKUZENyRGVSQiM3g2JlFgMlFbXSlGRDckRlokIjN1OyYqR2dXZyYzKUZENyRGaW4kIjMvPUNeLCo9Z2woRkQ3JEZebyQiM1xjdS87V1p6c0ZENyRGY28kIjMqNCFvPUBSRGVvRkQ3JEZobyQiMzF3aWNQWSYqUWtGRDckRl1wJCIzYl5CLi4sRFJnRkQ3JEZicCQiMzpSKSlcPiwlM28mRkQ3JEZncCQiM0UjXGdRVyVcaF9GRDckRlxxJCIzVi5uKHouUEoiXEZENyRGYXEkIjMhUVxMV1dac14lRkQ3JEZmcSQiM1p5OHZJP1Z3VEZENyRGW3IkIjNFPEghR08uViJRRkQ3JEZgciQiMyM0LUktJ3BFI1skRkQ3JEZlciQiM3g2YGM0IXkvOiRGRDckRmpyJCIzKFJIPEkoPTFnR0ZENyRGX3MkIjMrVSMpR2BtO2pERkQ3JEZkcyQiM2RnQTcjZjlQRiNGRDckRmlzJCIzRTtzJVFVcSFRP0ZENyRGXnQkIjNsPztka04xLD1GRDckRmN0JCIzOiZwI3Arel52OkZENyRGaHQkIjNsTl5PXislUVAiRkQ3JEZddSQiMzstZiIqKnpQaj4iRkQ3JEZidSQiMzZRQCMpUidbIT41RkQ3JEZndSQiM1pjKVE5J1tmbygpRltbbDckRlx2JCIzYHl5NmlOVj91RltbbDckRmF2JCIzbGlrLVBJcVRqRltbbDckRmZ2JCIzNCFmbGchUTcyYEZbW2w3JEZbdyQiMytYQDg7KDM2WSVGW1tsNyRGYHckIjNjXiNlW2Yuc3AkRltbbDckRmV3JCIzKCopSCd6Mk8peTAkRltbbDckRmp3JCIzSSF6emBLXS1cI0ZbW2w3JEZfeCQiM0d5L1tFWFtJP0ZbW2w3JEZkeCQiMzAhcGdUV1h1aiJGW1tsNyRGaXgkIjMiRydcXTdDWDk4RltbbDckRl55JCIzKnlwSEohKik+bzVGW1tsNyRGY3kkIjNhIkg4LVJwck8pISM+NyRGaHkkIjNXcUhDJlxoaG8nRmRobjckRl16JCIzSV5CT1Axa0tfRmRobjckRmJ6JCIzX0p6b2R3bTlURmRobjckRmd6JCIzZlNHRkYjKVJlSkZkaG43VTckRiokIiM3Riw3JEYxJCIzXm1tbSJwMGs6IkYvNyRGNiQiMzVMTDM8WFo9NkYvNyRGOyQiM2ttbTtXcCJlMiJGLzckRkAkIjNgbW07NG0oRy4iRi83JEZGJCIzOU5MJDNpLjkhKipGRDckRkskIjMjUm1tVCFSPTAmKkZENyRGUCQiMzMoKioqXFA4I1w0KkZENyRGVSQiM2xqbTsvc2lxJylGRDckRlokIjNnKCoqKlwoeSRwWiMpRkQ3JEZpbiQiMzNITEwkeWFFInlGRDckRl5vJCIzJClvbW0iPnMlSHVGRDckRmNvJCIzJW8qKioqKlwkKjQpKnBGRDckRmhvJCIzLSgqKioqKlxfJlxjJ0ZENyRGXXAkIjNHJyoqKioqXDFhWmhGRDckRmJwJCIzVGxtOy8jKVtvZEZENyRGZ3AkIjMxS0xMJD1leEomRkQ3JEZccSQiM0BPTExMMiRmJFxGRDckRmFxJCIzcysrXVBZeCJcJUZENyRGZnEkIjNzS0xMTDdpKTQlRkQ3JEZbciQiM24pKioqXFAncHNtJEZENyRGYHIkIjMpNCsrRCI0X2NLRkQ3JEZlciQiM0pMTEwzeCV6I0dGRDckRmpyJCIzQEpMJDNzJFFNQ0ZENyRGX3MkIjNqb21tInpyKTQ/RkQ3JEZkcyQiM2tsbTsvSyMqbzpGRDckRmlzJCIzOScqKipcaWgyJj0iRkQ3JEZedCQiMyNHam1tVDNecShGW1tsNyRGY3QkIjNRISkqKioqKipIQ0FNRltbbDckJCIzOCsrRDZ0cyc9JkYvJCIzRCgpKioqXCgpb3NLIkZbW2w3JEZodCRGLEYsNyQkIjNFK11QL3ElekEmRi9GX19vNyRGXXVGX19vNyRGYnVGX19vNyRGZ3VGX19vNyRGXHZGX19vNyRGYXZGX19vNyRGZnZGX19vNyRGW3dGX19vNyRGYHdGX19vNyRGZXdGX19vNyRGandGX19vNyRGX3hGX19vNyRGZHhGX19vNyRGaXhGX19vNyRGXnlGX19vNyRGY3lGX19vNyRGaHlGX19vNyRGXXpGX19vNyRGYnpGX19vNyRGZ3pGX19vLSUmQ09MT1JHNjglJFJHQkckIiM1ISIiJEYsRlxhb0ZdYW9GamBvRl1hb0ZdYW9GamBvRl1hb0ZdYW9GamBvRl1hb0ZdYW9GamBvRl1hb0ZdYW9GamBvRl1hb0ZdYW9GamBvRl1hb0ZdYW8tJStBWEVTTEFCRUxTRzYnUSJTNiJRIUZiYW8tJSVGT05URzYkJSpIRUxWRVRJQ0FHRlthbyUrSE9SSVpPTlRBTEdGaGFvLSUlVklFV0c2JDtGKkZnejskISNDISIjJCIyLSsrKysrU0EiISM6 NiMsKiMhI2IiIiMiIiItSSRlcmZHNiRJKnByb3RlY3RlZEdGK0koX3N5c2xpYkc2IjYjLCQqJiwmLUkjbG5HRio2IyMiI2IiI19GJyQiKysrXWlOISM2RidGJ0YmI0YnRiYkIitDclVHRyEiKiNGNkYmKiYtSSRleHBHRio2IyQhKysrKys/RjpGJywmRjtGJy1GKTYjLCQqJiwmRjJGJyQiKysrK3ZWISM3RidGJ0YmRjtGPEY7RichI19GQUY3 NiMkIitYMlFjPyEiKg==