TUZOV3RLVWI8b2I8Uj1NRExDZE5WWlpKOnROPkg6eFhWRXJwczo7Qk5TRE9FVGxNWGxnd2dpVzttRFtVVVVXVXNLaXRVZl1XZnZfaXZtaXhvWUtFVmNzSXl1eXZheXZVSXZfaW9peG9PV2tneHdpeXdPdmVDSHdnSXhpSXhteXFBWXNdSXdnWXRVaXVJWHBDSUZpU0lhQkFBc2E7R2JZeXZjaXhxeXhlWXdleXVZeXVXZE1XVHVVWXV5eXl5QTs6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6OjpaamlmRHFFdGtdYE5cXEBOZFxcUWdxeEhgandoU1dEUVZ5UHhQTEFJWFVgd3l5eVNVdW5gcltETlpdV21ValB1SlpdWV9sU0xxcVdpb1N4d3d0TEVRbEBVTkdpT0NAWHlRalhMWWJJdk48eHdhTG5BdD11T1pkUW5BdEU8U0lkUW5RSkxZUklkcTpgeEpZcnlxSkJoeU5Gdkw/Xl5Zb09BW3lZZWxvZmlHYnQ/d1t3W1BoZEs/Z1NPXkRHcExZZUpwXXQ/ZmpIb1xcSV86eW87SF1cXGBcXDpHb0RGXWBocUVodD13W0ZfYWxTPXdVVG9UdE9IUHdDYm9yWVt3Oj1FcFlkUllyWU1DaEtkRT9CRG1pZEtHPVFzQ19ZUm1IblFCTFlyP1FlRV9YX2tyaWdlOltpQlljZl9ERGFHZVNzXFxlVFBPYl93WXJ3c1hpcmRJdmlHYk53R107VFllVEttZ3l3dkpHQnN5Q3ldVmxtRmV5RVF3Y1g9amp5eDpgc1FNUF5cXFlQaG9fVGs+eE1zbXRzSVBNaEttWUx3TVh3SVdYcU14cVBJVWtFUVQ/bW9EaHRIRW9fbFlAbVFIUXBaRHlMVXJZSHBuPHlSdXRuSFV2PGxweEtZUFd3SVhSX3BgSWBwWGZXT3l5PmVNeV9KV3U9cWFSPnBwVnhPXmZ1bnI/R2BIdl5RaWFddnV1b2NKcHdVUWRUZ2RgX21leF1UdmZcXHhmcmhkYlh2cEllX0hzW0lpSD5uVW9udkBiS3BpWkh0WGBpYmhmS09gSkZkUFBrSXF2eV5xPD9tQHZ1dkFba2BmRGhia1lkTnF4al9jPl9mT2Z2X3dkeF5fRT91WVh5UUBvbEZxWUlmO19lXUl5UFZxbm9zZlB5SkFePWFzdXFbamBaUj9rRV55akh0SFFnT0h4U25cXHdZb0loYFRZXFxGZ2BSeGBJcVtWd3E6QF1UeXliUXh2QF1rPmtpdmRhWVxcdWlcXGRXaXJuW1BxclRncFBZYnhedHZGZmtXWmJpaGxZYT5eYktAd1R3c1Fodk95YkA/XWdxaHdvbW5nXz5vZz0+d3BHYXJBY11oaWJBeVhAZUxvZ1FuaGx5a0Q/czxfY1xcPmJAUXV2QV5reG1ecHBBWHZqVlpzRl5BRm9eX25WVmZsaXhyaWZoYXFpP2JISVxcSmZfXU9dc15gbHlzc0FzcF9iPUlaXWFrZFBtSm5pQXZeUG5hTkF3OkdpPlZxbWZ2Ukl1eUZfW05tcFFqYz9wSXFeUFdqaUZkVVlyYz5nbFBxaFBbQj9qTE5xS0F3eXhuVmhxXFxhallRXlpGVlF4az9lO19mQFViSVNzPz9UPGFCdz1mS11VeVl5W29SQU15Uj1Id2l3RVVIZm1SUFN0eV1Uc1N0YkFIeFN1WU1zXnlHS1V1PUlCPVF4ZW1VQT1ycndJO2FJWD1CSj9iXnNzW19UWEVZVENlRWt1R2dDTmdlRUtZOnl4RUtCTFdibXVCSGt2ak9ndmFjSV9XPV9kR2t0UmVnWXdyXVdGUT95VEtCQlV3SVtIVFlyQnlHanlGXFxXYndndnddU3hhd2FhV3M7eUF3VENBU155eGQ/WGQ9c0JneVJhRFc9RGpzVDo9aFxcS2dtTUdbYXZcXEtkXXNUSkVjdltkVjtmdmNoO3dTOl9Ea1l1XVF3T0NkTztzZz15b2V5dFNHYGtJbXNGeW9nXj94RU9CTENGVmlESWdJQEd5XW90XmlyUDtISz9oWk9zamdTXFxvSD9FVVN1REdNVUF1RkpJSGlfRktTV1V3UlRbaG89U3VjY147SXNfVlRVRT1JQ29TSXN3Q1dxUlpRRzxfaVVhY3JDZWhPY2FJUld1c3BxUmZZVEBjY2ZNdWhzeUNXclltSVBLSWJRaGRDZWhxeFxcc3Q/XURHXUVxTUlGWWZXXXJnVUNicXZJR1Nnb2ZMV2dgYUhKS2RsdUVxRWV1PWl4a3dRU3RyU1d0V2djZ3dKU0lHa3Veb3hnS1Z5UVdaRXReZ0JlS0daS3hjZWQ9SWRUT2hKRWZSW3hyTUJrS2debUdKXUhjW3RyT1Q6X1I/ZUZkX0ZWQ1haQ0Q/UUNxU1hdWWV0R0Y8RXVRZVVmY0NMTWhqR3ZWS3NfU1RrVXdeXUNFVUVsW2Y8P2hORXdkb0g/TVdmXUZiZXNQS1Vda2dIXWJTRVM7UVZWP2hIcWRUPWNlP2JwX2g9R0d1cUdEW3lAU1U9SUZLRXhFZVVXQWhOTVhed2RSWUZJTWV2S2VIWVdTc0NsW0hHYXVbQUVaaVJfaVRKVURTXVlja1hzb1Y+XUdCcWJAO1ZNPURsdVZIZ1Z1UWVEcXhMVUVdXVdTQVJfb0I/b3hMZ3I9PXZxa1I/TWNQQUVHXVdCS1ZQW0hWT0k+SXJFdUJrVWNxU2NrQ3dwc0ZvX1JjP2VCW2hoQ1hZU3JGQ2hWQVN0W1VVV1dzXWNlWUJoeUQ+YVVUTVdaO3ZEb1I8TWlnUUR1X1R0Q1V1UWVUcVhMT0k+UVZfQ2lJXXdfQ3J1RUhvc1J3b0ZmP0VjUWlKP2JoPHJUdVhbWG0+UU4/WXROZHBQUU1TeFVNPExxPXFASU5CQUtFVEpCaHhTdExzRXE6XFxWbVlNY0VKdkxNYFxcam9BV0tsdkxgb1RFeGJxUmBgdVJxSzs9UFg8TEF1c0NoTz9AbU5FamVhUF1JU1docEB5V2xcXFdjPXk8UWxQWEpRdVN3bFc9eHRZeXZKSE9UdEs7VFc+bE9JRE9EVEpaeU5vVVBSTHhIbHdQZWxLeFQ7dG9SRXZcXEFsY11rYnBwZmB5b2x5dlB2T01reERLPl11XFxFVkNdTkFhbkFZYz1GXktfdWRnZFtRXlZpXkRyPlt0Uj1IXFxhR18/R1RbcnRTcnVbWEJ1R0RzVUthZz9RVUVHRUtDaWdjR01lWW9HQjxVUkJJYltlYnZZRkFLYkd5R0s9Q01RQ1FdQ15bVWtVVEZjWFZFaD1dZzxbVkRvQkFJZ095WENnc1FzZGBDRmM9dWpRSEtdWWM7eE9PaUBZeGxPRlhFYnhPR2VDczxraGFJUlZJZ09tcz1lVE9JeVB5cmZ5QnRxdFZ1eVJFeTpvclBjZVxcSWdxa2JWTVVaQVg+c0hzVXVPa1lxQWdDPXN5b1lBc3ZgS0NoQVheV1JdX3h2Y0Zfa1JnQWNeSWNQW1NJX0RbVWY9TWVyb2ZjUUdvWUJmY2NhW1RpRVR2YWU7PUhHY3RhcXZ1V0hkO0liQWlKRVlkeXRHP2hCb3JkVVRYV0M7ZWJjaXNMW1V4WURDW2drUXJITWdIcWVibXJpa3ZqP0hyU2lQeXJja3hrQ3dRQURDb0llRXZiVWJHQWJvUWhYRWhbO2RcXEtISG1iX09GdFd1X3liX1VjUk90bndiUVVIakV1TD1VcF9SYj1VWUFnVU1FPmdDQWdDaXlTRW9zRVFVR3FkV01XcT1jP0VyS01XSUdGd09DZUd3Xz9jQFlCTT1zYFFlcz9VYEd2RElHdV1EaF9VXFxhRUNRQ2tpZ2BLWV5JZDxVckZTR2RpZENRZD93dmpzZ2pvY2Bhdj9BQlVjQ3FrRGJnVVFtWWRXVHlVSElFSV4/dk89eHJJWUBfSVhLeHlleTpXeV1ZUnJ1eERpR2lTdj91aUhHYlZRbDpES1tIbXJwUEhQcV5QbEVcXGtBTWt2bUxMeWxGQW9rbGpjZXY9bHFpPFlXdFJMZXdxSVFQXW51VGpTcXZvXXhndHI6YFREZW9zYHFzWG9VdHNePFFWQUtkPWxIRXdSUU9WWXFFeVRvYFk6YVlOUEtoXFxWPkFzTlF4O1R4cmRXXllKXnRqYVxcdkhkbmxVa1Jla29ZSnZYT1Zlc09xbFVNTkBtUG5VUG9YbVRAanRtVWRwS29IeGhtdUQ9UUhld2tcXG5CbE9odXFXWG93eXM8XFxWUGRrWkFKZ0VSb2BSQGV2W2V2cFRxYGFTeDxOVUF2eVVyYV1udnRSaUh1QlRRSVRzXFxxVlxceUxsXXJhWHRcXEBQQ0hTQHRyO1xcWG1EUztYbUZwVlJ5WHVIaklNS0JgbVppdlc9TkhUU0FcXHNycHRnTW1JQU5xZXVZQHFKTU9GaHJ4RUxpXXZvbXJQQGtnXVRFTVNORVhyZWxteXJveGtVXFxZeU1sbT1LYEF2dmFYV21LUXFtQTxRVElVXklRaG13XklZSFFxXlxcc1dsbEVdc2tscz1RWXdBUnRwVVBIVldpbj5US3llcWBETGtZUEQ9Vll4T1VpdTxRVG89dT1QVWNZWEx5a0VNbUJIWXd1T1Nkc211dV9kUm1dV2xwTElcXHhLbHF5QEteQU86SUpcXGFvO1lzZEh4UkhwT0B5REBMP0lwTHhyZFVwX0h2Y3B2Z111RVFWS1h3dmRucGBWTnFWVkB0W2xMX2lvO3FPSVVOd0xTZmlKOm10PnlWVHFOZ01Wb2FvUF1STmlWRFFPYEBWRGlzZEh5d3RZO0BWZnRMcVlsc3RyRTx2aG1yQmltVU1yPkVrSkF1R3h2WWlZSm1XeHhZcWRqR3hLbF1UQFFQVllPWWBMSmBtYGFqTlxcTUJNVlFteXNMTkRZc3E8b3BZeVNEbWBBdlBAcUJIbFBpT1xcQXg8cW9cXEBZZVhyU0hQUkBWZVlWR0FTcnhRWllQR3R0c1BrXWVMRWhXbzxQR0FQP1F4WkxYWDx1Y01TXFxsSnlkU1NNd0dAa1FMdmpBTVdUeVV0b3hVTGtVUFhUdTxQUVxcaHNhUGtkUEtOaHVIZGtBdHVDUVBaUU1LRFN2UVZQWXlwTFJUeHlNVFBWTVVVaHFzbW1EcG5jWWxYPU5xbHFreFJwZFBPZWtSeE1wQGtTbFVdSFc7eHQ/PVNfTG09QXRuYExxVVFFRVZCQVdCVW5uPXRCTFh4cHRGYE5TQVRkVU5HSFRFPFdOSU5QeFdOSVJHZXdKVE53SHU9TVlWQHVFQEs8TU1fZVNHRWtfRFBAYVZAbWxAPUxcXEV1Q1B2Y3l3U3BrYUB1W3RRaERwPGVTQEF2bUBVXFxNdj1BUVpEVztNTXdVa1E9bT1hVE1BWV1Mb3V1UzpNTlxceVBzXFxRWG1QVkl3dnFKb1BNVElwckFzXlFSdmxzU110TmRLQ0VsP3hUd21uW3Byal1XTXhLT0ROSUlMXkBzbjxZZkVrWEhRTmRRTXRXTGFQTVFMcVVUP015WnRXVGFSQ1Vsa0BYO0BtV2RLP21uVkRsRkB4dnRMVlFRc0ltXlRXczxvWD1LYUV1RWhZRUx0YF1RckB5VEhYUnR4Qlh1dkRyWmR0Xk1uVkhYSWFSeHFMS0xTR0FTTUh3XUBqZHlyTlRNPT1zZV1yYF1vRzxLXFw9VlBgWUtEalh1VHVJakVAd0NRU3hkTV5Ad1BQUz1Nc2I9az5MT1tcXG87dHNtXVVDSUxkRVZqPFM7UFRFaU5VTVZNbW9QdUpHTFRIVE5HcFhLUEtnRHBKbFVUSEtzdW88UGNKbl9jeHBARndhZ1pOWV9XcGVNYHFBV2ddaFxcZklzQWBiWl5hdGl0dXdgPkFpYXloW1ByRVFpZ3BiTU93eXZhSnZ4OkhnYllnO1htPE9yTW9nZFB3XFw+Xj8/a05GYVZYcUdQXmR5WndGcldHeEtuW2tmZ0w+YEdZblBZa2R3YktxYllYcHBoaE9Hcz55W1tGc3VWXUF2Okd1VktHWD9ydG1iVT1VQXlCWFFWSU93RHFFS0lFc29lOmFkXWtYSmF2UkdkZFtCd01jdUVZXFxlRURleF1bZHhPZV1BdVJJZEJLdlNbRD5DWGdBVkhdQndVVUdzQllpeEJ5ZlZ3dnJrU2FdQm1HV2djZnFfaWw7aGdpZzxBUnVhckh1aE5RZHFrSGtLV3FBZHBFY0dvR0NNSV5pd2FXY1d5RlNtU2xxc0lfV21nR2NxZVZpc21xYm9XaFlXR2RXU0ZHQ253dmtPSD9rRlBDc2pvaGFhRG95Z0I7Y2VzUlI9SFBBY0FbaXZhck5nckdReWY7Y2VbY1NhYj5vRjxlWEpdYkE/cmdJVUpLZXRtc2NxQ3FFVUxtVTxNZmVRd3J5dEtzWHRzdmd3ZlI9c3NlZV1lVnI9aVhnVHVzUkVLSXJhd1ZhZXNvRFtdeXFBZ1Q9VkU7WHdnVFg7Q3NpY21zaFpJeWhddGY7RD5hZ1lrU1o7Y0VlZU9VYnBnQktdQ15jY2g9WFtXY3c/VEtdRF9rd1pjcl5hSU1haXZdZ0M/dVJbaUFlR249Z1dJUkRZeGFBZ0x5aVVbY29zZlRNc2BXSEVDZVdpaFpTaVtTU0g/dkRdQ0NTeUpTZUxhd3A9ZEJfeE9PeEtxZ0BjeV4/YmJfd0pTQ2BVWVRbYk5FR2RdSFppSU9BYlZPckplRWNLSWo/Rm5hV2dBZWxTQ0FfV0tddFxcQVVBcVNNTURoU1ZSZ2VwP0JGTXVda1JEa3ZPS3dgXWl5XVNzQXhuVXQ+U2ZpO0lXdXV1R2V2PXJoV2RAcWdnWUZIdWNyPWZmb3ZMUUJfW2U+PWllQXdqQWRcXGloRFV0T1NIRmVXXndpQWVXbk1DY3NYc21oXW94Qz9nXz9mWWdybnNTdmd1dU1iW0lmP0NmYG1SUEV2WXVzbDtYUT1TcW1EUkFzdWN0UEd2c2d4SGdCS3VpXFw7RnZdR0xBRF5bQ29NRT4/U09XVWlnVEtVWVRZRFRjSVZDR3Q7YEZPXVVma1dWXkd3a2FXYHVma0tQWnQ/XFxEd15vT19cXF5pd0Bibl9tdk5eXT92S3FjYHZgdFZzPEB5Zk5fZWBcXER2c0ZnYGlZY0xIdHdfYkt3WmteaEtWbmhWZm1nXFxzSXdRUGZBR2BQbmh1aWtqd3c9YXRYUGJkd2VOQGpOV3VHYGFnaWJeVmlZZmVhUGlyP2c7X1s9T3I6eFxcdnFpPWlxWmlaVHZrQGFmWHZwZHFsU1dcXGlweWRhcllHXFxKRm09WV94Z2U/X1pCcWNESF9Od3JOWHhrUG49V1tsWW49YV89SFpbaW4/UG5BYXlUP3l4dmlpeGJoWHI7aWdkSGlvb1FdZURXZmJpRkBTeD1jdHZ3dVtRV0lpSV5zRV1XclBbRWJNdWdXUm5zQmRbeVJbZD9ZVT1NeWJjaGFkcVI9bG9NdndEU1dZUFd4UlxcPUxPZFdqPExJbXVMXXc8WE5BRVJPaHdLSEx4bUtmdFN4WVVmPEtFTXh4dExObFlLeHI+dVF0cFRWWVV0XFxtdWRUYVxcdF1lblpEUXBgWWtxU0lAb11xT3BgbG9dUEN4eHhBa3J4eF9Ec1RJbVFIU2F5VGhhUkVcXEpMdHB0cVJBdVF5WE9OZFV5VU47YXhbTVdzeHhOPXQ/PUpfeVU/UEs/UXg6SG5fWU5rPFNMSXlcXFVVRFlPaGh2Ok11Zm1Na0lPclVKcVhzbERyVUB3WXFRRWhMS0FrW1hua3F2X114Y1VLbD1UYUVuVUVOUG12SGlwVnBTSFl0Y1BRS3hUSVV3YFhYZUFxW2lYeWhxWGRXXXF1U2xrdVVOVFF3W2h5V2FsUkVTUWhLdWRNUHVSTkhRVlRwWjxWQUhuc0hsWXh2Xm1tVkhKdERZb1R3XWF0XFxYWF5ZeTxMakpga1ZAb110Uz1YSnFNUV1FcHFgc1ZJSl5wV3RQS2NMbWtkTklodVo9V1l5dmhUTnJwVFBpTXdUUWFxcVpFUm5YamU9bklgVXg9cFB1dHRdWF5Ad0dVS0x4als9cUpRblldVk1MVUZNcnFgd0lEc1hoUGNAUVxcQXlydFR1SE49SWtmZFVlcU9bbExxUXBIdHY9cUo+cU5ObVJVVUxeXXdSVExEQHFdYHNXaHhrVEs8aE5eZHJ4UW1rQVJGUXJWQXBKWFlybExXaFZLcXNgcHhSTFR3dVFqZHFmPHdSPWxPbXVfXU5cXGFYQ2lzZ21KdERYc2RVdFVMPFBrYFhwXjxVUlB1RVxcVHlAVUdES1hNS2xRTTpwdDxgbkh5dmc8a0V5bl5sVkZJUT5xUG51d0JleXJ1VG1IWW1YSl91ckRwS3FJUnBQTEx4UlZdSnRMU2t1anhtb2tFbHhNdXhBWE5ZV2NoUD1oUlh4VWpwdnFxbkdFbnY9WWpRclpUUF5lcFRxamFcXFBGRFludWNTUGRpZmxWZ1tqd1tId1xcal5fb3d1VlBkZmdfQ1ZnZFhuSGhoa1F3TVZzaGdaVHhvbFliaGBvakhxd19gZVhaXFw+d1hPbmVgbT9nb0xfd09uYEJ3YGFfdmZ5eVhHdUpHdWdmc29gbWdpdnRIbVheY1NwbVFhZlxcXl1ueWg9b1p4X3dQWG5PaXRyaWI6WHdZT2hwV3ldcWRsV3Z1YGNYPl1KZ3NtPnRkcXNzbl9GP2FuZnlOaFpLZ2dcXE5nRHlwO0FoOl9saEFzW3Z0REZcXE13aFxcZ3dCQWxbeWJNWF4/V2VdWWRFblp3aHk6UXc+YXV0QF9sT2w6PmhmZ2FveHVGUWJLbmJZSHBIUW9id15DX25XQHFEbnBjUXFFR2F3VlxcYEBybnBjbGhjaz5eWEdkTkBxZEF1W0ZmVUledV9cXDpgcWZ2cT9fc29HXFxVZ3VBQVxcQW5da1BsRk5kQkBzS1ZwZE50SF5nQWZvaXBkYUdkRUdsUHdiSlB0W09zUW5eVU5gbUZaWnZsbm9iPnlnTF53V1ltXFxWaGVWZU1HalBockpIZW5JYnBAeFxcd2VdWG9jYGBocGVgeHA6dnVYd2VNWWdbUHFUcG5pSGBvb1tKZ110P3NpQGBwdm9mSXRzblxcXklkYG92VmFnQXFsYUl4VkBdalZdZHZhUUZhbF9oYm93QU94RGBfYVlqSmhsb3FrV1lsSl5mQWZiaT5sTVBgUU5mW2dyWD5yQF9uSF9qX2FeVE52b3hpSlZyc15ldVBjb0BcXFFPW09fcEU+Z1lQbUBfbW9QXlVRX0JwZkVOY0hgak1uWmlZdG14YFZPZ3h2XFxmT3Fob2RAeW9XQW9ITmteV2JDWWRzT2hyeWdKbmRLdnFWWGJSQF1pPmpBSHlXXl1oP11meGdDSWNObl1Jb15sTndIRmY+QGdZQWtRVmNEQGlCP1xcVUdyVFZeaGZqRGlmZ155dEF5SXZcXFFgbXlWeGB2X0RRWlxcSHh0YF5RcV9zUW1AaGRDbnRUQGM9eGZnQGBVWW9cXFl4eGZwZ1lqSElgZGdnWW9fcWB0aEleV2BhPUdyQmhlVVZvUHdraHh5ZEdaU15ucD95Rl5tR2hodmhdVElePHFod3FeSEZcXHNRcFZHdG9vW0dhYklWXFxmQGZCeXdDP2pPd29HRlxcY0Z5cW5tbU5oZXduOndrZnhvYU9pcGhvOndfXnddR1hpQF54aVFlcUZpT25fZ0Feb1ZwVVluPE54RWdsP0lpZ2lgWlFobEd1V292QV94bmFcXFhOczp5Yl9QcHJYXkdpdjtGaHF4ZzxJdGVeZERGYWpIZlN2b1FZaT9XeFpQZGNJX05HbT1pWl5JdmA+ZFk/cD1xaG1QcD0+XU9gYklRd05nZWxRZD9WYllAaV9PYFxcSWJESWVaZnJtcGJsdmxaZlp5XnN2bnNuSWhtTmhbYXBqVmJtVmZVZlpbQXQ9YGZCZ3ZLZmdXeGtiP2Nmb2pkR3ZyaGlMZnZgWV9DX2RpcGdYd29DWHRzSGxcXG5dTlBabU9geVdcXGVeaFRfeERGbGhgW1BJXlpucFdwbURnWj1fY0dmY2NWdlpubkpZa1ZvZmdeaGxXdz5wYVxcYGxNcGZIUGpDUGo+R2ZuT2BUPmljdl1bZmpAdmt0UHJvbnltUFRnZGJQcD55TmRUbnBkb3JRbVRheTpEb0N4cmBpUD5RTGNoTkBEWERUcnlYeUldakc8dUVoTEBUdUE9bGV5T2ZcXFhyUHRacE9pbVlpcVFnUXJxWHZybFdxdXF0TFB2UXlwXXI/UVlOdXJgdU5XZVJgeHFeSGpGaXBnVVlYRHNBWXhOVHJtQE9KcVBJWW46ZW9GTVhZRXRjTFBScUxHSHdLbG5hVU1wSG9jTXdOPXlaVmZFXl5fSXFcXEZuQ19jVEhobldzVz9vTkBuYlBdXWhydkdicz9vcW5tQl9hW3h2bj5mY0BfRU9pPlhoTmZwdVZhQHhoTkljXm9ybUlxZmZvRl9tZkhjZ3lkZ05dX192OkhbcEBfUl5eS2BdZWFvTWlqUldcXFpPeWBoaEFwZVVwbUJoYEBAbVFgbltmaE1xYlRfb09maEVwaHRBbGA/bG5pZGJAdmh3YnFeeFZHY21vYG5odUphdmtOWkleYDxnZ2hIb0Q/YGhfbUlReUlIcUZWZ1txX1Rfbl1YYl5IXFxTZl5dbnRgd2ZXb19WbnZMSGZTbmJzWXl1aWBBUWdjcV1EPnNLQGZxZltDaGlTX1p2ZmZbSHRCUmxnRmA/c09TU3JxdmRBSUg/eHRtVV11V11HR0RhUmlTRj9NVHBLdm9raGJZRj5RaGhLVU4/RHk9VmZ3dj1FVkdXaXZvZUtfdUxhZ0RdclZDdE0/ZD89cz1HUlFvWWZJWTs9VXFvZXJNRHNnSUxvUz5FWDxtQ2tfR0ZBWVxcaUV2SURnZ0VBcWVocWdkV2hTYVNPVWRBVUJxZXlNTWNvO2dmXXdPX3Nsc0dFYUlsVXRvZWdnP0VhSVRrU3VgUXJ3b3Z4S3lea2VMXXJdR3ZqW1dKc2hPX2lnQURCd2ddYHZTbUxMSXNSSG9DZExsSE1GUHhmcGtMUHJgXWptUE9GUXh2SWtkTG4+eFZReXNEdHNFTXB0eHZLYU5gYHRobFJsYHFJYFhpaW1TZUo6XFxTaUxVPG1SXUhyYl14QW11ZFRrV2VQUWl1UTxsb3Vydl10TXRXVWFqaXlvOnBuVz1QRT1vTERRXnlMPmBvR2BqWHRMcGx4c2BtamB0W0hPRFBua0RPdGV0dEhxVkhLVT1UZ2xOYWRSXnhVS1lSO1huOzxZR2BQdFFZWFRPbFB0U0RRTkVsaU13XmRzdklwZnRQVWlrSGFsPkxWcXhrWHhSU2hVQ2BUaXB2YUFOO3Bzc1BPRVFwdGxMXm1ySz1NS3lUQ2B1Yz1tRWRwUmB1PUFxZXhqckBSa1VQcVhVR0FMbl1sPkh2Y1hLREB5Y0l2O1FXS3RTVW1VX1VSQmBrUmBNQkRsWFRuYkxucT1ZWFVuVlB0dXR5T2BReGBKZm1VYm1RR3hLbFFtVFFsTFVNaT1SWlhtRmVOYGxYVTx1UWF0eERZWHFRTGFWVVhyZHlVS015QWhrTVFxVERqdXFTVHhwSlRLQnF2XlF5bmxLXWRYbEBzWHBvZ3hSXnFPWGR2R3htWVlXSU1uZlxcT1l1dWVMTT9pUm5lUnd5d2NRUnJQbkZMUEc9VlJQT15pS108UnFwTVlscXVBeFdwWXdxS211ZnlmaXE/XmNAX0RObm5Zd0BvXWs+W15pdUJmaFFJeXJOdUBpbmxwW19XZnVeX1lwXUVBYWxfeV52ZV9gYlBeYWdnYT1BZVJ3YUhBalFvZURPeVZeZV1vcVVHYHk/Y2h3XT1OeHhPXFx3VndabmRrP2JWPnBQVmpQWWpEbmdAeGM9cXBRX2NPSFxcUT5cXF9mYEdmbV5vZG5RYD5YWmRXZVxcR21GRnZzUHVhcHA9YGxGVnRoaVtCeWE9TnFxRnhZWWNtcVxccGZzYUd4O0FcXERRbGRQd0JQcU1mXFxVb25zVnVoSGhYXmZjTm9Rb2NlT2toSXN0aV47cXRaR3RPbnhYP152cWRZUHBqb11oUGw8cWNAVmtUaGVHT3dFb2tuQFtOYWl5Xm1CUWpcXElyY09tT25rPG9rXXF0PXFyPVZvSVZfZF9qWHZ1cWlvd1h3W25jeEhfTFlddGdwdW5kUE5lR3Bpaz52W2l2UkZ0X0d1S25xRkl3XFxuW0RwaHRgdmdJa1phdl9RXU5xdmJBa21IcUFxZ1VPcV9IXlpYakU+Y0BPYklYanRub2ZHbz9xcV5nZWs+dXRAaXVgdHA/bFxcSXBPaHdDYXNcXGFvbmZpSXZtP1BkcVFnTU5oVlFqRmdlREFkU0ZmXFxWZURBZ2hxbFJHXFxVRnFNWWVwcF54dmxQUHZeR15lQGRGP15MQG1LSXVyUV5eeXhLRmlLbm1YR2o6SXZSeW12PmxnT2xyP3E8aHRKSV1sd1tecXVMQF9cXEdmdWdwPD94R3dbbF93P2llZ19pam5cXEReXFxlSW5Cd1s8Z1puSFxcb1FfRkFiVm5hRmFhUmBfPkdgckhuYV5za1BbQEZqaGloPmBiQnBiZ3Z2ckFhQnhwPUFxb3FiQHdfal5xO3FhXFxWdUZxW0B4W1ZAWz9BZ0peW2t4ZXJQZW5ocWZXYWZXeUBxZm1xdmN3WmtuZ1thd3FAb25hXkVHXUtncmJZXV9gXU54XmZHWnJPeGVBYFJGaDpAd0ZOamZoWl1OY2dpcT1ZZGRXXFxiXm1AQGVyRmNncVtMUW1SQXBlQVtxQXhBPnJUR3JCSG1KcHhveGF5cXFpYWdWT2BjSW5bPlteWHJASWhLPmI6WXdtXmtvTnhxZ1tYYGc8UXlDX2dNYFt4X2Z5PnRcXEBmYnhiO0FgaF5oZFFpS1h1ZkZuRkhdTGhxY2hqXmZ5T2VDW2Zjb3dZc3ZScUk7a1dbO3RlR1ZeP1RJVXdXW2V0R0RPbWhsVWZbdWJmY3lpXUlNb2dpd0heVWg6b3dBd2ZXYVJwbXRDQ3U9SUZkQWVeT2hzS3J1ZWVaRVNaVUJeY2luTUlmd0VfbWZza3hKc0VuP2ZPYWNAP0dKPXVYa0RJaUI6cVlPV0JVS1I/PXU/cWdOYXRoP3ZkR0ZQP3VgVXl3XWZDY2NVT3R2VUdZY2k/Q3U9Z1hwXXhJY0l5T3JBbXl4WXhJY1VOS0ZwUURZWVI7Q3dCPXNQcWlGP2VGP0M+VVRZZWY+PVJKTVRcXEVCQ0dVRGloaEtVUl12Z0FlblV2WmlDZWliV2tDY1VCb3Fkb29jd3lDQ19yQGlmO1NzcWtzd19EVm9CdXlZO3lyZHFmdVNjaEFZQVVncVN5U2tnWnVXbWtEPVFjdll3WVtia3NySnNFYUtXaVd1TT11X2djOkVnX1lmTVl0bnFIb2dSaXNFaT14cHNYcV1ncklocldGaUVDeVNVOm1SZV94Z0VCc1tDXFx5Um8/YzpzaWZDZ25bWEhlRFVnV0tvaWljRVBZSXFxSHhrV3dBZ1ljaHFnRGpnSTthRnBLRFdRdl9PQlpZdD5Xd2RBVVl1Ukw7Y295cU51VmhdVXFwdGRpaztwWGR5a3hUUUdsVFF0dlhgWVY9cTpAVm1AcmZUcklUbFR5all0eG9Aam5VU1NIU3NdcV91eXFdc0NlV05NdUhtSz5cXGpXPWx3ZE5aZG9CYFdzcG9QRExvSW1CVUphYVhYbVhQYUxTZW9uPVV1YW1RcW1wRVVXSHNMVHl2RVhUZU9paVl5aXF5PXFxXVhIQVJIYE1fSVFjPFlBZW5IeVFUSW55TWxSWUtLeGtmRVh5QW1ZeVFUPXB0aHlbZHV4QW1DQW1BeVl3VW9zZHlBYFBrUWtVdE95eHJgSVJedG1NTHZlUEw7RFlaZW5yeXZ1PFlHdU5yYHhrUHFQeE1ybW1DPVNAXFx2YXlKa014bnV0Y2F2PkBqVkR2YGxxZG1MO2l3ZlRrQmVMcWxVV1FQZlVYUVh5WnVMZ21zT3lRPUhWQHVOeGVRdXVtckBRT015a0RNY0xtcz1PSEBzUUlVPUhRbV1yeHVqbHV5R3FweFl1TnFVRF1ZeGFyO2F2Q2BtSFxcVTt0UDtsV2t1UXNJSm14TUdNTHZNTmFZUFNgbXJZT2hxWUpAU3JlUXhdVGJIT3JIbGNJbT9FeXBJam5obGtFd2dZbkZ0WG9oSlNVc2xMcnc8a3FgWGNYdUhtUzpJUFBNV15cXG5lTE5yRVlURUxraHlXRFc6cHZAUFJOQVNBcW9xYFlMeFV5dFhaTXZyYFRqbG1rVXE9aFFQUVJ2ZHRTRVhUWVV1VUxxQFl4QVQ/eVNMbVVWZE5vUU5WXFxvS1F3b1VYPW15SkB2XnlyZklzSmRPUj1YPj1XWEBVUkhqeUlWSXVYSUxMU2Bua1xcTHNkTEp4VlhwbWpQcnB5SjpZcWdNU1V1ckZYcWVhV0R4cFdgWTo9UnBQTExhbz9UUU1ETmxgUGRZcDxtWD55WFhFTEVgd1FtVD5dUWV0eVpFVUV4a3I9Uj95cmFJWFRJWWp1dWZpbnJpS1xcQHlTcUtFeFheXXhdYHJYPVI6TFZocFBCbFBlWVlcXGRMYlVOVEBQWGVTXmF3Wml3QXVralhXPDxtZWRqd1BgcXdlWE5ya29oX1lvPWhia1ddd2B3RUBlTldibHdraFhcXGheck5vYk1Rcnl3aUtZa1l5YHlXdF1pZmh5bHReaHlBeFxceFxceUhrVFlrQnhza195akFoPWh2W19eeWZ1YV9kXFw/a3R4YXI+am93Yll5WmxRYVBWbWpeaXdQcT9gW0FIXl5RYHBBeWB3ZFpgY01weXZpXjtGcm14cXl3dFBwd1dOYHFxYVtRZGxYeEpfaHdnWkxfckh4a2lHYVpXcVxceHFNSHdgeHZOQWhJX3Bqd3NMQHBNV3J4Xm9ocWlaP29oeXdeeHg8T3RRV1tjcWxyRmNab2FcXE5dU1diS3dxP2h5a3dqPXldYF9mZU5zZl9aW0BpXnh2eFhgaXZed15mWFdkUXZgdD9icVl2ZlZecXBoa2dubV9oT0lsZm5cXGZ5bm1weFA/eWZ2cHVpZXBuaHlmcnV5aT8+cVlvXnJ4dXhBdTtoeGVpYV93eVlmW0l2XWBIaXE/bkhXeER5XklQWnlobT9RYVZvaldhcFBnbW5JX3hgbl95aW1AYGp5cj9PeXlXXFxgeHFBQHVYTmFIV1xcP3BkQXlrdXZlSUBjO0FxbUl5Q3dcXGVBZj9Pb1dhcFBHeElIcHh5alJmanlha05gZ1VGeGFPY1N2eVlBbD9xZWZudTtmdklgeHRvd0hIcnBXaF5pY1lPaXRxeTxwZHlPbVM+eWluYFxceWE9aG50dm9bZ3RtcHlZVmpreWlPP2J0RmZ1YXhlX2V0eWhsWFp5eVo/eHBBeHlcXD5rWUZ2O1BbdG9gYkl2bEBuX2lmTUlhXVBgSXl5YFBjd154anlxRW5rc19iOm9qZElpPGZ4cFZ5cVF0S1lyeXlhYkhzS3ZcXEU/bFlndmpvdXU+ZHdeZXRIeWpYYHM/eVRxWlVZdVJeXlhgWj95a3hvaXQ/cW4/c0xgXztobGo/XV54ZHNoZ3N5cW1Bb053W3deeXlOW0J2XzxBZ09mdGloYFNJYT1oX2ZwYnhAdUlBZEl2bEhWXmY/YmRZZkBod1NfXzp2aF5Jd19Jb3lQb1ZZYkl2XFw9Vl9KPnA7Rm1oWWVKPnhhbl1iWW94b2l0Q1ZbYkJLZXFxcm9lR0JDST1NSGxReWNJd1tRdjs/Y0FheEo9UnhheV9FZlhLWXJ5eTpveHZjZHI9VFZZQ0N1Q3c/WDpJWDtDSXJJaUFvRWhTRXRpV2txRXQ9dz90S3dcXHg8dmpYbml1XnlBdl1BWWNOaWVkUGdqRDo7al5QTmFMTlFFTmpENkJcIlx7XH0=Solving Differential Equations
Laplace TransformsC. Madigan
Nova Scotia Agricultural College, Truro, NS, Canada B2N5E3 E Mail cmadigan@nsac.ca
<Text-field style="Heading 1" layout="Heading 1">Introduction</Text-field>Laplace transforms is one of the techniques for solving differential equations that is covered in an Intoductory Differential Equations course. This worksheet is a quick introduction to the Laplace transform, the Inverse Laplace and various applications for use by students in an introductory Differential Equations course.The second section of the worksheet on further applications shows how to handle solving differential equations with a periodic driving force. Included is a technique for graphically displaying the driving force together with the solution to the differential equation. As far as I could determine there is no simple way to deal directly with this situation in MAPLE. Two problems occur. First to adequately define the periodic function in question and then to determine the Laplace Inverse when it is a periodic function. periodic fcns
<Text-field style="Heading 1" layout="Heading 1">Initialization</Text-field>with(inttrans):with(DEtools):with(plots):with(student):
<Text-field style="Heading 1" layout="Heading 1">Introduction to the Laplace transform</Text-field>By the definition of the Laplace Transform, L{f(t)} = NiMtJSRpbnRHNiQqJi0lJGV4cEc2IywkKiYlInNHIiIiJSJ0R0YtISIiRi0tJSJmRzYjRi5GLS9GLjsiIiElKWluZmluaXR5Rw==, where f(t) is defined for t \342\200\26045 0.
<Text-field style="Heading 2" layout="Heading 2">Example 1. Calculate Laplace transform of f(t) = t</Text-field>f := t -> t;
F:=int(exp(-s*t)*f(t),t=0..infinity) assuming s>0;
laplace(f(t),t,s);
<Text-field style="Heading 2" layout="Heading 2">Example 2. Evaluate L{3t - 5 sin 2t}</Text-field>Solution:By linearity, L{3 t - 5 sin 2t} = 3 L{t} - 5 L{sin 2t} soRHS:=3*int(exp(-s*t)*t,t=0..infinity)-5*int(exp(-s*t)*sin(2*t),t=0..infinity);
eval(%) assuming s > 0;laplace(3*t-5*sin(2*t),t,s);
<Text-field style="Heading 2" layout="Heading 2">Example 3 Transform of a Piecewise Defined Function</Text-field>Evaluate L{f(t)} for f(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkobWZlbmNlZEdGJDYmLUknbXRhYmxlR0YkNjYtSSRtdHJHRiQ2Jy1JJG10ZEdGJDYoLUYjNiMtSSNtbkdGJDYkUSIwRicvRjNRJ25vcm1hbEYnLyUpcm93YWxpZ25HRi4vJSxjb2x1bW5hbGlnbkdGLi8lK2dyb3VwYWxpZ25HRi4vJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGUS1GPzYoLUYjNiotSSNtb0dGJDYxUSJ+RicvJStmb3JlZ3JvdW5kR1EsWzIwMCwwLDIwMF1GJ0ZHLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0Zbby8lKXN0cmV0Y2h5R0Zbby8lKnN5bW1ldHJpY0dGW28vJShsYXJnZW9wR0Zbby8lLm1vdmFibGVsaW1pdHNHRltvLyUnYWNjZW50R0Zbby8lJWZvcm1HRi4vJSdsc3BhY2VHUSQwZW1GJy8lJ3JzcGFjZUdGXHAvJShtaW5zaXplR1EiMUYnLyUobWF4c2l6ZUdRKWluZmluaXR5RidGWEZYLUZENiVGRkZmbkZHLUZZNjFRIjxGJ0ZmbkZHRmluRlxvRl5vRmBvRmJvRmRvRmZvL0Zpb1EmaW5maXhGJy9GW3BRL3RoaWNrbWF0aHNwYWNlRicvRl5wRl1xRl9wRmJwLUYsNiZRInRGJ0YvRmZuRjJGZ3AtRkQ2JVEiM0YnRmZuRkdGSUZLRk1GT0ZSRklGS0ZNLUY8NictRj82KC1GIzYjLUZENiRRIjJGJ0ZHRklGS0ZNRk9GUi1GPzYoLUYjNiVGX3EtRlk2MVElJmdlO0YnRmZuRkdGaW5GXG9GXm9GYG9GYm9GZG9GZm9GanBGXHFGXnFGX3BGYnBGYnFGSUZLRk1GT0ZSRklGS0ZNLyUmYWxpZ25HUSVheGlzRicvRkpRKWJhc2VsaW5lRicvRkxRJ2NlbnRlckYnL0ZOUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmJzLyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRl10LyUmZnJhbWVHRl10LyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGW28vJS1lcXVhbGNvbHVtbnNHRltvLyUtZGlzcGxheXN0eWxlR0Zbby8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0Zqc0ZHLyUlb3BlbkdRInxmckYnLyUmY2xvc2VHRi4=Solution:First, define and plot this function using Maple,f:=t->piecewise(t>=0 and t<3,0,t>=3,2);
plot(f(t),t=0..5,discont=true);
Using the definition of the Laplace Transform, since f(t) is piecewise continuous, the Laplace Transform is expressed as the sum of two integrals.L1 := int(0*exp(-s*t),t = 0..3);
L2 := int(2*exp(-s*t),t = 3..infinity);
L:=L1+L2;
Again this limit is 0 provided s > 0. Hence L = NiMqKCIiIyIiIi0lJGV4cEc2IywkKiYiIiRGJSUic0dGJSEiIkYlRixGLQ==laplace(f(t),t,s);
Maple cannot find the laplace transform of a piecewise defined function directly. We can redefine our function in terms of the Heaviside step function and then Maple can determine the transform for us.f1:=convert(f(t),Heaviside);
laplace(f1,t,s);
<Text-field style="Heading 3" layout="Heading 3">Example 3b. Consider the following piecewise continuous function</Text-field>f:=t->piecewise(t<=1,0,t>=1,2*t-3);
plot(f(t),t=0..5,discont=true);
We can rewrite this function is terms of the Heaviside function The unapply command is used to convert the answer into a function of the variable t.f:=unapply(factor(convert(f(t),Heaviside)),t);
Similarly below the unapply command is used to write the laplace transform as a function of sF := unapply(laplace((2*t - 3)*Heaviside(t-1),t,s),s);
laplace(f(t),t,s);
<Text-field style="Heading 2" layout="Heading 2">Example 4. Calculate <Equation executable="false" style="Heading 2" input-equation="Laplace(t^2*sin(kt))">NiMtJShMYXBsYWNlRzYjKiYpJSJ0RyIiIyIiIi0lJHNpbkc2IyUja3RHRio=</Equation></Text-field>We have the rule that NiMvLSUoTGFwbGFjZUc2IyomKSUidEclIm5HIiIiLSUiZkc2I0YpRisqJiksJEYrISIiRipGKy0tJSNAQEc2JCUiREdGKjYjLSUiRkc2IyUic0dGKw==f:=t->sin(k*t);
Laplace(t^2*f(t))=(-1)^2*diff(laplace(f(t),t,s),s,s);
simplify(%);laplace(t^2*f(t),t,s);
<Text-field style="Heading 2" layout="Heading 2">Example 5. Calculate L{kt cos kt + sin kt}</Text-field>This problem can be done three ways.First method:Note that kt cos kt + sin kt is the derivative of t sin kt.Hence, L{kt cos kt + sin kt} = L{d/dt (t sin kt)} = sL{t sin kt} = s(- d/ds L{sin kt}) Second method: L{ d/dt(f(t)) } = s L{ f(t) } - f(0) So define f(t) and apply this formula.Third method : Use the laplace command directly.
<Text-field style="Heading 2" layout="Heading 2">Example 6. Dirac Delta function.</Text-field>The Dirac delta function is found in Maple in the form Dirac[t - LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMEYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYn] where the impulse is at t =LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMEYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYn. The Laplace transform package can work with it.f:='f':t:='t':Dirac(t-t[0]);
laplace(Dirac(t-2),t,s);
laplace(f(t)*Dirac(t-2),t,s);
invlaplace(exp(-3*s),s,t);
invlaplace(1/(s^2-5)*exp(-3*s),s,t);
simplify(convert(%,exp));
<Text-field style="Heading 2" layout="Heading 2">Exercises Use the various methods to calculate the Laplace transforms and use the command laplace to check your answers.</Text-field>1. Use the definition to find the Laplace transform of f(t) = -1 if 0 <= t < 1,1 if t >= 1.2. Use the definition to find the Laplace transform of f(t) = 2t + 1 if 0 <= t < 1,0 if t >= 1.3. Use the definition to find the Laplace transform of NiMvLSUiZkc2IyUidEctJSRleHBHNiMsJkYnIiIiIiIoISIi4. Use the definition to find the Laplace transform of f(t) = NiMqJiklInRHIiIjIiIiLSUkZXhwRzYjKiYiIiRGJ0YlRidGJw== .5. Find the Laplace transform of f(t) = NiMqJiIiIyIiIiokKSUidEciIiVGJUYl.6. Find the Laplace transform of f(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUklbXN1cEdGJDYlLUkobWZlbmNlZEdGJDYkLUYjNiUtSSNtaUdGJDYlUSJ0RicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2MFEiK0YnL0Y7USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGRS8lKXN0cmV0Y2h5R0ZFLyUqc3ltbWV0cmljR0ZFLyUobGFyZ2VvcEdGRS8lLm1vdmFibGVsaW1pdHNHRkUvJSdhY2NlbnRHRkUvJSVmb3JtR1EmaW5maXhGJy8lJ2xzcGFjZUdRMG1lZGl1bW1hdGhzcGFjZUYnLyUncnNwYWNlR0ZXLyUobWluc2l6ZUdRIjFGJy8lKG1heHNpemVHUSlpbmZpbml0eUYnLUkjbW5HRiQ2JEZmbkZBRkEtRiM2Iy1GW282JFEiM0YnRkEvJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYnLUY+NjBRIn5GJ0ZBRkNGRkZIRkpGTEZORlAvRlNRIUYnL0ZWUSQwZW1GJy9GWUZbcEZaRmdu = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYrLUklbXN1cEdGJDYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiM0YnL0Y2USdub3JtYWxGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYwUSIrRidGPi8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGSS8lKXN0cmV0Y2h5R0ZJLyUqc3ltbWV0cmljR0ZJLyUobGFyZ2VvcEdGSS8lLm1vdmFibGVsaW1pdHNHRkkvJSdhY2NlbnRHRkkvJSVmb3JtR1EmaW5maXhGJy8lJ2xzcGFjZUdRMG1lZGl1bW1hdGhzcGFjZUYnLyUncnNwYWNlR0Zlbi8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJ0Y6LUYsNiVGLi1GIzYjLUY7NiRRIjJGJ0Y+RkBGQ0Y6Ri5GQy1GOzYkRmpuRj4=7. Find the Laplace transform of f(t) = cos 5t + sin 2t.8. Evaluate (a) L{LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYkLUYsNiVRIUYnRi9GMi1JJW1zdXBHRiQ2JS1JI21vR0YkNjBRLyZFeHBvbmVudGlhbEU7RicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHUSZ1bnNldEYnLyUqc2VwYXJhdG9yR0ZFLyUpc3RyZXRjaHlHRkUvJSpzeW1tZXRyaWNHRkUvJShsYXJnZW9wR0ZFLyUubW92YWJsZWxpbWl0c0dGRS8lJ2FjY2VudEdGRS8lJWZvcm1HUSdwcmVmaXhGJy8lJ2xzcGFjZUdRJDBlbUYnLyUncnNwYWNlR1EydmVyeXRoaW5tYXRoc3BhY2VGJy8lKG1pbnNpemVHRjkvJShtYXhzaXplR0Y5LUYjNiQtSSNtbkdGJDYkUSIzRidGQUYrLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJ0Y3}; (c) L{t^2}sin(kt)} , and (d) L{t LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JC1GOTYwUSomdW1pbnVzMDtGJ0Y8L0Y/USZmYWxzZUYnL0ZCRmpuL0ZERmpuL0ZGRmpuL0ZIRmpuL0ZKRmpuL0ZMRmpuL0ZOUSZpbmZpeEYnL0ZRUTBtZWRpdW1tYXRoc3BhY2VGJy9GVEZkby9GV1EiMUYnL0ZZUSlpbmZpbml0eUYnLUYsNiVRInRGJ0YvRjIvJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYncos( t )}9. Find the Laplace transform of t LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JC1JI21uR0YkNiRRIzEwRidGPC1GLDYlUSJ0RidGL0YyLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJw==.,LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JS1GOTYwUSomdW1pbnVzMDtGJ0Y8L0Y/USZmYWxzZUYnL0ZCRmpuL0ZERmpuL0ZGRmpuL0ZIRmpuL0ZKRmpuL0ZMRmpuL0ZOUSZpbmZpeEYnL0ZRUTBtZWRpdW1tYXRoc3BhY2VGJy9GVEZkby9GV1EiMUYnL0ZZUSlpbmZpbml0eUYnLUkjbW5HRiQ2JFEiMkYnRjwtRiw2JVEidEYnRi9GMi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRic= cos(4t), t(LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2Iy1GLDYlUSJ0RidGL0YyLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJw== + LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRjcvJSlzdHJldGNoeUdGNy8lKnN5bW1ldHJpY0dGNy8lKGxhcmdlb3BHRjcvJS5tb3ZhYmxlbGltaXRzR0Y3LyUnYWNjZW50R0Y3LyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdRIUYnLyUobWF4c2l6ZUdGTy1GIzYkLUkjbW5HRiQ2JFEiMkYnRjItSSNtaUdGJDYlUSJ0RicvJSdpdGFsaWNHUSV0cnVlRicvRjNRJ2l0YWxpY0YnLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJy1GLDYlLUYvNjBRIilGJ0YyL0Y2RmhuL0Y5USZmYWxzZUYnL0Y7RmhuL0Y9RmVvL0Y/RmVvL0ZBRmVvL0ZDRmVvL0ZFUShwb3N0Zml4RicvRkhRLnRoaW5tYXRoc3BhY2VGJy9GS0ZecC9GTlEiMUYnL0ZRUSlpbmZpbml0eUYnLUYjNiNGVEZbbw==JSMlP0c= .10. Find the Laplace transform of (3t + 1)U(t - 3), (LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1cEdGJDYlLUkobWZlbmNlZEdGJDYkLUYjNiUtSSNtaUdGJDYlUSJ0RicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2MFEqJnVtaW51czA7RicvRjtRJ25vcm1hbEYnLyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZFLyUpc3RyZXRjaHlHRkUvJSpzeW1tZXRyaWNHRkUvJShsYXJnZW9wR0ZFLyUubW92YWJsZWxpbWl0c0dGRS8lJ2FjY2VudEdGRS8lJWZvcm1HUSZpbmZpeEYnLyUnbHNwYWNlR1EwbWVkaXVtbWF0aHNwYWNlRicvJSdyc3BhY2VHRlcvJShtaW5zaXplR1EiMUYnLyUobWF4c2l6ZUdRKWluZmluaXR5RictSSNtbkdGJDYkRmZuRkFGQS1GIzYjLUZbbzYkUSIzRidGQS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRic=LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRjcvJSlzdHJldGNoeUdGNy8lKnN5bW1ldHJpY0dGNy8lKGxhcmdlb3BHRjcvJS5tb3ZhYmxlbGltaXRzR0Y3LyUnYWNjZW50R0Y3LyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdRIUYnLyUobWF4c2l6ZUdGTy1GIzYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnL0YzUSdpdGFsaWNGJy1GLzYwUSomdW1pbnVzMDtGJ0YyL0Y2USZmYWxzZUYnL0Y5RltvL0Y7RltvL0Y9RltvL0Y/RltvL0ZBRltvL0ZDRltvL0ZFUSZpbmZpeEYnL0ZIUTBtZWRpdW1tYXRoc3BhY2VGJy9GS0Zlby9GTlEiMUYnL0ZRUSlpbmZpbml0eUYnLUkjbW5HRiQ2JEZob0YyLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJy1GLzYwUSIpRidGMi9GNkZaRlxvL0Y7RlpGXm9GX29GYG9GYW8vRkVRKHBvc3RmaXhGJy9GSFEudGhpbm1hdGhzcGFjZUYnL0ZLRmlwRmdvRmlvU(t - 1).11. Find the Laplace transform of (t LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JS1GLDYlUSJ0RidGL0YyLUY5NjBRKiZ1bWludXMwO0YnRjwvRj9RJmZhbHNlRicvRkJGXW8vRkRGXW8vRkZGXW8vRkhGXW8vRkpGXW8vRkxGXW8vRk5RJmluZml4RicvRlFRMG1lZGl1bW1hdGhzcGFjZUYnL0ZURmdvL0ZXUSIxRicvRllRKWluZmluaXR5RictSSNtbkdGJDYkUSI1RidGPC8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRic=)U(t - 5).12. a) Write f(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkobWZlbmNlZEdGJDYmLUknbXRhYmxlR0YkNjYtSSRtdHJHRiQ2Jy1JJG10ZEdGJDYoLUYjNiQtSSNtbkdGJDYkUSIyRicvRjNRJ25vcm1hbEYnLUkjbW9HRiQ2MFEifkYnRkcvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRk8vJSlzdHJldGNoeUdGTy8lKnN5bW1ldHJpY0dGTy8lKGxhcmdlb3BHRk8vJS5tb3ZhYmxlbGltaXRzR0ZPLyUnYWNjZW50R0ZPLyUlZm9ybUdGLi8lJ2xzcGFjZUdRJDBlbUYnLyUncnNwYWNlR0Zqbi8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJy8lKXJvd2FsaWduR0YuLyUsY29sdW1uYWxpZ25HRi4vJStncm91cGFsaWduR0YuLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRltwLUY/NigtRiM2Jy1GRDYlUSIwRicvJStmb3JlZ3JvdW5kR1EsWzIwMCwwLDIwMF1GJ0ZHLUZKNjFRJiZsZXE7RidGZXBGR0ZNRlBGUkZURlZGWEZaL0ZnblEmaW5maXhGJy9GaW5RL3RoaWNrbWF0aHNwYWNlRicvRlxvRl5xRl1vRmBvLUYsNiZRInRGJ0YvRmVwRjItRko2MVElJmx0O0YnRmVwRkdGTUZQRlJGVEZWRlhGWkZbcUZdcUZfcUZdb0Zgby1GRDYlUSIzRidGZXBGR0Zjb0Zlb0Znb0Zpb0ZccEZjb0Zlb0Znby1GPDYnLUY/NigtRiM2JC1GSjYwUSomdW1pbnVzMDtGJ0ZHRk1GUEZSRlRGVkZYRlpGW3EvRmluUTBtZWRpdW1tYXRoc3BhY2VGJy9GXG9GY3JGXW9GYG9GQ0Zjb0Zlb0Znb0Zpb0ZccC1GPzYoLUYjNiUtRiw2JUZicUYvRjItRko2MVElJmdlO0YnRmVwRkdGTUZQRlJGVEZWRlhGWkZbcUZdcUZfcUZdb0Zgb0ZmcUZjb0Zlb0Znb0Zpb0ZccEZjb0Zlb0Znby8lJmFsaWduR1ElYXhpc0YnL0Zkb1EpYmFzZWxpbmVGJy9GZm9RJ2NlbnRlckYnL0Zob1EnfGZybGVmdHxockYnLyUvYWxpZ25tZW50c2NvcGVHRjEvJSxjb2x1bW53aWR0aEdRJWF1dG9GJy8lJndpZHRoR0ZbdC8lK3Jvd3NwYWNpbmdHUSYxLjBleEYnLyUuY29sdW1uc3BhY2luZ0dRJjAuOGVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0ZmdC8lJmZyYW1lR0ZmdC8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRk8vJS1lcXVhbGNvbHVtbnNHRk8vJS1kaXNwbGF5c3R5bGVHRk8vJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dGY3RGRy8lJW9wZW5HUSJ8ZnJGJy8lJmNsb3NlR0Yu in terms of unit step functions. Find the Laplace transform of the given function. b) Write f(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2MFEifGZyRicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHRjEvJSpzZXBhcmF0b3JHUSZmYWxzZUYnLyUpc3RyZXRjaHlHRjEvJSpzeW1tZXRyaWNHRj8vJShsYXJnZW9wR0Y/LyUubW92YWJsZWxpbWl0c0dGPy8lJ2FjY2VudEdGPy8lJWZvcm1HUSdwcmVmaXhGJy8lJ2xzcGFjZUdRLnRoaW5tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGTy8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJy1JJ210YWJsZUdGJDY2LUkkbXRyR0YkNictSSRtdGRHRiQ2KC1GIzYjLUkjbW5HRiQ2JFEiMEYnRjkvJSlyb3dhbGlnbkdGLi8lLGNvbHVtbmFsaWduR0YuLyUrZ3JvdXBhbGlnbkdGLi8lKHJvd3NwYW5HUSIxRicvJStjb2x1bW5zcGFuR0Zpby1GaW42KC1GIzYnLUZebzYlRmBvLyUrZm9yZWdyb3VuZEdRLFsyMDAsMCwyMDBdRidGOS1GNjYxUSYmbGVxO0YnRmJwRjkvRjxGP0Y9L0ZBRj9GQkZERkZGSC9GS1EmaW5maXhGJy9GTlEvdGhpY2ttYXRoc3BhY2VGJy9GUUZdcUZSRlUtRiw2JlEidEYnRi9GYnBGMi1GNjYxUSI8RidGYnBGOUZocEY9RmlwRkJGREZGRkhGanBGXHFGXnFGUkZVLUkmbWZyYWNHRiQ2KS1GIzYkLUZebzYlUSIzRidGYnBGOS1GLDYmUSUmcGk7RicvRjBGP0ZicEY5LUYjNiMtRl5vNiVRIjJGJ0ZicEY5LyUubGluZXRoaWNrbmVzc0dGaW8vJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGanIvJSliZXZlbGxlZEdGPy8lK2ZvcmVncm91bmRHRmRwRmFvRmNvRmVvRmdvRmpvRmFvRmNvRmVvLUZmbjYnLUZpbjYoLUYjNiQtRiw2JVEkc2luRidGYHJGOS1JKG1mZW5jZWRHRiQ2JC1GIzYjLUYsNiVGYXFGL0YyRjlGYW9GY29GZW9GZ29Gam8tRmluNigtRiM2JS1GLDYmUSJ4RidGL0ZicEYyLUY2NjFRJSZnZTtGJ0ZicEY5L0Y8USZ1bnNldEYnL0Y+Rlx1L0ZBRlx1L0ZDRlx1L0ZFRlx1L0ZHRlx1L0ZJRlx1RmpwRlxxRl5xL0ZTRi4vRlZGLi1GZnE2KC1GIzYkLUZebzYkRlxyRjktRiw2JUZfckZgckY5LUYjNiMtRl5vNiRGZXJGOUZmckZockZbc0Zdc0Zhb0Zjb0Zlb0Znb0Zqb0Zhb0Zjb0Zlby8lJmFsaWduR1ElYXhpc0YnL0Zib1EpYmFzZWxpbmVGJy9GZG9GanIvRmZvUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRl13LyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EmMC44ZW1GJy8lKXJvd2xpbmVzR1Elbm9uZUYnLyUsY29sdW1ubGluZXNHRmh3LyUmZnJhbWVHRmh3LyUtZnJhbWVzcGFjaW5nR1EsMC40ZW1+MC41ZXhGJy8lKmVxdWFscm93c0dGPy8lLWVxdWFsY29sdW1uc0dGPy8lLWRpc3BsYXlzdHlsZUdGPy8lJXNpZGVHUSZyaWdodEYnLyUwbWlubGFiZWxzcGFjaW5nR0Zldy1GNjYwRi5GOUZocEY9RmlwRkJGREZGRkgvRktGLi9GTlEkMGVtRicvRlFGX3lGUkZV in terms of unit step functions. Find the Laplace transform of the given function.
<Text-field style="Heading 1" layout="Heading 1">Inverse Transform</Text-field>
<Text-field style="Heading 2" layout="Heading 2">Example 1. Termwise Division and Linearity</Text-field>Evaluate L^(-1){(3 s + 5)/(s^2 + 7)}.Solution:The given function of s can be written as two expressions by means of termwise division:F := (3*s + 5)/(s^2 + 7);
F := expand(F);
LI := 3*cos(sqrt(7)*t) + 5/sqrt(7)*sin(sqrt(7)*t);
This can be done directly using the invlaplace which recognizes the linearity of the inverse Laplace transform.ls := invlaplace((3*s + 5)/(s^2 + 7),s,t);
rs := 3*invlaplace(s/(s^2 + 7),s,t) + 5/sqrt(7)*invlaplace(sqrt(7)/(s^2 + 7),s,t);
ls = rs;
<Text-field style="Heading 2" layout="Heading 2">Example 2 Partial Fractions and Linearity</Text-field>Evaluate NiMqJiklIkxHLCQiIiIhIiJGJzwjKiYsJiUic0dGJ0YnRidGJyomKUYsIiIjRicpLCZGLEYnRi9GJyIiJEYnRihGJw==Solution:LetF := s -> (s+1)/(s^2*(s+2)^3);
F1:=convert(F(s),parfrac,s);
T1:=invlaplace(op(1,F1),s,t);
T2:=invlaplace(op(2,F1),s,t);T3:=invlaplace(op(3,F1),s,t);T4:=invlaplace(op(4,F1),s,t);
y:=T1+T2+T3+T4;
Maple can also calculate the inverse of F(s) directlyinvlaplace(F(s),s,t);
<Text-field style="Heading 2" layout="Heading 2">Example 3 Completing the Square to find L^(-1)</Text-field>Evaluate L^(-1){s/(s^2 + 6 s + 11)}.Solution:Let F := (s) -> s/(s^2 + 6*s + 11);
F:=unapply(completesquare(F(s)),s);
F:=s->(s+3)/((s+3)^2+2)-3/((s+3)^2+2);
Applying the !st translation theorem we haveInv(F(s)):=exp(-3*t)*Inv(s/(s^2+2))-3/sqrt(2)*exp(-3*t)*Inv(sqrt(2)/(s^2+2));
f:=t->exp(-3*t)*cos(sqrt(2)*t)-3/2*2^(1/2)*exp(-3*t)*sin(sqrt(2)*t);
f := unapply(invlaplace(F(s),s,t),t);
<Text-field style="Heading 2" layout="Heading 2">Example 4 Convolutions</Text-field>
<Text-field style="Heading 3" layout="Heading 3">Definition</Text-field>Given two functions f(t) and g(t) their convolution f * g is defined as f*g= NiMtJSRpbnRHNiQqJi0lImZHNiMlJHRhdUciIiItJSJnRzYjLCYlInRHRitGKiEiIkYrL0YqOyIiIUYwThe Laplace{ f*g } = L{ f(t) }L{ g(t) } = F(s) G(s)
<Text-field style="Heading 4" layout="Heading 4">example</Text-field>Given f(t) =exp(t) g(t) = sin(t)f:=t->exp(t);g:=t->sin(t);
fstarg:=int(f(tau)*g(t-tau),tau=0..t);
gstarf:=int(g(tau)*f(t-tau),tau=0..t);
laplace(fstarg,t,s)=laplace(f(t),t,s)*laplace(g(t),t,s);
simplify(lhs(%));
<Text-field style="Heading 3" layout="Heading 3">Application</Text-field>Since L { f( * g } = F(s) G(s) then NiMpJSJMRywkIiIiISIi(F(s)G(s)) = f*g
<Text-field style="Heading 4" layout="Heading 4">Example Inverse Transform as a Convolution</Text-field>Evaluate L^(-1){1/((s - 1)(s + 4))}F:=s->1/(s-1); G:=s->1/(s+4);
f := unapply(invlaplace(1/(s-1),s,t),t);
g := unapply(invlaplace(1/(s+4),s,t),t);
IL := int(f(tau)*g(t - tau), tau = 0..t);
invlaplace(F(s)*G(s),s,t);
<Text-field style="Heading 2" layout="Heading 2">Exercises Using various methods calculate the Laplace inverse. Use the invlaplace command to check your answer.</Text-field>2. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkmbWZyYWNHRiQ2KC1JJW1zdXBHRiQ2JS1JKG1mZW5jZWRHRiQ2JC1GIzYlLUkjbWlHRiQ2JVEic0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNjBRIitGJy9GPlEnbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRkgvJSlzdHJldGNoeUdGSC8lKnN5bW1ldHJpY0dGSC8lKGxhcmdlb3BHRkgvJS5tb3ZhYmxlbGltaXRzR0ZILyUnYWNjZW50R0ZILyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUTBtZWRpdW1tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGWi8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJy1JI21uR0YkNiRRIjJGJ0ZERkQtRiM2I0Zdby8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictRiM2Iy1GLzYlRjYtRiM2Iy1GXm82JFEiM0YnRkRGY28vJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmRwLyUpYmV2ZWxsZWRHRkg=.3. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUkjbW5HRiQ2JFEiNUYnL0YzUSdub3JtYWxGJy1GIzYlLUklbXN1cEdGJDYlLUYsNiVRInNGJ0YvRjItRiM2Iy1GOzYkUSIyRidGPi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYwUSIrRidGPi8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGVi8lKXN0cmV0Y2h5R0ZWLyUqc3ltbWV0cmljR0ZWLyUobGFyZ2VvcEdGVi8lLm1vdmFibGVsaW1pdHNHRlYvJSdhY2NlbnRHRlYvJSVmb3JtR1EmaW5maXhGJy8lJ2xzcGFjZUdRMG1lZGl1bW1hdGhzcGFjZUYnLyUncnNwYWNlR0Ziby8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJy1GOzYkUSM0OUYnRj4vJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmNwLyUpYmV2ZWxsZWRHRlZGKw==.4. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbW9HRiQ2MFEifkYnLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGNC8lKXN0cmV0Y2h5R0Y0LyUqc3ltbWV0cmljR0Y0LyUobGFyZ2VvcEdGNC8lLm1vdmFibGVsaW1pdHNHRjQvJSdhY2NlbnRHRjQvJSVmb3JtR1EhRicvJSdsc3BhY2VHUSQwZW1GJy8lJ3JzcGFjZUdGRi8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJ0YrLUkmbWZyYWNHRiQ2KC1GIzYkLUkjbW5HRiQ2JFEjMTBGJ0YvLUkjbWlHRiQ2JVEic0YnLyUnaXRhbGljR1EldHJ1ZUYnL0YwUSdpdGFsaWNGJy1GIzYlLUklbXN1cEdGJDYlRlgtRiM2Iy1GVTYkUSIyRidGLy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictRiw2MFEiK0YnRi9GMkY1RjdGOUY7Rj1GPy9GQlEmaW5maXhGJy9GRVEwbWVkaXVtbWF0aHNwYWNlRicvRkhGXnBGSUZMLUZVNiRRIzE2RidGLy8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGaHAvJSliZXZlbGxlZEdGNA==.5. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUkjbW5HRiQ2JFEiMUYnL0YzUSdub3JtYWxGJy1GIzYlLUklbXN1cEdGJDYlLUYsNiVRInNGJ0YvRjItRiM2Iy1GOzYkUSIyRidGPi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSSNtb0dGJDYwUSomdW1pbnVzMDtGJ0Y+LyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZWLyUpc3RyZXRjaHlHRlYvJSpzeW1tZXRyaWNHRlYvJShsYXJnZW9wR0ZWLyUubW92YWJsZWxpbWl0c0dGVi8lJ2FjY2VudEdGVi8lJWZvcm1HUSZpbmZpeEYnLyUnbHNwYWNlR1EwbWVkaXVtbWF0aHNwYWNlRicvJSdyc3BhY2VHRmJvLyUobWluc2l6ZUdGPS8lKG1heHNpemVHUSlpbmZpbml0eUYnLUY7NiRRIzE2RidGPi8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGYnAvJSliZXZlbGxlZEdGVkYr .6. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUYsNiVRInNGJ0YvRjItRiM2KC1JJW1zdXBHRiQ2JUY6LUYjNiMtSSNtbkdGJDYkUSIyRicvRjNRJ25vcm1hbEYnLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJy1JI21vR0YkNjBRIitGJ0ZILyUmZmVuY2VHUSZmYWxzZUYnLyUqc2VwYXJhdG9yR0ZTLyUpc3RyZXRjaHlHRlMvJSpzeW1tZXRyaWNHRlMvJShsYXJnZW9wR0ZTLyUubW92YWJsZWxpbWl0c0dGUy8lJ2FjY2VudEdGUy8lJWZvcm1HUSZpbmZpeEYnLyUnbHNwYWNlR1EwbWVkaXVtbWF0aHNwYWNlRicvJSdyc3BhY2VHRl9vLyUobWluc2l6ZUdRIjFGJy8lKG1heHNpemVHUSlpbmZpbml0eUYnRkRGOi1GTjYwUSomdW1pbnVzMDtGJ0ZIRlFGVEZWRlhGWkZmbkZobkZqbkZdb0Zgb0Zib0Zlby1GRTYkUSIzRidGSC8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGY3AvJSliZXZlbGxlZEdGU0Yr.7. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUYsNiVRInNGJ0YvRjItRiM2JS1JKG1mZW5jZWRHRiQ2JC1GIzYlRjotSSNtb0dGJDYwUSomdW1pbnVzMDtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRkwvJSlzdHJldGNoeUdGTC8lKnN5bW1ldHJpY0dGTC8lKGxhcmdlb3BHRkwvJS5tb3ZhYmxlbGltaXRzR0ZMLyUnYWNjZW50R0ZMLyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUTBtZWRpdW1tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGaG4vJShtaW5zaXplR1EiMUYnLyUobWF4c2l6ZUdRKWluZmluaXR5RictSSNtbkdGJDYkUSIyRidGSEZILUZANiQtRiM2JUY6RkQtRmJvNiRRIjNGJ0ZIRkgtRkA2JC1GIzYlRjpGRC1GYm82JFEiNkYnRkhGSC8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGaHAvJSliZXZlbGxlZEdGTEYr .8. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUkjbW5HRiQ2JFEiMUYnL0YzUSdub3JtYWxGJy1GIzYkLUklbXN1cEdGJDYlLUYsNiVRInNGJ0YvRjItRiM2Iy1GOzYkUSIyRidGPi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSShtZmVuY2VkR0YkNiQtRiM2JUZCLUkjbW9HRiQ2MFEiK0YnRj4vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRmVuLyUpc3RyZXRjaHlHRmVuLyUqc3ltbWV0cmljR0Zlbi8lKGxhcmdlb3BHRmVuLyUubW92YWJsZWxpbWl0c0dGZW4vJSdhY2NlbnRHRmVuLyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUTBtZWRpdW1tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGZ28vJShtaW5zaXplR0Y9LyUobWF4c2l6ZUdRKWluZmluaXR5RictRjs2JFEiNEYnRj5GPi8lLmxpbmV0aGlja25lc3NHUSIxRicvJStkZW5vbWFsaWduR1EnY2VudGVyRicvJSludW1hbGlnbkdGZ3AvJSliZXZlbGxlZEdGZW5GKw==. 9. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkmbWZyYWNHRiQ2KC1JKG1mZW5jZWRHRiQ2JC1GIzYlLUkjbWlHRiQ2JVEic0YnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JI21vR0YkNjBRKiZ1bWludXMwO0YnL0Y7USdub3JtYWxGJy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGRS8lKXN0cmV0Y2h5R0ZFLyUqc3ltbWV0cmljR0ZFLyUobGFyZ2VvcEdGRS8lLm1vdmFibGVsaW1pdHNHRkUvJSdhY2NlbnRHRkUvJSVmb3JtR1EmaW5maXhGJy8lJ2xzcGFjZUdRMG1lZGl1bW1hdGhzcGFjZUYnLyUncnNwYWNlR0ZXLyUobWluc2l6ZUdRIjFGJy8lKG1heHNpemVHUSlpbmZpbml0eUYnLUkjbW5HRiQ2JEZmbkZBRkEtRiM2JC1JJW1zdXBHRiQ2JUYzLUYjNiMtRltvNiRRIjJGJ0ZBLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGJy1GLzYkLUYjNiVGX28tRj42MFEiK0YnRkFGQ0ZGRkhGSkZMRk5GUEZSRlVGWEZaRmduRmpuRkEvJS5saW5ldGhpY2tuZXNzR1EiMUYnLyUrZGVub21hbGlnbkdRJ2NlbnRlckYnLyUpbnVtYWxpZ25HRmZwLyUpYmV2ZWxsZWRHRkU=.10. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUkjbW5HRiQ2JFEiMUYnL0YzUSdub3JtYWxGJy1GIzYkLUYsNiVRInNGJ0YvRjItSShtZmVuY2VkR0YkNiQtRiM2JUZCLUkjbW9HRiQ2MFEiK0YnRj4vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRlAvJSlzdHJldGNoeUdGUC8lKnN5bW1ldHJpY0dGUC8lKGxhcmdlb3BHRlAvJS5tb3ZhYmxlbGltaXRzR0ZQLyUnYWNjZW50R0ZQLyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUTBtZWRpdW1tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGXG8vJShtaW5zaXplR0Y9LyUobWF4c2l6ZUdRKWluZmluaXR5RidGOkY+LyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0Zpby8lKWJldmVsbGVkR0ZQRis=.11. Find the inverse transform of LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUkjbW5HRiQ2JFEiMUYnL0YzUSdub3JtYWxGJy1GIzYkLUklbXN1cEdGJDYlLUYsNiVRInNGJ0YvRjItRiM2Iy1GOzYkUSIzRidGPi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictSShtZmVuY2VkR0YkNiQtRiM2JUZFLUkjbW9HRiQ2MFEiK0YnRj4vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRmVuLyUpc3RyZXRjaHlHRmVuLyUqc3ltbWV0cmljR0Zlbi8lKGxhcmdlb3BHRmVuLyUubW92YWJsZWxpbWl0c0dGZW4vJSdhY2NlbnRHRmVuLyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUTBtZWRpdW1tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGZ28vJShtaW5zaXplR0Y9LyUobWF4c2l6ZUdRKWluZmluaXR5RidGOkY+LyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZkcC8lKWJldmVsbGVkR0ZlbkYr.
<Text-field style="Heading 1" layout="Heading 1">Applications</Text-field>
<Text-field style="Heading 2" layout="Heading 2">Example 1 An Initial-Value Problem</Text-field>s := 's' : t := 't' : y := 'y' :
deq := diff(y(t),t,t) - 6*diff(y(t),t) + 9*y(t) = t^2*exp(3*t);
eqn := laplace(deq,t,s);
eqn := subs( [y(0) = 2, D(y)(0) = 6], eqn);
Y :=unapply(convert(solve(eqn, laplace(y(t), t, s)),parfrac,s),s);
y := unapply(invlaplace(Y(s),s,t),t);
y := 'y' : dsolve({deq, y(0) =2, D(y)(0) = 6}, y(t), method=laplace);
<Text-field style="Heading 2" layout="Heading 2">Example 2. Using a Unit Step Function</Text-field>Solve x'' + 16 x = f(t) , x(0) = 0, x'(0) = 1 wheref(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYmLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkjbW9HRiQ2MFEifGZyRicvRjNRJ25vcm1hbEYnLyUmZmVuY2VHRjEvJSpzZXBhcmF0b3JHUSZmYWxzZUYnLyUpc3RyZXRjaHlHRjEvJSpzeW1tZXRyaWNHRj8vJShsYXJnZW9wR0Y/LyUubW92YWJsZWxpbWl0c0dGPy8lJ2FjY2VudEdGPy8lJWZvcm1HUSdwcmVmaXhGJy8lJ2xzcGFjZUdRLnRoaW5tYXRoc3BhY2VGJy8lJ3JzcGFjZUdGTy8lKG1pbnNpemVHUSIxRicvJShtYXhzaXplR1EpaW5maW5pdHlGJy1JJ210YWJsZUdGJDY2LUkkbXRyR0YkNictSSRtdGRHRiQ2KC1GIzYkLUYsNiVRJGNvc0YnL0YwRj9GOS1JKG1mZW5jZWRHRiQ2JC1GIzYkLUkjbW5HRiQ2JFEiNEYnRjktRiw2JVEidEYnRi9GMkY5LyUpcm93YWxpZ25HRi4vJSxjb2x1bW5hbGlnbkdGLi8lK2dyb3VwYWxpZ25HRi4vJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGZXAtRmluNigtRiM2KC1GZ282JVEiMEYnLyUrZm9yZWdyb3VuZEdRLFsyMDAsMCwyMDBdRidGOS1GNjYxUSYmbGVxO0YnRl9xRjkvRjxGP0Y9L0ZBRj9GQkZERkZGSC9GS1EmaW5maXhGJy9GTlEvdGhpY2ttYXRoc3BhY2VGJy9GUUZqcUZSRlUtRiw2JkZccEYvRl9xRjItRjY2MVEiPEYnRl9xRjlGZXFGPUZmcUZCRkRGRkZIRmdxRmlxRltyRlJGVS1GIzYjLUYsNiZRJSZwaTtGJ0Zgb0ZfcUY5RitGXXBGX3BGYXBGY3BGZnBGXXBGX3BGYXAtRmZuNictRmluNigtRiM2Iy1GZ282JEZecUY5Rl1wRl9wRmFwRmNwRmZwLUZpbjYoLUYjNiVGKy1GIzYlRlxyLUY2NjFRJSZnZTtGJ0ZfcUY5RmVxRj1GZnFGQkZERkZGSEZncUZpcUZbckZSRlUtRiw2JUZlckZgb0Y5RitGXXBGX3BGYXBGY3BGZnBGXXBGX3BGYXAvJSZhbGlnbkdRJWF4aXNGJy9GXnBRKWJhc2VsaW5lRicvRmBwUSdjZW50ZXJGJy9GYnBRJ3xmcmxlZnR8aHJGJy8lL2FsaWdubWVudHNjb3BlR0YxLyUsY29sdW1ud2lkdGhHUSVhdXRvRicvJSZ3aWR0aEdGZnQvJStyb3dzcGFjaW5nR1EmMS4wZXhGJy8lLmNvbHVtbnNwYWNpbmdHUSYwLjhlbUYnLyUpcm93bGluZXNHUSVub25lRicvJSxjb2x1bW5saW5lc0dGYXUvJSZmcmFtZUdGYXUvJS1mcmFtZXNwYWNpbmdHUSwwLjRlbX4wLjVleEYnLyUqZXF1YWxyb3dzR0Y/LyUtZXF1YWxjb2x1bW5zR0Y/LyUtZGlzcGxheXN0eWxlR0Y/LyUlc2lkZUdRJnJpZ2h0RicvJTBtaW5sYWJlbHNwYWNpbmdHRl51LUY2NjBGLkY5RmVxRj1GZnFGQkZERkZGSC9GS0YuL0ZOUSQwZW1GJy9GUUZodkZSRlU=The function f(t) can be interpreted as an external force that is acting on a mechanical system for only a short period of time and then is removed. Although this problem can be solved conventionally, it is not at all convenient when f(t) is defined in a piecewise manner.f := t -> cos(4*t) - cos(4*(t-Pi))*Heaviside(t - Pi);
plot(f(t), t = 0..2*Pi, color = red);
x := 'x' :
deq := diff(x(t),t,t) + 16*x(t) = f(t);
eqn := laplace(deq,t,s) ;
eqn := subs([x(0) = 0, D(x)(0) = 1], eqn);
X := unapply(expand(solve(eqn, laplace(x(t),t,s))),s);
x := unapply(invlaplace(X(s),s,t),s);
plot(x(t), t = 0..3*Pi, color = blue);
<Text-field style="Heading 2" layout="Heading 2">Example 3. Solving an integral equation</Text-field>Solve NiMvLSUiZkc2IyUidEcsKComIiIkIiIiKiQpRiciIiNGK0YrRistJSRleHBHNiMsJEYnISIiRjMtJSRJbnRHNiQqJi1GJTYjJSR0YXVHRistRjA2IywkRjpGM0YrRidGMw== for f(t).f := 'f' :
eq1 := f(t) = 3*t^2 - exp(-t) - int(f(tau)*exp(t-tau), tau = 0..t);
eq := laplace(eq1,t,s);
F := unapply(convert(simplify(solve(eq,laplace(f(t),t,s))),parfrac,s),s);
f:= unapply(invlaplace(F(s),s,t),t);
<Text-field style="Heading 2" layout="Heading 2">Example 4. Using the Dirac Delta function</Text-field>Solve y'' + y = 4delta(t - 2 Pi), y(0) = 1, y'(0) = 0.y:='y':deq := diff(y(t),t,t) + y(t) = 4*Dirac(t - 2*Pi);
eqn := laplace(deq,t,s);
eqn := subs( [y(0) = 1, D(y)(0) = 0] , eqn);
Y := unapply(simplify(solve(eqn, laplace(y(t), t, s))),s);
y := unapply(invlaplace(Y(s),s,t),t);
plot(y(t), t = 0..4*Pi , color = red);
y := 'y' :
dsolve({deq, y(0) = 1, D(y)(0) = 0}, y(t),method= laplace);
<Text-field style="Heading 2" layout="Heading 2">Example 5. Solving a system of equations</Text-field>Solve LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic='' + 10LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic= - 4LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic= = 0 -4LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic= + LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic='' + 4LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic= = 0subject to LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic=(0) = 0, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMUYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic='(0) = 1, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic=(0) = 0, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1YkdGJDYmLUkjbWlHRiQ2JVEieEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y2USdub3JtYWxGJy8lL3N1YnNjcmlwdHNoaWZ0R1EiMEYnL0krbXNlbWFudGljc0dGJFEnYXRvbWljRic='(0) = -1.deq1 := diff(x1(t),t,t) + 10*x1(t) - 4*x2(t) = 0;
deq2 := -4*x1(t) + diff(x2(t),t,t) + 4*x2(t) = 0;
eqn1 := laplace(deq1,t,s);
eqn2 := laplace(deq2,t,s);
eqn1 := subs([x1(0) = 0, D(x1)(0) = 1, x2(0) = 0, D(x2)(0) = -1],eqn1);
eqn2 := subs([x1(0) = 0, D(x1)(0) = 1, x2(0) = 0, D(x2)(0) = -1],eqn2);
soln := solve({eqn1,eqn2}, {laplace(x1(t), t, s) ,laplace(x2(t), t, s)});
X1 := unapply(rhs(soln[2]),s);
X2 := unapply(rhs(soln[1]),s);
x1 := unapply(invlaplace(X1(s),s,t),t);
x2 := unapply(invlaplace(X2(s),s,t),t);
x1 := 'x1' : x2 := 'x2' :
dsolve({deq1,deq2,x1(0) = 0, D(x1)(0) = 1, x2(0) = 0, D(x2)(0) = -1},{x1(t),x2(t)});
(a)^(b);
<Text-field style="Heading 2" layout="Heading 2">Exercises</Text-field>Solve the following IVPs1. y' + 4y = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JS1GOTYwUSomdW1pbnVzMDtGJ0Y8L0Y/USZmYWxzZUYnL0ZCRmpuL0ZERmpuL0ZGRmpuL0ZIRmpuL0ZKRmpuL0ZMRmpuL0ZOUSZpbmZpeEYnL0ZRUTBtZWRpdW1tYXRoc3BhY2VGJy9GVEZkby9GV1EiMUYnL0ZZUSlpbmZpbml0eUYnLUkjbW5HRiQ2JFEiNEYnRjwtRiw2JVEidEYnRi9GMi8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRic= , y(0) = 2.2. y'' - 6y' + 13y = 0, y(0) = 0, y'(0) = -3.3. y'' - 4y' + 4y = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUklbXN1cEdGJDYlLUkjbWlHRiQ2JVEidEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiM0YnL0Y2USdub3JtYWxGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRic=LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUklbXN1cEdGJDYlLUkjbW9HRiQ2MFEvJkV4cG9uZW50aWFsRTtGJy9GM1Enbm9ybWFsRicvJSZmZW5jZUdRJnVuc2V0RicvJSpzZXBhcmF0b3JHRkAvJSlzdHJldGNoeUdGQC8lKnN5bW1ldHJpY0dGQC8lKGxhcmdlb3BHRkAvJS5tb3ZhYmxlbGltaXRzR0ZALyUnYWNjZW50R0ZALyUlZm9ybUdRJ3ByZWZpeEYnLyUnbHNwYWNlR1EkMGVtRicvJSdyc3BhY2VHUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYnLyUobWluc2l6ZUdGLi8lKG1heHNpemVHRi4tRiM2JC1JI21uR0YkNiRRIjJGJ0Y8LUYsNiVRInRGJ0YvRjIvJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYn, y(0) = 0, y'(0) = 0.4. y'' + y = sin t, y(0) = 1, y'(0) = -1.5. y' + y = f(t) where f(t) = LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUkobWZlbmNlZEdGJDYmLUknbXRhYmxlR0YkNjYtSSRtdHJHRiQ2Jy1JJG10ZEdGJDYoLUYjNiMtSSNtbkdGJDYkUSIwRicvRjNRJ25vcm1hbEYnLyUpcm93YWxpZ25HRi4vJSxjb2x1bW5hbGlnbkdGLi8lK2dyb3VwYWxpZ25HRi4vJShyb3dzcGFuR1EiMUYnLyUrY29sdW1uc3BhbkdGUS1GPzYoLUYjNictRkQ2JUZGLyUrZm9yZWdyb3VuZEdRLFsyMDAsMCwyMDBdRidGRy1JI21vR0YkNjBRJiZsZXE7RidGRy8lJmZlbmNlR1EmZmFsc2VGJy8lKnNlcGFyYXRvckdGXW8vJSlzdHJldGNoeUdGXW8vJSpzeW1tZXRyaWNHRl1vLyUobGFyZ2VvcEdGXW8vJS5tb3ZhYmxlbGltaXRzR0Zdby8lJ2FjY2VudEdGXW8vJSVmb3JtR1EmaW5maXhGJy8lJ2xzcGFjZUdRL3RoaWNrbWF0aHNwYWNlRicvJSdyc3BhY2VHRl9wLyUobWluc2l6ZUdRIjFGJy8lKG1heHNpemVHUSlpbmZpbml0eUYnLUYsNiVRInRGJ0YvRjItRmhuNjFRJSZsdDtGJ0ZaRkdGW29GXm9GYG9GYm9GZG9GZm9GaG9Gam9GXXBGYHBGYnBGZXAtRkQ2JUZkcEZaRkdGSUZLRk1GT0ZSRklGS0ZNLUY8NictRj82KC1GIzYjLUZENiRRIjVGJ0ZHRklGS0ZNRk9GUi1GPzYoLUYjNiVGKy1GIzYlLUYsNiZGanBGL0ZaRjItRmhuNjFRJSZnZTtGJ0ZaRkdGW29GXm9GYG9GYm9GZG9GZm9GaG9Gam9GXXBGYHBGYnBGZXBGXnFGK0ZJRktGTUZPRlJGSUZLRk0vJSZhbGlnbkdRJWF4aXNGJy9GSlEpYmFzZWxpbmVGJy9GTFEnY2VudGVyRicvRk5RJ3xmcmxlZnR8aHJGJy8lL2FsaWdubWVudHNjb3BlR0YxLyUsY29sdW1ud2lkdGhHUSVhdXRvRicvJSZ3aWR0aEdGYXMvJStyb3dzcGFjaW5nR1EmMS4wZXhGJy8lLmNvbHVtbnNwYWNpbmdHUSYwLjhlbUYnLyUpcm93bGluZXNHUSVub25lRicvJSxjb2x1bW5saW5lc0dGXHQvJSZmcmFtZUdGXHQvJS1mcmFtZXNwYWNpbmdHUSwwLjRlbX4wLjVleEYnLyUqZXF1YWxyb3dzR0Zdby8lLWVxdWFsY29sdW1uc0dGXW8vJS1kaXNwbGF5c3R5bGVHRl1vLyUlc2lkZUdRJnJpZ2h0RicvJTBtaW5sYWJlbHNwYWNpbmdHRmlzRkcvJSVvcGVuR1EifGZyRicvJSZjbG9zZUdGLg== , y(0) = 0.6. y'' + 4y = sin( t ) U(t - 2x0), y(0) = 1, y'(0) =0.7. NiMsJi0lImZHNiMlInRHIiIiKiYiIiNGKC0lJEludEc2JComLSUkY29zRzYjLCZGJ0YoJSR0YXVHISIiRigtRiU2I0YzRigvRjM7IiIhRidGKEYo = NiMsJiomIiIlIiIiLSUkZXhwRzYjLCQlInRHISIiRiZGJi0lJHNpbkc2I0YrRiY=8. y'' + 16 y = delta(t - 2 x0), y(0) = 0, y'(0) = 0.9. y'' + y = delta(t - LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUYsNiVRJSZwaTtGJy9GMFEmZmFsc2VGJy9GM1Enbm9ybWFsRictRiM2Iy1JI21uR0YkNiRRIjJGJ0Y/LyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZMLyUpYmV2ZWxsZWRHRj5GKw==) + delta(t - 3LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYlLUkjbWlHRiQ2JVEhRicvJSdpdGFsaWNHUSV0cnVlRicvJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YnLUYjNiMtSSZtZnJhY0dGJDYoLUYsNiVRJSZwaTtGJy9GMFEmZmFsc2VGJy9GM1Enbm9ybWFsRictRiM2Iy1JI21uR0YkNiRRIjJGJ0Y/LyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZMLyUpYmV2ZWxsZWRHRj5GKw==), y(0) = 0, y'(0) = 0.10. dx/dt = -x + y, dy/dt = 2x, x(0) = 0, y(0) = 1.11. LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYkLUkmbWZyYWNHRiQ2KC1GIzYkLUklbXN1cEdGJDYlLUkjbWlHRiQ2JVEiZEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y7USdub3JtYWxGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictRjQ2JVEieEYnRjdGOi1GIzYjLUYxNiUtRjQ2JVEjZHRGJ0Y3RjpGPUZFLyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZXLyUpYmV2ZWxsZWRHUSZmYWxzZUYnLUY0NiVRIUYnRjdGOg== + x - y = 0, LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkmbWZyYWNHRiQ2KC1GIzYkLUklbXN1cEdGJDYlLUkjbWlHRiQ2JVEiZEYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1GIzYjLUkjbW5HRiQ2JFEiMkYnL0Y7USdub3JtYWxGJy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRictRjQ2JVEieUYnRjdGOi1GIzYjLUYxNiUtRjQ2JVEjZHRGJ0Y3RjpGPUZFLyUubGluZXRoaWNrbmVzc0dRIjFGJy8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGJy8lKW51bWFsaWduR0ZXLyUpYmV2ZWxsZWRHUSZmYWxzZUYn + y - x = 0, x(0) = 0, x'(0) = -2, y(0) = 1, y'(0) = 1.12. dx/dt - 4 x + 2 y = 2 U(t-1), dy/dt + y - 3x = U(t-1), x(0) = 0, y(0) = 1/2.
<Text-field bookmark="applic periodic fcn" style="Heading 1" layout="Heading 1">Further applications</Text-field>
<Text-field style="Heading 2" layout="Heading 2">Introduction</Text-field>This section contains a set of solutions to differential equations directly and using Maples dsolve.It also contains a method to define a periodic function using a piecewise function definition so that one can obtain a visualization of the periodic function and then using this structure solve a differential equation having a periodic driving force.
<Text-field style="Heading 2" layout="Heading 2">Initialization</Text-field>with(inttrans):with(DEtools):with(plots):with(student):
<Text-field style="Heading 2" layout="Heading 2">Applications</Text-field>If you need assistance with Laplace transforms click here Intro to LPT
<Text-field style="Heading 3" layout="Heading 3">Example 1 Solving an IVP.</Text-field>Solve NiMvLCgtLS0lI0BARzYkJSJERyIiIzYjJSJ5RzYjJSJ0RyIiIiomIiIlRjAtLUYqRixGLkYwRjAqJiIiJ0YwLUYtRi5GMEYwLCZGMEYwLSUkZXhwRzYjLCRGLyEiIkYw, y(0)=D(y)(0)=0Define the eqn and take the laplace transform of it.y:='y':t:='t':Y:='Y':eqn:=(D@@2)(y)(t)+4*D(y)(t)+6*y(t)=1+exp(-t);eqn1:=laplace(eqn,t,s);Substituting in the initial values and solving for L{y(t)}eqn2:=subs({y(0)=0,D(y)(0)=0,laplace(y(t),t,s)=Y(s)},eqn1);solve(eqn2,Y(s));invlaplace(%,s,t);Solving the eqn directly using dsolve with the option laplacedsolve({eqn,y(0)=0,D(y)(0)=0},y(t),method=laplace);Solving the eqn with dsolve but without the laplace option.dsolve({eqn,y(0)=0,D(y)(0)=0},y(t));simplify(%);
<Text-field style="Heading 3" layout="Heading 3">Example 2. Working with a piecewise continuous or periodic function</Text-field>
<Text-field bookmark="piecewise fcn" style="Heading 4" layout="Heading 4">Defining a piecewise continuous function</Text-field>Defining a piecewise functionLUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYoLUkjbWlHRiQ2JVEiZkYnLyUnaXRhbGljR1EldHJ1ZUYnLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGJy1JKG1mZW5jZWRHRiQ2JC1GIzYjLUYsNiVRInRGJ0YvRjIvRjNRJ25vcm1hbEYnLUkjbW9HRiQ2MFEiPUYnRj0vJSZmZW5jZUdRJmZhbHNlRicvJSpzZXBhcmF0b3JHRkUvJSlzdHJldGNoeUdGRS8lKnN5bW1ldHJpY0dGRS8lKGxhcmdlb3BHRkUvJS5tb3ZhYmxlbGltaXRzR0ZFLyUnYWNjZW50R0ZFLyUlZm9ybUdRJmluZml4RicvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYnLyUncnNwYWNlR0ZXLyUobWluc2l6ZUdRIjFGJy8lKG1heHNpemVHUSlpbmZpbml0eUYnLUZANjBRIn5GJ0Y9RkNGRkZIRkpGTEZORlAvRlNRIUYnL0ZWUSQwZW1GJy9GWUZgb0ZaRmduLUYjNiUtRkA2MVEifGZyRicvJTBmb250X3N0eWxlX25hbWVHUSgyRH5NYXRoRidGPS9GREYxRkYvRklGMUZKRkxGTkZQL0ZTUSdwcmVmaXhGJy9GVlEudGhpbm1hdGhzcGFjZUYnL0ZZRl9wRlpGZ24tSSdtdGFibGVHRiQ2Ny1JJG10ckdGJDYnLUkkbXRkR0YkNigtRiw2JkY8Ri9GZ29GMi8lKXJvd2FsaWduR0Zeby8lLGNvbHVtbmFsaWduR0Zeby8lK2dyb3VwYWxpZ25HRl5vLyUocm93c3BhbkdRIjFGJy8lK2NvbHVtbnNwYW5HRmRxLUZocDYoLUYjNictRiw2JUZeb0YvRjItRiM2JS1JI21uR0YkNiVRIjBGJ0Znb0Y9LUZANjFRJSZsZTtGJ0Znb0Y9RkNGRkZIRkpGTEZORlBGXW9GX29GYW9GWkZnbkZqcC1GQDYzUSZ+YW5kfkYnLyUlYm9sZEdGMUZnby9GM1ElYm9sZEYnLyUrZm9udHdlaWdodEdGXHNGQ0ZGRkhGSkZMRk5GUEZdby9GVlEwbWVkaXVtbWF0aHNwYWNlRicvRllGYHNGWkZnbi1GIzYlRmpwLUZANjFRIjxGJ0Znb0Y9RkNGRkZIRkpGTEZORlBGUkZVRlhGWkZnbi1GYHI2JVEiMkYnRmdvRj1GW3JGXHFGXnFGYHFGYnFGZXFGXHFGXnFGYHEtRmVwNictRmhwNigtRiM2JS1GYHI2JVEiNEYnRmdvRj0tRkA2MVEoJm1pbnVzO0YnRmdvRj1GQ0ZGRkhGSkZMRk5GUEZSRl9zRmFzRlpGZ25GanBGXHFGXnFGYHFGYnFGZXEtRmhwNigtRiM2J0Zbci1GIzYlRmdzRmRzRmpwRmZyLUYjNiVGanBGY3JGYHRGW3JGXHFGXnFGYHFGYnFGZXFGXHFGXnFGYHEtRmVwNictRmhwNihGX3JGXHFGXnFGYHFGYnFGZXEtRmhwNigtRiw2JlEqb3RoZXJ3aXNlRidGL0Znb0YyRlxxRl5xRmBxRmJxRmVxRlxxRl5xRmBxLyUmYWxpZ25HUSVheGlzRicvRl1xUSliYXNlbGluZUYnL0ZfcVEnY2VudGVyRicvRmFxUSd8ZnJsZWZ0fGhyRicvJS9hbGlnbm1lbnRzY29wZUdGMS8lLGNvbHVtbndpZHRoR1ElYXV0b0YnLyUmd2lkdGhHRmR2LyUrcm93c3BhY2luZ0dRJjEuMGV4RicvJS5jb2x1bW5zcGFjaW5nR1EkMmVtRicvJSlyb3dsaW5lc0dRJW5vbmVGJy8lLGNvbHVtbmxpbmVzR0Zfdy8lJmZyYW1lR0Zfdy8lLWZyYW1lc3BhY2luZ0dRLDAuNGVtfjAuNWV4RicvJSplcXVhbHJvd3NHRkUvJS1lcXVhbGNvbHVtbnNHRkUvJS1kaXNwbGF5c3R5bGVHRkUvJSVzaWRlR1EmcmlnaHRGJy8lMG1pbmxhYmVsc3BhY2luZ0dRJjAuOGVtRidGW3JGW3I=f:=t->piecewise(t>=0 and t<2,t,t>2 and t<=4,4-t,0);plot(f(t),t=0..6);
<Text-field style="Heading 4" layout="Heading 4">Defining a periodic function</Text-field>f:=t->piecewise(t>=0 and t<2,t,t>=2 and t<4,4-t,0);f(t);Maple does not have a mechanism to handle periodic functions in general.We can make a piecewise continuous function periodic by defining another function g(k) = f(t-4k) which has a period of 4. and g(k) represents the kth period.Making f (t) defined in the previous section into a periodic function.g:=k->f(t-4*k);g(1);Plotting the first 4 periods of our periodic function.plot({g(k) $k=0..3},t=0..20);
<Text-field style="Heading 4" layout="Heading 4">Calculating a Laplace transform for a piecewise or periodic function</Text-field>In order to calculate the laplace transform we first convert f(t) and g(t) using the Heaviside functionconvert(f(t),Heaviside);laplace(%,t,s);invlaplace(%,s,t);convert(%,piecewise);convert(g(k),Heaviside);gl:=laplace(%,t,s);glv:=invlaplace(gl,s,t);convert(glv,piecewise,t);Plotting the first 11 periods of g(k)plot({glv $k=0..10},t=0..50);
<Text-field style="Heading 4" layout="Heading 4">Solving an IVP with a piecewise continuous driving function</Text-field>
<Text-field style="_pstyle18" selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" font_style_name="_pstyle18" size="12" layout="_pstyle18" linebreak="space" spaceabove="0" rightmargin="0" bullet="none" firstindent="0" layout_style_name="_pstyle18" linespacing="0.0" pagebreak-before="false" leftmargin="0" alignment="left" initial="0" spacebelow="0"><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12">Example 1. Solve </Font><Equation executable="false" style="_pstyle18" input-equation="(diff(y(t), t, t))+4*(diff(y(t), t))+6*y(t) = f(t)">NiMvLCgtJSVkaWZmRzYlLSUieUc2IyUidEdGK0YrIiIiKiYiIiVGLC1GJjYkRihGK0YsRiwqJiIiJ0YsRihGLEYsLSUiZkdGKg==</Equation><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12"> where and y(0) = y '(0) = 0</Font></Text-field>Defining the DEQ using the fact that f(t) was defined in the previous section.eqn1:=(D@@2)(y)(t)+4*D(y)(t)+6*y(t)=convert(f(t),Heaviside);Applying the Laplace transformeqnlp:=laplace(eqn1,t,s);Substituting eqnlp1:=subs({y(0)=0,D(y)(0)=0,laplace(y(t),t,s)=Y(s)},eqnlp);solve(eqnlp1,Y(s));sol:=invlaplace(%,s,t);solp:=convert(sol,piecewise);plot({sol},t=0..20);
<Text-field style="_pstyle18" selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" font_style_name="_pstyle18" size="12" layout="_pstyle18" linebreak="space" spaceabove="0" rightmargin="0" bullet="none" firstindent="0" layout_style_name="_pstyle18" linespacing="0.0" pagebreak-before="false" leftmargin="0" alignment="left" initial="0" spacebelow="0"><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12">Example 2. Solve IVP </Font><Equation executable="false" style="_pstyle18" input-equation="(diff(y(t), t, t))+3*(diff(y(t), t))+y(t) = g(t)">NiMvLCgtJSVkaWZmRzYlLSUieUc2IyUidEdGK0YrIiIiKiYiIiRGLC1GJjYkRihGK0YsRixGKEYsLSUiZ0dGKg==</Equation><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12"> where and y(0) = y '(0)= 1.</Font></Text-field>Defining g(t):gg:=t->piecewise(t<0,0,t>=0 and t<1,1,t>=1 and t<2,-1,t>2,0);Defining the IVPeqn:=(D@@2)(y)(t)+3*D(y)(t)+y(t)=gg(t);eqn1:=simplify(convert(eqn,Heaviside));Applying the Laplace transformleqn:=laplace(eqn1,t,s);Substituting the initial values and solving for Y(s)Y(s):=solve(subs({y(0)=1,D(y)(0)=1},leqn),laplace(y(t),t,s));Applying the inverse Laplace transformsol:=invlaplace(Y(s),s,t);Plotting the input and solution curves together.plot({gg(t),sol},t=0..5);convert(sol,exp);convert(%,piecewise);simplify(%);
<Text-field style="_pstyle18" selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" font_style_name="_pstyle18" size="12" layout="_pstyle18" linebreak="space" spaceabove="0" rightmargin="0" bullet="none" firstindent="0" layout_style_name="_pstyle18" linespacing="0.0" pagebreak-before="false" leftmargin="0" alignment="left" initial="0" spacebelow="0"><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12">Example 3. Solve IVP </Font><Equation executable="false" style="_pstyle18" input-equation="(diff(y(t), t, t))+3*(diff(y(t), t))+4*y(t) = f(t)">NiMvLCgtJSVkaWZmRzYlLSUieUc2IyUidEdGK0YrIiIiKiYiIiRGLC1GJjYkRihGK0YsRiwqJiIiJUYsRihGLEYsLSUiZkdGKg==</Equation><Font selection-placeholder="false" executable="false" family="Times New Roman" opaque="false" foreground="[0,0,0]" superscript="false" placeholder="false" readonly="false" subscript="false" bold="false" italic="false" underline="false" background="[255,255,255]" style="_cstyle17" font_style_name="_cstyle17" size="12"> where y(o) = y '(0) = 0</Font></Text-field>f:=t->piecewise(t<1,0,t>=1 and t<2,cos(t),t>=2,sin(t));eqn:=(D@@2)(y)(t)+3*D(y)(t)+4*y(t)=f(t);eqn1:=simplify(convert(eqn,Heaviside));leqn:=laplace(eqn1,t,s);Y(s):=solve(subs({y(0)=0,D(y)(0)=0},leqn),laplace(y(t),t,s));sol:=invlaplace(Y(s),s,t);plot(sol,t=0..10);soln1:=dsolve({eqn1,y(0)=0,D(y)(0)=0},y(t));plot(rhs(soln1),t=0..20);
<Text-field style="Heading 4" layout="Heading 4">Solving an IVP with a periodic driving force.</Text-field>Solve NiMvLCgtJSVkaWZmRzYlLSUieUc2IyUidEdGK0YrIiIiKiYiIiVGLC1GJjYkRihGK0YsRiwqJiIiJ0YsRihGLEYsLSUiZ0dGKg== where g(t) is the periodic function defined by cycling the piecewise continuous function f(t). andg(k) = f( t - 4k ) . g has a period of 4 units JSMlP0c=f:=t->piecewise(t>=0 and t<2,t,t>2 and t<=4,4-t,0);Making f (t) into a periodic function.g:=k->f(t-4*k);defining the DEQ and using only 2 cycles of the input functioneqn2:=(D@@2)(y)(t)+4*D(y)(t)+6*y(t)=simplify(sum(convert(g(k),Heaviside),k=0..1));Applying the Laplace transformleqn2:=laplace(eqn2,t,s);Subs the initial values and solving for Y(s)Y(s):=solve(subs({y(0)=0,D(y)(0)=0},leqn2),laplace(y(t),t,s));Applying the inverse transform to obtain the solution.sol2:=invlaplace(Y(s),s,t);Plotting the input function and the solution curve on the same graph.plot({sol2,g(k) $k=0..1},t=0..10);convert(sol2,piecewise);plot(%, t = 0 .. 10);
Legal Notice: The copyright for this application is owned by the author(s). Neither Maplesoft nor the author are responsible for any errors contained within and are not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact the author for permission if you wish to use this application in for-profit activities.
TUZOV3RLVWI8b2I8Uj1NRExDZE5WWlpKOkBMPkg6VEtHeE1rSjo8T2BMb1xcbFF4bFFXZE1XcHNIcVNobVdoWW9lWE9QbVRQbVZgbXZxeXhxPVhqPXhYcXVYYXhuYVhjRVdjPVVSPVV3ZVl3RUxLRExxdFBxPFI6PXJeYXZedVJBdXJaQG5adFZhdVZiPVdiTVl0TXl2YXl2WXl1WVl4bVl4cXl4cVl5dVl5RVlzRVlwbVhweXl5eXlwcXhwPUo6Pjo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6Ojo6OjpkeTxUeXBDPnFVTENUSmNEWG9YdXNUPGF1cGtjZldNWEBKQ2VVYGROdVRtV3h5eXlwcHVQQ0RTU3VMQ2x1Pjx4VHBRbWxzYl1NaWhVT2BxVGVYU1FPO0BKeFZdd09sOkBzeUZ2PHdcXHRAdHNOblFuXFxWP3c8d1xcP0ZxSmlqWHluWlZ2bnlIRXJtaUJfX3RXaXRbTXl4WVJJSVh2V2d0U1M9O2dRTXdBSUNdSVlyR1hSb2djW0VwcVl0c3huPUJWU1VHdUVBW1d4S3JXYVNIc3NvWUJQa3luS2N0cWdteVVLQVlRWVV3X3JzPXdib1lUV1hJP0lRS3lvW1hAd3lkcXl0WVJHQXlgaXhzW1NseVhhU3lxdXk6bWVsPWRYcXlkSWZ2Z1JJZVNVa1VtVUJHd3VaaXRTO2VRP1M+QWRNYXNua3lTR2JEU3VpbWJTYWJqeXROQXlNdVhsYVRXYUNwO3k/YXQ7X3R4YVR3YXRoP2NqPUdiZ1lWR0NBW2VBa2heaWh5YUlHb1ZkR3h5V2VRYXRhbVZIWXg6U0VJZXd5YWNtY1NCQXZnT3l5c3NFeUJWV0N3UUZ0WVd4WWRNZ2NZX3leVXk/Z2NlW1dYUUNEY3dHdXdITXc/cXd4W2dhY3NjR3J3T3R1S0ZYS3NjW0ZaSUJPcUlySUlda3VJQ2ZSb3NNX3lUU0VXV2NLUXNfcUdIZUlpYVdCc3ZhQVhXb0ZzWVR5dUlZU2RXQ2V0W2ZacE9ZdHZbXFxYU012Tj1YaGx1eGVsXXlsdlVuO1BZc3F2a21tQ3hTRVFQc01PZVVwUUVLTmB5VkFxY3FSUXBZeEhyW3hVXFxBdGdQVmV4bUhIUVlEWHB0TDtleV9cXFhIeHlUcExRPXFKaEprbHFBPXdQeHFPdHBQbXdRPWtXZFNTWWp4aFF0PWxpPFg9UHJcXEhveE1LeHBwZFVQR3hsYDxSYWRXc0VNVWhuTWluYXF2eVxcdF1wSndcXFB0dHQ6bHdfaHk7UHh1RWxXcGZ5cGlReWc8SWJnSHFRP3dSd3ZGZ2NRbm10SV1sWFpvYXV2d1xcXVZpXFw/eXVJakdxeUFfXWpeY2lhXFxedmFZZm1YWXZWX2ZveWRfd1phP3lJUGZOWHBPaW1iSW53aWllUXlaQFtqZltwX2BzP1xcTkBxYXdbPGFfPXFwZEl1XT5nbkhwVWlcXF5hW0FHY1NfeV1wbkhnX29JaT1Ya01gYkteeVVXakZoaENwaWY/bGxoZWxoa0txaz1xZ0NxcUlva0phZFpAXUlPc3BIamdRZ1V2Xk1wXltha1hOb2t4Y0ZheE1YPkVmeD1HSnlZXT11S1dYdWVmY1lDVl9ETztYXW9lRHdJXVVyaElYaEtkdFlndj1zWU14eU1oRUFiZEtkRkVEO01CaW1VWWd2TnNmQnVEZ3F3XnNSWm9pZXlpWUVmRUFzWU9jVTt1Zl9DXjtnPkVJVW1XeV14WltIP1VUaXdoYXliPEVXVUFobWdoVWVlXU9ETHlma1lkT1FETk1zbGVnXW1IR2t5blVyclVoamdidnN0cklDc09pVT91cFVodE1FX2NWVWV5d1dyU2VTdkl3SHFzRVV2d2FTYG12X2tDRWdERUVWT295ZlNGWUdYaFt4ZTt3ZnN5YT9IYmN1X1NpSFVmclN0cXNnSUNVS21SO0lFR0dpRVV4U1Nld2tCUmNpYz9mW0dIc11XQkNlRlNYTWVjQHF3UVlpT0NGaTtiZF9lcGdoQ2NyU0liclVGZktYcE9FPkNkR1VWSF9zcz1HYUVGXFxNaF91REpjWGVXR1NrSUE9VGBbdWhPaUtPeTtJZG9fc0JRZ1BHYmlNeFpJeFs9Uk5RSENVd2xJaFZBcz5NeHY9dDtJZWtlY1tpVG9lQl1ZU1ZzSV1VR2tNZ0M9eE1fY3ZdckNrR2xPeUU9d1ZzeW1vUlBFUkdVV29Lcz4/ZE5HY3FPdkw9RGNnVVVpZD1TZEJZdGFjQmN5VDtzQz8/c1hzQkZFSVBLZHdVaWJVVXVvd3RDeExFUnhHVVBPYz1lZVdXREpfdEJJRmpbUk1XWG9hSW5pRkRZeXZJZkZZSDtFaWZhV0FBZGtRZ1N1SW9ZSFM/c1xcYVlua1ljQ1JYQXk7PXVyU3NVRUdYb3Zta2RVP2JJa3V2SWhmO2hIS1Jtc0lxa0drQ0lFR1NRaVV5P3JbY2h5XURXP1VKd2VvX0hJO0lbaVJQdVlDY2VdeUlRR1NSPVNGY1lASUhOYWJFeWhUO0hcXGdDW2lpRXViWElZWz9GaGtmQWFSeWNjUTtEPE1CTGtzVUd2TV1GT1NXWmFGbm1VVk9CXU1oYGd1XWV3OkNTWFtWVVtkXmlXQ0lUTWtpbmdWbWNZO0V1SWtGWmdldGFTbGtlRF9TbFVkP1NVW1doYF9JSGt1TmFJQkVZQEtoUVtJYlNmbF9DcGdWXUlCZ2NmOkNyT1dXbGlWUFNETXVFa3dCWVFiZ0t4R2lXZmNkZ19jQ29YRHlGb0FGPENZZF9mWlNVS09YbVVFcm12cFdnYVFJZVdHeU1pdU9maGVGWVtVV2dkR3dlWztYQFloPG93c2tUd1VnallkdkVoblRQYExKYXRVbXlvXXhsa1VwZ1BTSG1TT2lTWHRNP0hzSGhXZ2xudT15cE1vc21QV1F0WG1sTERSXmVyYXBwQVBxQFR3dVxcbWY8eXRNb190TlFEbXd1VUJhbFtUS01dVVpcXFZzVVBnXFxPaFhVXWl3PmxUPlR0b2xZVWVNXFxgcTppTkZRa01ldUI8WV55cVtUcXdMeHlZa15tUERoVVRFTFtteGRZVHJVd0hZcHBgUl10c3lobTxcXHJkaE5cXF1WR2VqRXlUQkxsWGhVaWRTa2xWY0lta3VKQVxcT0ZBSnhYVEpcXG9ScFVyXFxxbkVVZjxQT2FvY2lvWHhZVVRSeGhtS0hub1V1QmF2dnh0XUBvcmR5cUlsYHR5Y0V5Zz1TdDxWO0xZYERvREVsQ2hXWWRrcElrU01vcGhuaHFrZU1XPFFYXmRvZ0VtTTxreEFZTT1tcFBLbVRUTW1YZVFMbnVLP0hNZUlVYGBUcU1TZGVOcW14SGVMSz1PVXB4XkBraVlwYHhYVmRvVUBMPVBwckFQSXVSW1FwQFlsdlBXd1FUb01wR2BqT1h5Rmh4QUVUaWVSQURLZ2lvVlBPeVhVbFhUOkl3YzxOZ2VNTnVwXFxYV3JkUUZQUXZsUD1Ub3Nlbz5xWGJpV09cXHlFPVBVaVBBQVNnTHR4WExHPVNUQVNBeGo9QFdpeHdYYFhPQXRIbG9JZW9IaUx2eXVvdU10THRUeUpzQXhCWHJAVHFXWE9zRUtvcHVBRVU8dXlPXFxMVHlQQVhtPXRPVVFuZWFORF1LT1l5THlYYnR4dWhtY1lyWE1raFxceWxMb19lcWB0U2VBT0hdbHFVd2lQbmtQd2xIUGdIcmVoWV5wS2hQd0dQSjs8TzxgcVU9dE14VVVFUFdAUmRJVGZZamphb3dUcU1RalhISlNcXE08RXZhcHBUQG1XTUpAaU9WaHlMUUtxXVQ9RXljPVVocU5hXVBKXFxYXFxMdVtEc1FAT1tYUnc8UmJgUGB0U3VlamNlWVhAVU49ckZleHVIbURta11YUkxhWUVsUm1JUF1QZWNoYHJ4bWE/YXJhYUN4dldRW1xcYVpgeWlGQWo/Z3ZWVmReQG1HeVtoaGp4UXZqSXdNVndQR3lYV19FcGpETm5zeV5FaHZFX2Q6UG5rT2FEQV5DbnhFQW9DaF9ld2M7cGJbSVtad2NVP2twR3d4dmNWVlxcT1dhWUdaV3FiR0dealZrQVFdbVhja2Z3VFZmb3ZaVm5aTHdmb0llUz5lQEh0Y3ZzZ1BuPFlxRE94Y3FiZE5tUHh0cXdoc2ZhZz5teU9lZGhxQ0ZrTldxc3B5XUBfVlFySUl1XW5jTEliPl94ZFFeW3l3XmBeWXFiU3hleWdhPk9rVkBmcFZmZU5obXhlU3duXj9fR09rbGZgUXFnS195Sz95akBweHZ3Ykh0SWB5WWFpP0h2Sl53dlF2WW5nQVZvPVhod2NSZUJJTWZsS1RVX2JgcXJGUUM8VUdSV1k9a1ZXQWl2XVg8Q1N5TXljeXdlb0U+P3R0a3NWZ0JUbXRHSVh2S0RUO0RgYXRwYUdRRVZBPWVmb0hAXVRnc3dzQ2ZXR0ViQ0NMSVl0U3dHO3RSYUM/XWhpW1Rmd1NQVWNTUVlaQ3Vsb0VbS1RuT1NUdURQcWZwUVVfWXhbP1VaPWJgeUN1RVRVZWN0Y3JzYVdJR2hQVVZkQ1hvW0RuO0dUb2Y9QVZCY1lSR2dhYVlic3Z0PVVCdVZJT2VaS2dHbWhIUXJdXXVtc2lmeVRQV3RuZXlaS3lkbUhqb1dSQXNTUUhld0RTPUhqXUM+cWRIW1hISWdrd1RHdXZJX3NnWURnYWJTc2lMWXJiXUljW3VaVXVDZUdOXUlueXlqaVZuTXVKaWJxXUU+PXNIW3RoUURYZ1RcXHFoTndUVm1HZG9TaUtzRF1ERF1VT2tzTz1mWDtYdklkYlV3Umlpc0NFdj90RUFTP2VIW0VIaU95W21jRT9oWTtld0tDclt4O0VDcFVFYUl0Uk1VZU1JQHdGPUd1cUlkcmlYbUFpSG91Ql1VRWt2Ym9EYF1iRGV1XlVIT3N4d0tTb2dWRV9HTlFiQkFkdU1ZUTtZX11YYnFCZVtGRllHRj10WGd4cnlZcEFGRG9pZElSSGdVZj91WEdnXVdndUdpZ11VUlFycDt1PU1IWUlYeGNJYW1zcUVsPHVSPFBNd3R3Tk1xTllNQj9cXGFJaXF2Ym94aGtud0RPdl1ecjphXFxbV2hFeHNuX2NkUW9ATmddb3JMUG5DcHRFP3dKcWk6YWRgP2dqWFxcQm9sOkBkSmlzW3ZlbF5wSz5dVHBjSUhob1Nab1hKT2h3W1dnc2VzdUJmRWddPXV1VVk9cVhaV1ZZTVNaRUNIV0hxZVg8U3VeRXV2WVg7QUZRUUNdXUZsXVNOcUlPPUlMUXdoSXdab2VxRW9PcVZZQFRUcHJXQU5xWXN1eE5BQFdqbHB1YVh5dG1YTVJrZHBJXUtcXExUQD1QZFxcU3hISlNYTmh1bEZZUW10d0poV0k8UXN1UlVwd21cXHJRREx5dU1nTXY+QHBTQHBmdFJpVW5pVFY6dVJSaWw8bFJZPHdsdFNWaUxoSEtEQHZWaVNgRE9mYVR2QXN5TXVLbVFVaHZxbFF1TFdAcWxyYFJkZFJLSW1eUVlBYVh4ZFBcXFR1VmxrdE1ZbXlQQWB4Uml2UlVvTHhLbUFOYWxMYHFWYGVURElPO01ZXFxIb1FpWW5Na0hMTnFoeWxVSlxcdFNedUtKSU1LQVlbcXVmTXJ4QVhmeEp5WHhlYFJQcXhPaW9ybEpXXVhFSFh3XFxsSnFyPVh3TjxUPmBuRlBrbEh2XkxUZF1rdml1Oll3bGhXa1R5RHBMU1VWVXFRQ0F1VFRsaVBvcHVvVEhOU1F5UnRzPklxS1lLaFROUU1zZUFqb2FsclF2YklzbE1wPVxcb2pMVU1EdURReW1hb2lRdWxtUE1FTHdocHVwbG5JdnlwUGBYbENETT5MWUBgcmRxdG95bkBNTEZUVVVQb1xcVVdSXFxXTWV0T0FvRWV3TElVY3RSd0B0XUVSR0BYdHFLdUhRV3FqV0xxWmBMVFVPVHVzbUhQY1lrP0ROPXVUXFxhWFNlTE51S3J0dGZAa0l1blVUWENNdFl5UlVRcGxYd2BYdj1pWHBwdUxtUlVxd1RNbVtdcXhoTEVsdD5sTmlAcVE9UV9sUkw8TmdlcmhoWHdBcnlBTD1pd11JeFlUVXloajtwb3FYUG1VZ0hHXFxnYW5mV2ZGPmhyQXd0d3lbWXM8VnVHWGhTR3hlUGpNXmV4blxcdmFiSE5qVGZmRll3RE5yZUBxb2hlSFdtb1dgXVBcXGdmcV1Ja3h4XFw/dmtubmNcXGdpdXBvdkloTWFaT0lraklkVnF0dj9lZm5oZWBpPU9peFZ1ZVZvcHhqSk91TllgW1dcXGpYXFxTTmtlcXJRX3BVZ2hqTmlOUXRwR1xcQ0llX0lhYllzQHd3QndcXExgeE8/cmBxWmk/Y0BXc1dgXkBmam9nZXBwamtJcG5Ya0tQbmRHYWRHaWRvY0U+bT9GamZfYllmXFxcXD9wXUhpZU5xV2dnZUl1Q0FuaGlad2FlcFlua2dlRnlqdk9odV9bR1FrcGlvU05hP25kaXByVUZqY1ZcXHBRbmd3XVI/XVdGZVd4YD5pX0hAdEF3ZGJueTx4X19PYEZ5Z2dxdWpBdEpoYWlBblNBcz14d3RwXmFZbmxvbG4/ZVlRdEFebUp2d0Q/a1xcUWxdeHFNUGNgX3NqVl1ndnJlT3NJT2twUF5WeV5bVndgT1tnd21McWldTm1aQGhCQXJpUF1PPltASGRtWVp5aXJbTm48WXBlTmZvbnNvXl1kbmZJWXVYd2tFQWNVeW5eQWBdVmV5WXVsUG9nQW47P1xcSz9tdF5ncF5qWEd4Zj55c2Zac2d1PWBzZWJfYUlFU1NKY1dld3RtQ3JFQ2ZnRVJhcUVOQ2hCO2ZeSXZ4WUw9UFNdPXlLWG1HZU1ZTG1yVFNCcExfYFVBbG1YbVhsVVRYRW5eRXNTbW1meVJFWHNERXdlbHZRcWxRYVhAQHRqPHBrVFlrRFNOcXhQUWpsdXNpVEpFTFhRXFxSd2BzUGFTVVlKd1BqZGVzX1FzS2BqQElqX0R1Rm1KbVBMbWxsaDxTU1BLVjxXW2VPYWFUTkB3TGx0dj1xZEBPT0hyYzxLPmh1aFBQPUFwU1VSUF1tYklWU3VybERMcXBLdWFWbGlWPklvT3hKeEx5R1hPaHF0PVFQQlFWSXRSamRWP11QRlBQQ3l2c11ZQl1SWEFzUEx5c1FUXk11TFVPRE11ZURQPVVQcEhzRlV4OlhKYGhObEVZS3lrcVFMUUhTRXVyXmFYX1hKSF1VeXh0Z01SQ1h0anVvP0VRV01MW2FSU2lraWRvZUxzVWR1V0VNdGhZWnlRW3F3eEhUW3RPdTxWR3hxYmBxcDxPUUFXT2VZSUl3XlR2YEhyTnlQO0VLaERMaVRxY1hMcTxOWGVqc0VLc2VUO01ZQTxvc211ZkBVQHR4VU1KWWFNRnV2VmFqVWVsdl14WGBuY3VUaFR4Qlxcd3h0dkNpdUBIc1FVUTptc0p5VVZYTE9lVUFMbWRhWV1UTW91cUVFeFdgeEs9UVFMeUdBeWlIUFxceE9mXXRHPmNKd2BneHdeZl1tSWRKd2dYaXliWF1fXlxcXXhdd1hvb3ZmSmB2Z1FrbFdyaHFgc3hxVGhkX0F1WEhvdGF1eHF2VlBzPmZYUUVHX1lHeXVqR1dxYUNPeUU+V1hbd3VFd3lzTUhzQUNhd1lmc0lpcXZXaVdwV0dvR1ltcXdBZWg7X1hxR1N5W1lRVVc8a0ZhVUdtdWhxZVlFO3hkd2JEVURkV1Y8T1lqbXdjXXJMP1RwdXdGX3NuV3VtaWlhQUlueUJbYVVieXhcXHl5YGNTTG1IeHNJbndZTHdmPW9iX2t0eGdVSldUQl1UdEl2S2tERE1JQ01WWkNIPFdXRjt2WGV1T0dlXlFlTHdpa11Ia0NmclVYdV9EZ29DW09JeXVoX0l5YltlRWhxcnlRP013VGV4SXVOYnVtdjxzT2l3eV11Tz5pZT9vTlhwbkZiXWl5a3l2QHBuTT9eYlFiY09wXUBwTV93T0laXFxpXXRWcEdJdT1QZGJIZk14Y3hYYXQ/YVdQWnN3dz54YUR2djx3cVF2eWtecGlBcl9AZmRZeWZveHNhY3RXX3V2Z0JQbXF2bUtfWk1BclpXWnlBdkNQbXVZZFxcQWJacF1aTmdYd3J5WGF4dmE+d2ZZcGNaZ2VtPnV4aXVbR2lZbnV3UXU8YWlKbnM/XFxVTnBxSGdqZndocVtiYWhiQHhDR2JIVmtrX25UUGVpb2JmeWNVZmBYbmF4aWRsd2lUSGptaGVGP3N3PnFXWHhUV3lnUWJ1cFp0WXBncXBrd3dmV3ZjSFpjQXdbaXVNaXliXm1FZnloX3l5WHNJSW9zWGRKZnh2cV0+eWFSX1pWeHlcXGJTP0ViQXdzXXddd3ZjT0ZvTWh3U1VSYWd5Q1lkaVR3QUJ1QUVHV0Z1U0lHb0VrS1lJR0ZZVVlddXdgdXdYb0d1QUZWV2tHd3F5ZmJAcXJyaWZqP3NZcHU9QF9db249Z1tRQGx0UWJRTlpEZlxcRldlXFx5cXV3WzxwdV4+bHZReFxcWXc8d1xcPFZ4UlBuPXl4aU5bQ05nQl5pck9wd0duRWZ5eVdudHF3Omd3RWZaU3BpX0dcXDw/YFFueFY/d3lnbTxOWl5xeWFHcHh4aU1wa19PaHFZcld4XFx0QHQ/QHZBQVxcZXFfclFxdj51eUB0eWFgV3l5Onh2bXlzWHd5WWZbTVd4b1dtSWd2b0U6O0I6TVRLV0RLV2dKO2VaMjpcIlx7XH0=