Xreferat.com » Рефераты по информатике и программированию » Экспериментальное исследование свойств методов Рунге-Кутты

Экспериментальное исследование свойств методов Рунге-Кутты

МИНИСТЕРСТВО НАУКИ И ОБРАЗОВАНИЯ УКРАИНЫ

Севастопольский национальный технический университет

Кафедра технической кибернетики


КУРСОВАЯ РАБОТА ПО ДИСЦИПЛИНЕ

«Вычислительные методы» на тему:

«Экспериментальное исследование свойств методов Рунге-Кутты»


Выполнила: студентка гр. А-31д

Воротилова Я.М.


Проверил: Мирянов В.И.


Севастополь

2004

СОДЕРЖАНИЕ


ВВЕДЕНИЕ

1 ПОСТАНОВКА ЗАДАЧИ

1.1 Приведение к нормальной форме Коши

1.2 Метод Рунге-Кутты второго порядка

2 ОПИСАНИЕ ПРОГРАММНЫХ МОДУЛЕЙ

2.1 Основная программа

2.2 Функция вычисления точного решения

2.3 Процедура вычисления правых частей системы уравнений в нормальной форме Коши

2.4 Процедура RK2

2.5 Процедура RK4

3 ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ МЕТОДОВ РУНГЕ-КУТТЫ

3.1 Анализ влияния величины шага на точность интегрирования методами Рунге-Кутты второго и четвертого порядка

3.2 Проверка гипотезы Рунге

3.3 Исследование поведение ошибки интегрирования как функции независимой переменной для обоих методов Рунге-Кутты при различных значениях шага

3.4 Сравнительный анализ эффективности методов Рунге-Кутты при различных требованиях к точности вычисления

ЗАКЛЮЧЕНИЕ

БИБЛИОГРАФИЯ

ПРИЛОЖЕНИЕ А

ПРИЛОЖЕНИЕ Б

ПРИЛОЖЕНИЕ В

ВВЕДЕНИЕ


Настоящая курсовая работа посвящена опытному исследованию свойств методов Рунге-Кутты и реализации на персональных компьютерах численных методов приближенного интегрирования ОДУ, наиболее часто применяющихся в практике моделирования и проектирования СА и У. Экспериментальные исследования проводятся с помощью составленных и отлаженных программ интегрирования обыкновенных дифференциальных уравнений на ЭВМ.

Задание предполагает:

закрепление теоретических навыков и знаний в вопросе о проблематике интегрирования ОДУ и численного решения задачи Коши методом Рунге-Кутты, изучение их основных свойств (точность, эффективность, устойчивость) и основных характеристик данных свойств (локальная и глобальная алгоритмические ошибки, порядок метода, ошибка вычисления и т.п.) ;

приобретение основных навыков составления и отладки процедур и функций интегрирования на основе методов Рунге-Кутты и программ интегрирования систем дифференциальных уравнений с использованием все тех же процедур и функций;

проведение опытных исследований зависимости точности, эффективности и устойчивости алгоритмов интегрирования от величины шага интегрирования и порядка метода Рунге-Кутты на ЭВМ.

В различных сферах технических и даже экономических отраслей приходится достаточно часто сталкиваться с математическими задачами, для которых не представляется возможным описать точное решение классическими методами или сие решение выражено крайне неудобочитаемыми соотношениями, которые представляют из себя неприемлемую для мозга пищу, не говоря уже об использовании или реализации на практике.

Разрабатываемые вычислительной математикой численные методы носят в основном ориентировочный характер, однако они позволяют получить итоговый числовой результат со сносной для практических нужд точностью. Численные методы представляют собой алгоритмы вычисления приблизительных значений искомого решения на определенной сетке значений аргумента. При определенных условиях значения аргумента могут являться точными.

Численные методы не позволяют найти общее решение: полученное решение является частным. Но одним из многочисленных плюсов данных методов можно назвать высокую степень применимости к обширным классам уравнений и всем типам вопросов и заданий к ним. Посему с появлением электронных вычислительных машин численные методы стали одними из основных технологий решения определенных практических задач решения ОДУ.

Большую значимость имеет вопрос о верности вычислений на ЭВМ, поскольку при практической реализации имеет место обширный объем обрабатываемой подсчитываемой информации и погрешности могут достаточно сильно исковеркать конечный результат, принимаемый нами за действительный с «поправками на ветер». Кроме сказанного оценка точности численного метода немаловажна и потому, что увеличить точность в некоторых пределах можно за счет увеличения объемов вычислений, а уменьшить временные затраты при решении задачи - за счет снижения точности получаемого результата.

