Sine function
l Image(I) ???
l subplot(222)
image(64*(I+1)/2)
colormap gray
  t = 0:pi/50:10*pi
I = sin(t)
subplot(221)
plot(t,I)