Array indices must be positive integers or logical values.

Question

Alexa Shumaker 2019-3-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/452643-array-indices-must-be-positive-integers-or-logical-values

编辑： Walter Roberson 2019-3-26

Writing function for Runge Kutta 4th order. ERROR --> Array indices must be positive integers or logical values.

command window call is

dydx = @(x,y) (1+4*x)*sqrt(y)
 RungeKutta(dydx, 0, 1,0.25,1)

code below

function RungeKutta(f, L, U,h,y0)
%f = function 
% L = lower bound
% U = upper bound
% y0 = inital condition y value
%h = step size 
t = L:h:U;
n = length(t);
 % need for loop 
 k1 = zeros (1,n);
 k2 = zeros (1,n);
 k3 = zeros (1,n);
 k4 = zeros (1,n);
 y = zeros(1, n);
 for i = 0: n
        k1(i) = f(t(i),y(i));
        k2(i)  = f(t(i)+0.5*h,y(i)+0.5*k1(i)*h);
        k3(i) = f(t(i) + 0.5*h, y(i)+0.5*k2(i)*h);  
        k4(i) = f(t(i) + h, y(i) + k3(i)*h); 
        y(i+1) =  y(i) + (1/6)*(k1(i)+2*k2(i) + 2*k3(i)+k4(i))*h;
 
 end %end of fo loop 
 hold on
 title('Runge-Kutta Method')
 xlabel('X axis')
 ylabel('Y axis')
 plot(t, f(t),'b', x, y, 'r')
 legend('Original Plot', 'True Plot')
end % end of function