Time-averaged moire fringes.

I'm trying to create time-averaged moire fringes. Please help me out.

%harmonic oscillation
 y = linspace(0,5,500);
 f = 10*pi*1000;  %(lambda= 0.2)
 wave = (0.5*sin(f*y));
 wave2= imresize(wave, [1 201]);
 
 % zero order bessel function of first kind
 z= 0:0.01:2;
 J = besselj(0,10*pi*z); %lambda= 0.2
 plot(z,J);
 subplot(2,2,1);
 title('bessel');
 
 %time-averaged moire fringes
 moire= 0.5 + 0.5*(wave2 .* J);
 onematrix= ones(size(moire));
 moire2= (onematrix.*moire);
 colormap(gray);
 subplot(2,2,2);
 imagesc(moire2);
 title('final moire');