clc
clear

k=1/100;
gamma=0.2;
n=9;

%Quitar los porcentajes para tomar los datos de un problema distinto al
%ejemplo hecho en clase

%n=input('Ingrese dimension de la matriz: ')
%k=input('ingrese el valor de k: ')
%gamma=input('ingrese el valor de gamma: ')



A=zeros(n) %Crea una matriz cuadrada nxn con ceros en todas las entradas

for i=2:n-1
    %
    A(i,i)=gamma;
    A(i,i+1)=1-3*gamma;
    A(i,i-1)=1-3*gamma;
end

A(1,1)=gamma;
A(1,2)=1-3*gamma;
A(n,n)=gamma;
A(n,n-1)=1-3*gamma;

v=zeros(n,1)

for j=1:n
    v(j)=(j*k)^2+1;
end

w=zeros(n+2,1);
w(1)=1;
w(n+2)=2;
w(2:n+1)=v;

plot(w)

hold on

for t=1:100
v = A*v+b
w(2:10,1)=v
plot(w)
end

    
    