N = 5;% Change this to 50000 for your case
a = 3;% The diagonal value
b = -1;% The off-diagonal value
% Create the matrix
M = diag(a*ones(1,N)) + diag(b*ones(1,N-1),1) + diag(b*ones(1,N-1),-1);
% Create the vector
v = [2; ones(N-2,1); 2];
% Solve [See //www.tatmou.com/help/matlab/ref/mldivide.html for details]
x = M\v
