function y=mult(a,b,n)

[fa,ea]=log2(a);
[fb,eb]=log2(b);

a=a*2^(-ea);
b=b*2^(-eb);

y=round(round(a*10^n)*round(b*10^n)*10^(-n))*10^(-n);

y=y*2^(ea+eb);