Solving the two mass system with linear algebra LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEocmVzdGFydEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNi1RIjtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRjEvJSlzdHJldGNoeUdGPS8lKnN5bW1ldHJpY0dGPS8lKGxhcmdlb3BHRj0vJS5tb3ZhYmxlbGltaXRzR0Y9LyUnYWNjZW50R0Y9LyUnbHNwYWNlR1EmMC4wZW1GJy8lJ3JzcGFjZUdRLDAuMjc3Nzc3OGVtRic= Here is the system of equations: JSFH QyQ+SSRlcTFHNiIvLUklZGlmZkclKnByb3RlY3RlZEc2JS1JI3UxR0YlNiNJInRHRiVGLkYuLCZGKyEiJi1JI3UyR0YlRi0iIiMiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtRiQ2JC1JI3UxRzYiNiNJInRHRitGLUYtLCZGKSEiJi1JI3UyR0YrRiwiIiM=

QyQ+SSRlcTJHNiIvLUklZGlmZkclKnByb3RlY3RlZEc2JS1JI3UyR0YlNiNJInRHRiVGLkYuLCYtSSN1MUdGJUYtIiIjRishIiMiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtRiQ2JC1JI3UyRzYiNiNJInRHRitGLUYtLCYtSSN1MUdGK0YsIiIjRikhIiM=

First we'll just have Maple solve the system: QyQtSSdkc29sdmVHNiI2IzwkSSRlcTFHRiVJJGVxMkdGJSIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

PCQvLUkjdTFHNiI2I0kidEdGJiwqKiZJJF9DMUdGJiIiIi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJjYjKiYiIicjRiwiIiNGKEYsRixGLComSSRfQzJHRiZGLC1JJGNvc0dGL0YyRixGLComSSRfQzNHRiZGLC1GLkYnRixGLComSSRfQzRHRiZGLC1GOkYnRixGLC8tSSN1MkdGJkYnLCpGKiMhIiJGNkY3RkVGO0Y2Rj5GNg==

QyQtSSV3aXRoRzYiNiNJLkxpbmVhckFsZ2VicmFHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

N2VySSMmeEc2IkkkQWRkR0YkSShBZGpvaW50R0YkSTNCYWNrd2FyZFN1YnN0aXR1dGVHRiRJK0JhbmRNYXRyaXhHRiRJJkJhc2lzR0YkSS1CZXpvdXRNYXRyaXhHRiRJL0JpZGlhZ29uYWxGb3JtR0YkSS1CaWxpbmVhckZvcm1HRiRJNUNoYXJhY3RlcmlzdGljTWF0cml4R0YkSTlDaGFyYWN0ZXJpc3RpY1BvbHlub21pYWxHRiRJJ0NvbHVtbkdGJEkwQ29sdW1uRGltZW5zaW9uR0YkSTBDb2x1bW5PcGVyYXRpb25HRiRJLENvbHVtblNwYWNlR0YkSTBDb21wYW5pb25NYXRyaXhHRiRJMENvbmRpdGlvbk51bWJlckdGJEkvQ29uc3RhbnRNYXRyaXhHRiRJL0NvbnN0YW50VmVjdG9yR0YkSSVDb3B5R0YkSTJDcmVhdGVQZXJtdXRhdGlvbkdGJEktQ3Jvc3NQcm9kdWN0R0YkSS1EZWxldGVDb2x1bW5HRiRJKkRlbGV0ZVJvd0dGJEksRGV0ZXJtaW5hbnRHRiRJKURpYWdvbmFsR0YkSS9EaWFnb25hbE1hdHJpeEdGJEkqRGltZW5zaW9uR0YkSStEaW1lbnNpb25zR0YkSStEb3RQcm9kdWN0R0YkSTZFaWdlbkNvbmRpdGlvbk51bWJlcnNHRiRJLEVpZ2VudmFsdWVzR0YkSS1FaWdlbnZlY3RvcnNHRiRJJkVxdWFsR0YkSTJGb3J3YXJkU3Vic3RpdHV0ZUdGJEkuRnJvYmVuaXVzRm9ybUdGJEk0R2F1c3NpYW5FbGltaW5hdGlvbkdGJEkyR2VuZXJhdGVFcXVhdGlvbnNHRiRJL0dlbmVyYXRlTWF0cml4R0YkSShHZW5lcmljR0YkSTJHZXRSZXN1bHREYXRhVHlwZUdGJEkvR2V0UmVzdWx0U2hhcGVHRiRJNUdpdmVuc1JvdGF0aW9uTWF0cml4R0YkSSxHcmFtU2NobWlkdEdGJEktSGFua2VsTWF0cml4R0YkSSxIZXJtaXRlRm9ybUdGJEkzSGVybWl0aWFuVHJhbnNwb3NlR0YkSS9IZXNzZW5iZXJnRm9ybUdGJEkuSGlsYmVydE1hdHJpeEdGJEkySG91c2Vob2xkZXJNYXRyaXhHRiRJL0lkZW50aXR5TWF0cml4R0YkSTJJbnRlcnNlY3Rpb25CYXNpc0dGJEkrSXNEZWZpbml0ZUdGJEktSXNPcnRob2dvbmFsR0YkSSpJc1NpbWlsYXJHRiRJKklzVW5pdGFyeUdGJEkySm9yZGFuQmxvY2tNYXRyaXhHRiRJK0pvcmRhbkZvcm1HRiRJMUtyb25lY2tlclByb2R1Y3RHRiRJKExBX01haW5HRiRJMExVRGVjb21wb3NpdGlvbkdGJEktTGVhc3RTcXVhcmVzR0YkSSxMaW5lYXJTb2x2ZUdGJEkkTWFwR0YkSSVNYXAyR0YkSSpNYXRyaXhBZGRHRiRJMk1hdHJpeEV4cG9uZW50aWFsR0YkSS9NYXRyaXhGdW5jdGlvbkdGJEkuTWF0cml4SW52ZXJzZUdGJEk1TWF0cml4TWF0cml4TXVsdGlwbHlHRiRJK01hdHJpeE5vcm1HRiRJLE1hdHJpeFBvd2VyR0YkSTVNYXRyaXhTY2FsYXJNdWx0aXBseUdGJEk1TWF0cml4VmVjdG9yTXVsdGlwbHlHRiRJMk1pbmltYWxQb2x5bm9taWFsR0YkSSZNaW5vckdGJEkoTW9kdWxhckdGJEkpTXVsdGlwbHlHRiRJLE5vVXNlclZhbHVlR0YkSSVOb3JtR0YkSSpOb3JtYWxpemVHRiRJKk51bGxTcGFjZUdGJEkzT3V0ZXJQcm9kdWN0TWF0cml4R0YkSSpQZXJtYW5lbnRHRiRJJlBpdm90R0YkSSpQb3BvdkZvcm1HRiRJMFFSRGVjb21wb3NpdGlvbkdGJEktUmFuZG9tTWF0cml4R0YkSS1SYW5kb21WZWN0b3JHRiRJJVJhbmtHRiRJNlJhdGlvbmFsQ2Fub25pY2FsRm9ybUdGJEk2UmVkdWNlZFJvd0VjaGVsb25Gb3JtR0YkSSRSb3dHRiRJLVJvd0RpbWVuc2lvbkdGJEktUm93T3BlcmF0aW9uR0YkSSlSb3dTcGFjZUdGJEktU2NhbGFyTWF0cml4R0YkSS9TY2FsYXJNdWx0aXBseUdGJEktU2NhbGFyVmVjdG9yR0YkSSpTY2h1ckZvcm1HRiRJL1Npbmd1bGFyVmFsdWVzR0YkSSpTbWl0aEZvcm1HRiRJOFN0cm9uZ2x5Q29ubmVjdGVkQmxvY2tzR0YkSSpTdWJNYXRyaXhHRiRJKlN1YlZlY3RvckdGJEkpU3VtQmFzaXNHRiRJMFN5bHZlc3Rlck1hdHJpeEdGJEkvVG9lcGxpdHpNYXRyaXhHRiRJJlRyYWNlR0YkSSpUcmFuc3Bvc2VHRiRJMFRyaWRpYWdvbmFsRm9ybUdGJEkrVW5pdFZlY3RvckdGJEkyVmFuZGVybW9uZGVNYXRyaXhHRiRJKlZlY3RvckFkZEdGJEksVmVjdG9yQW5nbGVHRiRJNVZlY3Rvck1hdHJpeE11bHRpcGx5R0YkSStWZWN0b3JOb3JtR0YkSTVWZWN0b3JTY2FsYXJNdWx0aXBseUdGJEkrWmVyb01hdHJpeEdGJEkrWmVyb1ZlY3RvckdGJEkkWmlwR0Yk