Для понижения погрешности методов интегрирования ОДУ, использующего разложения искомого решения в ряд Тейлора, необходимо принимать во внимание большее количество членов ряда. При всем при этом появляется потребность аппроксимации производных правых частей ОДУ. Ключевая идея методов Рунге-Кутты заключается в том, что производные аппроксимируются через значения функции Экспериментальное исследование свойств методов Рунге-Кутты в точках на интервале Экспериментальное исследование свойств методов Рунге-Кутты, которые выбираются из условия наибольшей близости алгоритма к ряду Тейлора. В зависимости от старшей степени Экспериментальное исследование свойств методов Рунге-Кутты, с коей учитываются члены ряда, построены всевозможные вычислительные схемы Рунге-Кутты разных порядков точности.

Среди достоинств схем Рунге-Кутты не следует обходить во внимании:

удобоваримую точность;

одноступенчатость, то есть дабы найти Экспериментальное исследование свойств методов Рунге-Кутты, необходима информация лишь о предыдущей точке Экспериментальное исследование свойств методов Рунге-Кутты;

координирование с рядом Тейлора вплоть до членов порядка Экспериментальное исследование свойств методов Рунге-Кутты, где степень Экспериментальное исследование свойств методов Рунге-Кутты неодинакова для различных методов и именуется порядком метода;

отсутствие необходимости вычисления производных от Экспериментальное исследование свойств методов Рунге-Кутты, причем накладывается требование вычисления всего-навсего самой функции.

Собственно благодаря вышеуказанному свойству c) методы Рунге-Кутты предпочтительней рядов Тейлора для реализации на практике. Тем не менее поводов для веселья мало, ибо перед нами стоит нелегкая задача неоднократного вычисления функции Экспериментальное исследование свойств методов Рунге-Кутты при неодинаковых значениях Экспериментальное исследование свойств методов Рунге-Кутты и Экспериментальное исследование свойств методов Рунге-Кутты для вычисления последующей точки решения. Это Богом дарованное наказание за преподнесенную нам численным методом поблажку, заключающуюся в отсутствии какой бы то ни было надобности вычисления иной раз весьма громоздких производных, но трудностей боятся кто угодно, только не мы.

1 ПОСТАНОВКА ЗАДАЧИ


1.1 Приведение к нормальной форме Коши


Нормальной формой Коши принято называть общую форму записи ОДУ, то есть представление в виде системы Экспериментальное исследование свойств методов Рунге-Кутты дифференциальных уравнений первого порядка:

Экспериментальное исследование свойств методов Рунге-Кутты (1)

ДУ второго порядка, заданное согласно варианту №3 имеет вид:

Экспериментальное исследование свойств методов Рунге-Кутты (2)

Задание предполагает нахождение решения Экспериментальное исследование свойств методов Рунге-Кутты на интервале Экспериментальное исследование свойств методов Рунге-Кутты при следующих начальных условиях:

Экспериментальное исследование свойств методов Рунге-Кутты (3)

Для решения ДУ его просто необходимо представить согласно нормальной формы Коши. Для этого руководствуемся следующими обозначениями:

Экспериментальное исследование свойств методов Рунге-Кутты (4)

В итоге имеется система ДУ первого порядка вида:

Экспериментальное исследование свойств методов Рунге-Кутты (5)

Произведя все вышеописанные манипуляции над заданным в варианте уравнением, получим следующую систему:

Экспериментальное исследование свойств методов Рунге-Кутты (6)

Система (6) есть решение уравнения (2).

1.2 Метод Рунге-Кутты второго порядка


В методах Рунге-Кутты интеграл заменяется линейной комбинацией значений подынтегральной функции, вычисленных при разных значениях аргумента:

Экспериментальное исследование свойств методов Рунге-Кутты (7)

Метод Рунге-Кутты представим в виде:

Экспериментальное исследование свойств методов Рунге-Кутты


(8)


Из вышеуказанных общих формул (8) получают формулы метода Рунге-Кутты 2-ого порядка m=2;

Экспериментальное исследование свойств методов Рунге-Кутты (9)

