1. Бабешко В.А Динамика неоднородных линейно-упругих сред/
В.А. Бабешко, Е.В. Глушков, Ж.Ф. Зинченко.– М. Наука. 1989.
2. Гантмахер Ф.Р. Теория матриц – М. Наука. 1988.
3. Фаддеев Д.К. Вычислительные методы линейной алгебры/
Д.К. Фаддеев, В.Н. Фаддеева – М. 1960.
4. Демидович Б.П. Основы вычислительной математики/
Б.П. Демидович, И.А. Марон.– М.Наука.1970.
![](data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/wAALCAI7A8EBAREA/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/9oACAEBAAA/APf6KKKK5a08aSahb/abHwzrVzbF2VJkEAD7WKkgNKDjIPUVN/wk1/8A9Cjrn523/wAeo/4Sa/8A+hR1z87b/wCPUf8ACTX/AP0KOufnbf8Ax6j/AISa/wD+hR1z87b/AOPUf8JNf/8AQo65+dt/8eo/4Sa//wChR1z87b/49R/wk1//ANCjrn523/x6j/hJr/8A6FHXPztv/j1H/CTX/wD0KOufnbf/AB6j/hJr/wD6FHXPztv/AI9R/wAJNf8A/Qo65+dt/wDHqP8AhJr/AP6FHXPztv8A49R/wk1//wBCjrn523/x6j/hJr//AKFHXPztv/j1H/CTX/8A0KOufnbf/HqP+Emv/wDoUdc/O2/+PUf8JNf/APQo65+dt/8AHqP+Emv/APoUdc/O2/8Aj1H/AAk1/wD9Cjrn523/AMeo/wCEmv8A/oUdc/O2/wDj1H/CTX//AEKOufnbf/HqP+Emv/8AoUdc/O2/+PUf8JNf/wDQo65+dt/8eo/4Sa//AOhR1z87b/49R/wk1/8A9Cjrn523/wAeo/4Sa/8A+hR1z87b/wCPUf8ACTX/AP0KOufnbf8Ax6j/AISa/wD+hR1z87b/AOPUf8JNf/8AQo65+dt/8eo/4Sa//wChR1z87b/49R/wk1//ANCjrn523/x6j/hJr/8A6FHXPztv/j1H/CTX/wD0KOufnbf/AB6j/hJr/wD6FHXPztv/AI9R/wAJNf8A/Qo65+dt/wDHqP8AhJr/AP6FHXPztv8A49R/wk1//wBCjrn523/x6j/hJr//AKFHXPztv/j1H/CTX/8A0KOufnbf/HqP+Emv/wDoUdc/O2/+PUf8JNf/APQo65+dt/8AHqP+Emv/APoUdc/O2/8Aj1H/AAk1/wD9Cjrn523/AMeo/wCEmv8A/oUdc/O2/wDj1H/CTX//AEKOufnbf/HqP+Emv/8AoUdc/O2/+PUf8JNf/wDQo65+dt/8eo/4Sa//AOhR1z87b/49R/wk1/8A9Cjrn523/wAeo/4Sa/8A+hR1z87b/wCPUf8ACTX/AP0KOufnbf8Ax6j/AISa/wD+hR1z87b/AOPUf8JNf/8AQo65+dt/8eo/4Sa//wChR1z87b/49R/wk1//ANCjrn523/x6j/hJr/8A6FHXPztv/j1H/CTX/wD0KOufnbf/AB6j/hJr/wD6FHXPztv/AI9R/wAJNf8A/Qo65+dt/wDHqP8AhJr/AP6FHXPztv8A49R/wk1//wBCjrn523/x6j/hJr//AKFHXPztv/j1H/CTX/8A0KOufnbf/HqP+Emv/wDoUdc/O2/+PUf8JNf/APQo65+dt/8AHqP+Emv/APoUdc/O2/8Aj1H/AAk1/wD9Cjrn523/AMeo/wCEmv8A/oUdc/O2/wDj1H/CTX//AEKOufnbf/HqP+Emv/8AoUdc/O2/+PUf8JNf/wDQo65+dt/8eo/4Sa//AOhR1z87b/49R/wk1/8A9Cjrn523/wAeo/4Sa/8A+hR1z87b/wCPUf8ACTX/AP0KOufnbf8Ax6j/AISa/wD+hR1z87b/AOPUf8JNf/8AQo65+dt/8eo/4Sa//wChR1z87b/49R/wk1//ANCjrn523/x6j/hJr/8A6FHXPztv/j1H/CTX/wD0KOufnbf/AB6j/hJr/wD6FHXPztv/AI9R/wAJNf8A/Qo65+dt/wDHqP8AhJr/AP6FHXPztv8A49R/wk1//wBCjrn523/x6j/hJr//AKFHXPztv/j1H/CTX/8A0KOufnbf/HqP+Emv/wDoUdc/O2/+PUf8JNf/APQo65+dt/8AHqP+Emv/APoUdc/O2/8Aj1H/AAk1/wD9Cjrn523/AMeo/wCEmv8A/oUdc/O2/wDj1H/CTX//AEKOufnbf/HqP+Emv/8AoUdc/O2/+PUf8JNf/wDQo65+dt/8eo/4Sa//AOhR1z87b/49R/wk1/8A9Cjrn523/wAeo/4Sa/8A+hR1z87b/wCPUf8ACTX/AP0KOufnbf8Ax6j/AISa/wD+hR1z87b/AOPUf8JNf/8AQo65+dt/8eo/4Sa//wChR1z87b/49R/wk1//ANCjrn523/x6j/hJr/8A6FHXPztv/j1H/CTX/wD0KOufnbf/AB6j/hJr/wD6FHXPztv/AI9R/wAJNf8A/Qo65+dt/wDHqP8AhJr/AP6FHXPztv8A49R/wk1//wBCjrn523/x6j/hJr//AKFHXPztv/j1H/CTX/8A0KOufnbf/HqP+Emv/wDoUdc/O2/+PUf8JNf/APQo65+dt/8AHqda+K/N1mz0y70TU7CW8D+S9yIijFF3EfJIxHHtXQ1ieKfEkXhjR5b57Se7dVLLDCvUDqS3RQB6/hmpdT1DVYJo4tL0cXpKb3kluRCi+wO1iT7YHHese38Z3OqSrbaRo4uLxI2e5jnuvKWPbI0ZCttbd8yN2HGD3xRa+NLnVnWDRtHFzdJGZLmOe68lY8O0eFba275o37DjB71u6Hq0WuaRDqEUbRrIXVkbqrIxRh7jKnB71kfDz/kSLH/rpcf+jnrqKKKKKKKKKKKKKKKKKKKKKKKKKKKKz9V1m10hI/OEsssh/dwQJvkcD7xCjsByT/UgG1aXdvf2kd1azJNBKNyOhyCKmoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omuprB8bWdzqPgrV7Ozhaa4mt2WONerH0qr4o1fVbEW1jpumahL565mvLWFZPIX0Clhlz27DrWRYK+hawmq2ugaqbG4sFtxCkQadZVlkd2kBb+LcDnJySaq+H7PU/Ct8+oXekXtzHeQFVjs0EjxkzyyYcEjHyyL6859K2dD8IWdz4WsbXX9NimmhknlRJOTGJJWcDjvgjPvU3w4RY/AmnxooVEedVA7ATOAK6qiiiiiiiiiiiiiiiiiiiiiiiiiiiqGqaRZ6xAsV3Gx2NuR0Yq6+uGHIyMg+oNW7e3htbeOC3iWKGNdqIgwFFSUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/wB67/8ARNdTRRRRRXL/AA7/AORHsf8ArpP/AOjnrqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/wCia6miiiiq97am9s5bcXM9sZBjzYGCuv0JB/lXO/DhdngTT03M215xuY5J/fPyfeuqoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/wCR68Jf713/AOia6miiiiiuX+Hf/Ij2P/XSf/0c9dRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu/wD0TXU0UUUVBeJdSWki2U0UNyR8kksZkVT7qCM/mK5z4cBh4E08OQzh59xAwCfOfJA7V1VFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/ACPXhL/eu/8A0TXU0UUUUVy/w7/5Eex/66T/APo566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/AJHrwl/vXf8A6JrqaKKKKr3stzDZySWdstzcKPkhaTywx/3sHH5Vzvw4LHwJp5dQrl5yyg5wfOfIz3rqqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/6JrqaKKKKK5f4d/8AIj2P/XSf/wBHPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/3rv/0TXU0UUUVXvbsWNnLctDPMIxny4Iy7t9FHJrnfhw2/wJp77WXc85wwwRmZ+COxrqqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/wCia6miiiiiuX+Hf/Ij2P8A10n/APRz11FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/ACPXhL/eu/8A0TXU1z3jTWtR0Hw5dX+m2cc8sUbOzzNiOMDuQOST2A79SKj1DUtWuvEi6Lpdxa2bJZLePPPbmbducptChlxjbnOT16Vk6X4k1zxFeHTrOazsZ7aJ2mnkt2mWUrM8XyrvXb/qy3U/ex2yTSfEuueJbk2VlNZ2E1tCXmmlt2mEpEskXyrvXb/qiep+9jtk7Ok+KrSXwvb6vrFxbWCvLJAzySBELo7JwT67Ccdveq/w4dZPAmnujBkZ52VgeCDM+DXVUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/3rv8A9E11NY/irS59b8K6lplq0az3MBjQyEhQT64qlqGm6ta+JF1rS7e1u2eyWzeCecwhdrl9wYK2c7sYwOnWsrS/DWueHbw6jZRWV9PcxOs8Etw0KxFpnl+VgjbuZCvQfdz3wDSfDWueGrn7bYxWV/NcQlJopbhoRETLJL8rBG3f60joPu574HReGtHfRvD9vYXDpLIrySOQOAXdnIHrjdjPfHas/wCHf/Ij2P8A10n/APRz11FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/ACPXhL/eu/8A0TXU0UUUUVy/w7/5Eex/66T/APo566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/AJHrwl/vXf8A6JrqaKKKKZKnmxPGWZd6ldyHBGe4PY1zHw4XZ4E09AWba84yxyTiZ+Se5rqqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/6JrqaKKKKK5f4d/wDIj2P/AF0n/wDRz11FFFFFFFFFFFFFFFFFFFFFFFFFFFFNaREKh3VS52rk4yfQflTqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/wCia6miiiiq97FczWcsdncrbXDD5Jmj8wKf93Iz+dc78OAw8CaeHbc4ecMwGMnznycdq6qiiiiiiiiiiiiiiiiiiiiiiiiiiiua8ReGZtZvIZ4rtVAXyysoJ8nnPmRYxiTjGT2Ptg9FEnlQpGXZ9ihd7nLNjufen0UUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/3rv/0TXU0UUUUVy/w7/wCRHsf+uk//AKOeuooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooqC8e6SzkayhimuQP3ccshjVj7sASPyNc58OCx8CaeXAVy8+4KcgHznzg966qiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6miiiiiuX+Hf8AyI9j/wBdJ/8A0c9dRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/AHrv/wBE11NFFFFV725NnZy3C209yUGfJgUF2+gJA/Wud+HDb/AmnttZdzznaw5H75+D711VFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv8Aeu//AETXU0UUUUVy/wAO/wDkR7H/AK6T/wDo566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/8AomuporjfFWv3Ona1DZHVoNGszbiUXktuJvNk3EGMDIxgAHPvVPS9c8QeIruOwgvYdPeK1+0G5FsJftKmV41baSNoIQN3+9jtSaRr+veJ7o2NtfwadJawl5ZhbCXziJpIvukjb/qs8Z+9jtWppnjWwTw7ZX2uXcFpLO80QznDmORkJHpnAOO2aX4cOsngPT5EYMjvOykdwZnINdVRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu/wD0TXU0Vzmp6Vq0fiQa3pH2KWZ7RbSSG7do1Chy4YFVYk5bGCKytO8K61oF19v02TT7m6uI3W4iuXeNELSvKdjKpLcyEcgcAHvRpvhXW/Dtx9s0ySwurieEpcR3TvGiEyyS/IVViRmUjkDhQe9dH4b0htD0KCwklEsiNJI7KMDc7s5A9gWx+FZnw8/5Eix/66XH/o566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6miiiikJCjJIAHc1zHw7/AORHsf8ArpP/AOjnrqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/wCia6miiiioLuztr+1ktby3juLeQYeKVQysPcHrXOfDlFj8CaeiKFRXnVVA4AEz4FdVRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/wAj14S/3rv/ANE11NFFFFFcv8O/+RHsf+uk/wD6Oeuooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/wCR68Jf713/AOia6miiiiq97ai+s5bZpp4RIMeZBIUdfow5Fc78OF2eBNPQFm2vOMsck4mfqe5rqqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv8AkevCX+9d/wDomupooooorl/h3/yI9j/10n/9HPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/3rv8A9E11NFFFFV72K5ms5I7O5W2uGHyTNH5gU/7uRn86534cBh4E08O25w84ZgMZPnPk47V1VFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv967/9E11NFFFFFcv8O/8AkR7H/rpP/wCjnrqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/6JrqaKKKKgvHuks5GsoYprkD93HLIY1Y+7AEj8jXOfDgsfAmnlwFcvPuCnIB8584Peuqoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooorl/h3/AMiPY/8AXSf/ANHPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/wB67/8ARNdTRXCeJ9Qd/GK6VLJq7QLp63EVvpLtHK8hkZWJYEcBQOCfpUOmufEWtrpo1rVGsbawWVWjnMEwk82RGSRkwWK7AvOeQTznNVfDc994p1CSw1DVL9IbO3JjNrcNC7H7RNHl2Ugsdsa9e+T3rf8AD2vXY8F2t/c217qU/nzQlbZFaQqkrqpOSM8KMnPJp3w4bf4E09trLuec7WHI/fPwfeuqoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/wCR68Jf713/AOia6misLVfD9xd6quqadqkmnXvkC3kcQrKHjDFgNrcA5J5rPi8FSaeyT6RrE1leMjJcTmFJfOy7SE7W4X53Y8euO1EPgqTTikujaxNY3RjMdxKYUl8752fO1uF+Z3PHrjsK3tG0qHRNKhsIGZkjLMWbqzMxZj+LMTjtWN8O/wDkR7H/AK6T/wDo566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/8AomupoooopksiQxPLIwWNFLMx7Ada5j4cOsngTT5EYMjvOykdwZnINdVRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu/wD0TXU0UUUUVy/w8/5Eix/66XH/AKOeuooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooorl/h3/AMiRY/8AXSf/ANHPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/wB67/8ARNdTRRRRUF5aR31pJbStKscgwTDK0bj6MpBH4Guc+HCCPwJp6AkhXnUbjk8TP1Peuqoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooorl/h3/AMiPY/8AXSf/ANHPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8j14S/wB67/8ARNdTRRRRVe9hnuLOWK1umtZmGEnVA5T3w3B/Gud+HAZfAmnqz72DzgtjG4+c/OK6qiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6miiiiiuX+Hf8AyI9j/wBdJ/8A0c9dRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/AHrv/wBE11NFFFFV703gs5TYLA11j92J2KoT7kAn9K534cb/APhBNP8AMCh98+4L0z5z5x7V1VFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv967/APRNdTRRRRRXL/Dv/kR7H/rpP/6Oeuooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupormfGmq3Nppqadp0/k6nqCukMoxmFVXLy4PB2jHHGc8VgaEj+Ip9IsdRvL14E8PWl3+6u5ImaVywZmZGBYkKOpNV/DBm8T6lLYare3slvZ27GERXUkLZ+0TR5ZkILfLGg+YnpnqTW54d13UV8J2M7abfatK0txGzQNECqpKyruLuueABnk8c1P8OGLeBNPYqULPOSrdV/fPwcV1VFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv967/APRNdTRWRrHhjRNfO/U9MtbqURmNJZYlZkB9CRx1rMtvA9vplvZR6LqFzpkltaraNLDHGxlRSSNwdSudxYkgDJNOj8FxWJjk0nVLzT7jyzHNNGschmG9nyQ6sB8zuflA+96AVt6Tpdvo2mRWFru8uPccsclmZizE/ViT6c1ifDv/AJEex/66T/8Ao566iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6miiiiq97ewadZy3d0zJDEMuyoWI7dFBJ/Cud+HDrJ4D0+RDlXedlOMZBmfFdVRRRRRRRRRRRRRRRRRRRRRRRRRRRTXdI42kkZURQSzMcAAdzSo6yIrowZGGVZTkEeopaKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/6JrqaKKKKK5f4ef8AIkWP/XS4/wDRz11FFFFFFFFFFFFFFFFFFFFFFFFFFFFcb4p1a0vi+l7pUSGaNZrgjMSynlI5E6urd8cDj8N/QdJOjaWtq0xlYu0jYGEQnkqi/wAKjsK06KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/wCia6miiiimSyxwxtJLIsca9Wc4A/GuZ+HRB8DWJBBBknII7/vnrqaKKKKKKKKKKKKKKKKKKKKKKKKKKKqSaZYzahDqElrE13CpWOYr8yg/5/DJ9TVuiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/AJHrwl/vXf8A6JrqaKKKKiubaC8t3t7qCOeBxh45UDK31B4Nc38OVVPAtgiqFVXnAAGAB5z8Cupoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooorl/h3/wAiRY/9dJ//AEc9dRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu//RNdTRRRRUF5aR31pJbStKscgwTDK0bj6MpBH4Guc+HCCPwJp6AkhXnUbjk8TP1Peuqoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooooorl/h3/wAiPY/9dJ//AEc9dRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu//RNdTRRRRVe9hnuLOWK1umtZmGEmVA5T3w3B/Gud+HAZfAmnqzb2DzgtjG4+c/OK6qiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6mioLm8tbNVa6uYYFY4UyyBQfzpJr+zt4o5Z7uCKOQZRnkChvoT1omv7O2jjknu4IkkGUaSQKG+hPWp1ZXUMpDKRkEHIIrmPh3/yI9j/10n/9HPXUUUUUUUUUUUUUUUUUUUUUUUUUUUUVg2viuxutTW1VJRDLI0MFyR8ksq/fQDqMepGDz6c71FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv8Aeu//AETXU0VheILfQowt/q1hFeTbRDDG8Xms5zkKiHI3E+g+vFcp4c0BV8VS6drVrb3ITThcR20yiVLfzLiYhQDkAhcLxx8oxwBVTwNZ2uqa3d2+oW0N5Db2zeVFcIJFi/0u4X5Q2QvCgcdlA7Vt+FZNc/4Qyw/shLCUie6V/tsrrtAncKF2qeABjHbAq98ON/8Awgmn+YFD759wU8Z85849q6qiiiiiiiiiiiiiiiiszWNah0kQp5MtzdTkiG3hxvfHLHngAD19vWrWn39tqlhDe2knmQSjKtjHsQR2IIIqzRRRRRRRRWRf6Ijre3emLDaatcRbBd+Xkj/9fTPsOuAKr+FtP1GwtJlvmaONnzDbPL5zRf3iZO+4847fjgb9FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/I9eEv967/APRNdTRWPq3hbSNcvILzULeV7i3UpFJHcyxFAeuNjDrUUvg/RJoIIXt58QgqjreTLJgnOC4fcwyTgEnGeKdceEdDuYYIms2jEC7ENvPJCxHXDMjAtzk8k8knqTWpaWlvYWkdraxLFBGMKijgf/X7571z3w7/AORHsf8ArpP/AOjnrqKKKKKKKKKKKKKKKKKqajplnq1r9mvrdZodwYK2RgjuCOR+HqR3qzHGkMSRRIqRoAqqowFA6AD0p1FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFctr3/ACPXhL/eu/8A0TXU0UUUU2RikTuqNIyqSEXGW9hniuX+HDFvAmnsVKFnnJVuq/vn4OO9dVRRRRRRRRRRRRRRRRRRRRRRRRRRUF5eW+n2ct3dzLDBENzu3QVV0rWrbVhIsaSwXEWPMt502SID91iPQjkH6jqCBo0UUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVy2vf8AI9eEv967/wDRNdTRRRRRXL/Dv/kR7H/rpP8A+jnrqKKKKKKKKKKKKKKKKK5vxJ4kudFuoIobRWDLv3zbsTHOPKj2g/vDnPPp35K9DDIZYY5DG8ZdQxR8blz2OMjIp9FFFFFFFMlijnheGaNZI3BVkYZDA9iKpWWnafoNjKLWExwrmSQjc7HA/EnAAAHoABVHQ/FEGtXLQC3eFmj8+Alg4eInAY7fuNkH5T+vIG9RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/AHrv/wBE11NFFFFI7rGjO7BUUZZmOAB6muX+HLK/gWwZWDKzzkEHII85+a6miiiiiiiiiiiiiiiiue1DxVZWeq/ZJLeSSKGREnnwP3Mj/cAU/M2eclRx7846GiiiiiiiiszX9Vk0fSnuo4DK24ICeEjz/G5HIUd8f/XpNA1WTWNLW5lh8tg5QsnMcuP44yeSp7Z/XqbltY2lk0zW1vFCZ38yUooG9vU+9WKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK5bXv+R68Jf713/6JrqaKKKKRlV1KsAykYII4Irl/h0AvgaxAAAEk4AA4H7566miiiiiiiiiiiiiiikZgilmICgZJJ4AryyfxDcN4mfUjbrcy2zFI5pYwtr9nY4WSJvvHDZDtyB9Cpr0DQNWbWNMFy8PlurtGzKcxyEcboz3U9jWpRRRRRRRWF4m1t9HtoVjgiZrhinm3LbYEAGSHb1I4A71oaReLqGk2t2ttJbLLGCIZFwUHpj09PbFYun+K5L3X2sGsmWNmZFUHM0RXq0q/wAKNkYOf58dPRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRXLa9/yPXhL/eu/wD0TXU0UUUUVy/w8/5Eix/66XH/AKOeuoooooooooooooooorhfiH44i8NfZtMjJFxdo0k0ip5htYBwZCnVueB+PpVnwbowjsxqdzCiQz26LbwSsZGijxknc33Q2c7ccd/QdPp93ZXtmsmnzRS2ykopiI2jbxgY9KtUUUUUUUVx/i7xFNp2pWelpHbs1yhkiimjMhu3Uj9yAPunvuPpx0OOthZ3gjeSMxSMoLRkg7TjkZHBxWHN4l0211iSEwttDrb3F4qDakh+5Gx6nqeegyPWt+iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/8Aomupooooqve2UGoWclpcqzQyjDBHZD+akEfga534cIsfgPT41BCo86gE54Ez966qiiiiiiiiiiiiiisrxBrsGg6RcXzoZ5IxiO3jI3yvjIRfc1474Y0U/FTxIviq5ldbFJV+22syZHmR/djTPKgZySpweCeSQPa9VtLW+0m6tr1itq8Z81g5XCjnOa5zwfPp6XE6LLcPd3SiSOaaIRLcwqNqsiLxgdD3PB6FQOvooooooo6DJrzKLHiDxmzR/NYzzfv7WI5ceUSFn87oEJ3KVQ9V453ivTa5HxDqmnaV4ggupNLjlu4It73EjBD5ZO3Eef8AWOCeg6AkZBbB6yOQSxJIoYB1DAMMHn1FOoooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooornPF/iF9AtrPbLb2wupvJF1cfMkbYyF2jklsYHYdz0B0dC1ZdY0m3uyqRSyKS8QcMVIOD+o78jvg1pUUUUVy2vf8AI9eEv967/wDRNdTRRRRRXL/Dv/kR7H/rpP8A+jnrqKKKKKKKKKKKKKoanrelaJEsuqalaWSOdqtcTKgY+gyeTXj/AIw/aBsrf/Q/CVu17clwDcTxFY8dwq8MT25AH1ra0bwXq+uw2+papvshdHzJ47g7rgAgdRjGT056ADIxhF9NsrK2060jtbSJYoUHCj9ST1JJ5JPJPWrFeXQXEXhPxXeXUsqW6SXiRXj5wsiPnyljhOCoDbslSQSWPzZOz1GiiiiiiuH8d+KEg0+60fT7kpeyxlJJ0AYW4I9DwWwRx2B5IyM5/wAJvBd94d0p7vVwFuZPktbZXZo7aLHVVflGc5Zhn0yARgekVznja1im8OyXMq7/ALGwuBH5e5ZCONjY5VSCQWBGAcnK5Bh8CXMj6PLZ3DE3VrJiRUfzIY9w3BI37qBxgnI9gRXU0UUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUm9c43DPpmk8xP76/nRvTGd649c0nmx/31/Ose48Y+GLR1S58Q6XCzIHVZLtFJU9DyehpIPGfhe5LiDxFpUhjQyPsu4ztUcknnoPWov+E88If9DPo/8A4Gx/41KfGfhcWq3R8RaULdnMYl+1x7SwAJGc9cEfnWJc/F7wJa6gbKTxBA0oYLujR3j5x/GoK45654qCb40+AIZpIm10FkYqSltKykj0IXBHuKcvxm8AtavcDXk2IwUqYJAxJz0XbkjjqBxx61F/wu74f/8AQcP/AICTf/E1k6j+0B4T0+6+zi31C5YKCzW6xsqn0zvwT9MioYP2ifCErsslnqsICMwZ4UIJAyF4Y8np6euKiP7R3hXHGmaxn/rnH/8AF1iS/tLIIITF4YZpSD5oe7wo542nac8deBzVK6/aR1S5jWLTfDcEVwzjBkuGlBGeQFCryfrVC7+OHj3UNRm0u00yG0nllMSRw2jyXEXPQBjhmHuv4Cug8MeCfF3jO4vJvGLXY026DKX1BU84rkBfKix/o7jbkt0ORwetUNV+DnjDRdSe90S7a/USMiywX72168RUBVd2+XapVeB1/lT1e++Nfh+9FjGL+bEas01nbfaUlOOpYq2G4wQMDjOOcmLTvip8UPDt15OtaVcXrXCO8Ud5ZNE/yqSSu0DIHBPB4HUVX/4aF8a/8+Ok/wDgNJ/8XXrHgj4m3Wv2umrq2ktbXNyNpaLJDMfusE5IUjJ65HX7oLD0muW17/kevCX+9d/+ia6miiiioLyCW5s5IYbqW1kcYWaIKWT3AYEfmK5z4cKV8CaerOXIecFj1b98/JxXVUUUUUUUUVDdXdtY2z3N3cRW9vGMvLM4RV7ck8CuavPiZ4KsZpop/EdjvijEreW5kG0nHBUEE5PQZPfFcFrf7RGl6fLCth4fv7hJE377phb7lP3WUYbcp554rltZ+K/jzxZdJo/hnTbi3E0anzILYrLIjcb8kkRqc8Nnj+9VvRPgv4n1/WJNU8TahLaQuCnl3cq3l0U6Y3kbVOOjDJHoK9X8K/Dbw/4UupL63jmvNTl3eZf3snmStk5PsPTIAz3rr6KK4f4meH7rU9C/tLSoi+q2Ct5YWNXd4nGJUTdkKxXo2CRjA61k/CTxdHeaePC1/chtXsI96pvMjCDPCu/QuuQDjA5HuB6dRRRWNr/izQfC8Il1rVLezDYKo5y7AnGQgyxHuBXlHin45STaxbaN4Z0u9kE52vM8RSaQN90wqfUYIJHIPQda7rwp4KFmY9R1VA91nfHAeREeuW65fJJ6nBJ5JJY9tRTJolngkhcsFkUqSrFTgjHBHIPuK8E03WtS8B/FgadqM7zQS4tXjRQA1sctFMqKAqhCGDMzE4JOK99VldFdGDKwyCDkEUtFFFFFFFFFFFFFFFFFFFFFFFFFFFFZreIdFS7ltG1jT1uYgxkhNym9AoycrnIwOTVb/hMfC/8A0Mmj/wDgdF/8VUv/AAlPh77N9p/t7S/I3+X5v2yPbuxnbnOM45xUR8Y+GsfJr2nTN2jguVkdvoqkkn6CsOf4v+Are4khk8RQ742KtshkYZHoQpB+opU+LvgN7eWdfEUGyIqGzFIG5zjClcnp2HHeov8Ahcvw/wD+hij/APAeb/4ioNX+NfgjSbSGdNSe/MwDLFZxlnCnPJ3YA6dCc8jisq3/AGhPBc0ux49UgXaTvkt1xwCccMTk9Bx1NRf8NFeDf+fTWP8AvxH/APF1K/7QngtLeKUR6o7Sbt0a267o8euWxz2wT+FZOpftJaNDLGNM0O9uoyPnaeRYSp9gN2fzFVLj9pa0W4kW38MzSQhiEd7sIzD1I2HH5mlj/aWszbzNJ4amWYbfKRbsFW55ydvGB7HPtUX/AA0xF/0Kr/8AgcP/AI3WNe/tI+IHu5GsNG02G2ONiT75HHHOWDKDz7Cl0/8AaI8VXF6kbaHp9yCG/dW8cgdsAnj5j06njoDVb/ho7xX/ANAzR/8Av3J/8XVmX9onxUthby/2Hp6F2cGZ0k8uTGOF+bqM88nqOlQWnxK+Lur6jBd6fYXL2txKpihj039wwJ6byudp9d341FL48+NAmcG21RcMflXRwQPYHy6nt/HPxpayvGFhqDhVXMj6UFeP5h90bBuz0PBwOeOtVf8AhPfjR/zw1b/wTD/43V7U4PjV4lNteJ9vgIgVWhtZvsu3jOXBKgs3U4zjOOMYEdj4X+OC/acXeqR5gYHz9SVt3T5V+Y4Y9jx9RVX/AIRb43/89dc/8Gq//HKtP4X+N/8AZkK/a9UIErny11JRKvC8s27lT2GTjB4GeWWnwQ+IOqXcGqahqdvbXUkivJJPdu88ZB+8SAQWAGR83pzSXH7PnjWS5lf+09Ml3OT5j3EgZ+ep+Q8n61MvwA8ajS3g/tmwGZlf7MJ5PLbAI3k7fvDOOnQnmqo/Z68a5/4/9J/8CZP/AIiunvv2cZb24Fw3iyR5nRTM89t5jNJj5iDvHGegP50yD9mlVMnn+J2YGNgnl2m3D44Jy5yPUcZ9RUP/AAzPL/0NSf8AgCf/AI5Ux/ZpX7IqjxO32gOSW+yfJtwMDbvznOec/hWtp37OHh6GCE6jquoXM6nMnlbI4356YwSOOOtWJv2dPCUk0jpf6tErMSI1lQhR6DKE4HvTl/Z28IC1eI3mqtIzBhMZk3KBnIA24wc9xniov+GcfCn/AEE9Y/7+R/8AxFdLJ8GPAUyQrJoa5ijEYZZ5ELAd22sMn1PWlg+DPgG2dnTQUYsjIRJPI4wRg8Fuvoeo7VEPgl8Pwc/2Gf8AwKm/+KrpX8F+F5IIYJPDulPFACIka0QhATk444yealsPCvh7SroXWn6Hp1pcAECWC2RGAPUZArW2LnO0Z9cUtFFIQCckA4pPLj/uL+VV4dNsbe5a5gs4I523bpEjAY7jlufcjNWq5bXv+R68Jf713/6JrqaKKKqX2q6dpaodQv7W0EhIQ3Eyx7sdcZIzVlHSWNZI2V0YBlZTkEHoQa5n4d/8iPY/9dJ//Rz11FFVrrUbGxDG7vLe3CrvPmyqmFzjPJ6ZIH1NYmp/EDwlpESSXniCwVXbavlyiU5+iZIrI1P4xeCNM06K9/tdbsS/citkLSEA4JwcY/HHtXM3X7RfhaOWAW1hqM8bMRKzIqGMeoGTu+nFczqf7Sl26x/2V4ehibJ8w3UxkBHbG3bjvWZP8aviLr0KSaLpawxxsVkks7Jpgx44JbcBj29ahin+NHi6VLpBqaw7hA6ELbRH13R8bhzycHI4rQu/gP4z1XXHbUNbs3tpSDJcb2JHy9BHjoCAAM9BXR6T+zho8NjMmravdXN04wkluojWPkHIBzk8Ec8c9K7jS/hL4J0uFo10K3uiwALXmZyMem7OPoMV2SQxRACONECqFG1cYA6D6U+iiiiivNvEPgVdHuJtY8OQ+QXl8+eKAYdX/vrjr3yO2TjILK0fhb4w6VqKvFrObFElFvFfyjZFdSdCAOqkd+3rt4FekvPDFGskksaRsQAzMACScAZ9yRj615x4t+N/hfw1JJbWrPqt9GxVorc4RSDgguePXoD0rzq5+JvxK8W6PG2g6VPBHdTm3M9pFkBiMBUY8r6lieD3HSrvhb4GeIL/AFBNb8U63Pa3LN5u2CTfPuIzkucgHJOev1r2Hwz4G8PeERMdHsBFJMQXlkdpHOOg3MSQPauiooorC8R+GLTX4NxCw30a4iuAOR32nuVyAfYgEYIBryS/+I2ufDDUo9LvtDnuNOB+ZpDtQEg4EL4xjuQR9AOc+geH/i54N8Q+UkWqpaXMmcQXg8thgZ6/d/Wu4oJAGTwBWQ3ivw6ilm17TMAZP+lx/wCNYkvxX8CwxwyP4jtdsy70wrk4yRyAOOQeDVC6+NfgS2mgjXVzOJd2XhhYrHgZ+bIHXoMZrJ/4aF8Ff889V/8AAZf/AIqpW+P/AIIW2jmDagzOzAxC3G5MY5PzYwc8YJ6Gov8AhoXwV/zz1X/wGX/4qpbn4/8Agi3uHiV9QnCnAkitxtb6ZYH9Kfa/Hrwbd71iXUzMNojhFrl5STjCAE5POe1c5dftJ6dFdSxweHLqWJXIR3uAjMPUrtOPpmkn/aV01HUQeHbqRSili1yq4bHI+6cgHjPeltv2lNMe4Rbnw7dRQk/M6XKuV/DaM/nUD/tL22xtnhiYPg7d12MZ9/lrNu/2hvEFtbRldG0uVi5BuIzKYW+VTtXJB3DPPPp9TXt/2jvEcl1EjaHpsis4BSISb256L8x59ODST/tH+I0uJUXQ9NjVXICSCTcvPQ/MOfwpI/i/8TdeivL3R9JX7KkQUm2sjIsTBhlgTnc3OMc8HOO9Uv8AhZfxg/5977/wUD/4irV98R/jDG0GdOuYN0CNhNMDb8j7x+U4J9OMegqTSfiB8YL27aA2suGifMlxpvlpFx9/IXOR1xzn0NUf7G+OJ583Xv8AwPH/AMXVo6N8cv7KQfaNWx57fIL1fN+6OSd33fQZ654qGDRvjl9oixPrY+ccyXylRz3G7pVjUfhx8UfFt0w1G4uFa2LEtfXwMTuzHJhCg7VwF4I/wDbT4FePoo7pf7XtIPMhK7UunIl+YHY3AwOM9+gqr/woXx9/z+WH/gY//wATVq9+Bfj6YW3/ABN7S52QKn7y6ceVjPyDg5A9eOvStHTv2b76SNZdT8QxxzPG2+OCEttcg4+YnkZ68DIzR/wzRef9DPB/4Bn/AOLqb/hmqX7Fs/4SZftHmZz9lOzbj03ZznvnpS2X7NkkV9DJdeJVe3VwZFhtijle4B3HB98Vuz/s4+F5LiSSLU9UhjZiVjVkIQegJXJx70qfs5eFltpY21HVHkcqUlLoCmM5AG3Bzx19Ki/4Zu8N/wDQY1b84/8A4mtG0/Z98G288ck5vrpFh2NHJPtDvnO8lcHOOMA49q1Lf4JeAbabzF0ZnO0riS4kdeQR0J688ehqL/hRfgD/AKBM3/gZL/8AFVNJ8E/AUlvDCdGIEW7DLcSBmz/eOcn2z0q7pfwm8D6SjrD4ftZt7Bs3QMxH0L5x+FaVx4C8I3dxJcXHhvS5ZpGLO7WqEsfU8UsfgTwlFbzW8fhvS1hm2+YgtUw+DkZ47VF/wrrwX/0K2kf+Aif4Vuw6fZW8McMNpAkUahERYwAoHAAqRbaBG3LDGpHcIKT7Jbf8+8X/AHwKU20BUKYY9o6DYMCpFUKoVQAo6AClooooooooooooooooooooooooooooooorlte/5Hrwl/vXf/omupooork0t4Lv4l6lHcQxzImlWxVZFDBSZJuQD0zj9KpeFP7cbwZY/wBjSaeCLi7En21XbgTuAF2kYxz+lX/hxv8A+EE0/wAzbv3z7tvTPnPnHtXVV4N4+8AfELxd4ia4hvVEUI2pH5xiiUHvHjrnHOeR6ms3TPgP4rv9Ku7bWteS2+ZXgiVzMjMOpbpjjOPerumfs1qJX/tXxCWj2/ILWHac++7PFdDafs8eE4JLdri5v7lY1YSK0gUSkng8dMe3Wum0P4R+C9AvPtVrpKyzcbWuXMuz3AbpW3png7w3ozSNp2h2Ns0oAcpCOcdP51sRQQ26bIYkjUnOEUAZ/CpKKKKKKKKKKKK4Dx98KNG8cRxzb2sdQiVgk8IG1sg4DL3GcHI5rzTRPhD47TVrC2vNZ8qxtJG8ws5kjVG6hFPD7gMEHGOK7vQPgZ4V0y7e8v7YX0zMWEMhJhjyCMBT1HP8Wa9Mihit4xHDGkcY6KigAfgKfRRRRRRVPVNJ0/W7FrLU7OG7tXILRTIGUkHINeWav8BdNZSdC1GS0QTJOlndIJ4C44JYHk5GeM1w3/Cq/iN4Y1M39neXFxDFJtDadd7ZnjJxlVbgcdj0qdPhR8R9T0qWWW+FrcXQWS4Sa+dnuSXLfOOilQeg61Mv7NWo7hu8R2u3POLds4/Oukl/Zw8OvHCserajGyLh2+Q+Yck5xjjggcelaGh/ADwnpVy81691qmRhEuGCqvBB4XGevfpirn/ChvAf/Phdf+Bb/wCNSv8AA3wG9tHD/Zkw8tmbeLhwzZxwTnkDHH1NRf8AChvAf/Phdf8AgW/+NS3PwN8B3Nw8x0yWPec7Irh1UfQZ4p0HwQ8DW6SLHp0+5wMSG5fcmDnKnPB4rVt/hb4It7aOEeG7GTYoXfLHvdvck8k1Zn+Hng65ZWm8N6a7IixqTAOFUYA/ACltvh54Ps7hLi38N6bHKhyrCAcVGnw28Fxurp4Z00MpyD5A4NdCbCzMYjNpBsByF8sYB/KkXTrFGDLZ26sDkERLkH8qG06xZizWVuWJySYlyf0qaKGKBNkMaRpnO1FAH6U+iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiuW17/kevCX+9d/+ia6miiisvVvDmka40b6jZJO0YIVtzKQPTKkZq/bW0FnbR21tEsUMY2oijAArhfCutz6D4eh0y88P66Z4JJtxisS6nMrsCCDyMEVs/wDCZJ/0L3iH/wAFzf41Db+PrO7877Nouvy+TIYpNunt8jjqp9xmpv8AhMk/6F7xD/4Lm/xo/wCEyT/oXvEP/gub/Go4fHVvcKzQ6F4hcI5RiNOfhgcEfgak/wCEyT/oXvEP/gub/Gj/AITJP+he8Q/+C5v8ajj8dW8zypHoXiFmibZIP7Of5TgHH5EH8ak/4TJP+he8Q/8Agub/ABo/4TJP+he8Q/8Agub/ABqNPHVs9xJAuheITLEFLr/Zz/KDnH54P5VJ/wAJkn/QveIf/Bc3+NH/AAmSf9C94h/8Fzf41B/wn9n9v+w/2Lr32vyvO8n+z23bM43fTPFT/wDCZJ/0L3iH/wAFzf40f8Jkn/QveIf/AAXN/jUMnj6ziuobWTRdfW4nDGKM6e2XC43Y+mR+dTf8Jkn/AEL3iH/wXN/jR/wmSf8AQveIf/Bc3+NQz+PrO1kgjn0XX43uJPKhDae3zvgnaPfCk/gam/4TJP8AoXvEP/gub/Gj/hMk/wChe8Q/+C5v8ahufH1nZRpJdaLr8KPIsSl9Pbl2ICr9SSB+NTf8Jkn/AEL3iH/wXN/jR/wmSf8AQveIf/Bc3+NQ3nj+z0+0ku7zRdegt4hl5H09gqj3qb/hMk/6F7xD/wCC5v8AGj/hMk/6F7xD/wCC5v8AGmT+Obe2t5J59C8QRwxKXd209sKoGST+FEPjiC4gjnh0LxBJFIodHXT2wykZBFP/AOEyT/oXvEP/AILm/wAaP+EyT/oXvEP/AILm/wAahtPH1pf2kd1aaJr81vKu6ORNPbDD1FTf8Jkn/QveIf8AwXN/jR/wmSf9C94h/wDBc3+NQ2vj6zvYmktdF1+ZFdo2ZNPbAZSVYfUEEfhU3/CZJ/0L3iH/AMFzf40f8Jkn/QveIf8AwXN/jUMHj6zunnS30XX5Gt5PKlC6e3yPgHaffDA/jU3/AAmSf9C94h/8Fzf40f8ACZJ/0L3iH/wXN/jUSeO7aS4lgTQvEJliwZF/s5/lz0/lUv8AwmSf9C94h/8ABc3+NH/CZJ/0L3iH/wAFzf41H/wndt9q+y/2F4h88J5mz+zn+7nGfTrUn/CZJ/0L3iH/AMFzf40f8Jkn/QveIf8AwXN/jUL+PbSO7htX0TX1uJlZo4zp7ZcLjcR9Nw/Opv8AhMk/6F7xD/4Lm/xo/wCEyT/oXvEP/gub/Goz46txcrbnQvEPnMhcJ/Zz8qCAT+ZH51J/wmSf9C94h/8ABc3+NH/CZJ/0L3iH/wAFzf41HJ47top4oX0LxCskxIjX+zn+bAyf0qT/AITJP+he8Q/+C5v8aP8AhMk/6F7xD/4Lm/xqKbx3bW5jE2heIU81xGmdOc7mPQcfSpf+EyT/AKF7xD/4Lm/xo/4TJP8AoXvEP/gub/Gobrx9Z2NuZ7vRdfhhDKpd9PbGWIVR+JIH41N/wmSf9C94h/8ABc3+NH/CZJ/0L3iH/wAFzf41Fc+PLWytZbq60TX4beJS8kj6e2FUdSalHjSNgCPD/iEg8j/iXN/jR/wmSf8AQveIf/Bc3+NNk8bQxRtJJoHiFUQFmY6c2AB1PWm2/ju2u7aK5t9D1+WCZBJHIunth1IyCPYipP8AhMk/6F7xD/4Lm/xo/wCEyT/oXvEP/gub/GoLPx/Z6hapdWei69PA+dsiaexBwSD+oIqf/hMk/wChe8Q/+C5v8aP+EyT/AKF7xD/4Lm/xqC28f2d6JTbaLr0wilaGQpp7fK6nDKfcGp/+EyT/AKF7xD/4Lm/xo/4TJP8AoXvEP/gub/GoYfH1ncT3EEOi6+8tuwSZF09sxsVDAH8CD+NTf8Jkn/QveIf/AAXN/jR/wmSf9C94h/8ABc3+NQp49tJLuW0TRNfa4iVXkjGntlVbO0n67T+Rqb/hMk/6F7xD/wCC5v8AGj/hMk/6F7xD/wCC5v8AGoT4/sxfCxOi699qaIzCH+z23FAQC30yQPxqb/hMk/6F7xD/AOC5v8aP+EyT/oXvEP8A4Lm/xqCXx/ZwXdvay6LryXFxu8mM6e2ZNoy2PoKn/wCEyT/oXvEP/gub/Gj/AITJP+he8Q/+C5v8ajl8dW0BjEuheIVMr+WmdOf5mwTj8gak/wCEyT/oXvEP/gub/Gj/AITJP+he8Q/+C5v8ajm8dW9tGJJtC8QohYLk6c55JwBx6kgVJ/wmSf8AQveIf/Bc3+NH/CZJ/wBC94h/8Fzf41FdePbSxtZbq70TX4beJd0kj6e2FHqal/4TNP8AoXvEX/gub/Gj/hMk/wChe8Q/+C5v8aiuPHdtaQNPcaF4hjiXG5jprnGTjtUv/CZJ/wBC94h/8Fzf40f8Jkn/AEL3iH/wXN/jTJfHEEELzS6D4hSONSzsdObgDkmnJ41ikRXTw/4hKsMg/wBnNyKX/hMk/wChe8Q/+C5v8aP+EzQDP/CPeIv/AAXNTIPHNvcwJPBoXiF4pFDIw058EHvT/wDhMk/6F7xD/wCC5v8AGj/hMk/6F7xD/wCC5v8AGobXx9Z30Hn2mi6/NEWZN6ae2MqxVh+BBH4VN/wmSf8AQveIf/Bc3+NH/CZJ/wBC94h/8Fzf41Db+PrO7Mwt9F1+UwSGKXbp7fI4wSp9+R+dTf8ACZJ/0L3iH/wXN/jR/wAJkn/QveIf/Bc3+NQReP7Oe6uLaLRdeee3KiaMae2Y9wyufqOarNfz67408Pzw6PqtvBZ/aWmlu7UxKN0eByepzXbUUUUUUVja/rFxpn2G2s4Ipby/n8iDzmKxqwRnJYgE4wp6A84qn4XuWXUNY0y4sba1voZUublrV2aKVpgTuBYA5+Xnj0rpaKKKKKKKKKKZ5MXn+f5aedt2eZtG7bnOM+mafRXMXmpxt4+0mwfSWLCK48q+lO3adqFlQfxAgrk9OMdc109FMeGKVkaSNHaNtyFlB2nGMj0OCafRTJIYplCyxpIoYMAyg4IOQfqDT6KZLDFPE0U0aSRsMMjqCD9QafRTXRZEZHUMjAhlYZBB7GlRFjRURQqKMKqjAA9KWsvVptbV44tIs7STILPLdzFEH+yAoLE9+mODVTwfdQzaI1pDZx2g0+eSyeGJi0YZDg7Ceduemea36KZFDFCpWKNI1LFiEUAEk5J+pNPopiQxRFzHGiGRtzlVA3Hpk+p4FPoooooopsjMsTsiF2CkhAcbj6c1xr6nfw+I9Ln1jTdN8yWU2lv9luGkuLfzBuO8YC4wgBIJ5xiu0oooooooooopksMU8eyaNJEyDtdQRkHI/Wn0Vj+KLqxtfDd82o2/2q3eMxm1B+a4LcCNf9pjwPc0/wAM6o2teGNM1NoBbtdW6SmENnZkfdz7Vq0jKGUqwBUjBBHBpI40ijWONFSNAFVVGAAOgAp1FMihigiEUMaRxjoqKAB36Cn0UyOGKEMIo0TexdtqgZY9Sfc0+imJDFG8jpGitIcuyqAWOMZPrxT6o6zqUej6Nd6hKUCwRFgHbaGb+Fc+5wPxrE8Mtql7fyapf2ekQSXNtGJjaXDyTKRkqjgqAMbm/GupopnkxeeJvLTzQuwPtG7b1xn0zT6KY0MTyxytGjSR52OVBK564PbNPooooooorjtb1O71a81TQ7ez0ySztwkN0NRuGjEpdA4CbQex745rrLZ5JLWJ5YDBIyAtEWDFD6ZHBx7VLRRRRRRRRRRTIoYoE2QxpGmSdqKAMk5J/On1XvjeCykOnrA11j92LhiqE57kAnpntXMeErvUZNb1q0m0/TILeKUST3FlNI4luWA3D5lHIAGfqtdfRTFhiSSSRI0WSTG9goBbHAye9Pooooooorn/ABNYXs82kX9lbm6bTrs3D26uqvIDG6YUsQucuDyRwDTPD1nfnXNY1m9s2shfLAkdvI6vIvlhgSShK4O7jB7V0dFFFFFFFFFFFFFYWpaddXHi/Qb+KMNbWkV0sz7gNpcIF46nO01u0UUUUUUUUUUUUVz/AInu9dijgttF0ue4WYkT3MMsSvAv+yJGALHseg689Db8PQLbaRHGNMm04hjuhmkSR2PdmZWYMT1znNatFFFFFFFFFFFFNk3iNjGFMmDtDHAJ7ZNcW+nahqevafdjw2ulXUF0s13qAkhP2hFUqYwyHewJKkBgBhRnB4rtqKKKKKKKKKKKKK5zxJoGpazeWFxZarDarZs0ggmtPOSVyCAW+demcj0ODUngnTdQ0fwbpen6oyG7ggCOEUALj+HgkHHr3rfoooooooooooorN1+wTU9DubR7GG+WQAm2mcosmGBwSPp9PXiuf0PSJl8Sw6hbeHl0CzitpIpoAYV+0OxUq2IiQdoVvvf3uO9dlRRRRRRRRRRRRRXBap4ekTxNrF/N4bGtR6gIjA2+IrblYwhysrABsjOVBOO/aus0CyutN8PafZXs/n3VvbpHNLuLb2AAJyeTz61o0UUUUUUUUUUUVW1FrtNMu2sER7xYXNuj/daTadoPtnFU/Dml/wBj6FbWzBvPK+ZcM7BmaVuWJbvzxn0ArVooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooor/9k=)
4
ПРИЛОЖЕНИЕ 2
#include <iostream.h>
#include <math.h>
double alpha(double x,double b)
double talpha=sqrt(x*x-b*b);
return talpha;
double alpha1(double d)
double talpha=sqrt(fabs(d));
return talpha;
double s(double x,double b)
double ss=0.5*x-b*b;
return ss;
double f(double b,double a31,double a32,double s)
double fff=a32*a32*b*b*(a31*a31+s)+a31*a32*(b*b*b*b+s*a32*a32);
return fff;
double ff(double b,double a12,double a1,double a2,double s)
double kk=a2*b*b*(a1+s)+a12*(b*b*b*b+s*a2);
return kk;
double compl(double x,double b)
double g=x*x-b*b;
return g;
void koren(double b,double w)
cout<<"При w="<<w<<"есть корень, равный "<<b<<endl;
void main()
{
double p,tt,t,k,r=0.0,bb=0,e=0.001,amega=10.0,betta=0.0,c1=5.6,c2=3.5,c,x,x1=0,x2=0,s1=0;
double h_betta=0.1,h_amega=0.1,func1,func2;
double f(double,double,double,double);
double compl(double,double);
double ff(double,double,double,double,double);
int v=0;
amega=h_amega;
do {
amega+=h_amega;
betta=h_betta;
x1=amega/c1;
x2=amega/c2;
do {
c=betta;
do {
betta+=h_betta;
if (betta>x1) betta=x2+h_betta;
else
{
func1=f(betta, alpha(x1,betta), alpha(x2,betta), s(x2,betta));
if (func1==0)
{koren(betta,amega);v=1;}
}
bb=betta+h_betta;
if (bb>x1)
bb=x2+h_betta;
else
{
func2=f(bb, alpha(x1,bb), alpha(x2,bb), s(x2,bb));
if (func2==0)
{koren(bb,amega);v=1;}
}
if (func1*func2<0)
{x=(betta+bb)/2;
r=f(x,alpha(x1,x),alpha(x2,x),s(x2,x));
if (r*func1<0)
bb=x;betta=bb;
if (r*func2<0) betta=x;
if (r*func1==0)
{koren(betta,amega);v=1;}
}
if (betta>x2&&betta<2*x2)
{
t=compl(x1,betta);
tt=compl(x2,betta);
func1=ff(betta, alpha1(t*tt), alpha1(t*t), alpha1(tt*tt), s(x2,betta));
if (func1==0)
{koren(betta,amega);v=1;}
bb+=h_betta;
if (bb<2*x2)
{
t=compl(x1,bb);
tt=compl(x2,bb);
func2=ff(bb, alpha1(t*tt), alpha1(t*t), alpha1(tt*tt), s(x2,bb););
if (func2==0) {koren(bb,amega);v=1;}
if (func1*func2<0)
{
x=(betta+bb)/2;
r=f(x,alpha1(compl(x1,x)),alpha1(compl(x2,x)),s(x2,x));
if (r*func1<0) bb=x;betta=bb;
if (r*func2<0) betta=x;
if (r*func1==0) {koren(x,amega);v=1;}
}
while (r>e|v==0|betta<2*x2);
betta=c;
}
while (betta<2*x1);
}
while (amega<10);
if (v==0) cout<<"в данном уравнении нет вещественных корней!"<<endl;
}
ПРИЛОЖЕНИЕ 3
#include <fstream.h>
#include <iostream.h>
#include <math.h>
#include <complex.h>
#include <io.h>
#include <process.h>
double alpha31(double x,double a1,double a2)
{ double talpha;
talpha=sqrt(fabs(x*x-a1*a1-a2*a2));
return talpha;}
double s(double x,double a1,double a2)
{ double ss=0.5*x-a1*a1-a2*a2;
return ss;}
double cabs(double x,double y)
{ double g=sqrt(x*x+y*y);
return g;}
complex det(double a1,double a2, double a31,double a32,double s, double \
h1,double h2,double h3)
{ complex det1[9],det2;
complex ii(0,1);
complex x12=exp(ii*ii*a31*(h2-h1));
complex x13=exp(ii*ii*a31*(h3-h1)),x21=exp(ii*ii*a31*(h2-h1)),\
x32=exp(ii*ii*a31*(h2-h3));
complex x23=exp(ii*ii*a31*(h3-h2)),x31=exp(ii*ii*a31*(h1-h3));
complex y12=exp(ii*ii*a32*(h2-h1)),y13=exp(ii*ii*a32*(h3-h1)),\
y21=exp(ii*ii*a32*(h2-h1)),y32=exp(ii*ii*a32*(h2-h3));
complex y23=exp(ii*ii*a32*(h3-h2)),y31=exp(ii*ii*a32*(h1-h3));
complex A[9][9];
complex delta1=s*a31+a32*(a1*a1+a2*a2);
complex delta2=a32*a32*(a1*a1+a2*a2)*(a31*a31+s)+a31*a32*(a1*a1+\
a2*a2+s*a32*a32);
complex M[9][9]={{(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*\
(a31-a32)/delta2,0,(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*\
(a31-a32)/delta2,0,(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*\
(a31-a32)/delta2,0,},\
{a2*(2*a1*a1-a32*a32)/(a1*delta1),2*(a1*a1*a31+a32*(a2*a2+s))/delta2,0,\
a2*(2*a1*a1-a32*a32)/(a1*delta1),2*(a1*a1*a31+a32*(a2*a2+s))/delta2,0,\
a2*(2*a1*a1-a32*a32)/(a1*delta1),2*(a1*a1*a31+a32*(a2*a2+s))/delta2,0},\
{0,0,1/delta1,0,0,1/delta1,0,0,1/delta1},\
{(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0,\
(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0,\
(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0},\
{a2*(2*a1*a1-a32*a32)/(a1*delta1),2*(a1*a1*a31+a32*(a2*a2+s))/delta2,0},\
{0,0,1/delta1,0,0,1/delta1,0,0,1/delta1},\
{(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0,\
(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0,\
(a31*(2*a2*a2+a32*a32)+a32*(a1*a1-a2*a2))/delta2,2*a1*a2*(a31-a32)/delta2,0},\
{a2*(2*a1*a1-a32*a32)/(a1*delta1),2*(a1*a1*a31+a32*(a2*a2+s))/delta2,0},\
{0,0,1/delta1,0,0,1/delta1,0,0,1/delta1}};
complex N[9][9]={{1,0,0,1/2*((a1*a1*x12)+(a2*a2+a31*a32)*y12),a1*a2*(x12-y12)/2,\
a1*a32*(x12-y12)/2,1/2*((a1*a1*x13)+(a2*a2+a31*a32)*y13),a1*a2*(x13-y13)/2,\
a1*a32*(x13-y13)/2},\
{0,1,0,a1*a2*(x12-y12)/2,1/2*(a2*a2*x12+(a1*a1+a31*a32)*y12),a2*a32*(x12-y12)/2,\
a1*a2*(x13-y13)/2,1/2*(a2*a2*x13+(a1*a1+a31*a32)*y13),a2*a32*(x13-y13)/2},\
{0,0,1,a1*a31*(x12-y12)/2,a2*a31*(x12-y12)/2,(a31*a32*x12+(a1*a1+a2*a2)*y12)/2,\
a1*a31*(x12-y12)/2,a2*a31*(x12-y12)/2,(a31*a32*x12+(a1*a1+a2*a2)*y12)/2},\
{(a1*a1*x21+(a2*a2+a31*a32)*y21)/2,a1*a2*(x21-y21)/2,a1*a32*(y21-x21)/2,1,0,0,\
1/2*((a1*a1*x23)+(a2*a2+a31*a32)*y23),a1*a2*(x23-y23)/2,a1*a32*(x12-y23)/2},\
{a1*a2*(x21-y21)/2,(a2*a2*x21+(a1*a1+a31*a32)*y21)/2,a2*a32*(y21-x21)/2,0,1,0,\
a1*a2*(x23-y23)/2,1/2*(a2*a2*x23+(a1*a1+a31*a32)*y23),a2*a32*(x23-y23)/2},\
{a1*a31*(y21-x21)/2,a2*a31*(y21-x21)/2,(a31*a32*x21+(a1*a1+a2*a2)*y21)/2,0,0,1,\
a1*a31*(x23-y23)/2,a2*a31*(x23-y23)/2,(a31*a32*x23+(a1*a1+a2*a2)*y23)/2},\
{(a1*a1*x31+(a2*a2+a31*a32)*y31)/2,a1*a2*(x31-y31)/2,a1*a32*(y31-x31)/2,\
(a1*a1*x32+(a2*a2+a31*a32)*y32)/2,a1*a2*(x32-y32)/2,a1*a32*(y32-x32)/2,1,0,0},\
{a1*a2*(x31-y31)/2,(a2*a2*x31+(a1*a1+a31*a32)*y31)/2,a2*a32*(y31-x31)/2,a1*a2*(x32-y32)/2,(a2*a2*x32+(a1*a1+a31*a32)*y32)/2 ,a2*a32*(y32-x32)/2,0,1,0},{a1*a31*(y31-x31)/2,a2*a31*(y31-x31)/2,(a31*a32*x31+(a1*a1+a2*a2)*y31)/2,\
a1*a31*(y32-x32)/2,a2*a31*(y32-x32)/2,(a31*a32*x32+(a1*a1+a2*a2)*y32)/2,0,0,1.0}};
complex C[9][9];
for (int t=0;t<9;t++){
for (int r=0;r<9;r++){
for (int q=0;q<9;++q){
C[t][r]+=M[t][q]*N[q][r];} } }
det1[0]=C[0][0];
int k=0;
int u=9;
int j;
int i;
for (k=0;k<u-1;k++)
{ for (i=k;i<u;i++)
{ for (j=k;j<u;j++)
{ A[i][j]=C[i][j]-C[i][k]*C[k][j]/C[k][k];}
det1[k+1]=A[k][k]*det1[k];
for (i=k;i<u;i++)
{ for (j=k;i<u;i++)
{ C[i][j]=A[i][j]; } } }
complex d1=det1[u];
cout<<"determinant="<<cabs(real(d1),imag(d1))<<endl;
return d1;
}
struct element
{float zn;};
ofstream& operator <<(ofstream& out,element el)
{ out<<" "<<el.zn<<endl;
return out;}
double main()
{ double h1,h2,h3,amega=20.0,alpha1=0.0,alpha2=0.0,c1=1.5,c2=5.6,x1=0.0,x2=0.0,\
s1=0.0,h_alpha1=0.01,h_alpha2=0.01,h_amega=0.01;
complex func1,func2;
complex f(double,double,double,double);
double s(double,double,double);
double alpha31(double,double,double);
double cabs(double,double);
char ans;
do{
ofstream file1("data.txt",ios::trunc);
cout<<endl<<"h1=";
cin>>h1;
cout<<endl;
cout<<"h2=";
cin>>h2;
cout<<endl<<"h3=";
cin>>h3;
alpha1=h_alpha1;
alpha2=h_alpha2;
x1=amega/c1;
x2=amega/c2;
do{
func1=det(alpha1,alpha2,alpha31(x1,alpha1,alpha2),\
alpha31(x2,alpha1,alpha2),s(x2,alpha2,alpha2),h1,h2,h3);
ofstream file1("data.txt",ios::app);
if (!file1)
{cerr<<"Error open file!";
return 1;}
file1<<alpha1<<" "<<cabs(real(func1),imag(func1))<<endl;
alpha1+=h_alpha1;
alpha2+=h_alpha2; }
while (alpha1<(2*x2));
file1.close();
cout<<" Press 'N' for QUIT ";
cin>>ans; }
while (ans!='N');
}