- •Цель и задачи
- •1. Постановка задачи.
- •2. Описание методов решения.
- •2. 1. Суть задачи.
- •2. 2. Геометрический смысл задачи.
- •2. 3. Численные методы решения задачи Коши.
- •2. 4. Метод Эйлера.
- •2. 5. Метод Эйлера модифицированный.
- •2. 7. Решение поставленной задачи методами Эйлера и Эйлера модифицированного.
- •2. 7. 1. Метод Эйлера.
- •2. 7. 2. Метод Эйлера модифицированный.
- •3. Алгоритм решения задачи.
- •3. 1. Алгоритм функции.
- •3. 2. Алгоритм программы.
- •4. Форма программы.
- •5. Листинг программы.
- •6. Решение задачи в MathCad.
2. 7. 2. Метод Эйлера модифицированный.
1. Строим оси координат;
2. Отмечаем А(1; 0) – первую точку интегральной кривой;
3. Ищем угол наклона касательной к графику в точке A:
4. Строим касательную l0 в точке А под углом α0;
5. Находим х1 по формуле: xi = х0 + ih, где h – шаг интегрирования
x1 = 0 + 1 · 0,2 = 1,2;
6. Отмечаем середину отрезка x0x1: x0 + h/2, проводим прямую из этой точки до прямой l0, отмечаем точку B(xB; yB);
7. Ищем координаты В:
xB = x0 + h/2 = 0 + 0,2/2 = 0.1
yB = y0 + h/2 · f(x0; y0) = -1.8 + 0,2/2 · 0 = -1.8
Следовательно, точка B имеет координаты (0.1; -1,8);
8. Ищем угол наклона касательной к графику в точке B:
αB = arctg(f(xB; yB)) = arctg((0.1 –( -1.8 )2+2*0.1*(-1.8)) = arctg(-0.036) = 2°
9. Строим касательную l1 в точке B под углом αB;
10. Проводим прямую x = x1 = 1,2 до пересечения с прямой l1, отмечаем точку C(x1; y1);
11. Ищем y точки C:
y1 = yB + h/2(f(xB;yB)) = -1.8+ 0,2/2 · (-0.036) = -1.8036
Следовательно, точка C имеет координаты (1,2; -1.8036).
3. Алгоритм решения задачи.
3. 1. Алгоритм функции.
3. 2. Алгоритм программы.
n = (kx-x0)/h
x(0):=x0
e(0):=y0
em(0):=y0
o (0):=
y0
i:=0…..N
x (i)=
(x(0)+i*h
e(i+1)=e(i)+h*f((x(i);e(i))
i:=0…..N
x (i)=x(0)+i*h
em(i+1)=em(i)+h*f((x(i)+h/2;em(i)+h/2*f((x(i),em(i))
i:=0…..N
x (i)=x(0)+i*h
o(i-1)=(-2)/(1+f(x(i),o(i)*(Exp(-x(i))^2)
i:=0…..N
z1 =
Round(600 + (x(i) - x0) * kx)
z2 =
Round(6480 - (e (i) - miny) * ky)
z3 =
Round(600 + (x(i+1) - x0) * kx)
z4 =
Round(6480- (e (i+1) - miny) * ky)
i:=0…..N
z1 =
Round(600 + (x(i) - x0) * kx)
z2 =
Round(6480 - (em (i) - miny) * ky)
z3 =
Round(600 + (x(i+1) - x0) * kx)
z4 =
Round(6480- (em(i+1) - miny) * ky)
i:=0…..N
z1 =
Round(600 + (x(i) - x0) * kx)
z2 =
Round(6480 - (o (i) - miny) * ky)
z3 =
Round(600 + (x(i+1) - x0) * kx)
z4 =
Round(6480- (o (i+1) - miny) * ky)