restart;
  s1:=proc(n,x) 
  local y,xx,i,j,zz::array(1..n,1..n);
  for i from 1 to n do for j from 1 to n do zz[i,j]:=x[i]*(1+x[j]^2);od:
  od:
  y:=array(1..2,[(1)=x, (2)=zz]):
  for j from 1 to n do xx[i]:=zz[i,i]/(add(zz[i,j],j=1..n));od:
  0;
  end proc:
  #showstat(s1);
  
  x0:=Vector(2,[1,1]);
  
  s1(2,x0);
  
  stopat(s1,4);
  
  s1(2,x0);
  
  next
  Y:=eval(y);
  quit
  eval(Y);
  
  unstopat(s1);