integral = 5.822549E–02. Лабораторная работа № 3.

integral = 5. 822549E–02

Лабораторная работа № 3.

Решить методом Эйлера задачу Коши на отрезке [1; 2] для уравнения

, если y(1) = 1

Текст программы на QBASIC.

CLS

PRINT " Reshenie differencialnogo uravneniyа metodom Eulera"

INPUT " x0=", x0

INPUT " xn=", xn

INPUT " y0=", y0

INPUT " h=", h

DIM g(100)

n = (xn – x0) / h

DEF fnf (X, Y) =2*X*Y / (3*X^2–Y^2)

g(1) = y0

X = x0

PRINT " Y("; X; " )="; y0

FOR i = 1 TO n

g(i + 1) = g(i) + fnf(x, g(i)) * h

X = X + h

PRINT " Y("; X; " )="; g(i + 1)

NEXT I

Вычисления по программе для h=0. 1 дают следующий результат:

Y(1)= 1

Y(1. 1)= 1. 1

Y(1. 2)= 1. 2

Y(1. 3)= 1. 3

Y(1. 4)= 1. 4

Y(1. 5)= 1. 5

Y(1. 6)= 1. 6

Y(1. 7)= 1. 7

Y(1. 8)= 1. 8

Y(1. 9)= 1. 9

Y(2)= 2