
Problemas con ciclo for
Publicado por Patricio (5 intervenciones) el 16/11/2017 22:11:37
Estimados, tengo un problema al ejecutar un programa para resolver edo's mediante el método de heun, pero visto como un método de Runge-Kutta-. Adjunto el código:
Al momento de ejecutar cada sentencia en la ventana de comandos de Matlab estas funcionan correctamente, pero al momento de hacer correr el código como programa me arroja un error, indicando que en la línea 11 el índice excede las dimensiones de la matriz.
Saludos y desde ya agradezco alguna respuesta.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
function [X,Y]=edoheun(f,y0,x0,xn,n)
if nargin(f)~= 2
disp('Error: ingrese una función anónima de dos variables');
return;
end
h = (xn-x0/n);
X = x0:h:xn;
Y = zeros(1,n+1);
Y(1) = y0;
for i=1:n
k1 = f(X(i),Y(i));
k2 = f(X(i)+h,Y(i)+k1*h);
Y(i+1) = Y(i)+(0.5*k1+0.5*k2)*h;
end
end
Al momento de ejecutar cada sentencia en la ventana de comandos de Matlab estas funcionan correctamente, pero al momento de hacer correr el código como programa me arroja un error, indicando que en la línea 11 el índice excede las dimensiones de la matriz.
Saludos y desde ya agradezco alguna respuesta.
Valora esta pregunta


0