Assuming H is a scalar:

Code:

%inserting values
R = 12000
C = 1.6e-9
n = 0
HH=[]; %empty array to hold the output
%the while loop
while n<5
n = n+1;
f=n;
RC= 1/(j*2*pi*f*C);
tRC=2/(j*2*pi*f*C);
%insert of matrix 1
M1 = [0 (2/R)+(1/tRC)_(-1/R);
(2/RC)+(2/R) 0 (1/RC);
(-1/RC) (-1/R) (2/RC)];
%insert of matrix 2
M2 = [10/R;
10/RC;
0];
%calculate unknown values in the matrix
M = M1\M2
Vin = 10
%grab the value i need or my calculation
Vout= M(3,1)
H = Vout/Vin
HH=[HH,H]; % stuff result into extended output array
end

CB