x = linspace( 0, 1000, 100 );
% Alternatively:
% h = ( 1000 - 0 ) / ( 100 - 1 ); x = 0 : h : 1000;
f = @( x )( ( x - log( x + 1 ) ).^4 );
fx = f( x );
figure( 1 );
plot( x, fx );
figure( 2 );
semilogx( x, fx );
figure( 3 );
semilogy( x, fx );
figure( 4 );
loglog( x, fx );