Here is the coefficient matrix for the system: JSFH QyQ+SSJBRzYiLUknTWF0cml4R0YlNiM3JDckISImIiIjNyRGLCEiIyIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKSM0Iyo+Jg==

JSFH This matrix is symmetric, so it can be diagonalized and it has real eigenvalues. JSFH QyQ+SSJURzYiLUktRWlnZW52ZWN0b3JzR0YlNiNJIkFHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

NiQtJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiNJJ2NvbHVtbkdGKTYjL0kkJWlkR0YpIiljSFFsLUknTWF0cml4R0YmNiMvRi4iKTs4UWw=

JSFH Let's call P the matrix whose columns are independent eigenvectors of A. P will diagonalize the matrix A. JSFH QyQ+SSJQRzYiJkkiVEdGJTYjIiIjIiIi

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKTs4UWw=

QyQ+SSJFRzYiLUkwZGVsYXlEb3RQcm9kdWN0RzYkJSpwcm90ZWN0ZWRHL0krbW9kdWxlbmFtZUdGJUksVHlwZXNldHRpbmdHNiRGKUkoX3N5c2xpYkdGJTYkLUYnNiQqJEkiUEdGJSEiIkkiQUdGJUYzIiIi

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKUdoU20=

QyQqJEkiUEc2IiEiIiIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKV9cJXkn

U will be the column vector of the original functions we are solving for. QyQ+SSJVRzYiLUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiM3JDcjLUkjdTFHRiU2I0kidEdGJTcjLUkjdTJHRiVGMCIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKU8pRyRv

V will be the column vector of the new variables. JSFH QyQ+SSJWRzYiLUkwZGVsYXlEb3RQcm9kdWN0RzYkJSpwcm90ZWN0ZWRHL0krbW9kdWxlbmFtZUdGJUksVHlwZXNldHRpbmdHNiRGKUkoX3N5c2xpYkdGJTYkKiRJIlBHRiUhIiJJIlVHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKTdLMlw=

JSFH The sysyem of equations for the new functions V corresponds to the diagonal matrix D=P^(-1)AP. QyQ+SSRlcTNHNiIvLUklZGlmZkclKnByb3RlY3RlZEc2JS1JI3YxR0YlNiNJInRHRiVGLkYuLCRGKyEiIiIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtRiQ2JC1JI3YxRzYiNiNJInRHRitGLUYtLCRGKSEiIg==

QyQ+SSRlcTRHNiIvLUklZGlmZkclKnByb3RlY3RlZEc2JS1JI3YyR0YlNiNJInRHRiVGLkYuLCRGKyEiJyIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtRiQ2JC1JI3YyRzYiNiNJInRHRitGLUYtLCRGKSEiJw==

QyQ+SSJURzYiLUknZHNvbHZlRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiM8JEkkZXEzR0YlSSRlcTRHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

PCQvLUkjdjFHNiI2I0kidEdGJiwmKiZJJF9DM0dGJiIiIi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJkYnRixGLComSSRfQzRHRiZGLC1JJGNvc0dGL0YnRixGLC8tSSN2MkdGJkYnLCYqJkkkX0MxR0YmRiwtRi42IyomIiInI0YsIiIjRihGLEYsRiwqJkkkX0MyR0YmRiwtRjVGPUYsRiw=

JSFH Let's check this by solving the system with some initial conditions. First we'll just have Maple solve for us. JSFH QyQ+SSNTMUc2Ii1JJ2Rzb2x2ZUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYkNyhJJGVxMUdGJUkkZXEyR0YlLy1JI3UxR0YlNiMiIiEiIiIvLUkjdTJHRiVGMiIiJC8tLUkiREdGKDYjRjFGMiMhIiIiIiMvLS1GPDYjRjdGMkY0PCQtRjE2I0kidEdGJS1GN0ZHRjQ=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

PCQvLUkjdTFHNiI2I0kidEdGJiwqKiYtSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiY2IyomIiInIyIiIiIiI0YoRjRGNEYyRjMjISIjIiM6LUkkY29zR0YtRjAjRjciIiYtRixGJyMiIiQiIzUtRjpGJyMiIihGPC8tSSN1MkdGJkYnLCpGKiNGNEY4RjkjRjRGPEY9I0Y/RjxGQSMiIzlGPA==

JSFH Here are the solutions for u1 and u2 extracted as expressions. JSFH QyQ+SSNGMUc2Ii1JJXN1YnNHJSpwcm90ZWN0ZWRHNiQmSSNTMUdGJTYjIiIiLUkjdTFHRiU2I0kidEdGJUYt

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LCoqJi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjKiYiIicjIiIiIiIjSSJ0R0YpRi5GLkYsRi0jISIjIiM6LUkkY29zR0YmRiojRjIiIiYtRiU2I0YwIyIiJCIjNS1GNUY5IyIiKEY3

QyQ+SSNGMkc2Ii1JJXN1YnNHJSpwcm90ZWN0ZWRHNiQmSSNTMUdGJTYjIiIjLUkjdTJHRiU2I0kidEdGJSIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LCoqJi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjKiYiIicjIiIiIiIjSSJ0R0YpRi5GLkYsRi0jRi4iIzotSSRjb3NHRiZGKiNGLiIiJi1GJTYjRjAjIiIkRjYtRjRGOCMiIzlGNg==

