This does essentially what you want it to do. I didn't put it in a neat ‘function’ form because I'm not certain how you want to use it:
x = [pi; 10*pi; 100*pi; eps; -eps; 1/eps; realmax/10; realmin];
for k1 = 1:length(x)
if abs(x(k1)) >= 1
xpnt(k1) = fix(log10(abs(x(k1))))+1;
x1(k1) = x(k1)./(10.^(xpnt(k1)));
elseif abs(x(k1)) < 1
xpnt(k1) = fix(log10(abs(x(k1))));
x1(k1) = x(k1)./(10.^(xpnt(k1)));
end
fprintf(1,'\n%18.15fE%0+4d\n', x1(k1), xpnt(k1))
end
You can create formatted strings for each value of x by replacing the fprintf with sprintf and making the appropriate changes.
