For[num = 2, 
      num < 10, num++, 
eq1 =  
        Sum[X^i Y^j, {i, 1, num}, {j, 1, num}] == 
          1; 
eq2 =  
        Sum[i/j X^i Y^j, {i, 1, num}, {j, 1, num}] == 
          1/2; 
sol = 
        FindRoot[{eq1, 
            eq2}, {{X, 0, 1}, {Y, 0, 
              1}}] ; 
Print["N=", 
        num, "\n"]; 
Print[
        sol];
]