Курсовая работа: Численные методы при решении задач

double y[3], xs, xe;

int i;

y[0] = 1.; y[1] = 0.1; y[2] = 0.; // Начальные условия

xs = .0; xe = .1; // Начало интегрирования

printf ("x = %5.3lg, y(%4.2lg) = %10.3lg\n", xs, xs, y[0]);

for (i = 0; i < 20; i++)

{

Adams (func, y, 3, xs, xe, 10, 1.e-3);

xs += 0.1; xe += 0.1;

printf ("x = %5.3lg, y(%4.2lg) = %10.3lg\n", xs, xs, y[0]);

}

return 0;

}

Результат решения задачи 17 на ЭВМ

Для работы программу необходимо скомпилировать в модели не ниже SMALL. Использовался компилятор Micro$oftC 6.00 из одноимённого пакета. После запуска программа выводит следующее:

Программа численного интегрирования системы дифференциальных

уравнений экстраполяционным методом Адамса

Разработчик: студент гр. ПС-146

Нечаев Леонид Владимирович

17.03.2004

Дифференциальное уравнение имеет вид y''' + 2y'' + 3y' + y = x^2 + 5

Итак, зависимость y[x]:

x = 0, y( 0) = 1

x = 0.1, y(0.1) = 1.01

x = 0.2, y(0.2) = 1.02

x = 0.3, y(0.3) = 1.04

x = 0.4, y(0.4) = 1.07

x = 0.5, y(0.5) = 1.11

x = 0.6, y(0.6) = 1.16

x = 0.7, y(0.7) = 1.22