x=1:30; c=abs( cos(1./x)-1 ); s = abs( sin(1./x)-1./x); figure(1), plot( x, s, 'o', x, c, 'x') figure(2), loglog( x, s, 'o', x, c, 'x')