Resonance (Example 9.4.2)
unit:= (t,a,b) -> Heaviside(t-a) - Heaviside(t-b);
Zio2JUkidEc2IkkiYUdGJUkiYkdGJUYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCYtSSpIZWF2aXNpZGVHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IywmOSQiIiI5JSEiIkY0LUYtNiMsJkYzRjQ5JkY2RjZGJUYlRiU=
Square wave function
F1 := t-> piecewise(0<t and t<=Pi, 10, Pi<t and t<=2*Pi,-10);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkqcGllY2V3aXNlRyUqcHJvdGVjdGVkRzYmMzIiIiE5JDFGMEkjUGlHRisiIzUzMkYyRjAxRjAsJEYyIiIjISM1RiVGJUYl
use Heaviside functions to get a few periods (at least for t>0, which is what we care for)
F1per := t-> sum(Heaviside(t-2*Pi*n)*F1(t-2*Pi*n),n=0..5);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc3VtRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQqJi1JKkhlYXZpc2lkZUdGKzYjLCY5JCIiIiomSSNQaUdGLEY1SSJuR0YlRjUhIiNGNS1JI0YxR0YlRjJGNS9GODsiIiEiIiZGJUYlRiU=
plot(F1per(t),t=0..6*Pi);
NictJSdDVVJWRVNHNiQ3XnI3JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMkBfcTZPZVJHIiEjPSQiJCsiISIiNyQkIjJVL1RCczt6YyMhIz0kIiQrIiEiIjckJCIxbTpeJDN2PSZRISM8JCIkKyIhIiI3JCQiMiUpMyNvV00kZTgmISM9JCIkKyIhIiI3JCQiMUtKLW4sdi54ISM8JCIkKyIhIiI3JCQiMnhUTyopb21yLSIhIzwkIiQrIiEiIjckJCIybGkvTS5dMmEiISM8JCIkKyIhIiI3JCQiMmEkRyh5UExWMCMhIzwkIiQrIiEiIjckJCIxYCM0bzErOjMkISM7JCIkKyIhIiI3JCQiMXJjdWJubTNUISM7JCIkKyIhIiI3JCQiMVc+UFY3OCcqZSEjOyQiJCsiISIiNyQkIjE7IykqNHQmZiRvKCEjOyQiJCsiISIiNyQkIjJgJT5RRHhScTYhIzskIiQrIiEiIjckJCIyKikpKTR2JTM1djohIzskIiQrIiEiIjckJCIyUWhValshKXkoPiEjOyQiJCsiISIiNyQkIjI8PUI0eDM4TiMhIzskIiQrIiEiIjckJCIyMidmekhAKHp0IyEjOyQiJCsiISIiNyQkIjIwIUhCSGcieiRIISM7JCIkKyIhIiI3JCQiMi4lKXAnRypmeTgkISM7JCIkKyIhIiI3JCQiMCMzKGZNOy46JCEjOSQhJCsiISIiNyQkIi89RmpGeGlKISM4JCEkKyIhIiI3JCQiMih6RmQhPUhfPCQhIzskISQrIiEiIjckJCIyJmZQKHlmJm8oPSQhIzskISQrIiEiIjckJCIxPmRaSyUpZjdLISM6JCEkKyIhIiI3JCQiMih5dzJuN15QSyEjOyQhJCsiISIiNyQkIjEpZiJHT3BMKEckISM6JCEkKyIhIiI3JCQiMnZeJlswRTtQTCEjOyQhJCsiISIiNyQkIjJqTiQqUSVSIm9WJCEjOyQhJCsiISIiNyQkIjJZPixCR2xrYCQhIzskISQrIiEiIjckJCIwND5HZXM5dSQhIzkkISQrIiEiIjckJCIyWClwTCQpKXprJVIhIzskISQrIiEiIjckJCIxM1w1bjRpMlYhIzokISQrIiEiIjckJCIxUFhSIkghPTlaISM6JCEkKyIhIiI3JCQiMTo2TCgpKjNDNyYhIzokISQrIiEiIjckJCIxL1IhKXlIImVeJiEjOiQhJCsiISIiNyQkIjE5MEtOc1YlcCYhIzokISQrIiEiIjckJCIxQnIkPVxoSShlISM6JCEkKyIhIiI3JCQiMiZSLXpdQUV6ZiEjOyQhJCsiISIiNyQkIjFjTHU0SVkmMychIzokISQrIiEiIjckJCIyWCIqPiMqUWomUWghIzskISQrIiEiIjckJCIxdGtwb1BtIj4nISM6JCEkKyIhIiI3JCQiMV9aVmVSQD1pISM6JCEkKyIhIiI3JCQiMUpJPFtUd1dpISM6JCEkKyIhIiI3JCQiMjA8VUlDUiFlaSEjOyQhJCsiISIiNyQkIjBKNnpMOThGJyEjOSQhJCsiISIiNyQkIjBYIXlLJSplJUcnISM5JCIkKyIhIiI3JCQiMiYqZVx3X2t5SCchIzskIiQrIiEiIjckJCIxVCJwU1onenhrISM6JCIkKyIhIiI3JCQiMjpwKVs/JUd4bCchIzskIiQrIiEiIjckJCIxYXonem5OajIoISM6JCIkKyIhIiI3JCQiMTRzNUpTKG9XKCEjOiQiJCsiISIiNyQkIjEuMm5KTFRgeSEjOiQiJCsiISIiNyQkIjE3cihSKVtgUyMpISM6JCIkKyIhIiI3JCQiMSc+IjQ3Y1hXJykhIzokIiQrIiEiIjckJCIxJnkvNCQ0UTohKiEjOiQiJCsiISIiNyQkIjE7K0BrdVU6IyohIzokIiQrIiEiIjckJCIxW19eKCpSWjolKiEjOiQiJCsiISIiNyQkIjElUSRSMDVZRyUqISM6JCEkKyIhIiI3JCQiMUE6RjghWzlXKiEjOiQhJCsiISIiNyQkIjBtXDYtTldYKiEjOSQhJCsiISIiNyQkIjEpekYhSD9VbiUqISM6JCEkKyIhIiI3JCQiMXNTeVdnUiRcKiEjOiQhJCsiISIiNyQkIjFZLmFnK1A+JiohIzokISQrIiEiIjckJCIxKCpHMCM0PThkKiEjOiQhJCsiISIiNyQkIjFaYWNCaEVCJyohIzokISQrIiEiIjckJCIxWjBmJz1pcnMqISM6JCEkKyIhIiI3JCQiMVtjaFwjZTUkKSohIzokISQrIiEiIjckJCIyJ3klPTE+JT4sNSEjOiQhJCsiISIiNyQkIjJDUnZpYiNHPjUhIzokISQrIiEiIjckJCIyb0UpPTxXTmU1ISM6JCEkKyIhIiI3JCQiMnJnJnBzJT4oKTQiISM6JCEkKyIhIiI3JCQiMiM+LGEvKTMjUTYhIzokISQrIiEiIjckJCIybkguT3c7azwiISM6JCEkKyIhIiI3JCQiMjd5TGRaR3c+IiEjOiQhJCsiISIiNyQkIjJjRWt5PVMpPTchIzokISQrIiEiIjckJCIyMFJgNy1xJEc3ISM6JCEkKyIhIiI3JCQiMmBeVVkmKSoqeUIiISM6JCEkKyIhIiI3JCQiMngyUDh4a0VDIiEjOiQhJCsiISIiNyQkIjIta0ohKW9IdUMiISM6JCEkKyIhIiI3JCQiMjkjKnlqOTcpXDchIzokISQrIiEiIjckJCIyRj9FWmclPl83ISM6JCEkKyIhIiI3JCQiMkwlKSoqUSRlUWA3ISM6JCEkKyIhIiI3JCQiMSVbdEkxeFhEIiEjOSQhJCsiISIiNyQkIjJXN1pBSG9kRCIhIzokISQrIiEiIjckJCIxbDJVQCZmcEQiISM5JCIkKyIhIiI3JCQiMlAnNCJSeTR0RiIhIzokIiQrIiEiIjckJCIyQjssay9nd0giISM6JCIkKyIhIiI3JCQiMjMlM0tyMmFNOCEjOiQiJCsiISIiNyQkIjFjUzhYNidbUCIhIzkkIiQrIiEiIjckJCIyRzNfWVoqejc5ISM6JCIkKyIhIiI3JCQiMkQ9Si8nZlhfOSEjOiQiJCsiISIiNyQkIjJuImYjUVVGN1wiISM6JCIkKyIhIiI3JCQiMkQmcEwpW0E6XiIhIzokIiQrIiEiIjckJCIxKSl6JUdiPD1gIiEjOSQiJCsiISIiNyQkIjIob0daSzNmVDohIzokIiQrIiEiIjckJCIyJFx4NDdUT146ISM6JCIkKyIhIiI3JCQiMicqPTU+dl1pYiIhIzokIiQrIiEiIjckJCIyKkhFcyJSUDZjIiEjOiQiJCsiISIiNyQkIjAmKUc7ciFlajohIzgkIiQrIiEiIjckJCIyLTJOOi5DZ2MiISM6JCIkKyIhIiI3JCQiMi49KVsicFhzYyIhIzokIiQrIiEiIjckJCIyLkhUOU5uJW86ISM6JCIkKyIhIiI3JCQiMi9TJVI2ISpvcDohIzokIiQrIiEiIjckJCIyMF5aOG41NGQiISM6JCEkKyIhIiI3JCQiMSlcXzgmKSoqM2YiISM5JCEkKyIhIiI3JCQiMmFbZDguKikzaCIhIzokISQrIiEiIjckJCIyI1I+T11qYF07ISM6JCEkKyIhIiI3JCQiMidwZCIqPSJvcG8iISM6JCEkKyIhIiI3JCQiMjlWLHEuQihHPCEjOiQhJCsiISIiNyQkIjImUjddTCJwZ3ciISM6JCEkKyIhIiI3JCQiMUF2c2h5KWUhPSEjOSQhJCsiISIiNyQkIjE4YFJsJSoqUiU9ISM5JCEkKyIhIiI3JCQiLVIpZWJcKT0hIzUkISQrIiEiIi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIrI2ZiXCk9ISIpJShERUZBVUxURy0lKkFYRVNTVFlMRUc2IyUnTk9STUFMRy0lKFNDQUxJTkdHNiMlLlVOQ09OU1RSQUlORURHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkIVIhIiItJSlCT1VORFNfWUc2IyQiJD8iISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlZ04hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlXVAhIiItJSlDSElMRFJFTkc2Ig==
b1 := n-> (1/Pi)*int(F1(t)*sin(n*t),t=0..2*Pi);
Zio2I0kibkc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiZJI1BpRyUqcHJvdGVjdGVkRyEiIi1JJGludEc2JEYrSShfc3lzbGliR0YlNiQqJi1JI0YxR0YlNiNJInRHRiUiIiItSSRzaW5HRi82IyomOSRGN0Y2RjdGNy9GNjsiIiEsJEYqIiIjRjdGJUYlRiU=
b1(n) assuming integer;
LCQqKEkjUGlHJSpwcm90ZWN0ZWRHISIiLCZGJiIiIilGJkkibkc2IkYoRihGKkYmISM/
Sawtooth wave function
F2 := t-> 10*t * unit(t,-Pi,Pi);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLCQqJjkkIiIiLUkldW5pdEdGJTYlRissJEkjUGlHJSpwcm90ZWN0ZWRHISIiRjFGLCIjNUYlRiVGJQ==
Use Heaviside functions to get a few periods
F2per := t-> sum(Heaviside(t-(2*n+1)*Pi)*F2(t-(2*n+2)*Pi),n=-1..5);
Zio2I0kidEc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlLUkkc3VtRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiQqJi1JKkhlYXZpc2lkZUdGKzYjLCY5JCIiIiomLCZJIm5HRiUiIiNGNUY1RjVJI1BpR0YsRjUhIiJGNS1JI0YyR0YlNiMsJkY0RjUqJiwmRjhGOUY5RjVGNUY6RjVGO0Y1L0Y4O0Y7IiImRiVGJUYl
plot(F2per(t),t=0..6*Pi);
NictJSdDVVJWRVNHNiQ3XnA3JCQiIiEhIiIkIiIhISIiNyQkIjFyY3Vibm0zVCEjOyQiMXJjdWJubTNUISM6NyQkIjE7IykqNHQmZiRvKCEjOyQiMTsjKSo0dCZmJG8oISM6NyQkIjJgJT5RRHhScTYhIzskIjJgJT5RRHhScTYhIzo3JCQiMiopKSk0diUzNXY6ISM7JCIyKikpKTR2JTM1djohIzo3JCQiMlFoVWpbISl5KD4hIzskIjE5RU0nWyEpeSg+ISM5NyQkIjI8PUI0eDM4TiMhIzskIjEjPUI0eDM4TiMhIzk3JCQiMjInZnpIQCh6dCMhIzskIjIwJ2Z6SEAoenQjISM6NyQkIjIwIUhCSGcieiRIISM7JCIyLyFIQkhnInokSCEjOjckJCIyLiUpcCdHKmZ5OCQhIzskIjIuJSlwJ0cqZnk4JCEjOjckJCIwIzMoZk07LjokISM5JCExbSgzNycqb0c4JCEjOTckJCIvPUZqRnhpSiEjOCQhMmp5MlJhNy83JCEjOjckJCIyKHpGZCE9SF88JCEjOyQhMmwhb2dFaCZ6NSQhIzo3JCQiMiZmUCh5ZiZvKD0kISM7JCEybiNlSTQoKlwmNCQhIzo3JCQiMT5kWkslKWY3SyEjOiQhMW5RcXVvZXFJISM5NyQkIjIoeXcybjdeUEshIzskITJ3IT41U1NuWEkhIzo3JCQiMSlmIkdPcEwoRyQhIzokITIlKSl6KjNQW2UqSCEjOjckJCIydl4mWzBFO1BMISM7JCExcFNwLEYtWUghIzk3JCQiMmpOJCpRJVIib1YkISM7JCEwQidHajhQWUchIzg3JCQiMlk+LEJHbGtgJCEjOyQhMj1SeVstP251IyEjOjckJCIwND5HZXM5dSQhIzkkITJtXGdWczc8YSMhIzo3JCQiMlgpcEwkKSl6ayVSISM7JCEyPWdVUVUwbkwjISM6NyQkIjEzXDVuNGkyViEjOiQhMXlZMlNWY3Y+ISM5NyQkIjFQWFIiSCE9OVohIzokITImXF15Ol0rcDohIzo3JCQiMTo2TCgpKjNDNyYhIzokITI6WlspPmp4ZzYhIzo3JCQiMS9SISl5SCJlXiYhIzokITI5I292JEdCUG4oISM7NyQkIjFCciQ9XGhJKGUhIzokITFMWVVgIlE3NSUhIzo3JCQiMiYqZVx3X2t5SCchIzskIjFELitaPyN6WSIhIzs3JCQiMjpwKVs/JUd4bCchIzskIjFgNTRMNlZYUCEjOjckJCIxYXonem5OajIoISM6JCIxeU8peXEuOiR6ISM6NyQkIjE0czVKUyhvVyghIzokIjJCaUZScylvajYhIzo3JCQiMS4ybkpMVGB5ISM6JCIycDYiXEMhRy1kIiEjOjckJCIxN3IoUilbYFMjKSEjOiQiMUV2encmXHQmPiEjOTckJCIxJz4iNDdjWFcnKSEjOiQiMGg2XElxN08jISM4NyQkIjEmeS80JDRROiEqISM6JCIyJik+RFBpJj5LRiEjOjckJCIxOytAa3VVOiMqISM6JCIwVklxOlVBJEghIzg3JCQiMVtfXigqUlo6JSohIzokIjI6bU4ucClHS0ohIzo3JCQiMSVRJFIwNVlHJSohIzokITImKXlsKjMnNHo4JCEjOjckJCIxQTpGOCFbOVcqISM6JCEyMGwoMyxFI1w3JCEjOjckJCIwbVw2LU5XWCohIzkkITJEXjRLZk4+NiQhIzo3JCQiMSl6RiFIP1VuJSohIzokITJYUEpgZVsqKTQkISM6NyQkIjFzU3lXZ1IkXCohIzokITIuNXYmcFgoSDIkISM6NyQkIjFZLmFnK1A+JiohIzokITFFKT1RYitxLyQhIzk3JCQiMSgqRzAjND04ZCohIzokITF3aUlBRDAmKkghIzk3JCQiMVphY0JoRUInKiEjOiQhMmNzJHohXC9KJUghIzo3JCQiMVowZic9aXJzKiEjOiQhMUQnb3hVMyNSRyEjOTckJCIxW2NoXCNlNSQpKiEjOiQhMlpfVlpPN2B0IyEjOjckJCIyJ3klPTE+JT4sNSEjOiQhMmtRdSIzKEdXYiMhIzo3JCQiMkNSdmliI0c+NSEjOiQhMVtfZ15dYXRCISM5NyQkIjJvRSk9PFdOZTUhIzokITJYXXdDV0VHKT4hIzo3JCQiMnJnJnBzJT4oKTQiISM6JCEyOTUudSllPHo6ISM6NyQkIjIjPixhLykzI1E2ISM6JCEyMyl6JipvREclPSIhIzo3JCQiMm5ILk93O2s8IiEjOiQhMV8/RSR5SD8tKSEjOjckJCIyY0VreT1TKT03ISM6JCEyQztsck4ob3pQISM7NyQkIjFsMlVAJmZwRCIhIzkkIjElKnlaWykqZkNLISM8NyQkIjJCOyxrL2d3SCIhIzokIjEzWF8nXClILVQhIzo3JCQiMjMlM0tyMmFNOCEjOiQiMV9CXCkpNFAheighIzo3JCQiMWNTOFg2J1tQIiEjOSQiMnZRIilwJDNDIz0iISM6NyQkIjJHM19ZWip6NzkhIzokIjFjOztLVGloOiEjOTckJCIyRD1KLydmWF85ISM6JCIxYEUmKiopKik9ZT4hIzk3JCQiMm4iZiNRVUY3XCIhIzokIjJZKioqKlFpLmZNIyEjOjckJCIyRCZwTClbQTpeIiEjOiQiMV8uLHBVJilbRCEjOTckJCIxKSl6JUdiPD1gIiEjOSQiMTMyNzlcIT12IyEjOTckJCIyKG9HWkszZlQ6ISM6JCIxOSZwLnJQJlxHISM5NyQkIjIkXHg0N1RPXjohIzokIjItSz1tXXFzJUghIzo3JCQiMicqPTU+dl1pYiIhIzokIjJMc1VaIXA4JypIISM6NyQkIjIqSEVzIlJQNmMiISM6JCIyazduR0kuXS8kISM6NyQkIjAmKUc7ciFlajohIzgkIjFHJEg+XU8lcEkhIzk3JCQiMi0yTjouQ2djIiEjOiQiMiZIOio0cXBRNCQhIzo3JCQiMi49KVsicFhzYyIhIzokIjItakEwSSczMUohIzo3JCQiMi5IVDlObiVvOiEjOiQiMUpQMCtIST1KISM5NyQkIjIvUyVSNiEqb3A6ISM6JCIyPSRbZSpcPjA4JCEjOjckJCIyMF5aOG41NGQiISM6JCEyUGxqITMjXC85JCEjOjckJCIxKVxfOCYpKiozZiIhIzkkITIkelAsM3ViU0ghIzo3JCQiMmFbZDguKikzaCIhIzokITEwUid6Z2wxdSMhIzk3JCQiMiNSPk9damBdOyEjOiQhMmpPPnpUIz5XQiEjOjckJCIyJ3BkIio9Im9wbyIhIzokITFqNVFLWigpej4hIzk3JCQiMjlWLHEuQihHPCEjOiQhMl9XQzpiREJjIiEjOjckJCIyJlI3XUwicGd3IiEjOiQhMWtqXydlaykpPSIhIzk3JCQiMUF2c2h5KWUhPSEjOSQhMSRSTkUvdG4heiEjOjckJCIxOGBSbCUqKlIlPSEjOSQhMS1qJmVuN2M0JSEjOjckJCItUillYlwpPSEjNSQhMiUpKmZYbWYoUXckISNCLSUmQ09MT1JHNiYlJFJHQkckIiM1ISIiJCIiISEiIiQiIiEhIiItJSVWSUVXRzYkOyQiIiEhIiIkIisjZmJcKT0hIiklKERFRkFVTFRHLSUqQVhFU1NUWUxFRzYjJSdOT1JNQUxHLSUoU0NBTElOR0c2IyUuVU5DT05TVFJBSU5FREctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQhUiEiIi0lKUJPVU5EU19ZRzYjJCIjISohIiItJS1CT1VORFNfV0lEVEhHNiMkIiVnTiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiU1USEiIi0lKUNISUxEUkVORzYi
b2 := n-> (1/Pi)*int(F2(t)*sin(n*t),t=-Pi..Pi);
Zio2I0kibkc2IkYlNiRJKW9wZXJhdG9yR0YlSSZhcnJvd0dGJUYlKiZJI1BpRyUqcHJvdGVjdGVkRyEiIi1JJGludEc2JEYrSShfc3lzbGliR0YlNiQqJi1JI0YyR0YlNiNJInRHRiUiIiItSSRzaW5HRi82IyomOSRGN0Y2RjdGNy9GNjssJEYqRixGKkY3RiVGJUYl
b2(n) assuming integer;
LCQqJkkibkc2IiEiIilGJkYkIiIiISM/
Convolution solutions obtained using Laplace transform
x1(t) := (1/8)*int(sin(4*(t-tau))*F1per(tau),tau=0..t):
x2(t) := (1/8)*int(sin(4*(t-tau))*F2per(tau),tau=0..t):
The displacement is periodic
plot(x1(t),t=0..6*Pi,color=black);
NigtJSdDVVJWRVNHNiM3Z3g3JCQiIiEhIiIlKnVuZGVmaW5lZEc3JCQiMkBfcTZPZVJHIiEjPSQiMXlnLTtQJDRDKSEjPjckJCIyVS9UQnM7emMjISM9JCIyJnB5eHEtPyVIJCEjPjckJCIxbTpeJDN2PSZRISM8JCIxKSpHaCJ5MVFTKCEjPTckJCIyJSkzI29XTSRlOCYhIz0kIjJAWHZsYTJVSiIhIz03JCQiMUtKLW4sdi54ISM8JCIxVCw0ejkpUiVIISM8NyQkIjJ4VE8qKW9tci0iISM8JCIwYSczZjtjLF8hIzs3JCQiMmxpL00uXTJhIiEjPCQiMk5gKVIwImUpXDYhIzw3JCQiMmEkRyh5UExWMCMhIzwkIjIzNztRb1dTKj4hIzw3JCQiMlYvVEJzO3pjIyEjPCQiMVZdRFBSQDxJISM7NyQkIjFgIzRvMSs6MyQhIzskIjImM3pkXG9Md1QhIzw3JCQiMWl1RjZNMyZmJCEjOyQiMWhPRUFwbUFhISM7NyQkIjFyY3Vibm0zVCEjOyQiMXhuIzQ+IXoubiEjOzckJCIxMillJioqKilSLV0hIzskIjE6TFJTK1BjKSkhIzs3JCQiMVc+UFY3OCcqZSEjOyQiMiMqXGckPSZRejEiISM7NyQkIjF6XT0oW2opKnknISM7JCIyJHp1JSoqR0tVPiIhIzs3JCQiMTsjKSo0dCZmJG8oISM7JCIxNGgqUioqWyZbNyEjOjckJCIxRU5xNjZwKW8pISM7JCIyaWcpPU5hWzo3ISM7NyQkIjFOKTNDXCd5JHAqISM7JCIyb3RyPGlkIykzIiEjOzckJCIyVlQ2dD0pKSlwNSEjOyQiMUUpXF91VzooKSkhIzs3JCQiMmAlPlFEeFJxNiEjOyQiMWAuWWt2WFVrISM7NyQkIjI4JFIiZl10OkYiISM7JCIwV0NUKGY7blIhIzo3JCQiMTxmVydHXEZQIiEjOiQiMidcU08pNCNwZz0hIzw3JCQiMS56KHAxRFJaIiEjOiQiMUpUQy5uZ0xZISM8NyQkIjIqKSkpNHYlMzV2OiEjOyQiMS9dNDovQWsjKiEjPzckJCIyLTI9c3YmenY7ISM7JCIway16OyQpPVYmISM7NyQkIjI6REVwbSFcdzwhIzskIjJiM00xW0QpKSo+ISM8NyQkIjJEVk1tZCY9eD0hIzskIjEoRysvPGpbOCUhIzs3JCQiMlFoVWpbISl5KD4hIzskIjFsV29KJkclNG0hIzs3JCQiMWN4W2R2QnI/ISM6JCIxJUgobz9SJTQnKSkhIzs3JCQiMSkqR2pHWWZrQCEjOiQiMkdyNz4heUR2NSEjOzckJCIwL3koKnBeekQjISM5JCIxYipSNmhmQj8iISM6NyQkIjI8PUI0eDM4TiMhIzskIjIkZllgXzEpKVw3ISM7NyQkIjJrUFQxaHV6VyMhIzskIjFUITRLKEhOMzchIzo3JCQiMXImZi5YU1lhIyEjOiQiMjtZJFFsLXAhMyIhIzs3JCQiMW14MiFIMThrIyEjOiQiMVRgIzN3RHUmKSkhIzs3JCQiMjInZnpIQCh6dCMhIzskIjEkbygpW0BcSF8nISM7NyQkIjFKV156UyV6JEchIzokIjFrI3ovYTwwMiUhIzs3JCQiMjAhSEJIZyJ6JEghIzskIjJDdS9BeigpPic+ISM8NyQkIjBQXip5eil5LiQhIzkkIjElKnp3RUFtK2AhIzw3JCQiMi4lKXAnRypmeTgkISM7JCIxXkpqVCMpW21wISM/NyQkIjIoeXcybjdeUEshIzskITJsTVkhW10uV1ghIz03JCQiMnZeJlswRTtQTCEjOyQhMSh5TDoiRyNvIj0hIzs3JCQiMmpOJCpRJVIib1YkISM7JCEwKjROUWZMdVEhIzo3JCQiMlk+LEJHbGtgJCEjOyQhMSNvXzIqekwvaiEjOzckJCIxVSxjSypvKlFPISM6JCExJCpmPD4/OCJ6KSEjOzckJCIwND5HZXM5dSQhIzkkITIjNFdnZyNvYzMiISM7NyQkIjJ1LnlJQndSJVEhIzskITJASk5UTWVlQCIhIzs3JCQiMlgpcEwkKSl6ayVSISM7JCEyYyNIL1Y9NVs3ISM7NyQkIjJiJyp5VTpsbi4lISM7JCEyOm1ZZCYpNDI+IiEjOzckJCIxWTRBRC8wRlQhIzokITE5PDkheU0uMSIhIzo3JCQiMUZIOydwTnRAJSEjOiQhMWE3KSlHXSd6dCkhIzs3JCQiMTNcNW40aTJWISM6JCExb3JfMWhnXmwhIzs3JCQiMTp0PCl6ZyM0VyEjOiQhMTg7TUQiSCRlUyEjOzckJCIxQShcI0gxITReJSEjOiQhMCNlTD4/TkE+ISM6NyQkIjA4QS5ZU0RoJSEjOSQhMHYrbTxeKj1cISM7NyQkIjFQWFIiSCE9OVohIzokITJPdCF5LTJSLzshI0A3JCQiMSJveWBZUGkiWyEjOiQhMUp3R1k5OzpgISM8NyQkIjFER09SWUg9XCEjOiQhMiw3XmpoKnAtPyEjPDckJCIxKSpbTkVLS3BcISM6JCEyRXFPJWZcRT9JISM8NyQkIjAocE04PU4/XSEjOSQhMiY0K0ZeSiM+PCUhIzw3JCQiMVUhUi5TIVFyXSEjOiQhMXklM1RSZik0YSEjOzckJCIxOjZMKCkqM0M3JiEjOiQhMC5SInpubiNvJyEjOjckJCIxNyQqPiYpKmYyQSYhIzokITFpR1ApZTwoUiEqISM7NyQkIjE0djEkKTQ2PmAhIzokITJ3X3U0R2FxNCIhIzs3JCQiMTJkJDQpPlk8YSEjOiQhMlppZzYjZiw9NyEjOzckJCIxL1IhKXlIImVeJiEjOiQhMl5SLiwwdyRbNyEjOzckJCIxNEExMl43MGMhIzokITIxMl1eTXNLPiIhIzs3JCQiMTkwS05zViVwJiEjOiQhMjNkRWE2NWsxIiEjOzckJCIxPil5Tk9cUHkmISM6JCExOl5WVigpM1EpKSEjOzckJCIxQnIkPVxoSShlISM6JCEwQF9gc2xfbychIzo3JCQiMSJvODgoPTtFZiEjOiQhMUwkKjRjRzhoYCEjOzckJCIyJlItel1BRXpmISM7JCExVSpcYjtdcDIlISM7NyQkIjEpem0tamlCLichIzokITEkekgnNGBWISpHISM7NyQkIjFjTHU0SVkmMychIzokITJuUSwmZWkiXCY9ISM8NyQkIjF0a3BvUG0iPichIzokITE0Yih5c0k6OSUhIzw3JCQiMiYqZVx3X2t5SCchIzskIjIjbzMhMzEpM3g1ISM+NyQkIjFtJGYzXUl5USchIzokIjFSZlI+KEhlUiYhIzw3JCQiMVQicFNaJ3p4ayEjOiQiMSVmcGlJZioqeiIhIzs3JCQiMTsqeXNXaXhjJyEjOiQiMmN4MDtaSy9qJCEjPDckJCIyOnApWz8lR3hsJyEjOyQiMT4/anonKlwneiYhIzs3JCQiL09uRlYwNW4hIzgkIjEwYSUqKj4iZixyISM7NyQkIjEyJmVbQiFRaW4hIzokIjFgImYtPzcmcCQpISM7NyQkIjJYVFY/OTFaIm8hIzskIjEpPm1uMT9cYSohIzs3JCQiMUIkRyNcPy5ubyEjOiQiMk1FJ3pBNWxkNSEjOzckJCIxUSIpZmpRb3JwISM6JCIyMCdHLV5yai43ISM7NyQkIjFheid6bk5qMighIzokIjEyJnovYSpwXDchIzo3JCQiMW5GRG0tKCpvciEjOiQiMiJRNzJwZUcrNyEjOzckJCIxImVQWCZbZ2hzISM6JCIyTyhvKzRleXM1ISM7NyQkIjEmUkFHV1JVTighIzokIjF3J1JNc2JdJSkpISM7NyQkIjE0czVKUyhvVyghIzokIjEmXEtKazgraCchIzs3JCQiMVh3blY5cChcKCEjOiQiMSg0Y3UqPSpITSYhIzs3JCQiMSMzW2kmKTMmW3YhIzokIjFGXiQ0YjxMNiUhIzs3JCQiMT4mPSlvaUsqZighIzokIjJhbDdBRUM7KEghIzw3JCQiMWMqKVEib1YsbCghIzokIjJLSVYpUkgjXCc+ISM8NyQkIjImSClIbF15PHYoISM7JCIwPjhIJlxTXV4hIzs3JCQiMS4ybkpMVGB5ISM6JCIyJz0jM3AlbyRbaCIhI0E3JCQiMTB0dT5QPl16ISM6JCIyOkxqYScpUTpkJSEjPTckJCIxMlIjeTV1cC8pISM6JCIyZEVrXnNQO3giISM8NyQkIjBeK2ZcYVA5KSEjOSQiMUIpWyhcSCMqW1AhIzs3JCQiMTdyKFIpW2BTIykhIzokIjFMXyVHb3BqNCchIzs3JCQiMU1jK21dXlQkKSEjOiQiMTkqMzBnMVdjKSEjOzckJCIxYVQuW19cVSUpISM6JCIyXz5bbSZcKmYxIiEjOzckJCIxdUUxSWFaViYpISM6JCIyYTQ+MClIZC83ISM7NyQkIjEnPiI0N2NYVycpISM6JCIyaiEpXFw9cilcNyEjOzckJCIxJWYlelRwPVAoKSEjOiQiMjM/UFZ6dkY/IiEjOzckJCIxIip6XHIjPSpIKSkhIzokIjJUJzRsU1M0eDUhIzs3JCQiMSpRLDdnXEUjKikhIzokIjFBRy0nbzAjKiopKSEjOzckJCIxJnkvNCQ0UTohKiEjOiQiMS1DRiFSUXJtJyEjOzckJCIvdWIoPi9hNiohIzgkIjEkKSkpPilbQV8/JSEjOzckJCIxOytAa3VVOiMqISM6JCIyZUsiKlJBIUdtPyEjPDckJCIxTEUnM3RdYUoqISM6JCIxLmcmZTAnZSIpZSEjPDckJCIxW19eKCpSWjolKiEjOiQiMiYzQ3ojMyhvRlYhIz83JCQiMVkuYWcrUD4mKiEjOiQhMmxWVlIiSHI/VyEjPTckJCIxWmFjQmhFQicqISM6JCEyTk9yQnJcJm89ISM8NyQkIjEoKnoyYlRAdicqISM6JCEyJEg9Mlg9SCMpRyEjPDckJCIxWjBmJz1pcnMqISM6JCEyY18xWig0IzQvJSEjPDckJCIxKTQuIj0tNnooKiEjOiQhMWwlKXBkK2YlSCYhIzs3JCQiMVtjaFwjZTUkKSohIzokITFYISlbK0tPKmUnISM7NyQkIjE7LSF6Mis6IyoqISM6JCExdjhKUz89dygpISM7NyQkIjIneSU9MT4lPiw1ISM6JCEyJT1AbVR2Zmo1ISM7NyQkIjJjJHBXdCRRLSwiISM6JCEyTl1eS1speiM+IiEjOzckJCIyQ1J2aWIjRz41ISM6JCEyTVhyVFImXFs3ISM7NyQkIjE2T11AMDBINSEjOSQhMmNPLyo+Ryw9NyEjOzckJCIyJ0g9dCdbPSlRNSEjOiQhMmJ4YTw6WSIpNCIhIzs3JCQiMiRbKyc+WCdlWzUhIzokITEoKXBuYHVpcCEqISM7NyQkIjJvRSk9PFdOZTUhIzokITE0a08yV3pLbiEjOzckJCIxLV4xIm9YJW81ISM5JCExdygpUTNkWVlVISM7NyQkIjFQPiVcJXBgeTUhIzkkITIlXC5DQD88IzMjISM8NyQkIjFzKD0pMyNHJykzIiEjOSQhMThVZUV4KHkoZSEjPDckJCIycmcmcHMlPigpNCIhIzokITI5U2hoNkIwXiQhIz83JCQiMl50Y2MhPWYzNiEjOiQhMV5TKD0wem4uJSEjPDckJCIySid5aFFUWT02ISM6JCEyO2V0czZSUHEiISM8NyQkIjI3Kip5OlpPJEc2ISM6JCEyKFIjUmdfI2UuUCEjPDckJCIyIz4sYS8pMyNRNiEjOiQhMSMqPiUqeUVRJjQnISM7NyQkIjJPImZJJXpneDkiISM6JCExTT1RP3A4TyUpISM7NyQkIjEzPDIleTd0OiIhIzkkITIkUTt1LTc7WTUhIzs3JCQiMkJdUFF4a287IiEjOiQhMVoqeTAiKnB6PSIhIzo3JCQiMm5ILk93O2s8IiEjOiQhMkUleWNBJillWzchIzs3JCQiMVImbyc+RS0oPSIhIzkkITJQIlFBQixmNTchIzs3JCQiMjd5TGRaR3c+IiEjOiQhMVJ5Mz06eG81ISM6NyQkIjJBU21QU0pIPyIhIzokITEkb1BBXy08bSohIzs3JCQiMkwtKnpKVkIzNyEjOiQhMFU7ZWBmRlspISM6NyQkIjJXa0opZnNgODchIzokITEzSVAvRXQucyEjOzckJCIyY0VreT1TKT03ISM6JCExYjonZVJhPillISM7NyQkIjIwUmA3LXEkRzchIzokITIneWh4YydIc2UkISM8NyQkIjJgXlVZJikqKnlCIiEjOiQhMi1sRT8zbFpuIiEjPDckJCIyLWtKISlvSHVDIiEjOiQhMmwxJFJNJiozIj4lISM9NyQkIjFsMlVAJmZwRCIhIzkkIjEsSEQ3KyYqKT4mISM/NyQkIjJXJ2VtX1k4bjchIzokIjEkSGZDKmZbSGEhIzw3JCQiMlAnNCJSeTR0RiIhIzokIjJjPENmSCF6PD8hIzw3JCQiMkxeTCZcdFIjRyIhIzokIjI7J3lDJlsoZU1JISM8NyQkIjFqZzo6XFsoRyIhIzkkIjJWYmh3bCczJT0lISM8NyQkIjJGaHkyW3NESCIhIzokIjFYJSoqKipIWik9YSEjOzckJCIyQjssay9nd0giISM6JCIxKFJBdzU1em8nISM7NyQkIjEjM0p3QSEpb0kiISM5JCIxKD52S1VyZyEqKSEjOzckJCIyOitoKTMvNTs4ISM6JCIydzVeRU8ubjIiISM7NyQkIjI3IzQ0IWY/YEsiISM6JCIxaSg+dXFfPz8iISM6NyQkIjIzJTNLcjJhTTghIzokIjJ3JFI/JnAoelw3ISM7NyQkIjJaOXVaJzNpVzghIzokIjItUio+XEwkZj8iISM7NyQkIjIlW3VBZTRxYTghIzokIjJ2Vl5JUCgqKW81ISM7NyQkIjJAdiFvXjV5azghIzokIjFEaCp5ZnptZykhIzs3JCQiMWNTOFg2J1tQIiEjOSQiMXN1bFhcVFloISM7NyQkIjJGYzh2c1hWUSIhIzokIjEmbyJSRmFVUlEhIzs3JCQiMiVwSSopNC4kUVIiISM6JCIyWCQqZTlPX18oPSEjPDckJCIyaGRzQSpbSi45ISM6JCIxeXMlXFR5QUwmISM8NyQkIjJHM19ZWip6NzkhIzokIjFFSy4kKmUpZj8lISM+NyQkIjJ5Jm80J2Y4RlUiISM6JCIxLnIpKillP1ArJSEjPDckJCIyRWpUdnJGRVYiISM6JCIxJ3prKCo9eldxIiEjOzckJCIyd1MnKSpRPWFVOSEjOiQiMWhtIlJRN1RyJCEjOzckJCIyRD1KLydmWF85ISM6JCIxJnkoWyhcXHQ2JyEjOzckJCIxbSl6aSMpW0BZIiEjOSQiMTZjdG9OaSpbKSEjOzckJCIyJlwmR0BwVD1aIiEjOiQiMkQlZjUjeVVIMCIhIzs3JCQiMUxzKHpiTTpbIiEjOSQiMVJqeVNiUiQ+IiEjOjckJCIybiJmI1FVRjdcIiEjOiQiMmRtPV5jcSVcNyEjOzckJCIyWlYiM2NcUCw6ISM6JCIyLkBSLCNlKSkzNyEjOzckJCIyRCZwTClbQTpeIiEjOiQiMkI0QngrO00yIiEjOzckJCIyLlojZj8rbkA6ISM6JCIxUVtsRSJlMmwpISM7NyQkIjEpKXolR2I8PWAiISM5JCIxak5MKjM9czwnISM7NyQkIjIob0daSzNmVDohIzokIjFfLDEmPTY4IVEhIzs3JCQiMiRceDQ3VE9eOiEjOiQiMV9vZ0soKSlbeiIhIzs3JCQiMipIRXMiUlA2YyIhIzokIi8kKXp5RywyWSEjOjckJCIyMF5aOG41NGQiISM6JCExQGR3ZnBLUGwhI0A3JCQiMlYrXThFMDRlIiEjOiQhMk4zdEMjb0dTXSEjPTckJCIxKVxfOCYpKiozZiIhIzkkITEnKXBrU2o9OT4hIzs3JCQiMjsqXE5UVyozZyIhIzokITFqTztvSTkzUyEjOzckJCIyYVtkOC4qKTNoIiEjOiQhMSN6ISlcJz5tYmshIzs3JCQiMigpZjM2JzMhM2kiISM6JCExOyZ6YldbPiYpKSEjOzckJCIyQnJmM3A3MmoiISM6JCEybzRkNDxjVzMiISM7NyQkIjFFM2g/WGlTOyEjOSQhMlcvVkNDT0NAIiEjOzckJCIyI1I+T11qYF07ISM6JCEyKGUxaiopKXojXDchIzs3JCQiMShSK0RIVydmOyEjOSQhMjh2U0FcPnc+IiEjOzckJCIyVyYpUVlCXyhvOyEjOiQhMng6LlF0JHpxNSEjOzckJCIyPUp4bjxneW4iISM6JCExcE1NLF1dYSkpISM7NyQkIjIncGQiKj0ib3BvIiEjOiQhMUojSHo/IT5mbSEjOzckJCIyc1p3T1soPSNwIiEjOiQhMWIkKVwnRyx4TiYhIzs3JCQiMSY9UCVbb1MocCIhIzkkITJFWDRieGpcNCUhIzw3JCQiMkQqeT44aWktPCEjOiQhMkRHJkh3biJlI0ghIzw3JCQiMjBnZXpkWHlxIiEjOiQhMmFCPD8oWy4sPiEjPDckJCIxOytbMlZHPTwhIzkkITI5X0UzOFJRYSUhIz03JCQiMjlWLHEuQihHPCEjOiQhMnczLVhYZXRlJCEjPzckJCIyT1FFNmNmIVE8ISM6JCEyWWhMb0okKVI2JiEjPTckJCIyYUxeXzMnUlo8ISM6JCEyc2ApKXpMSS4iPSEjPDckJCIydEd3JDRFdGM8ISM6JCExR3YpeT44OHMkISM7NyQkIjImUjddTCJwZ3ciISM6JCExbCNweG4qKjMpZiEjOzckJCIwInliOlEtdzwhIzgkITE2bSJHeklKVSkhIzs3JCQiMjJROXdceWZ5IiEjOiQhMFwjKTRRTEQwIiEjOTckJCIyOCY0bnpKJGZ6IiEjOiQhMUooNFgneSZlPiIhIzo3JCQiMUF2c2h5KWUhPSEjOSQhMSpcTyZlMScpXDchIzo3JCQiMihwV2tpZFQ6PSEjOiQhMVE+eUljJCo0NyEjOjckJCIydVRoTm1WXCM9ISM6JCExOCRvR1hzZzMiISM6NyQkIjFsJHlXY3JXJD0hIzkkITF6I3ojcFlYZyopISM7NyQkIjE4YFJsJSoqUiU9ISM5JCExXj96KDNON24nISM7NyQkIjInWydwMikqPSJcPSEjOiQhMCc+cWw3UCVSJiEjOjckJCIyWilSOSdcUVUmPSEjOiQhMkwmZl0xSkVgVCEjPDckJCIyM0s9OixlJGY9ISM6JCEyJXAsJil6XHgqKkghIzw3JCQiMm1sIypvX3hXJz0hIzokITI9P2xIczRAKT4hIzw3JCQiMUc4a2RscnU9ISM5JCExQipwVHNFIXBeISM8NyQkIi1SKWViXCk9ISM1JCExLTVnZHRVMHIhI0ktJSVWSUVXRzYkOyQiIiEhIiIkIisjZmJcKT0hIiklKERFRkFVTFRHLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpBWEVTU1RZTEVHNiMlJ05PUk1BTEctJShTQ0FMSU5HRzYjJS5VTkNPTlNUUkFJTkVERy0lJVJPT1RHNictJSlCT1VORFNfWEc2IyQiJD8lISIiLSUpQk9VTkRTX1lHNiMkIiMhKiEiIi0lLUJPVU5EU19XSURUSEc2IyQiJUlOISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJTVRISIiLSUpQ0hJTERSRU5HNiI=
The amplitude of the displacement increases with time (linearly) we have resonance.
plot(x2(t),t=0..6*Pi,color=black);
NigtJSdDVVJWRVNHNiM3Z3o3JCQiIiEhIiIkIjFFMV0pNCMqM1clISNKNyQkIjJhJEcoeVBMVjAjISM8JCIyKXpiO3osKHBSIiEjPTckJCIxcmN1Ym5tM1QhIzskIjIuIVxLdjFhNDUhIzw3JCQiMVc+UFY3OCcqZSEjOyQiMW1ZK2xfcyNlIyEjOzckJCIxOyMpKjR0JmYkbyghIzskIjF3JSlmbCZRZXAlISM7NyQkIjFOKTNDXCd5JHAqISM7JCIxTi5HeHJaMnIhIzs3JCQiMmAlPlFEeFJxNiEjOyQiMTc3bTRbdXcpKSEjOzckJCIxPGZXJ0dcRlAiISM6JCIxZGsqeiJRLCNwKiEjOzckJCIyKikpKTR2JTM1djohIzskIjFuKVs8UHl1IikqISM7NyQkIjI6REVwbSFcdzwhIzskIjFrK19NXnBkKiohIzs3JCQiMlFoVWpbISl5KD4hIzskIjJKKWVvJipRPSEzIiEjOzckJCIxKSpHakdZZmtAISM6JCIyOl80SDUwWEMiISM7NyQkIjI8PUI0eDM4TiMhIzskIjJqSV87ZzlsWSIhIzs3JCQiMXImZi5YU1lhIyEjOiQiMjhUeHE7UXRwIiEjOzckJCIyMidmekhAKHp0IyEjOyQiMiMqPmh1XUx0Jz0hIzs3JCQiMjAhSEJIZyJ6JEghIzskIjI+MyRRJT10KVw+ISM7NyQkIjIuJSlwJ0cqZnk4JCEjOyQiMmstOD1LJlxqPiEjOzckJCIyKHl3Mm43XlBLISM7JCIyKSpmQzZIWCV6OyEjOzckJCIydl4mWzBFO1BMISM7JCIwPW8oUSNHL00pISM6NyQkIjJqTiQqUSVSIm9WJCEjOyQhMSY+MGJtMSQzViEjOzckJCIyWT4sQkdsa2AkISM7JCEyZCJbcndjMzI+ISM7NyQkIjFVLGNLKm8qUU8hIzokITEpPkFeJ3ktI1IkISM6NyQkIjA0Pkdlczl1JCEjOSQhMWArPFk0aillJSEjOjckJCIyT2NbelNDRnokISM7JCExPDAhKT1CZzVdISM6NyQkIjJ1LnlJQndSJVEhIzskITF1JnlbLDl6RyYhIzo3JCQiMlZ4VWM5LSdwUSEjOyQhMWg7NyswLW9gISM6NyQkIjE2dj9lIUdfKlEhIzokITF1TG8kUjV3UyYhIzo3JCQiMndDczIoUiYzI1IhIzskITImZVQnUSYqKjQxYSEjOzckJCIyWClwTCQpKXprJVIhIzskITFYdzYtcExqYCEjOjckJCIyYicqeVU6bG4uJSEjOyQhMmxDUy8vWj8qWyEjOzckJCIxWTRBRC8wRlQhIzokITEkKXllYDB3cVIhIzo3JCQiMUZIOydwTnRAJSEjOiQhMicpPnA5NmE1ciMhIzs3JCQiMTNcNW40aTJWISM6JCEyLzEleWNXO283ISM7NyQkIjE6dDwpemcjNFchIzokIjA9Q2NXRj5fJCEjOjckJCIxQShcI0gxITReJSEjOiQiMlc0UmhbJz5DPCEjOzckJCIwOEEuWVNEaCUhIzkkIjE0UF50LmBNRSEjOjckJCIxUFhSIkghPTlaISM6JCIxPCV6YUpVXiVIISM6NyQkIjEib3lgWVBpIlshIzokIjJPJm8iR0FLSmgjISM7NyQkIjFER09SWUg9XCEjOiQiMilvJjM5I1IoNHEiISM7NyQkIjAocE04PU4/XSEjOSQiMSo9NGtdSDFwJCEjOzckJCIxOjZMKCkqM0M3JiEjOiQhMksnZSgqcCI+SzoiISM7NyQkIjE3JCo+JikqZjJBJiEjOiQhMUJfKUdNb21iIyEjOjckJCIxNHYxJCk0Nj5gISM6JCEyRCFlbFdqJjRuJCEjOzckJCIxMzsrI1snR29gISM6JCExdyFSJHBuKTQxJSEjOjckJCIxMmQkNCk+WTxhISM6JCExOm09JiopXGtKJSEjOjckJCIxY0ZTSShcP1cmISM6JCExb3IrJGVgK1IlISM6NyQkIjExKXApenVqbWEhIzokITFGIWZ3IXpJRVchIzo3JCQiMWJvTEhfQSJcJiEjOiQhMTYpNCRSUnJDVyEjOjckJCIxL1IhKXlIImVeJiEjOiQhMSkzemxiUV5RJSEjOjckJCIxNEExMl43MGMhIzokITJaR04helNJSFIhIzs3JCQiMTkwS05zViVwJiEjOiQhMid6Ky07NCIzLiQhIzs3JCQiMT4peU5PXFB5JiEjOiQhMjFVdVBebmd6IiEjOzckJCIxQnIkPVxoSShlISM6JCExKSo0dSd5QSVSUCEjOzckJCIyJlItel1BRXpmISM7JCIyKmV3OnlyIj5LIiEjOzckJCIxY0x1NElZJjMnISM6JCIxLndDbmwtXEYhIzo3JCQiMlgiKj4jKlFqJlFoISM7JCIxYjZVS1YyJEckISM6NyQkIjF0a3BvUG0iPichIzokIjJgVHU2aixibSQhIzs3JCQiMV9aVmVSQD1pISM6JCIybCpRXSNva1l6JCEjOzckJCIxSkk8W1R3V2khIzokIjIuZGlQKD1oISlRISM7NyQkIjBKNnpMOThGJyEjOSQiMmBUIVJKPGNBUiEjOzckJCIyJiplXHdfa3lIJyEjOyQiMWhlKVJhRy0jUiEjOjckJCIxVCJwU1onenhrISM6JCIxLiY+KipcZ3ohRyEjOjckJCIyOnApWz8lR3hsJyEjOyQiMk9UY3FvUj5qJCEjPDckJCIxMiZlW0IhUWluISM6JCEyKSpRLVwnSEF6NiEjOzckJCIxQiRHI1w/Lm5vISM6JCEybGVSZmgjR21DISM7NyQkIjFKS1Rjek4+cCEjOiQhMiVRKikqKXpsOlFIISM7NyQkIjFRIilmalFvcnAhIzokITFwMmBbK1ZrSyEjOjckJCIxI2YhPjxvJXkqcCEjOiQhMihveW47NiF5TyQhIzs3JCQiMVlJeXEoNFMtKCEjOiQhMiRlJlxTRk0lSE0hIzs3JCQiL2JQQ0Y8XXEhIzgkITJtSUl6P3clW00hIzs3JCQiMWF6J3puTmoyKCEjOiQhMSJvOCNSJ1JYVSQhIzo3JCQiMW5GRG0tKCpvciEjOiQhMWEkek0iSCUqKipIISM6NyQkIjEiZVBYJltnaHMhIzokITE+ZmJpVCskNCMhIzo3JCQiMSZSQUdXUlVOKCEjOiQhMTk8aUoxbyo9KSEjOzckJCIxNHM1SlMob1coISM6JCIxbCVmOHVGNGQnISM7NyQkIjEjM1tpJikzJlt2ISM6JCIxRjJwXCp6LEcjISM6NyQkIjFjKilRIm9WLGwoISM6JCIyTyQ+XkVyW2dPISM7NyQkIjImSClIbF15PHYoISM7JCIwPFVhZGtMZSUhIzk3JCQiMS4ybkpMVGB5ISM6JCIxJyoqM1lQRygzXCEjOjckJCIxMHR1PlA+XXohIzokIjFvQChRbXRIaSUhIzo3JCQiMTJSI3k1dXAvKSEjOiQiMjxwYlw2OXMhUSEjOzckJCIwXitmXGFQOSkhIzkkIjJsJFxda104ImYjISM7NyQkIjE3cihSKVtgUyMpISM6JCIyaz1HQlpwTzsiISM7NyQkIjFNYyttXV5UJCkhIzokITEsRCZvdTEoR0ohIzs3JCQiMWFULltfXFUlKSEjOiQhMi8jPjNWZitLOiEjOzckJCIxOCVbIVJgKUhcKSEjOiQhMlZGNSlIdiZ6KD4hIzs3JCQiMXVFMUlhWlYmKSEjOiQhMkNwTjAnUk8oRyMhIzs3JCQiMTApcGJaPyhvJikhIzokITIuSl40aShIJ1EjISM7NyQkIjFOcDJAYidSZikhIzokITJ4WXFzRGxqVyMhIzs3JCQiMW1TZW0wQD4nKSEjOiQhMnVyaXI3JXptQyEjOzckJCIxJz4iNDdjWFcnKSEjOiQhMnhLUDhNOnNXIyEjOzckJCIxJWYlelRwPVAoKSEjOiQhMkNVIno5JVJwLiMhIzs3JCQiMSIqelxyIz0qSCkpISM6JCEyJ1Iqb1QnWyo0OSIhIzs3JCQiMSpRLDdnXEUjKikhIzokIjJXXiopZnIiUm03ISM8NyQkIjEmeS80JDRROiEqISM6JCIybyFbMCgzRzlnIiEjOzckJCIvdWIoPi9hNiohIzgkIjEnXFQsZGZEPyQhIzo3JCQiMTsrQGt1VTojKiEjOiQiMUs8RGMiUnVkJSEjOjckJCIxTEUnM3RdYUoqISM6JCIxXmNpRCJ5KD1iISM6NyQkIjFbX14oKlJaOiUqISM6JCIxM1xGPGR3KCllISM6NyQkIjFZLmFnK1A+JiohIzokIjFtMmFDT09PYCEjOjckJCIxWmFjQmhFQicqISM6JCIyTUEheiEpZi1iTiEjOzckJCIxKCp6MmJUQHYnKiEjOiQiMS5bT2ByUSRIIyEjOjckJCIxWjBmJz1pcnMqISM6JCIxQVIkKXAnM01iKSEjOzckJCIxKTQuIj0tNnooKiEjOiQhMSpmQUsoelplcCEjOzckJCIxW2NoXCNlNSQpKiEjOiQhMUt6Pik0WT9IIyEjOjckJCIxS3p2aiJ6aSgpKiEjOiQhMiZwMXFlPWZqTyEjOzckJCIxOy0hejIrOiMqKiEjOiQhMVgzbyg0Yy8oXCEjOjckJCIvRC8jKjRzbSoqISM4JCExJ2Z6JSoqUTJwaCEjOjckJCIyJ3klPTE+JT4sNSEjOiQhMVZbSDFDVD5zISM6NyQkIjJjJHBXdCRRLSwiISM6JCExSjhuKyxZUygpISM6NyQkIjJDUnZpYiNHPjUhIzokITEmekNWRHYkSCQqISM6NyQkIjE2T11AMDBINSEjOSQhMEFqY18pM0QpKSEjOTckJCIyJ0g9dCdbPSlRNSEjOiQhMSlmMD0vMV0/KCEjOjckJCIyKlFmTXBDcVY1ISM6JCExWWpoTS9BV2chIzo3JCQiMiRbKyc+WCdlWzUhIzokITImZkZuJykpPVRxJSEjOzckJCIydjp1WFZxTTAiISM6JCExUUhfckhfTUshIzo3JCQiMm9FKT08V05lNSEjOiQhMmlUKkhdMj0hcCIhIzs3JCQiMlZvRSJcK1NqNSEjOiQhMSY9XE1KRCZReCEjPDckJCIxLV4xIm9YJW81ISM5JCIyRi1WYVZCJipbIiEjOzckJCIyJz5OKzg4XHQ1ISM6JCIydSFcNnNpQFtIISM7NyQkIjFQPiVcJXBgeTUhIzkkIjFabEh2c3ZTVSEjOjckJCIxcyg9KTMjRycpMyIhIzkkIjE2SDNUSFdKaCEjOjckJCIycmcmcHMlPigpNCIhIzokIjFZTW8mXEB5J28hIzo3JCQiMl50Y2MhPWYzNiEjOiQiMXVxRSpHeWhPJyEjOjckJCIySid5aFFUWT02ISM6JCIwRydlWSYpPlVaISM5NyQkIjFGJSk0MC5TQjYhIzkkIjIuMCozWS5nJWUkISM7NyQkIjI3Kip5OlpPJEc2ISM6JCIxT04oelldYEQjISM6NyQkIjJfYmYhUUVGTDYhIzokIjFEaFp3bCt0ISkhIzs3JCQiMiM+LGEvKTMjUTYhIzokITFRbG4jSDQzLSghIzs3JCQiMmssQiUqeiUpSDkiISM6JCEyKHovN01ZXmtAISM7NyQkIjJPImZJJXpneDkiISM6JCEwQ0BNQFBTZCQhIzk3JCQiMjMiKSk9KnlPRDoiISM6JCExTFJ5MW1GeVshIzo3JCQiMTM8MiV5N3Q6IiEjOSQhMXhDUz06cEdnISM6NyQkIjJCXVBReGtvOyIhIzokITFvR1RRK1YueCEjOjckJCIybkguT3c7azwiISM6JCExKVwtJltrNlskKSEjOjckJCIxUiZvJz5FLSg9IiEjOSQhME5UUW9FI1J4ISM5NyQkIjI3eUxkWkd3PiIhIzokITEvcEdQXFBOZSEjOjckJCIyQVNtUFNKSD8iISM6JCExd3Zybl8sI1wlISM6NyQkIjJMLSp6SlZCMzchIzokITFeQyUzMVlCJ0ghIzo3JCQiMldrSilmc2A4NyEjOiQhMl9gc1pVXE1KIiEjOzckJCIyY0VreT1TKT03ISM6JCIybUwwQk8iW0FRISM8NyQkIjFHKWVYNTBPQSIhIzkkIjJDKipwXiQpKlIlKT0hIzs3JCQiMjBSYDctcSRHNyEjOiQiMXdMdzRuIzNKJCEjOjckJCIxYHoleiRcOEw3ISM5JCIxTCdmYTJZNGglISM6NyQkIjJgXlVZJikqKnlCIiEjOiQiMTsvXmUoXChRZCEjOjckJCIyLWtKISlvSHVDIiEjOiQiMUB5aVs7LUV0ISM6NyQkIjFsMlVAJmZwRCIhIzkkIjF3KkgwUEdMJnkhIzo3JCQiMlcnZW1fWThuNyEjOiQiMlh2I3lBTWd0ciEjOzckJCIyUCc0IlJ5NHRGIiEjOiQiMUNfXzQzZUtgISM6NyQkIjJMXkwmXHRSI0ciISM6JCIxSCVwYkZBdzElISM6NyQkIjFqZzo6XFsoRyIhIzkkIjEtJD1kPk45ayMhIzo3JCQiMkZoeTJbc0RIIiEjOiQiMmw5WSxYIj45NiEjOzckJCIyQjssay9nd0giISM6JCEwQ0pSN2t3XCUhIzo3JCQiMkE3O3E4cUFJIiEjOiQhMiZwQTl5RmFWPSEjOzckJCIxIzNKd0EhKW9JIiEjOSQhMSYqZlljbDhsSiEjOjckJCIyPC9ZIz0uXDY4ISM6JCExMUIkKVx2dm9WISM6NyQkIjI6K2gpMy81OzghIzokITFHU15NZWg3YSEjOjckJCIyNyM0NCFmP2BLIiEjOiQhMVFVd1pXPyMpbyEjOjckJCIyMyUzS3IyYU04ISM6JCExXFpfWWteb3QhIzo3JCQiMlo5dVonM2lXOCEjOiQhMEtfdCNScCNwJyEjOTckJCIyJVt1QWU0cWE4ISM6JCEydmAwbz4mR2JbISM7NyQkIjItNWFcK1QoZjghIzokITI7LWcocDlyJGUkISM7NyQkIjJAdiFvXjV5azghIzokITFcKnlmYHUzOSMhIzo3JCQiMlRTMiUpNEApcDghIzokITEoM1tXdzwlUmUhIzs3JCQiMWNTOFg2J1tQIiEjOSQiMiUzLUAhKXlIRDUhIzs3JCQiMiQ0UUtPTWd6OCEjOiQiMk1qJFxSVS1JRCEjOzckJCIyRmM4dnNYVlEiISM6JCIyRFxUbTJkOihSISM7NyQkIjJpSi4oPSEpMypRIiEjOiQiMlhNemRcVyMqSCYhIzs3JCQiMiVwSSopNC4kUVIiISM6JCIyMTg7RyMpPWxZJyEjOzckJCIyaGRzQSpbSi45ISM6JCIxJGUwPSplemoiKSEjOjckJCIyRzNfWVoqejc5ISM6JCIxYUpaYUVXSSkpISM6NyQkIjJ5Jm80J2Y4RlUiISM6JCIxUFIiNClHIlFMKSEjOjckJCIyRWpUdnJGRVYiISM6JCIxJltXTElyWnEnISM6NyQkIjAtayN5WmVQOSEjOCQiMWJbYHA2cVRiISM6NyQkIjJ3UycpKlE9YVU5ISM6JCIxXiZRQWwleTBVISM6NyQkIjJfejMoKiopKVxaOSEjOiQiMid6IjMpcG5lXUYhIzs3JCQiMkQ9Si8nZlhfOSEjOiQiMXY1bl9vTU03ISM6NyQkIjJWX2JMUi10WCIhIzokITIkcEwpb2cwKilbIyEjPDckJCIxbSl6aSMpW0BZIiEjOSQhMipHP2EiMz5lbiIhIzs3JCQiMnk/LyNmXypwWSIhIzokITJ0VWl5dlk9KkghIzs3JCQiMiZcJkdAcFQ9WiIhIzokITEmUWcnR2pgWVQhIzo3JCQiMUxzKHpiTTpbIiEjOSQhMVRxaV1qLi1lISM6NyQkIjJuImYjUVVGN1wiISM6JCExJWU8JHAob3ZRJyEjOjckJCIyWlYiM2NcUCw6ISM6JCEwbyVvNih5PnYmISM5NyQkIjJEJnBMKVtBOl4iISM6JCExJ1wpb0JdM0xSISM6NyQkIjI5cmtXRCdmOzohIzokITF4WXVHc2hpRSEjOjckJCIyLlojZj8rbkA6ISM6JCExVjokM3RpaEAiISM6NyQkIjIiSC1zJ3lWbl8iISM6JCIxIkhWcSFSMSNbJCEjOzckJCIxKSl6JUdiPD1gIiEjOSQiMml3KSpHI3pkbj4hIzs3JCQiMiRHL20jPi9uYCIhIzokIjJgLCNvPj1QPU4hIzs3JCQiMihvR1pLM2ZUOiEjOiQiMWNpJlJySz0rJiEjOjckJCIxNGBHc3VaWTohIzkkIjFBbihlc0RFTychIzo3JCQiMiRceDQ3VE9eOiEjOiQiMkQmM1c+eTNddiEjOzckJCIyKkhFcyJSUDZjIiEjOiQiMSRHTV8sX3FCKiEjOjckJCIyMF5aOG41NGQiISM6JCIwcFwjPlFOPCkqISM5NyQkIjJ1dltqJ3ohZmQiISM6JCIxJVskPj5zSHMmKiEjOjckJCIyVitdOEUwNGUiISM6JCIwKXkheTE4InApKSEjOTckJCIxXjdOY0QhZmUiISM5JCIxdUNUKHBbcXQoISM6NyQkIjEpXF84JikqKjNmIiEjOSQiMSdRVHg1SkNBJyEjOjckJCIyW3VgajkoKmVmIiEjOiQiMVh0ImUvQW9RJSEjOjckJCIyOypcTlRXKjNnIiEjOiQiMCVbS3pqYy9CISM5NyQkIjIoUWlOTzwqZWciISM6JCIyVy5FRm5SSilmISM9NyQkIjJhW2Q4LiopM2giISM6JCEyTndoJGZddWNBISM7NyQkIjFVSUJZXCVlaCIhIzkkITFLXXZGWStMWCEjOjckJCIyKClmMzYnMyEzaSIhIzokITFvPkRwLWYocCchIzo3JCQiMmQ6JSlmeGNkaSIhIzokITAqZmgqZi5YbSkhIzk3JCQiMkJyZjNwNzJqIiEjOiQhMXYzYiEqKVtiLiIhIzk3JCQiMUUzaD9YaVM7ISM5JCEyc2FkLFJMYEUiISM6NyQkIjIjUj5PXWpgXTshIzokITFSdXVgP11BOCEjOTckJCIxKFIrREhXJ2Y7ISM5JCEyMyUzQyJIP1JAIiEjOjckJCIyVyYpUVlCXyhvOyEjOiQhMXUjKkdKNHhXJyohIzo3JCQiMSQzM2Q/MUxuIiEjOSQhMVtITidlIzRqeiEjOjckJCIyPUp4bjxneW4iISM6JCExSE4vb2ZvaGchIzo3JCQiMjJhWXk5OUNvIiEjOiQhMU8hNFE7RkQrJSEjOjckJCIyJ3BkIio9Im9wbyIhIzokITFlNEk2QiJHJj0hIzo3JCQiMnNad09bKD0jcCIhIzokIjFjIzQmKXpAekonISM7NyQkIjEmPVAlW29TKHAiISM5JCIyY0F5UnZObS4kISM7NyQkIjJEKnk+OGlpLTwhIzokIjFsanBxZHBlXyEjOjckJCIyMGdlemRYeXEiISM6JCIxLj5hXHcpRz8oISM6NyQkIjIlMyQ+RiVcMTg8ISM6JCIwcm1PRiU+J3kpISM5NyQkIjE7K1syVkc9PCEjOSQiMUU5eD04RlQqKiEjOjckJCIyUHNTQW4uTnMiISM6JCIyc3ZQSz5QPjEiISM6NyQkIjI5VixxLkIoRzwhIzokIjIiUShIW1BZI3o1ISM6NyQkIjJPUUU2Y2YhUTwhIzokIjEnKUdKPnMrUCkqISM6NyQkIjJhTF5fMydSWjwhIzokIjF3XiE+V2MpKVIoISM6NyQkIjI6IlFKWlYxXzwhIzokIjFFdSlmN25XcyYhIzo3JCQiMnRHdyQ0RXRjPCEjOiQiMkNqSm9MeUAjUSEjOzckJCIyTXdROSgzU2g8ISM6JCIyRGxfbVpUImY8ISM7NyQkIjImUjddTCJwZ3ciISM6JCExVUkpPTdMIT5SISM7NyQkIjJbX0hYWlk1eCIhIzokITIleTgqXC5gcHAjISM7NyQkIjAieWI6US13PCEjOCQhMSNwQEddek4jXCEjOjckJCIyYTQnZWM2KyJ5IiEjOiQhMTg3TjpfZSMpcCEjOjckJCIyMlE5d1x5ZnkiISM6JCExOipvQSZcUiJ6KSEjOjckJCIxbUVrUWUmNHoiISM5JCEyIlIpRzoiPXRGNSEjOjckJCIyOCY0bnpKJGZ6IiEjOiQhMjlCRmAwVSFRNiEjOjckJCIybUIqcD8wIjQhPSEjOiQhMj0zejZDKmUwNyEjOjckJCIxQXZzaHkpZSE9ISM5JCEyJWVhI2YlPWVGNyEjOjckJCIyKHBXa2lkVDo9ISM6JCEydUtXWCFlYFM2ISM6NyQkIjJ1VGhObVZcIz0hIzokITEvITNNUypmZyopISM6NyQkIjEiKik+U2gyKEg9ISM5JCEybGl2ZiU9OVpzISM7NyQkIjFsJHlXY3JXJD0hIzkkITFiNCU+Zz5RRyYhIzo3JCQiMVJvJFxeTiNSPSEjOSQhMSJbXnQ4SzE5JCEjOjckJCIxOGBSbCUqKlIlPSEjOSQhMT1BMCYzMDQlKikhIzs3JCQiMidbJ3AyKSo9Ilw9ISM6JCIyZXFFKT5BaVY6ISM7NyQkIjJaKVI5J1xRVSY9ISM6JCIxYyEzVzRwdSFSISM6NyQkIjIzSz06LGUkZj0hIzokIjEkKnAhPSswKzUnISM6NyQkIjJtbCMqb194Vyc9ISM6JCIxJyp6SCxQJTQuKSEjOjckJCIyQipwRVVxZnA9ISM6JCIxbiczW3hJNGkqISM6NyQkIjFHOGtkbHJ1PSEjOSQiMWBrJlwoZlshMyIhIzk3JCQiMlFtOkkyTyl6PSEjOiQiMkNOeTtHZk06IiEjOjckJCItUillYlwpPSEjNSQiMWZoNFhzNHk2ISM5LSUlVklFV0c2JDskIiIhISIiJCIrI2ZiXCk9ISIpJShERUZBVUxURy0lJkNPTE9SRzYmJSRSR0JHJCIiISEiIiQiIiEhIiIkIiIhISIiLSUqQVhFU1NUWUxFRzYjJSdOT1JNQUxHLSUoU0NBTElOR0c2IyUuVU5DT05TVFJBSU5FREctJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQhUiEiIi0lKUJPVU5EU19ZRzYjJCIjISohIiItJS1CT1VORFNfV0lEVEhHNiMkIiVnTiEiIi0lLkJPVU5EU19IRUlHSFRHNiMkIiU1USEiIi0lKUNISUxEUkVORzYi