Please, how can I solve runKut4 in line 6?

t = 0.0;

tStop = 6.0;

y = [1.5 0];

h = 0.1;

[T,Y] = runKut4(@dEqs,t,y,tStop,h);

plot(T,Y(:,1),'ko-')

xlabel('time (s)'); ylabel('h/r')

grid on

function F = dEqs(x,y)

p = sqrt(2*y(1) - y(1)^2);

q =1 - y(1);

F = zeros(1,2);

F(1) = y(2);

F(2) = 2/pi*(atan(q/p) + q*p);

end % function dEqs

end

Answer 1

[Alan Stevens]

Try

t = 0.0;
tStop = 6.0;
y = [1.5 0];
h = 0.1;
[T,Y] = ode45(@dEqs,[0 tStop],y);
plot(T,Y(:,1),'ko-')
xlabel('time (s)'); ylabel('h/r')
grid on
function F = dEqs(x,y)
p = sqrt(2*y(1) - y(1)^2);
q =1 - y(1);
F = zeros(2,1);
F(1) = y(2);
F(2) = 2/pi*(atan(q/p) + q*p);
end % function dEqs

[Stephen Akosah]

Thank you Steven, I appreciate your kind support.