To clarify i have already wrote some code for it but i keep getting errors
%% Forward Euler Method
clc
a = 0;
b = 1000;
n = 100;
h = (b-a)/n;
y(1,:)=[2 7];%Initial condition
i=0;
for x = a:h:b
i = i +1;
y(i+1,:) = y(i,:) + h * Fun (x,y(i,:));
xx(i) = x;
end
plot(xx,y(1:n+1,1),xx,y(1:n+1,2))
legend('first Eq','second Eq')
grid on
% main function
function [f]=Fun(x,y)
f1_dot = y*beta(y)*(1-(y/max(y))-(x/max(x)));
f2_dot = x*beta(x)*(1-(y/max(y))-(x/max(x)));
f = [f1_dot f2_dot];
end