Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Курсач.doc
Скачиваний:
3
Добавлен:
08.09.2019
Размер:
236.03 Кб
Скачать

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)