mysqrt[aa_]:=Block[{x,i,y,a},
  a = N[aa];
  x = a;
  For[i=0, i<5, i++,
    y = (x+a/x)/2;
    x = y;
  ];
  Return[y];
]

(* Type in <<sqrt2.m  to load *)
