doublepi(longn){doubles=.;longi;for(i=;i<=n;i++)s=s+?;return(sqrt(s));}