Éder,

A lattice::contourplot() só interpola se o grid fornecido for igualmente espaçado em todas as direções. Ao fazer o log de x você causa deformação do espaçamento. Então você deve iniciar com a escala log de x igualmente espaçada, algo que tentei reproduzir aprimorando seu CMR

logx <- seq(log10(min(x[1])), log(x[100]), l=100)
grid <- expand.grid(logx=logx, y=y)
r <- as.vector(sqrt(outer(10^(2*logx), y^2, "+")))
grid$z <- cos(r^2) * exp(-r/(pi^3))
contourplot(z~logx*y,
            grid,col.regions=colorRampPalette(c("green","yellow","red")),region =TRUE)

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================