{
TASK:atom
LANG:pascal
}
var a:array[-1..110,1..2] of longint;
    n,i,j,k,z:longint;
    a1,b1,c1,a2,b2,c2:extended;
    cosG,SinG:extended;
    dy,otg:extended;
    p,q,r,t,t1,x,y,x1,y1,otgx,otgy:extended;
begin
   readln(n);
   for i:=1 to n do readln(a[i,1],a[i,2]);
   for i:=1 to n-2 do
      for j:=i+1 to n-1 do
         for k:=j+1 to n do
            begin
               p:=sqr(a[i,1]-a[j,1]) + sqr(a[i,2]-a[j,2]);
               q:=sqr(a[i,1]-a[k,1]) + sqr(a[i,2]-a[k,2]);
               r:=sqr(a[k,1]-a[j,1]) + sqr(a[k,2]-a[j,2]);
               if (p+q=r) then
                  begin
                     t:=r/4;
                     x:=(a[k,1]+a[j,1])/2;
                     y:=(a[k,2]+a[j,2])/2;
                     for z:=1 to n do
                        if (z<>i)and(z<>k)and(z<>j) then
                           begin
                              r:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                              if r<t then t:=r;
                           end;
                  end
               else
                  if (p+r=q) then
                  begin
                     t:=q/4;
                     x:=(a[k,1]+a[i,1])/2;
                     y:=(a[k,2]+a[i,2])/2;
                     for z:=1 to n do
                        if (z<>i)and(z<>k)and(z<>j) then
                           begin
                              q:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                              if q<t then t:=q;
                           end;
                  end
               else
                  if (r+q=p) then
                  begin
                     t:=p/4;
                     x:=(a[i,1]+a[j,1])/2;
                     y:=(a[i,2]+a[j,2])/2;
                     for z:=1 to n do
                        if (z<>i)and(z<>k)and(z<>j) then
                           begin
                              p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                              if p<t then t:=p;
                           end;
                  end
               else if (p+q<r)and(p+r<q)and(r+q<p) then
                  begin
                     t:=sqrt(q);
                     t1:=sqrt(r);
                     cosG:=(r+q-p)/(2*t*t1);
                     sinG:=sqrt(1-cosG*cosG);
                     if sinG <>0 then
                        begin
                           t:=p/(2*SinG);
                           Dy:=t*cosG;
                           a1:=a[j,2]-a[i,2];
                           b1:=a[i,1]-a[j,1];
                           c1:=a[i,1]*a[j,2]-a[j,1]*a[i,2];
                           a2:=b1;
                           b2:=-a1;
                           c2:=(a1-b1)*((a[i,1]+a[j,1])/2) + (a1+b1)*(a[i,2]+a[j,2])/2;
                           y:=Dy + ((a1*c2-c1*b1)/(a1*a1+b1*b1))/b1;
                           x:=(a1*y-c2)/b1;
                           for z:=1 to n do
                              if (z<>i)and(z<>j)and(z<>k) then
                                 begin
                                    p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                                    if p<t then t:=p;
                                 end;
                        end
                     else t:=-1;
                  end
               else
                  begin
                     t:=sqrt(q);
                     t1:=sqrt(r);
                     cosG:=(r+q-p)/(2*t*t1);
                     sinG:=sqrt(1-cosG*cosG);
                     if sinG <>0 then
                        begin
                           t:=p/(2*SinG);
                           Dy:=t*cosG;
                           a1:=a[j,2]-a[i,2];
                           b1:=a[i,1]-a[j,1];
                           c1:=a[i,1]*a[j,2]-a[j,1]*a[i,2];
                           a2:=b1;
                           b2:=-a1;
                           c2:=(a1-b1)*((a[i,1]+a[j,1])/2) + (a1+b1)*(a[i,2]+a[j,2])/2;
                           y:=Dy + ((a1*c2-c1*b1)/(a1*a1+b1*b1))/b1;
                           x:=(a1*y-c2)/b1;
                           for z:=1 to n do
                              if (z<>i)and(z<>j)and(z<>k) then
                                 begin
                                    p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                                    if p<t then t:=p;
                                 end;
                     a1:=a[j,2]-a[i,2];
                     b1:=a[i,1]-a[j,1];
                     c1:=a[i,1]*a[j,2]-a[j,1]*a[i,2];
                     a2:=y-a[k,2];
                     b2:=a[k,1]-x;
                     c2:=a[k,1]*y-a[k,2]*x;
                     x:=(a[i,1]+a[j,1])/2;
                     y:=(a[i,2]+a[j,2])/2;
                     t:=100000000;
                     for z:=1 to n do
                        begin
                           p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                           if t>p then
                              begin
                                 t:=p;
                                 x1:=x;
                                 y1:=y;
                              end;
                        end;
                     if otg<t then
                        begin
                           otg:=t;
                           otgx:=x1;
                           otgy:=y1;
                        end;
                                             t:=100000000;
                     a1:=a[k,2]-a[i,2];
                     b1:=a[i,1]-a[k,1];
                     c1:=a[i,1]*a[k,2]-a[k,1]*a[i,2];
                     a2:=y-a[j,2];
                     b2:=a[j,1]-x;
                     c2:=a[j,1]*y-a[j,2]*x;
                     x:=(a[i,1]+a[k,1])/2;
                     y:=(a[i,2]+a[k,2])/2;
                     for z:=1 to n do
                        begin
                           p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                           if t>p then
                              begin
                                 t:=p;
                                 x1:=x;
                                 y1:=y;
                              end;

                        end;
                                             if otg<t then
                        begin
                           otg:=t;
                           otgx:=x1;
                           otgy:=y1;
                        end;
                                             t:=100000000;
                     a1:=a[j,2]-a[k,2];
                     b1:=a[k,1]-a[j,1];
                     c1:=a[k,1]*a[j,2]-a[j,1]*a[k,2];
                     a2:=y-a[i,2];
                     b2:=a[i,1]-x;
                     c2:=a[i,1]*y-a[i,2]*x;
                     x:=(a[k,1]+a[j,1])/2;
                     y:=(a[k,2]+a[j,2])/2;

                     for z:=1 to n do
                        begin
                           p:=sqr(x-a[z,1])+sqr(y-a[z,2]);
                           if t>p then
                              begin
                                 t:=p;
                                 x1:=x;
                                 y1:=y;
                              end;

                        end;
                     x:=x1;
                     y:=y1;
                  end;
                  end;
            if t>otg then
                begin
                   otg:=t;
                   otgX:=x;
                   otgY:=y;
                end;
            end;
      otg:=sqrt(otg);
      otg:=otg*100;
      otg:=trunc(otg);
      otg:=otg/100;
      writeln(otg:0:2);
      otgx:=otgx*100;
      otgx:=trunc(otgx);
      otgx:=otgx/100;
      otgy:=otgy*100;otgy:=trunc(otgy);otgy:=otgy/100;
      writeln(otgx:0:2,' ',otgy:0:2);
end.