Now let's do it ourselves. First we get the initial conditions in ters of the new variables V. Remember that V=P^(-1)*U. First we compute the column vector giving v1(0) and v2(0), and then we compute the column vector giving v1'(0), v2'(0). JSFH JSFH QyQtSTBkZWxheURvdFByb2R1Y3RHNiQlKnByb3RlY3RlZEcvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0c2JEYmSShfc3lzbGliR0YpNiQqJEkiUEdGKSEiIi1JJ01hdHJpeEdGKzYjNyQ3IyIiIjcjIiIkRjY=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKSVbNGcn

QyQtSTBkZWxheURvdFByb2R1Y3RHNiQlKnByb3RlY3RlZEcvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0c2JEYmSShfc3lzbGliR0YpNiQqJEkiUEdGKSEiIi1JJ01hdHJpeEdGKzYjNyQ3IyNGMCIiIzcjIiIiRjk=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKVs3XGw=

Now we solve the diagonal system with these initial conditions for the new variables. JSFH QyQ+SSNNMUc2Ii1JJ2Rzb2x2ZUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYjNyVJJGVxM0dGJS8tSSN2MUdGJTYjIiIhIyIjOSIiJi8tLUkiREdGKDYjRjBGMSMiIiRGNSIiIg==

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JI3YxRzYiNiNJInRHRiUsJi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJUYmIyIiJCIiJi1JJGNvc0dGK0YmIyIjOUYw

QyQ+SSNNMkc2Ii1JJ2Rzb2x2ZUc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYjNyVJJGVxNEdGJS8tSSN2MkdGJTYjIiIhIyIiIiIiJi8tLUkiREdGKDYjRjBGMSMiIiNGNUY0

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

Ly1JI3YyRzYiNiNJInRHRiUsJiomLUkkc2luRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiMqJiIiJyMiIiIiIiNGJ0YzRjNGMUYyI0YzIiM6LUkkY29zR0YsRi8jRjMiIiY=

QyQ+SSN2MUc2Ii1JJXN1YnNHJSpwcm90ZWN0ZWRHNiRJI00xR0YlLUYkNiNJInRHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LCYtSSRzaW5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2I0kidEdGKCMiIiQiIiYtSSRjb3NHRiVGKSMiIzlGLQ==

QyQ+SSN2Mkc2Ii1JJXN1YnNHJSpwcm90ZWN0ZWRHNiRJI00yR0YlLUYkNiNJInRHRiUiIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LCYqJi1JJHNpbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjKiYiIicjIiIiIiIjSSJ0R0YpRi5GLkYsRi0jRi4iIzotSSRjb3NHRiZGKiNGLiIiJg==

JSFH Finally, we convert back to the original functions u1, u2 and see that we get the same solution. JSFH QyQtSTBkZWxheURvdFByb2R1Y3RHNiQlKnByb3RlY3RlZEcvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0c2JEYmSShfc3lzbGliR0YpNiRJIlBHRiktSSdNYXRyaXhHRis2IzckNyNJI3YxR0YpNyNJI3YyR0YpIiIi

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

print(); # input placeholder

LUknTWF0cml4RzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMvSSQlaWRHRiciKVt6KEgn

Here is a plot of the solutions. JSFH QyQtSSVwbG90RzYiNiQ3JEkjRjFHRiVJI0YyR0YlL0kidEdGJTsiIiEiIz8iIiI=

LUkjbWlHNiMvSSttb2R1bGVuYW1lRzYiSSxUeXBlc2V0dGluZ0dJKF9zeXNsaWJHRic2JVE1b3V0cHV0fnJlZGlyZWN0ZWQuLi5GJy8lJ2l0YWxpY0dRJXRydWVGJy8lLG1hdGh2YXJpYW50R1EnaXRhbGljRic=

%;