Для определения метода необходимо найти значения вещественных коэффициентов: Экспериментальное исследование свойств методов Рунге-Кутты. Для этого интеграл, заменяемый линейной комбинацией значений подынтегральной функции, вычисленных при разных значениях аргумента, можно представить как:

Экспериментальное исследование свойств методов Рунге-Кутты (10)

А его, в свою очередь, можно представить рядом Тейлора:

Экспериментальное исследование свойств методов Рунге-Кутты (11)

где Экспериментальное исследование свойств методов Рунге-Кутты- сумма элементов ряда Тейлора, степень которых не ниже 3.

Осталось найти неизвестные значения Экспериментальное исследование свойств методов Рунге-Кутты

Экспериментальное исследование свойств методов Рунге-Кутты (12)

В результате таких бесхитростных манипуляций получаем искомый ряд Тейлора:

Экспериментальное исследование свойств методов Рунге-Кутты (13)

Приравняем коэффициенты при одинаковых степенях Экспериментальное исследование свойств методов Рунге-Кутты в выражениях

(11) и (13). В итоге получим систему уравнений вида:

Экспериментальное исследование свойств методов Рунге-Кутты (14)

Из свойств системы (14) следует отметить, что она не обладает единственным решением. При Экспериментальное исследование свойств методов Рунге-Кутты значение Экспериментальное исследование свойств методов Рунге-Кутты, значение Экспериментальное исследование свойств методов Рунге-Кутты, а Экспериментальное исследование свойств методов Рунге-Кутты (15)

Подставив полученные коэффициенты в соотношение (8), получаем следующие формулы метода Рунге-Кутты 2-ого порядка:

Экспериментальное исследование свойств методов Рунге-Кутты

Экспериментальное исследование свойств методов Рунге-Кутты (16)

2 ОПИСАНИЕ ПРОГРАММНЫХ МОДУЛЕЙ


Составленная в ходе курсовой работы программа вычисляет решения дифференциального уравнения, с предварительно заданными начальными условиями. Интегрирование происходит согласно двум методам: Рунге-Кутты второго и четвертого порядков.

Программа состоит из следующих модулей:

Основная программа;

Процедура вычисления точного решения ДУ;

Процедура вычисления правых частей;

Процедура выполняющая шаг интегрирования методом Рунге-Кутты 2-ого порядка;

Процедура выполняющая шаг интегрирования методом Рунге-Кутты 4-ого порядка.


2.1 Основная программа


Блок программы осуществляет следующие операции:

запрашивает у нерадивого пользователя величину шага интегрирования и шаг вывода на экран;

вычисляет количество шагов;

с заданным шагом вызывает процедуры интегрирования методом Рунге-Кутты 2-ого и 4-ого порядков на отрезке интегрирования;

вычисляет погрешность и оценку погрешности интегрирования;

выводит замечательные результаты работы программы с заданным шагом вывода на экран.

Для простоты понимания укажем следующие переменные, содержащиеся в программе:

h - шаг интегрирования. Вводится нерадивым пользователем с клавиатуры;

n – число шагов интегрирования;

h_screen - шаг вывода результатов на экран. Вводится нерадивым пользователем с клавиатуры;

i_screen – счётчик вывода результатов на экран. Когда i_screen > h_screen, то происходит вывод результатов и обнуление i_screen;

i, j – переменные, используемые циклом;

e2, e4– ошибки интегрирования для методов Рунге-Кутты 2-ого и 4-ого порядков соответственно. Подсчитываются из соотношения(1):

Экспериментальное исследование свойств методов Рунге-Кутты (1)

e2max, e4max – оценки погрешностей интегрирования для методов Рунге-Кутты 2-ого и 4-ого порядков соответственно. Подсчитываются из соотношения(2):

Экспериментальное исследование свойств методов Рунге-Кутты (2)


t – значения независимой переменной;

t0, tf – пределы интегрирования

y2, y4 – вектора решения для методов Рунге-Кутты 2-ого и 4-ого порядка соответственно в узле tk ;

outfile– переменная файлового типа. Определена для вывода результатов в текстовой файл;

name – переменная строкового типа. Используется для передачи имени файла.

Текст основной программы приведён в приложении А, схема в приложении Б.

2.2 Функция вычисления точного решения


function clearsolve (t: real): real

Функция предназначена для вычисления точного решения для дифференциального уравнения по формуле (3):

Экспериментальное исследование свойств методов Рунге-Кутты (3)

Текст функции приведен в приложении 2, схема в приложение 7.


