T(1.2): индексы должны быть либо целыми числами от 1 до (2^63)-1, либо логическими.
Итак, эта проблема возникает, когда я пытаюсь запустить свой код в 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);
Может ли кто-нибудь помочь решить мою проблему?