NictJSdDVVJWRVNHNiQ3YXY3JCQiIiEhIiIkIjEqKioqKioqKioqKioqKioqISM7NyQkIjJMTEwkM3gmKSozIiEjPCQiMSl5SE4nM1BCJiohIzs3JCQiMm1tbW1UOih6QCEjPCQiMSU9MD9mNEBAKiEjOzckJCIxTEwkM0ZXWXMjISM7JCIxIltaZSpHbUYiKiEjOzckJCIrREpkcEshIzUkIjFKM09QZWsjNCohIzs3JCQiMW5tO3o+XTlRISM7JCIxMVkqWzFlbTUqISM7NyQkIjFMTExMM1ZmViEjOyQiMTwkUTo+LHk7KiEjOzckJCIsRGMqKWZEJyEjNiQiMSg9IT0mXCE+KW8qISM7NyQkIjFubW0iSFtEOikhIzskIjF2RypwYCFIWTUhIzo3JCQiMU1MZTl3KSo9IyohIzskIjBBNyRcT1YnMyIhIzk3JCQiLXYkcFUmRzUhIzYkIjI3LVdqd2FUNiIhIzs3JCQiMkskZSopZlknPTMiISM7JCIyPU4tZSpHKTM3IiEjOzckJCIybG1UZ2knPU42ISM7JCIxKTRTImV6Z0A2ISM6NyQkIi92PSNmMyYpPSIhIzgkIjE8M1NvMVo6NiEjOjckJCIyTExMJGUwJD1DIiEjOyQiMmgpZlJQUXAsNiEjOzckJCIyTkxMTEJLbFgiISM7JCIxWUdOKTQhKipvJiohIzs3JCQiMk9MTCQzUkJyOyEjOyQiMURLQC8wamhsISM7NyQkIjJvbVRnX3UheTwhIzskIjJMMlFxeCRlN1ghIzw3JCQiMi0rXVA5OlwpPSEjOyQiMDlTa0pIJGZAISM6NyQkIjJNTGU5d2I8Kj4hIzskITJXMD1VcXpYQyUhIz03JCQiMm1tbSJ6amYpNCMhIzskITEpejVdJj5YW0ohIzs3JCQiMk9MTCQzOGwoPiMhIzskITFVTT5rInAwciYhIzs3JCQiLHZCMW5IIyEjNSQhMGtQJCk0Qi1BKSEjOjckJCIybG1tbTtoZFIjISM7JCEyUG9OWj4+I2Y1ISM7NyQkIjJNTExlNDtbXCMhIzskITIlKj1MYSZbZnU3ISM7NyQkIjFubTt6dCUqKnAjISM6JCEyUFAyQFlYPGkiISM7NyQkIixEbXldIUghIzUkITIyIioqKikzKHk5IT0hIzs3JCQiMk9MTGUvKGZKSCEjOyQhMkxUZF1JTTkiPSEjOzckJCIxbm07SGE2ZUghIzokITJ6eG8ieUBNPT0hIzs3JCQiMi4rK0QiUWolKUghIzskITIqeWYjeTZFQSM9ISM7NyQkIjJPTExlPl82LCQhIzskITJLWWk5Uz5KIz0hIzs3JCQiMW5tO3owblBJISM6JCEyMiFRMj5ZMUA9ISM7NyQkIjIuKytEJyopPWtJISM7JCEyTSYzKTRZOGgiPSEjOzckJCIyT0xMZU0yMjQkISM7JCEyJ29OYnBrSzM9ISM7NyQkIjFubTtIZEE8SiEjOiQhMk8oUSdmdnN4eiIhIzs3JCQiMi4rK0RFKkhCSyEjOyQhMSVmJjNrVGJIPCEjOjckJCIyUExMZXpzJEhMISM7JCEyYCZbIjRMI2VDOyEjOzckJCIybm1tVF1SM2EkISM7JCEyMD0nSEo6Q1A4ISM7NyQkIixEQDFCdiQhIzUkITFHIz4mcEJUMDUhIzo3JCQiMU1MZTlkIylwUiEjOiQhMSQ9PmZ1PjsqbyEjOzckJCIxbm1tO19NKD0lISM6JCEyWHZmW1IzYV0lISM8NyQkIjEsK103bCQqeVYhIzokITJrMDlaLjtVPyQhIzw3JCQiMU1MTDN5X3FYISM6JCEyOXNPNio0YVVEISM8NyQkIjFubTtIVUAneSUhIzokITJ5JCo+JSopKT5QPSMhIzw3JCQiMjArKytsKz4rJiEjOyQhMlBdVUp5JylHbiIhIzw3JCQiKUZaPV8hIigkITFaMiQqKWVeI1xXISM8NyQkIip2V11WJiEiKSQiMmVIKDMzSDRiPSEjPDckJCIpJSlSUmIhIigkIjFBUyxJIip6bEwhIzs3JCQiKjBfUGsmISIpJCIxLzlHRDJJMF4hIzs3JCQiKWQ1W2QhIigkIjEpUUxyInkpUi0oISM7NyQkIipOZkMmZSEiKSQiMUIib19fWl0wKiEjOzckJCIxbm0iSGQmKT4vJyEjOiQiMjojNE9JKSlSdzchIzs3JCQiMU1MJGV6NjpCJyEjOiQiMVBnWiJ5KWUuOyEjOjckJCIxbjsvLFY+V2ohIzokIjItYlJmcnRhdiIhIzs3JCQiLUQxbyhvWCchIzYkIjI4akk8d3pPJz0hIzs3JCQiMW5UJillIT1LXichIzokIjJCPGtTISlSISoqPSEjOzckJCIxTSRlOUpmJnBsISM6JCIyJlIkWyJRbic0Iz4hIzs3JCQiL0QxazAhZmknISM4JCIyPDJ2LlYkKSpHPiEjOzckJCIxbm1tOz1DI28nISM6JCIyRVtAPlswRyM+ISM7NyQkIjFvbW1UYjp0byEjOiQiMm9kU1ZoInkneiIhIzs3JCQiMW5tbW0jcFMxKCEjOiQiMnAhM0s2IVJBXyIhIzs3JCQiLkQxSDNePCghIzckIjIpKikqbyEpeSoqekkiISM7NyQkIjFMTGU5dDknRyghIzokIjJLelZYRylvbDUhIzs3JCQiMW47YVFqPShSKCEjOiQiMSNIVFIhRyllMSkhIzs3JCQiLERPRCMzdiEjNSQiMS9QdDdoZkNhISM7NyQkIjFtO2E4UF4xdyEjOiQiL1hLTkQjbzgkISM5NyQkIjFMTGVrPyFbcSghIzokIjF4JXA3Jyo+X3MqISM8NyQkIi5EY1QhNC55ISM3JCEyZDlTdWdUYSsiISM8NyQkIjFtbW1tKHk4IXohIzokITJ2OCozZ0tIWkYhIzw3JCQiMU1MZWtYMDwiKSEjOiQhMS1tY0FbbyhmJiEjOzckJCIxLCtdaS50SyQpISM6JCExclguIyopWyw0KCEjOzckJCIuREomUjIlUSkhIzckITFpJnosLF4pcHMhIzs3JCQiLXZWdlROJSkhIzYkITFzaT9vQiJRUighIzs3JCQiMSxdUE02dydbKSEjOiQhMSdvSGpieiJwdSEjOzckJCIrRFo1USYpISIqJCExdWttPnFvLnYhIzs3JCQiMSwrRDE+elMnKSEjOiQhMUtGcyc+PkpbKCEjOzckJCIsdjN6TXUpISM1JCExX0wxSk8xK3UhIzs3JCQiMU0zLVEyMCh6KSEjOiQhMXVrPypITVlOKCEjOzckJCIxbzthKVFBMSYpKSEjOiQhMSczZnBdRik+dCEjOzckJCIxLEQxUlM+LyopISM6JCEwTENveDlVSSghIzo3JCQiMU5MZSpvbHgmKikhIzokITFYQDAqUUBjSighIzs3JCQiMW9UNVN0TDYhKiEjOiQhMVBKJFw5cTRPKCEjOzckJCIxLV1pISoqM1wxKiEjOiQhMSY+blchPjlZdSEjOzckJCIxTmU5VDFbPSIqISM6JCExLF0vYmAmZWQoISM7NyQkIjFvbW0iSF8/PCohIzokITFPVExMIm9OdighIzs3JCQiMW9tVCZHTSlvJCohIzokITEuW1RUOWhQKSkhIzs3JCQiMW5tO3ppaGwmKiEjOiQhMipbcShwQndAMCIhIzs3JCQiLXZWQSh5eCohIzYkITFKN2Q/YiFvRSIhIzo3JCQiMUxMTDMjRywqKiohIzokITJsOE54XjR3WCIhIzs3JCQiMkskM19ObC41NSEjOiQhMlQ/J1JoI29cXyIhIzs3JCQiMkxMMy1EZzUtIiEjOiQhMiYqW3BCblIhZjohIzs3JCQiMkwzX3Y1c2wtIiEjOiQhMi9RbzgnKXA3YyIhIzs3JCQiMk0kZSpbJ1IzSzUhIzokITIwSG9OaUdFYiIhIzs3JCQiMk1lUkEjZWZQNSEjOiQhMi1zIypvM1VEYCIhIzs3JCQiMk1MJGV6dzVWNSEjOiQhMiMzPm9CT2YrOiEjOzckJCIybm1tOy4rQjEiISM6JCExPVwqcCRHeSVIIiEjOjckJCItdiRRI1wiMyIhIzUkITE4KkdgZ3FFXSohIzs3JCQiMk0kM3g7bCY9NCIhIzokITE4OnVvaCFmOyghIzs3JCQiMm5tInpcMUEtNiEjOiQhMSM+LSVHZyFwZSUhIzs3JCQiL0QiR3klZTc2ISM3JCExSytkWCRIUyY9ISM7NyQkIjJOTExlIipbSDciISM6JCIxPlFnUyZbK08qISM8NyQkIjIsK3YkemdsTDYhIzokIjEpb25OKDNhc1AhIzs3JCQiMm9tO0hDalY5IiEjOiQiMSI0aWNFZ21YJyEjOzckJCIyTUxla1NxXToiISM6JCIxIjNkbFN0MCopKSEjOzckJCIqZHhkOyIhIigkIjJ5VVpsNFQiKjQiISM7NyQkIjIsK103SkZuPSIhIzokIjIkPlhLW3YjPVIiISM7NyQkIixEMHh3PyIhIiokIjElPj4nRypvUV4iISM6NyQkIjIsdm8vYVdGQCIhIzokIjI6YyEqPkQqcD06ISM7NyQkIjItXVAlRz8ieUAiISM6JCIydThfQSRHeDk6ISM7NyQkIjBEMWteekdBIiEjOCQiMUdbKVwvJXAtOiEjOjckJCIudlYrWnpBIiEjNiQiMkpEVishUTokWyIhIzs3JCQiL0RKISk+M1E3ISM3JCIxPHlTTEUiW1UiISM6NyQkIjIsK11pJnBAWzchIzokIjItdEpoYEhsTSIhIzs3JCQiLnY9R0IyRiIhIzYkIjJiOlwpPVRAUTYhIzs3JCQiLHZnSEtIIiEiKiQiMFFzcyhIJ3pWKiEjOjckJCIybG0ienU1TS44ISM6JCIxXCFIJ29UcG0oKSEjOzckJCIyTEwkM1VEWDg4ISM6JCIwZGNFSmRVRSkhIzo3JCQiMm47SGRGMyY9OCEjOiQiMFopKSk+ISpweiEpISM6NyQkIi52JDRTY0I4ISM2JCIxKnBAIlEoRyJSeiEjOzckJCIyTCQzLVYoPidHOCEjOiQiMVBmTjtbIzMleSEjOzckJCIybG1tbVp2T0wiISM6JCIxWi1JZyUqKj55KCEjOzckJCIyS0xlOWJ0IVI4ISM6JCIvL1cuPVlleCEjOTckJCIyKSoqKlxpaXJXTSIhIzokIjFDImVgUngocHghIzs3JCQiMm5tVDVxcClcOCEjOiQiMVo8KVF2TCE0eSEjOzckJCIyTExMZXhuX04iISM6JCIxLlYjUS9hJG95ISM7NyQkIjJvbTthI1IxbTghIzokIjFITiVvXEI9LCkhIzs3JCQiMiwrK10yZ29QIiEjOiQiMVQjeVpjbFU3KSEjOzckJCIybiJIMkw6diJRIiEjOiQiMWZmIylvd1pWIikhIzs3JCQiMkwkZTkiKkhrJ1EiISM6JCIxUiNSYmAuRDgpISM7NyQkIjB2PSNcV2AiUiIhIzgkIjFXLyEqPlQ9JTMpISM7NyQkIjJubSJIMmZVJ1IiISM6JCIxWW43WWlxIip6ISM7NyQkIi92VkIpM2lTIiEjNyQiMSdmdVVXUiVcdyEjOzckJCIyS0wkZVI8KmZUIiEjOiQiMSd6RS51VEIxKCEjOzckJCIybm0iSGlCUVA5ISM6JCIxRF1gKkh0XXolISM7NyQkIjItKytdKUh4ZTkhIzokIjJFdjMqUiE9KGY2ISM8NyQkIjJvO3olUmskKW85ISM6JCExX3prbVE6OiYqISM8NyQkIjJMTGVSKikqKil5OSEjOiQhMjEoXFtuJVtjQyQhIzw3JCQiL3ZWW0wnKilbIiEjNyQhMSIqPSUqeiRINWwmISM7NyQkIjJtbTtIIW8tKlwiISM6JCExJmUoKj1TcGgzKSEjOzckJCIyKioqKipcN2dhNDohIzokITJjRU9iIkcqbzAiISM7NyQkIjJNTCQzQV8xPzohIzokITJCNGdaMXUpKUciISM7NyQkIjJvbW07ViVlSTohIzokITJZWlFYcTVdXCIhIzs3JCQiMiwrXTdrLjZhIiEjOiQhMjtTLFs0L29tIiEjOzckJCIybzthODQpUV46ISM6JCEyY3A7IXkiKilbeiIhIzs3JCQiMk1MZTlhcztjIiEjOiQhMilRWTFTKEcpeT0hIzs3JCQiMm5UNWx3OW9jIiEjOiQhMls5YFlmU0whPiEjOzckJCIvRGMiKnAmPmQiISM3JCEweCdRVFErOz4hIzk3JCQiMkxlOW1AKjR4OiEjOiQhMnY/NT8lWyRvIj4hIzs3JCQiMm1tbTtXVEFlIiEjOiQhMj1lQyNlJSkpZiE+ISM7NyQkIjFMJGVSPHZQZyIhIzkkITI8NCczcEheVTwhIzs3JCQiLUQxKjNgaSIhIzUkITI4IilSKUdFO0U5ISM7NyQkIjFMMy15J3ljaiIhIzkkITIyRTcncGxdUjchIzs3JCQiMm5tInpcJVtnayIhIzokITJlWHB4SCVmVTUhIzs3JCQiMi5daTpBPWtsIiEjOiQhMSJHXjk/TkxXKSEjOzckJCIyTkxMTCp6eW07ISM6JCExIls3J0g6OUlsISM7NyQkIjJPTCQzc3Iqem8iISM6JCEyJmVTS1kzYD1KISM8NyQkIjJPTEwzTjEjNDwhIzokIS9DYyRSaF8hbyEjOjckJCIyLSt2byEqUi10IiEjOiQiMTA3TkQwP1RyISM8NyQkIjFubSJIWXQ3diIhIzkkIjInR0QxdiJHI1s4ISM8NyQkIjJQTGVrNiwxeCIhIzokIjFjL19CJEc3biIhIzs3JCQiMi0rKyt4RyoqeSIhIzokIjFUYCFRKkh3KD0jISM7NyQkIjJRTCQzVS8zNz0hIzokIjFMXk0lNEBRWyQhIzs3JCQiMW5tOzlAQk09ISM5JCIxS3onKW9UbCN6JiEjOzckJCIyLitdUCRbL2E9ISM6JCIxNDYhXDVrcWopISM7NyQkIjJPTExMYmRRKD0hIzokIjFLInpOZyopcD0iISM6NyQkIjJvbSJ6Vz8pXCo9ISM6JCIyRj5rOCRSQjo6ISM7NyQkIi1ET2w1Oz4hIzUkIjJaQ3clSFExYzwhIzs3JCQiL3ZvVzc7QD4hIzckIjI6SCE0JjRhTHoiISM7NyQkIi5ESiZmQEU+ISM2JCIxMUopcE9WNyM9ISM6NyQkIjIpXFBNMkx1Rz4hIzokIjInM1JoSFhaSj0hIzs3JCQiL0RjaDFGSj4hIzckIjI9Wi5tIUc5Uj0hIzs3JCQiMi5EInk6ISl6TD4hIzokIjFfL0VCOT5XPSEjOjckJCIqUERqJD4hIigkIjJGRyxVIjNkWT0hIzs3JCQiMi1dUCV5K1FUPiEjOiQiMjtPR0QkNDtWPSEjOzckJCIyLyt2b3lNayU+ISM6JCIybS1rSTVvJ0c9ISM7NyQkIi9ESiZcKltePiEjNyQiMio9MEx4SihIIT0hIzs3JCQiMi0rXVA/V2wmPiEjOiQiMmE7KUc3YzNtPCEjOzckJCIvRCJHOjN1Jz4hIzckIjJcYG06UVooXDshIzs3JCQiMi4rdj01cyN5PiEjOiQiMmRseHdWJCllWyIhIzs3JCQiL3YkNDBPIiopPiEjNyQiMmlZWmpJKFEhRyIhIzs3JCQiJCsjISIiJCIyJnpnTGMzTVQ1ISM7LSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSdDVVJWRVNHNiQ3ZXQ3JCQiIiEhIiIkIiNJISIiNyQkIjJrbW07YSlHXGEhIz0kIjJNLjw2d3clW0khIzs3JCQiMkxMTCQzeCYpKjMiISM8JCIxVz5gNUxrJTMkISM6NyQkIixEYyd5TTshIzYkIjJjdTRSYGEiM0ohIzs3JCQiMm1tbW1UOih6QCEjPCQiMWhLKVw2bCg9SiEjOjckJCIxTEwkM0ZXWXMjISM7JCIxV0VPJ1FMajYkISM6NyQkIitESmRwSyEjNSQiMjkuSityQjM1JCEjOzckJCIxbm07ej5dOVEhIzskIjBYImZYdElzSSEjOTckJCIxTExMTDNWZlYhIzskIjJPOEIsemg0LiQhIzs3JCQiLERjKilmRCchIzYkIjFYUTojXEA+eiMhIzo3JCQiMW5tbSJIW0Q6KSEjOyQiMk8jKm9qQnFFVSMhIzs3JCQiLXYkcFUmRzUhIzYkIjJHR01zQ3E6Kj0hIzs3JCQiMkxMTCRlMCQ9QyIhIzskIjFmSmc1OG0qRyIhIzo3JCQiMk1MTGVSIj1cOCEjOyQiMV9CaFxyL3AoKiEjOzckJCIyTkxMTEJLbFgiISM7JCIxR2lbWE45ZG0hIzs3JCQiMk9MTDMyJClRYyIhIzskIjJ2ZyVIJVFpOWgkISM8NyQkIjJPTEwkM1JCcjshIzskIjE6MyI0WCZ6Pm4hIzw3JCQiMi0rXVA5OlwpPSEjOyQhMXNKbE9hKjN3JSEjOzckJCIybW1tInpqZik0IyEjOyQhMUA8Q3VeZXQmKiEjOzckJCIsdkIxbkgjISM1JCEydD54OClRVV44ISM7NyQkIjJNTExlNDtbXCMhIzskITIpUkQlKXBpJVFxIiEjOzckJCIxbm07enQlKipwIyEjOiQhMV9DS2cvdEs/ISM6NyQkIixEbXldIUghIzUkITIwMz8pKnknM0VCISM7NyQkIjFubTtIZEE8SiEjOiQhMilvMkZKNyE9ZSMhIzs3JCQiMlBMTGV6cyRITCEjOyQhMWptcz05Pm1GISM6NyQkIjFubSJIWlJBUSQhIzokITJZOipbIT1cdXojISM7NyQkIio6MV5WJCEiKSQhMkwlKno4PSZ6QEchIzs3JCQiMk1MJDNGRyh6WyQhIzskITImXGFkXk5xUUchIzs3JCQiMm5tbVRdUjNhJCEjOyQhMiYqcCQ+OzxtWkchIzs3JCQiLUQiPTFQZiQhIzYkITElKm8qKT4neiJbRyEjOjckJCIyS0xMJGVHZFlPISM7JCEyWShbWSQpKSl6UkchIzs3JCQiMmxtO2FgUiUqcCQhIzskITE6JlxqdCxAI0chIzo3JCQiLERAMUJ2JCEjNSQhMWBfNnMmPlp6IyEjOjckJCIxTUxlOWQjKXBSISM6JCExWiJlKnl0ZnVEISM6NyQkIjFubW07X00oPSUhIzokITAuKWZWVDN5QCEjOTckJCIxLCtdN2wkKnlWISM6JCEyTD9oIT1KMCVwIiEjOzckJCIxTUxMM3lfcVghIzokITJHLlNmSCo9NTYhIzs3JCQiMSwrdj01UHlZISM6JCExX14iNG4xU14oISM7NyQkIjFubTtIVUAneSUhIzokITFSaFxoUjwxUSEjOzckJCIxTExlUnUwJSpbISM6JCEwKW8mb0w7LDEnISM8NyQkIjIwKysrbCs+KyYhIzskIjFrImUlKW87MGskISM7NyQkIjIwKytdbic9NV4hIzskIjE7JypvRzUoSEIoISM7NyQkIilGWj1fISIoJCIycFdIL2kuSDEiISM7NyQkIitEKGVuSyYhIiokIjJpVShbRlRrdzghIzs3JCQiKnZXXVYmISIpJCIySyVSeDRHYGY7ISM7NyQkIiowX1BrJiEiKSQiMXJXaCpvKXAwQCEjOjckJCIqTmZDJmUhIikkIjEpKkg0Ink+V1QjISM6NyQkIjFNJGU5WUFzJWYhIzokIjImekcuKy5mNUQhIzs3JCQiMW5tIkhkJik+LychIzokIjElKSo0Jm8/TCJlIyEjOjckJCIudlZvW244JyEjNyQiMSQpSC8lemInR0UhIzo3JCQiMU1MJGV6NjpCJyEjOiQiMlldOlgqZnphRSEjOzckJCIvdlZbSSZ5RychIzgkIjEnKj44ZSZSN20jISM6NyQkIjFuOy8sVj5XaiEjOiQiMjF3RXQmKXo5bSMhIzs3JCQiMU1la2BiYCtrISM6JCIxIz04VFlsZmwjISM6NyQkIi1EMW8ob1gnISM2JCIyYT9KPS43XmsjISM7NyQkIjFNJGU5SmYmcGwhIzokIjInXDFEeVIkKTNFISM7NyQkIjFubW07PUMjbychIzokIjJ3ZzxXNjJdYiMhIzs3JCQiMW9tbVRiOnRvISM6JCIyJ3kjZm11NnNVIyEjOzckJCIxbm1tbSNwUzEoISM6JCIxI3olM0o8JEdEIyEjOjckJCIxTExlOXQ5J0coISM6JCIxOFAobyVHeSIpPiEjOjckJCIsRE9EIzN2ISM1JCIybSVwMW5QYz87ISM7NyQkIjFMTGVrPyFbcSghIzokIjJXcDlsJHpCODchIzs3JCQiMW1tbW0oeTgheiEjOiQiMUUsdXlfJXA/KCEjOzckJCIuRGNtOyM0ISkhIzckIjF1eXptVGFyVCEjOzckJCIxTUxla1gwPCIpISM6JCIxKT1FdTNQNFIqISM8NyQkIjFvO2FqQypbQSkhIzokITJ1KkglKWZQVVdDISM8NyQkIjEsK11pLnRLJCkhIzokITF5SCEpZiNbOiNmISM7NyQkIi12VnZUTiUpISM2JCExQiFcWUEpUWUjKiEjOzckJCIrRFo1USYpISIqJCEyKHpuRVlMY2I3ISM7NyQkIjEsK0QxPnpTJykhIzokITFmOURCJUhZZCIhIzo3JCQiLHYzek11KSEjNSQhMnMvMVZtKUd3PSEjOzckJCIxTkxlKm9seCYqKSEjOiQhMS5XdiJbazxVIyEjOjckJCIxb21tIkhfPzwqISM6JCEyJlEjUj1JNE8iRyEjOzckJCIxbjthKUdWL0YqISM6JCEyOWBiY1dfRCRIISM7NyQkIjFvbVQmR00pbyQqISM6JCEyWS0mPTshMzIsJCEjOzckJCIxb1QmUXlIIT0lKiEjOiQhMlhaXiVwR05NSSEjOzckJCIxbztII0dEc1kqISM6JCEydCFlQHdReFpJISM7NyQkIjFuIkgyeT9rXiohIzokITImNHVRb0s1XkkhIzs3JCQiMW5tO3ppaGwmKiEjOiQhMm4hXE8scmBXSSEjOzckJCItdlZBKHl4KiEjNiQhMVxzbGJYLTVIISM6NyQkIjFMTEwzI0csKioqISM6JCEyVmA3PF4jUUdFISM7NyQkIjJMTDMtRGc1LSIhIzokITEmcFYtZ2FbQSMhIzo3JCQiMk1MJGV6dzVWNSEjOiQhMkMrLzotO2V2IiEjOzckJCIybm1tOy4rQjEiISM6JCEyYiZvVU9GP0Q4ISM7NyQkIi12JFEjXCIzIiEjNSQhMT5KV1YiXF0iKikhIzs3JCQiMm5tInpcMUEtNiEjOiQhMUBiJmY9cT1GJSEjOzckJCIyTkxMZSIqW0g3IiEjOiQiMiZ5MTxNUjgpWyQhIz03JCQiMm9tO0hDalY5IiEjOiQiMU8iPXhpZ1g8JiEjOzckJCIqZHhkOyIhIigkIjJMPDYhRydSNCwiISM7NyQkIjIsK103SkZuPSIhIzokIjJvWWA/TUEvXSIhIzs3JCQiLEQweHc/IiEiKiQiMSIqNHR2VT16PiEjOjckJCIudlYrWnpBIiEjNiQiMSZcIj5EcDcuQyEjOjckJCIyLCtdaSZwQFs3ISM6JCIxVic9bD1mN3YjISM6NyQkIi9EMT4sWmY3ISM3JCIyWU8sOC92cypHISM7NyQkIi52PUdCMkYiISM2JCIyJ0hHPmhUXyxJISM7NyQkIjBEIkdqKVxqRiIhIzgkIjFBSi0kUXlpLiQhIzo3JCQiL3ZvV2soPkciISM3JCIxUXlMbyUqb2VJISM6NyQkIjIsdiQ0RUlnKEciISM6JCIyQTJwLUwnSG9JISM7NyQkIix2Z0hLSCIhIiokIjIvc15Pa1BaMSQhIzs3JCQiMkxMJDNVRFg4OCEjOiQiMXhhPTkpeS8lSCEjOjckJCIybG1tbVp2T0wiISM6JCIyeG42IlIlZlxrIyEjOzckJCIyKSoqKlxpaXJXTSIhIzokIjJ1NVJAbGlLVSMhIzs3JCQiMkxMTGV4bl9OIiEjOiQiMngxMDMnW0dqQCEjOzckJCIyb207YSNSMW04ISM6JCIwWD5XQyJIcj0hIzk3JCQiMiwrK10yZ29QIiEjOiQiMWZeJkhPJVJhOiEjOjckJCIyTCRlOSIqSGsnUSIhIzokIjEnZmwjejc1XzchIzo3JCQiMm5tIkgyZlUnUiIhIzokIjFBYCwnKnAoUlQqISM7NyQkIi92VkIpM2lTIiEjNyQiMWEqKjMqKilRJXppISM7NyQkIjJLTCRlUjwqZlQiISM6JCIyT00+KzAicHFKISM8NyQkIi92JDQwKG9FOSEjNyQhMlAuail6XWg9OSEjPTckJCIybm0iSGlCUVA5ISM6JCExJD0wKj1jWTdMISM7NyQkIjJOJGVrdHcyWzkhIzokITEtIykqR2gvJylIJyEjOzckJCIyLSsrXSlIeGU5ISM6JCExIz1WXilbKClwISohIzs3JCQiMkxMZVIqKSoqKXk5ISM6JCEydFNrQWdTWk8iISM7NyQkIjJtbTtIIW8tKlwiISM6JCEyWSFlOCIpZjNSPCEjOzckJCIyTUwkM0FfMT86ISM6JCExKVE1JWUwUV0/ISM6NyQkIjIsK103ay42YSIhIzokITIkZlIjNEgsU0gjISM7NyQkIjJNTGU5YXM7YyIhIzokITJRcFdaaiVweEMhIzs3JCQiMm1tbTtXVEFlIiEjOiQhMHh6QF1lKDRFISM5NyQkIjIpKlw3eUkzSWYiISM6JCEyREJVcyUqKXliRSEjOzckJCIxTCRlUjx2UGciISM5JCEyOyRvXkdQLCRvIyEjOzckJCIwREpxZ2UiNDshIzgkITJfK2UsT1cnKW8jISM7NyQkIjJsOy8sL1VYaCIhIzokITFwIj5lQislKW8jISM6NyQkIjJLM3hKWkQqPjshIzokITIkKSpccTNRJT1vIyEjOzckJCItRDEqM2BpIiEjNSQhMkRyK0U8SSZvRSEjOzckJCIybm0ielwlW2drIiEjOiQhMm5YMDt5LGlhIyEjOzckJCIyTkxMTCp6eW07ISM6JCExIypcKjNqNGFIIyEjOjckJCIyT0wkM3NyKnpvIiEjOiQhMXgqelVPeVUqPSEjOjckJCIyT0xMM04xIzQ8ISM6JCExOW9BdXhlXzghIzo3JCQiMW5UJilHSnM+PCEjOSQhMVU2VTF0RlI1ISM6NyQkIjItK3ZvISpSLXQiISM6JCExI0hnKikqKmZ3LSghIzs3JCQiMk8kZSpbb2MydSIhIzokITIkXF4kW05vPlwkISM8NyQkIjFubSJIWXQ3diIhIzkkIjJUNyJIPHFhWzkhIz03JCQiL3ZvKkdQNHciISM3JCIyOTEhSExtJj1eJCEjPDckJCIyUExlazYsMXgiISM6JCIxKHBIJypcZzolbyEjOzckJCIxbiJISyVcRSF5IiEjOSQiMndBa2UyI1wyNSEjOzckJCIyLSsrK3hHKip5IiEjOiQiMGZlVE5nYkoiISM5NyQkIjFuOy8xWSssPSEjOSQiMmJuQkZvTEtrIiEjOzckJCIyUUwkM1UvMzc9ISM6JCIxM2pqKEh4dCQ+ISM6NyQkIjIwK0QieWk6Qj0hIzokIjIlKSlcZjMib0U+IyEjOzckJCIxbm07OUBCTT0hIzkkIjFmempiWE8wQyEjOjckJCIyLitdUCRbL2E9ISM6JCIyTXUmM0AkNFRuIyEjOzckJCIyT0xMTGJkUSg9ISM6JCIxaVlYWiczPiFHISM6NyQkIjE8ej53J1Eiej0hIzkkIjJOQFEncCZHUSJHISM7NyQkIjItXWkhKno+Vyk9ISM6JCIyLjcnPnQkZXIiRyEjOzckJCIyTjNGPiM0cSopPSEjOiQiMWAqMz1uW0IiRyEjOjckJCIyb20ielc/KVwqPSEjOiQiMi9uWFdfeCkqeiMhIzs3JCQiMk0kM18hSFdiIT4hIzokIjI6JFxRbmQlUnYjISM7NyQkIi1ET2w1Oz4hIzUkIjIuR2YrKT5SJG8jISM7NyQkIipQRGokPiEiKCQiMSJcJT1SQWAkXCMhIzo3JCQiMi0rXVA/V2wmPiEjOiQiMUhpaWZRJypcQSEjOjckJCIyLit2PTVzI3k+ISM6JCIxVGRZSHdAVj4hIzo3JCQiJCsjISIiJCIyRmJDaXNEQmYiISM7LSUmQ09MT1JHNiYlJFJHQkckIjInXF8nZU0leWc+ISM8JCIxaWo5LmVAUiEpISM7JCIyJ1xfJ2VNJXlnPiEjPC0lJVZJRVdHNiQ7JCIiISEiIiQiJCsjISIiJShERUZBVUxURy0lK0FYRVNMQUJFTFNHNictSSNtaUc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkc2IjY1USJ0NiIvJSdmYW1pbHlHUSE2Ii8lJXNpemVHUSMxMDYiLyUlYm9sZEdRJmZhbHNlNiIvJSdpdGFsaWNHUSV0cnVlNiIvJSp1bmRlcmxpbmVHUSZmYWxzZTYiLyUqc3Vic2NyaXB0R1EmZmFsc2U2Ii8lLHN1cGVyc2NyaXB0R1EmZmFsc2U2Ii8lK2ZvcmVncm91bmRHUShbMCwwLDBdNiIvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XTYiLyUnb3BhcXVlR1EmZmFsc2U2Ii8lK2V4ZWN1dGFibGVHUSZmYWxzZTYiLyUpcmVhZG9ubHlHUSZmYWxzZTYiLyUpY29tcG9zZWRHUSZmYWxzZTYiLyUqY29udmVydGVkR1EmZmFsc2U2Ii8lK2ltc2VsZWN0ZWRHUSZmYWxzZTYiLyUscGxhY2Vob2xkZXJHUSZmYWxzZTYiLyU2c2VsZWN0aW9uLXBsYWNlaG9sZGVyR1EmZmFsc2U2Ii8lLG1hdGh2YXJpYW50R1EnaXRhbGljNiJRITYiLSUlRk9OVEc2JSUoREVGQVVMVEclKERFRkFVTFRHIiM1JStIT1JJWk9OVEFMRyUrSE9SSVpPTlRBTEctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiRJJSEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJWdNISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJV1RISIiLSUpQ0hJTERSRU5HNiI=