2.3 Процедура вычисления правых частей системы уравнений в нормальной форме Коши


procedure right(t: real; var x,f: vector_n);

Процедура вычисляет правые части системы однородных дифференциальных уравнений в нормальной форме Коши по формуле (4):


Экспериментальное исследование свойств методов Рунге-Кутты (4)


Текст процедуры приведен в приложении А, а схема в приложение Б.


2.4 Процедура RK2


procedure RK2(t: real; h: real; var x: vector_n);

Укажем формальные параметры:

t – независимая переменная ;

h – шаг интегрирования;

x – массив решений. При входе в процедуру решение в текущем узле интегрирования, при выходе в следующем.

Процедура RK2 выполняет шаг интегрирования системы ОДУ методом Рунге-Кутты 2-ого порядка из соотношения (5):

Экспериментальное исследование свойств методов Рунге-Кутты (5)

где

Экспериментальное исследование свойств методов Рунге-Кутты (6)

Процедура обращается к процедуре вычисления правых частей right с различными параметрами для вычисления Экспериментальное исследование свойств методов Рунге-Кутты и Экспериментальное исследование свойств методов Рунге-Кутты(6). Затем с Божьей помощью (5) считает значение Экспериментальное исследование свойств методов Рунге-Кутты.

Текст процедуры приведен в приложении А, схема в приложение Б.


2.5 Процедура RK4


procedure RK4(t: real; h: real; var x_4: vector_n);

Формальные параметры:

t – независимая переменная ;

h – шаг интегрирования;

x – массив решений. При входе в процедуру решение в текущем узле интегрирования, при выходе в следующем.

Процедура RK4 выполняет шаг интегрирования системы обыкновенных дифференциальных уравнений (1.1.2) методом Рунге-Кутты 4-ого порядка (7).

Экспериментальное исследование свойств методов Рунге-Кутты (7)

где

Экспериментальное исследование свойств методов Рунге-Кутты (8)

Процедура четыре раза обращается к процедуре вычисления правых частей right с разными параметрами для вычисления Экспериментальное исследование свойств методов Рунге-Кутты,Экспериментальное исследование свойств методов Рунге-Кутты,Экспериментальное исследование свойств методов Рунге-Кутты,Экспериментальное исследование свойств методов Рунге-Кутты(8). Затем с Божьей помощью (7) считает значение Экспериментальное исследование свойств методов Рунге-Кутты.

Текст процедуры приведен в приложении А, схема в приложение Б.

3 ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ МЕТОДОВ РУНГЕ-КУТТЫ


3.1 Анализ влияния величины шага на точность интегрирования методами Рунге-Кутты второго и четвертого порядков


Плоды деятельности на ПЭВМ приводятся в приложении В. Результат упорного труда программы представлены в графическом виде (рисунок 1).


Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 1 – Зависимость оценки e2 от шага интегрирования


Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 2 – Зависимость оценки e4 от шага интегрирования


На рисунках изображены зависимости оценки погрешности интегрирования от величины шага интегрирования для обоих методов. Из него видно, что практические результаты соответствуют теоретическим положениям, но не совсем.

При шаге от 0.1 до 1 бесспорно влияние погрешности интегрирования как в методе второго, так и в методе четвертого порядков. Ошибка на данном интервале начинает лавинообразно возрастать, что связано с нарушением устойчивости алгоритма.

С дальнейшим уменьшением шага до 0.1 – 0.001 величина погрешности уменьшается за компанию, и наблюдается достаточно большая точность вычислений.

Дальнейшее уменьшение шага (менее 0.001) вызывает увеличение полной ошибки, а также скверное её поведение. Это связано с возрастанием влияния ошибки вычислений из-за увеличения количества вычислений, необходимых для получения решения. На фоне уменьшения алгоритмической погрешности решающую роль играет погрешность вычислений, которая представляет собой сумму всех ошибок округления при реализации данного метода на конкретной ПЭВМ.

Из двух методов 2-го и 4-го порядков при одинаковых значениях шага точнее метод четвертого порядка, но при уменьшении шага, точность методов постепенно выравнивается.


3.2 Проверка гипотезы Рунге


Согласно гипотезе, суммарная погрешность алгоритма при интегрировании ДУ с постоянным шагом пропорциональна величине шага в степени, равной порядку метода.

Экспериментальное исследование свойств методов Рунге-Кутты (1)

