Цикл с начальными условиями

function my(N);
H = 1/N;
p1=(1/(H^2)) + 3;%первый и последний элемент
p2=(-1/(H^2));%верхняя и нижняя диагональ
p3=(2/(H^2)) + 3;
n1=N-1;
t=1/100;
A=diag(p3*ones(1,n1))+diag(p2*ones(1,n1-1),1)+diag(p2*ones(1,n1-1),-1);
A(1,1)=p1;
A(n1,n1)=p1;
A
L = ones(n1,1);
for i=1:n1
L(i)=cos(H*i*pi);%вектор-столбец
end
F=t*L;
F
S=t*A;
E=eye(n1);
B=E-S;
B
x0=F;%начальные условия
x0
if(norm(B,inf) < 1)
    for i = 0:1:n1-1
        
    end
else
    fprintf('Измените значение переменной t');
end

for i = 0:1:n1-1
        
    end
else
    fprintf('Измените значение переменной t');

x ( i0+1 ) = B * x ( i0 ) + F

x ( i0+1 ) = B.*x ( i0 ) + F

x0=F(1);%начальные условия

if(norm(B,inf) < 1)
    for i = 1:1:n1
        for j = 1:1:n1
        x(i)=x(j)*B(i,j);
        Res(i)=x(i)+F(i);
        end
    end
    
else
    fprintf('Измените значение переменной t');
end