Actual source code: ex61m.m
petsc-3.4.2 2013-07-02

      2: f = PetscBinaryRead('perfectmovie','cell');
      3: s = size(f); s = s(2);
      4: clear M
      5: for i=1:s,
      6:   cv = f{i};
      7:   cv = reshape(cv(2:5:5*96*96),96,96);
      8:   pcolor(cv);
      9:   caxis([0 1])
     10:   shading interp;
     11:   colorbar;
     12: M((i-1)/1 + 1) = getframe;
     13: end
     14: movie(M)