Результаты вычислений на ПЭВМ приводятся в приложении В. Результат работы программы наглядно представлены в графическом виде (рисунок 2).


Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 3 – Зависимость отношения оценки погрешности к величине шага интегрирования в степени, равной порядку метода от шага интегрирования


Гипотеза Рунге экспериментально подтверждается только при таких значениях шага (1-0.01), где погрешность вычисления ПЭВМ влияет на результат в благоразумных пределах. При дальнейшем уменьшении шага пропорциональность скрывается из виду. Это можно истолковать тем, что предположение Рунге не учитывает влияния ошибки вычисления ПЭВМ на полученный результат.


3.3 Исследование поведение ошибки интегрирования как функции независимой переменной для обоих методов Рунге-Кутты при различных значениях шага


Для проведения своих изысканий нам необходимо выбрать 3 значения шага в начале, середине и конце диапазона изменения шага и построить графические зависимости решений от независимой переменной и погрешностей интегрирования от независимой переменной.

Результаты вычислений на ЭВМ приводятся в приложении В. Результат работы программы наглядно представлены в графическом виде (рисунок 3 и 4).


Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 4 – График точного и приближенных решений (h=0.001, 0.1, 0.5) дифференциального уравнения


На рисунке 3 слегка, но все-таки заметны отклонения приближённых решений при разной величине шага от точного, что подтверждает высокую точность получаемого приближенного решения относительно точного.

Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 5 – График зависимости величины ошибки интегрирования от независимой переменной при h=0.001, 0.1, 0.5


При шаге, численно равном 0.5 из рисунка 5 очевидно, что различия между погрешностями методов 2-ого (E5) и 4-ого (E6) порядков довольно значительны, что объясняется более высокой точностью метода 4-ого порядка.

Разница ошибок интегрирования для обоих методов при шаге интегрирования 0.1 остается достаточно высокой, но следует подметить, что точность метода Рунге-Кутты 4-ого порядка при данном значении шага выше, чем при меньшем h=0.001 (это обусловлено тем, что алгоритмическая ошибка стремится к нулю, а вычислительная ещё не показывает свой отвратительный характер).

При величине шага интегрирования h=0.001 ошибки интегрирования (E1 и E2), равно как и при больших h возрастают с ростом независимой переменной, что слава Богу не противоречит теории, однако прослеживается скачкообразный характер графиков, что связано с внушительной лептой вычислительной ошибки в общую погрешность.

3.4 Сравнительный анализ эффективности методов Рунге-Кутты при различных требованиях к точности вычисления


Для анализа надлежит оценить затраты машинного времени на интегрирование ДУ в зависимости от величины погрешности интегрирования.

Так как для довольно непростых уравнений львиная доля времени уходит на вычисление правых частей уравнений, то в качестве оценки затрат машинного времени можно принять количество вычислений правых частей при нахождении решения на всём отрезке интегрирования, что есть произведение количества шагов интегрирования на порядок метода. Зависимость оценки ошибки интегрирования от количества шагов интегрирования определена.

Замечательные результаты вычислений на ПЭВМ приведены все в том же приложении В. Результат работы программы представлены графически (рисунок 5).


Экспериментальное исследование свойств методов Рунге-Кутты

Рисунок 6 – График зависимости величины ошибки интегрирования от количества вычислений правой части


Из близлежащего графика практически определяется метод для различных требованиях к точности и времени работы нашей замечательной программы.

График наглядно демонстрирует, что наличествует такой предел шага, ниже которого программа гоняет свои байты зря, так как оценка погрешности растёт из-за ошибки вычисления ПЭВМ. Стало быть, на практике следует выбирать определенный промежуток шага, желательно в котором алгоритм устойчив.

ЗАКЛЮЧЕНИЕ


При реализации на практике задания для курсовой работы, заключающегося в интегрирования ОДУ, была составлена и отлажена программа, приведенная в приложении А. С помощью данной программы проведена серия опытных исследований свойств методов Рунге-Кутты второго и четвёртого порядков.

При задании определенного интервала значений шага интегрирования ошибка интегрирования уменьшается с уменьшением шага. Подтверждение сего факта можно без труда найти в теории. Однако нельзя обойти вниманием тот труднооспоримый аргумент, что для довольно незначительных значений Экспериментальное исследование свойств методов Рунге-Кутты с уменьшением шага интегрирования ошибке свойственно увеличивается. Это связано с лавинообразным ростом числа требуемых для получения решения вычислений и c увеличением ошибки вычислений.