JSFH TTdSMApJNVJUQUJMRV9TQVZFLzUxOTkyMDkyWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjJSIjIiMhIiYiIiNGKCEiI0YmTTdSMApJNVJUQUJMRV9TQVZFLzY1MzgyOTU2WColKmFsZ2VicmFpY0c2IjYiW2dsISMlISEhIiMiIyEiIiEiJ0YmTTdSMApJNVJUQUJMRV9TQVZFLzY1MzgxMzE2WCwlKmFsZ2VicmFpY0c2IjYiW2dsISIlISEhIyUiIyIjIyIiIiIiI0YoISIjRihGJg==TTdSMApJNVJUQUJMRV9TQVZFLzY1MzgxMzE2WCwlKmFsZ2VicmFpY0c2IjYiW2dsISIlISEhIyUiIyIjIyIiIiIiI0YoISIjRihGJg==TTdSMApJNVJUQUJMRV9TQVZFLzY2NDA2MTI4WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjJSIjIiMhIiIiIiFGKCEiJ0YmTTdSMApJNVJUQUJMRV9TQVZFLzY3ODQ0OTUyWCwlKmFsZ2VicmFpY0c2IjYiW2dsISIlISEhIyUiIyIjIyIiIyIiJiMhIiNGKSMiIiVGKSMiIiJGKUYmTTdSMApJNVJUQUJMRV9TQVZFLzY4MzI4ODM2WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjIyIjIiItJSN1MUc2IyUidEctJSN1MkdGKUYmTTdSMApJNVJUQUJMRV9TQVZFLzQ5MDczMjEyWCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjIyIjIiIsJi0lI3UxRzYjJSJ0RyMiIiMiIiYtJSN1MkdGKiMiIiVGLiwmRigjISIjRi5GLyMiIiJGLkYmTTdSMApJNVJUQUJMRV9TQVZFLzY2MDA5NDg0WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjIyIjIiIjIiM5IiImIyIiIkYpRiY=TTdSMApJNVJUQUJMRV9TQVZFLzY1NDkxMjQ4WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjIyIjIiIjIiIkIiImIyIiI0YpRiY=TTdSMApJNVJUQUJMRV9TQVZFLzYyOTc3OTQ4WCwlKWFueXRoaW5nRzYiNiJbZ2whIiUhISEjIyIjIiIsKiomLSUkc2luRzYjKiYiIicjIiIiIiIjJSJ0R0YvRi9GLUYuIyEiIyIjOi0lJGNvc0dGKyNGMyIiJi1GKjYjRjEjIiIkIiM1LUY2RjojIiIoRjgsKkYoI0YvRjRGNSNGL0Y4RjkjRjxGOEY+IyIjOUY4RiY=