Математическая функция в СИ

[прямИСТ]

f = 1 + x^2/2!+x^4/4!+...
Рассчитывается в диапазоне заданном пользователем: start - начальное значение х, finish - конечное значение х, shag - шаг, на который увеличивается х, eps - точность вычисления ряда, ak – член ряда, sum – текущее значение ряда, k – счетчик итераций.
#include<stdio.h>
#include<locale.h>
#include<process.h>
#include<math.h>

void main()
{
double start, finish, shag, eps, ak, sum, x;
int k;
setlocale(LC_ALL, "RUS");
printf("введите точность расчета");
scanf_s("%lf", &eps);
printf("введите нижний порог x");
scanf_s("%lf", &start);
printf("введите верхний порог x");
scanf_s("%lf", &finish);
printf("введите шаг");
scanf_s("%lf", &shag);
for (x = start; x <= finish; x += shag)
{
for (ak = x, sum = ak, k = 0; ak < eps; k += 2, ak *= x*x / (k*(k - 1)))
{
sum += ak;
printf("%lf%lf", x, sum);
}
}
}

[olej.tsil]

Сообщение от прямИСТ

for (x = start; x <= finish; x += shag)
{
for (ak = x, sum = ak, k = 0; ak < eps; k += 2, ak *= x*x / (k*(k - 1)))
{
sum += ak;
printf("%lf%lf", x, sum);
}
}
}

В принципе, идея всё вычислить в заголовке цикла - не очень хорошая, но если очень хочется чтобы было так, то вот примерно:

   for( double x = start; x < finish + eps; x += shag ) {
      double sum = 0.;
      for( double ak = 1., k = 0.; ak > eps; sum += ak, k += 2, ak *= x * x / k / ( k - 1 ) );
      printf( "%f => %f\n", x, sum );
   }

[olej@dell float]$ ./iter 1 2 .1
1.000000 => 1.543081
1.100000 => 1.668519
1.200000 => 1.810656
1.300000 => 1.970914
1.400000 => 2.150898
1.500000 => 2.352410
1.600000 => 2.577464
1.700000 => 2.828315
1.800000 => 3.107473
1.900000 => 3.417731
2.000000 => 3.762196