0

因此，当我尝试在 Octave 中运行我的代码时，会出现这些问题：

*

T(1.2)：下标必须是整数 1 到 (2^63)-1 或逻辑

我尝试使用欧拉方法来研究，这是我的代码

function fxy = fa(x,y)

% fa merupakan fungsi persamaan diferensial

fxy = x^2 - 4*y;

end

function [T, Y] = euler(f, a, b, y0, n)

% fungsi EULER penyelesaian numerik persamaan diferensial dengan metode
% Euler explisit

h = (b-a)/n;
Y = zeros(n+1,1);
T = zeros(n+1,1);

%syarat awal
T(1) = a;
Y(1) = y0;

% jalankan sebanyak n langkah
for i = 1:n
    fi = f(T(i), Y(i));
    Y(i+1) = Y(i) +  h*fi;
    T(i+h) = T(i) + h;
end

    x0 = 0;   % batas bawah/ bayas awal
    xm = 1;   % batas atas / batas akhir
    y0 = 1;   %syarat batas awal
    n = 5;  % banyaknya segmen / langkah


    [x,y] = euler(@fa, x0, xm, y0, n);

有人可以帮助解决我的问题吗？

4

0 回答 0