Здравствуйте. В коде
Код:
clc
clear all
a=1;
T=5;
syms x1 t
a0=x1*5+5;
da0dx1=diff(a0,x1,1);
f=t+0.001*x1;
N=15;
h1=a/N;
h3=T/(20*N);
x(1)=0;
for i=1:N-1
x(i+1)=x(1)+i*h1;
end;
z(1)=0;
for j=1:20*N-1
z(i+1)=z(1)+j*h3;
end;
w=0;
for k=1:20*N-1
for i=1:N-1
f0=subs(f,[x1 t],[x(i) z(k)]);
w=w+1;
end;
end;
какая-то противоречащая разуму ошибка.
Код:
Index exceeds matrix dimensions.
Error in test (line 34)
f0=subs(f,[x1 t],[x(i) z(k)]);
Для интереса я завёл специально переменную

, которая фиксирует, на каком шаге произошла ошибка. И каждый раз эта ошибка вылетает на

шаге.
Помогите, пожалуйста, разобраться.