CHỨC NĂNG fRK4 (X, Y)
{$ fRK4
fRK4 là một chức năng người dùng cung cấp để đánh giá dY / dX. Điều này
chức năng được sử dụng với chức năng RK4 để giải quyết khác biệt
phương trình bằng phương pháp Runge-Kutta. Nhập fRK4 (X, Y)
chức năng trong cửa sổ phương trình để đánh giá dY / dY cho bạn
vấn đề. Xem các chức năng RK4 để biết thêm thông tin}.
FRK4: = (Y + X) ^ 2
END
CHỨC NĂNG RK4 (LowX, HighX, StepX, Y0)
{$ RK4
RK4 là một chức năng mục đích chung mà giải quyết một lệnh đầu tiên
phương trình vi phân của hình thức dY / dX = fRK4 (X, Y) bằng cách sử dụng
thuật toán để 4 Runge-Kutta. Các chức năng RK4 gọi hàm
fRK4 (X, Y) được cung cấp bởi người sử dụng để đánh giá dY / dX tại các giá trị cụ thể
của X và Y. Người sử dụng phải cung cấp các chức năng. FRK4
RK4 đòi hỏi bốn thông số đầu vào. LowX là ban đầu
giá trị của độc lập HighX biến X. là giá trị cuối cùng
của biến X độc lập và StepX là kích thước bước. Y0 là
giá trị của Y khi X bằng LowX}.
X: = LowX
Y: = Y0;
Tol: = 0,1 * StepX
10:
IF (X> HighX-Tol) THEN GOTO 20
k1: = fRK4 (X, Y ) * StepX
k2: = StepX * fRK4 (X + 0,5 * StepX, Y + 0.5 * k1)
k3: = StepX * fRK4 (X + 0,5 * StepX, Y + 0.5 * k2)
k4: = StepX * fRK4 (X + StepX, Y + k3)
Y: = Y + k1 / 6 + (k2 + k3) / 3 + k4 / 6
X: = X + StepX
GOTO 10;
20:
RK4: = Y
END
đang được dịch, vui lòng đợi..