На заданном интервале значений Экспериментальное исследование свойств методов Рунге-Кутты (при условии, что на интервале ошибка вычислений велика не до безобразия), при одинаковых значениях шага интегрирования метод Рунге-Кутты четвёртого порядка имеет достаточно малую ошибку вычисления относительно ошибки метода Рунге-Кутты второго порядка.

Следует отметить, что коэффициент Экспериментальное исследование свойств методов Рунге-Кутты (формула (1) раздела 3) для таких значений шага интегрирования, при которых можно скрипя сердцем считать его неизменным, имеет меньшее значение для метода Рунге-Кутты четвёртого, нежели второго порядка. Отсюда мораль: ошибка интегрирования с помощью метода Рунге-Кутты четвёртого порядка меньше ошибки интегрирования при использовании метода Рунге-Кутты второго порядка и обусловлено это не только тем, что она обратно пропорциональна величине шага в четвёртой, а не во второй степени, но и оттого, что коэффициент пропорциональности при этом существенно меньше.

Метод Рунге-Кутты четвёртого порядка имеет значительно меньшую ошибку вычисления, нежели метод Рунге-Кутты второго порядка (при условии идентичных машинных затрат). При неизменной требуемой точности для метода Рунге-Кутты четвёртого порядка необходимы меньшие затраты машинного времени относительно метода Рунге-Кутты второго порядка, отчего при корректном выборе шага интегрирования метод Рунге-Кутты четвёртого порядка значительно результативнее, чем метод Рунге-Кутты второго порядка.

БИБЛИОГРАФИЯ


Самарский А.А., Гулин А.В. Численные методы .М.: Наука, ГРФМЛ, 1989.- 432с.

Бахвалов Н.С. Численные методы. М.: Наука, ГРФМЛ,1987.-600 с.

Ляшко И.И., Макаров В.Л. Интегральные методы вычислений. Киев, 1977, 408с.

Мантуров О.В. Курс высшей математики М.:В.Ш.-1991.-448с.

Маликов В.Т. Вычислительные методы и применение Киев: В.Ш.-1989.-213 с.

ПРИЛОЖЕНИЕ А

(Справочное)


Основная программа


program DU(input,output);

uses crt;

type vector_n = array [1..2] of real;

var t0, tf, k, j, n: integer;

t, yt, h, h_screen, e2, e4, e2max, e4max, i_screen: real;

y2, y4: vector_n;

name: string;

outfile: text;


begin

clrscr;

writeln('Please enter file name');

readln(name);

writeln('Please enter h,h_screen');

readln(h,h_screen);

clrscr;

writeln;

assign(outfile,name);

rewrite(outfile);

t0:=0;

tf:=10;

n:=round((tf-t0)/h);

y2[1]:=1;

y2[2]:=0;

y4:=y2;

e2:=0;

e4:=0;

e2max:=0;

e4max:=0;

t:=t0;

i_screen:=h_screen;

for k:=0 to n do

begin

yt:=clearsolve(t);

e2:=abs(yt-y2[1]);

e4:=abs(yt-y4[1]);

if e2>e2max then e2max:=e2;

if e4>e4max then e4max:=e4;

if i_screen>h_screen-0.00001 then

begin


yt:=clearsolve (t);

e2:=abs(yt-y2[1]);

e4:=abs(yt-y4[1]);

if e2>e2max then e2max:=e2;

if e4>e4max then e4max:=e4;

if i_screen>h_screen-0.00001 then

begin

writeln(' t=',t:6:3,'; yt=',yt:9:3,'; y2=',y2[1]:9:3,'; y4=',y4[1]:9:3,'; e2=',e2:6:3,'; e4=',e4:8:6);

writeln(outfile,'

Если Вам нужна помощь с академической работой (курсовая, контрольная, диплом, реферат и т.д.), обратитесь к нашим специалистам. Более 90000 специалистов готовы Вам помочь.
Бесплатные корректировки и доработки. Бесплатная оценка стоимости работы.

Поможем написать работу на аналогичную тему

Получить выполненную работу или консультацию специалиста по вашему учебному проекту
Нужна помощь в написании работы?
Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Пишем статьи РИНЦ, ВАК, Scopus. Помогаем в публикации. Правки вносим бесплатно.

Похожие рефераты: