function y=f13per(i,v)  % f13per.m
  persistent  palex;

  if isempty(palex)
    palex= [10 20 30 40 50];
    fprintf('Init  : %s\n',mat2str(palex));
  end
  fprintf('Before: %s\n',mat2str(palex));
  if (nargin ==2)
    palex(i) = v;
  end
  fprintf('After:  %s\n',mat2str(palex));
end
