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