function S=scale(a,b)
% scale

S=[a 0 0; 0 b 0; 0 0 1];