Вычисление суммы ряда через функции

считает сумму ряда.
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
FILE *f;
void start (){
 f=fopen("tablica.txt", "a");
  }
const double pi = 3.1415926535897932384626433832795;
double S (double x, double E){
    int k=1;
    double y, sum=0;
    double T=1;
  while (fabs(T)>E)
   {    sum=sum+T;
        T=T*(1/(((2*k-1)*(2*k-1))+x*x));
        k=k+1;
    return sum;
   }
   y=(4*x/pi)*sum;
}
 
void finish (){
        fprintf(f,"-----------------------------------\n");   
        fclose (f);
}
int main(int argc, char *argv[]) {
    start ();
    double x, Xn, Xk, dx, E, y;
  printf("Vvedite nachal'noe znachenie argumenta: \n");
  scanf("%lf", &Xn);
  printf("Vvedite konechnoe znachenie argumenta: \n");
  scanf("%lf", &Xk);
  printf("Vvedite shag izmeneniya argumenta: \n");
  scanf("%lf", &dx);
  printf("Vvedite tochnost' vychislenij: \n");
  scanf("%lf", &E);
    for (x=Xn; x<Xk; x=x+dx) {y=S(x, E);};
    printf("y="); printf("%lf\n", y);
    finish ();
    system("PAUSE");    
  return 0;
}