function A1_makeall
  clear all;

  %define the interval, we want to see
  X = -1:0.1:1;
  Y = -1:0.1:1;

  %calc the Lyapunov Function
  for myi = 1:size(X, 2)
    for myj = 1:size(Y, 2)
      Lyapunov_Function(myi, myj) = ...
          X(myi)^2 - 2*X(myi)*Y(myj) + Y(myj)^2; %a)
          %X(myi)^2 - X(myi)*Y(myj) + Y(myj)^2;  %b)
    end
  end

  meshc(X, Y, Lyapunov_Function);           
  
  %%c)
  %for myi = 1:size(X, 2)
  %  Lyapunov_Function_c(myi) = exp(X(myi)^2);
  %end
  %plot(X, Lyapunov_Function_c);
  
end