a = [1 2 3 4; 5 6 7 8; 9 10 11 12];
b = [20 21 22 23];
c = [30 31 32 33];
for col = 1:size(a, 2)
x(:, col) = (a(:, col) - b(col)) / c(col); % aのn列目は、bとcのn番目を使って計算
disp([char(col+'0') '列目:']);
curcol = regexprep(num2str(a(:, col)'),'\s+',';');
disp(['x = (a - b) / c = ([' curcol '] - ' num2str(b(col)) ') / ' num2str(c(col))]);
x % 結果は4×4 (3x4?) の行列で表示
end