Курсовая работа: Решение прикладных задач методом дихотомии
Соответствующее уравнение касательной:y - yn -1 = k ( x - xn -1 )
Отсюда получаем значение х=хn , соответствующее точке: х n =х n -1 + h ,
А именно: yn - yn -1 = kn -1 ( xn -1 + h - xn -1 ), или
yn =yn-1 +h·kn-1
yn =yn-1 +h·f(xn-1, yn-1 )
Полученная формула является основной расчетной формулой метода Эйлера.
Процесс вычислений заканчивается, когда аргумент после очередного приращения выйдет за пределы исследуемого отрезка .
2. Дифференциальное уравнение:
x0 = 0 , y0 = 1, xmax =1, Δx = 0.01; 0.005; 0.001
3. Схема алгоритма:
5. Таблица идентификаторов:
Обозначение | Идентификатор | Тип |
s | s | int |
i | i | int |
x | x | double |
xmax | x_max | double |
x1 | x1 | double |
Δx | h[i] | double |
y | y | double |
d | d | double |
f(x) | f(x) | double |
k | k(x,y) | double |
6. Листинг программы:
#include<stdio.h>
#include<math.h>
double k(double x,double y )
{
return ((x/exp(x*x))-2.*x*y);
}
double f(double x)
{
return ((1./exp(x*x))*(1+x*x/2.));
}
int main(void)
{
int s,i;
double x,x1,x_max=1,y,d;
double h[3]={0.01,0.005,0.001};
FILE*file;