Решение задачи Коши для обыкновенных дифференциальных уравнений, страница 2

Номер шага i

Значение аргумента xi

Приближенное решение yi

Точное решение y(xi)

Погрешность  |y(x i) - yi|

0

0

0

0

0

1

0,03927

-0,03927

-0,03926

1,0092E-05

2

0,07854

-0,07851

-0,07846

5,0476E-05

3

0,11781

-0,11766

-0,11754

0,00012134

4

0,15708

-0,15666

-0,15643

0,0002231

5

0,19635

-0,19545

-0,19509

0,00035641

6

0,235619

-0,23397

-0,23345

0,00052219

7

0,274889

-0,27216

-0,27144

0,00072156

8

0,314159

-0,30997

-0,30902

0,00095595

9

0,353429

-0,34734

-0,34612

0,00122703

10

0,392699

-0,38422

-0,38268

0,00153675

11

0,431969

-0,42055

-0,41866

0,00188739

12

0,471239

-0,45627

-0,45399

0,00228153

13

0,510509

-0,49134

-0,48862

0,00272209

14

0,549779

-0,52571

-0,5225

0,00321236

15

0,589049

-0,55933

-0,55557

0,00375602

16

0,628319

-0,59214

-0,58779

0,00435718

17

0,667588

-0,62411

-0,61909

0,0050204

18

0,706858

-0,6552

-0,64945

0,00575073

19

0,746128

-0,68535

-0,6788

0,00655373

20

0,785398

-0,71454

-0,70711

0,00743556

21

0,824668

-0,74273

-0,73432

0,00840296

22

0,863938

-0,76987

-0,76041

0,00946335

23

0,903208

-0,79594

-0,78532

0,01062484

24

0,942478

-0,82091

-0,80902

0,01189629

25

0,981748

-0,84476

-0,83147

0,0132874

26

1,021018

-0,86745

-0,85264

0,01480868

27

1,060288

-0,88897

-0,8725

0,01647161

28

1,099557

-0,9093

-0,89101

0,01828858

29

1,138827

-0,92842

-0,90814

0,02027306

30

1,178097

-0,94632

-0,92388

0,02243955

31

1,217367

-0,963

-0,93819

0,02480371

32

1,256637

-0,97844

-0,95106

0,02738235

33

1,295907

-0,99265

-0,96246

0,03019348

34

1,335177

-1,00563

-0,97237

0,0332564

35

1,374447

-1,01738

-0,98079

0,03659163

36

1,413717

-1,02791

-0,98769

0,04022101

37

1,452987

-1,03724

-0,99307

0,04416765

38

1,492257

-1,04537

-0,99692

0,04845594

39

1,531526

-1,05234

-0,99923

0,05311151

40

1,570796

-1,05816

-1

0,0581612

                         Модифицированный метод Эйлера (предиктор-корректор)

В данном методе вычисление yi+1 состоит из двух этапов:

i+1 = yi + h*f(xi, yi),

yi+1 = yi + h/2*[f(xi, yi) + f(xi+1, ỹi+1)], i = 0, 1, 2,…..M.

Данный метод также называется методом предиктор-корректор. На первом этапе приближенное значение предсказывается с первым порядком точности, а на втором этапе это предсказание исправляется, так что результирующее значение имеет второй порядок точности.

M=10