¿Cómo generar una tabla con columnas de encabezado especificadas en Matlab?

Question

Mi técnica de crear un vector a partir de variables de "columna" en una tabla sólo funciona parte del tiempo. ¿Qué diferencia a K1 del índice en el siguiente código? Estoy depurando un método numérico y necesito el índice de columnas, X, K1, K2, K, Y. ¿Todo funciona bien hasta que agrego K1, K2 y K?

CÓDIGO DE FUNCIÓN:

function [index,X,K1,K2,K,Y] = impeulerT(x,y,x1,n)
% modified version of Improved Euler method found in
% Elementary Differential Equations by Edwards and Penney 
X=x;               % initial x
Y=y;               % initial y
x1 = x1;           % final x
n = n;             % number of subintervals
h = (x1-x)/n;      % step size
index = 0;         % initialize index
k1=0;
k2=0;
k=0;
for i=1:n;         % begin loop
k1=f(x,y);         % first slope
k2=f(x+h,y+h*k1);  % second slope
k=(k1+k2)/2;       % average slope
x=x+h;             % new x
y=y+h*k;           % new y
X=[X;x];           % update x-column       
Y=[Y;y];           % update y-column
index = [index;i]; % update index-column
K1=[K1;k1];   Line 22
K2=[K2;k2];
K= [K;k];
end                % end loop
ImprovedEulerTable=table(index,X,K1,K2,K,Y)
end

CÓDIGO DE LLAMADA:

[index,X,K1,K2,K,Y] = impeulerT(0,1,1,10);

REGISTRO:

>> [index,X,K1,K2,K,Y] = impeulerT(0,1,1,10);
Undefined function or variable 'K1'.

Error in impeulerT (line 22)
K1=[K1;k1];

22  K1=[K1;k1];
>> 

Answer 1

Este código es mejor. La primera línea de salida podría cambiarse para que coincida mejor con lo que se hace a mano. Básicamente, lo que eso significa es que las columnas K1, K2 y K podrían desplazarse una fila hacia arriba. El resto de columnas se podrían dejar igual. Sin embargo, esto coincide con Edwards y Penney. Observe cómo se definen aquí K1, K2 y K. ¡Salud! MM

Código de función:

function [index,X,Y,K1,K2,K] = impeulerT(x,y,x1,n)
% modified version of Improved Euler method found in
% Elementary Differential Equations by Edwards and Penney 
X=x;               % initial x
Y=y;               % initial y
x1 = x1;           % final x
n = n;             % number of subintervals
h = (x1-x)/n;      % step size
index = 0;         % initialize index
% Initialize the lower-case variables
k1=0;
k2=0;
k=0;
% Initialize the upper-case variables
K1=k1;
K2=k2;
K =k;
for i=1:n;         % begin loop
k1=f(x,y);         % first slope
k2=f(x+h,y+h*k1);  % second slope
k=(k1+k2)/2;       % average slope
x=x+h;             % new x
y=y+h*k;           % new y
X=[X;x];           % update x-column       
Y=[Y;y];           % update y-column
index = [index;i]; % update index-column
K1=[K1;k1];        % update K1 column
K2=[K2;k2];        % update K2 column
K= [K;k];          % update K column
end                % end loop
ImprovedEulerTable=table(index,X,Y,K1,K2,K)
end

Código de llamada:

% Improved Euler
[index,X,Y,K1,K2,K] = impeulerTX(0,1,1,10);

Answer 2

El error está en el hecho de que en la línea 22 K1=[K1;k1];tiene K1 entre corchetes, pero no estaba definido antes.

La solución es definir K1=[]; antes del bucle for.

Editar: lo mismo ocurre con todas las demás variables. Entonces este código funciona

function [index,X,K1,K2,K,Y] = impeulerT(x,y,x1,n)
% modified version of Improved Euler method found in
% Elementary Differential Equations by Edwards and Penney 
X=x;               % initial x
Y=y;               % initial y
x1 = x1;           % final x
n = n;             % number of subintervals
h = (x1-x)/n;      % step size
index = 0;         % initialize index
k1=0;
k2=0;
k=0;
% Initialize the upper-case variables
K1=[];
K2=[];
K=[];
for i=1:n;         % begin loop
  k1=f(x,y);         % first slope
  k2=f(x+h,y+h*k1);  % second slope
  k=(k1+k2)/2;       % average slope
  x=x+h;             % new x
  y=y+h*k;           % new y
  X=[X;x];           % update x-column       
  Y=[Y;y];           % update y-column
  index = [index;i]; % update index-column
  K1=[K1;k1];  % Line 22
  K2=[K2;k2];
  K= [K;k];
end                % end loop
ImprovedEulerTable=table(index,X,K1,K2,K,Y)
end