% FIGURE 4.8.1
      function X = UPDATE(X,Y,S,Z,N)
%                              CALCULATE X = Y + S*Z
      for I=2:N
      for J=2:N
         X(I,J) = Y(I,J) + S*Z(I,J);
      end
      end