Сумма ряда в С

[Eovin]

IDE Pelles C, задание :
Вычислить сумму ряда (приближенное значение функции arctg для |x|<1) с точностью eps. Значения x и eps вводятся с клавиатуры. Полученный результат и значение библиотечной функции arctg вывести на экран.
ряд.bmp - собсно ряд.

#include<stdio.h>
#include<math.h>
int main (int argc, char *argv[])
{      int i;
	double arctg,a,x,t,eps;
	printf("vvedite x,eps\n");
	scanf("%lf %lf", &x,&eps);
	if (x>1 || x<-1) 
               {puts("Enter -1<x<1");
		         return 0;
	        }
	
	else {   arctg=0;
		   i=0;
		   double _arctg=x;
                   t=x;
		   while (labs(arctg-_arctg)>eps) 
			{      arctg=_arctg;
				i=i+1;
				t=t*x*x;
				if(i&1) _arctg=_arctg-(t/(2*i+1));
				else _arctg=_arctg+(t/(2*i+1));
			}
	        }
	arctg=floor(arctg/eps)*eps;
	printf("Znachenie funkcii = %lf\n", arctg); 
	return 0;
}

Вопрос заключается в следующем: мне нужно, например, вывести результат 0,12345678 с точностью 0,00001, т.е. 0,12346. У меня же получается 0,123460. А если с точностью 0,0000001 - имеем 0,123457. Как-то так, криво. В общем помогите нормальный вывод сделать.
И еще - что за "значение библиотечной функции arctg "? Это просто вывести значение арктангенса введенного числа?

[Vago]

Вы, вообще, на результат, выводимый своей программой, смотрели? Он Вас не смутил?.. У меня, например, выдаётся 0, какое бы x я ни вводил... А i, заодно с результатом, распечатать не пробовали? То, что оно остаётся нулём, какую бы точность ни задали, не удивило?..

[Eovin]

Ого, нет не смотрела, признаюсь =/ и в чем тут проблема?

[Vago]

Значит, у кого-то из нас - "неправильный мёд"... У Вас, похоже, транслятор (long) double превращает в long double, а у меня - в long int. Чтобы не было недоразумений, labs() на fabs() поменяйте, от греха подальше.

Сообщение от Eovin

arctg=floor(arctg/eps)*eps;

Это, вообще-то, опасный фокус. И, чем меньше eps, тем опаснее. Просто, уберите эту строку вообще. Точность у Вас и так достигается.

Сообщение от Eovin

И еще - что за "значение библиотечной функции arctg "?

Библиотечная ф-я atan()

[Eovin]

0 там получается если ввести 1, значит не работает нормально проверка на ввод -1<x<1.
И насчет точности: вот если я задаю например 0,00000001, то все равно в ответе получаю число с точностью 0,000001.

[Vago]

Ещё раз. Медленно. По пунктам...

1. 0 (и i=0) я получаю всегда, кроме крайних значений (-1 и 1). Причина этого - использование Вами для принятия решения о выходе из итерационного цикла ф-ии labs(), которая ожидает аргументя типа long. О том, что будет конверсия с непредсказуемым р-том, M$VC 6 меня честно предупредил на этапе трансляции. Как только я ставлю туда fabs(), программа начинает считать нормально.

2. В подправленной программе для -1 и +1 у меня всё считается нормально (-0.785 и 0.785 соответственно).

Сообщение от Eovin

И насчет точности: вот если я задаю например 0,00000001, то все равно в ответе получаю число с точностью 0,000001.

3.а Ещё раз повторяю - уберите floor(),
3.b Укажите явно длины дробной части в полях форматов вывода, например "%.10g".

[Eovin]

Удалила, поменяла, причем с первого раза. Ну вроде как работает, так у меня условие -1<x<1, значит х не может равняться -1 и 1. Я или условия неправильно записала или уже запуталась. А явно размер я задавала - ну так тогда нули остаются, я ж заранее не знаю какую точность пользователь выберет.

[Vago]

Сообщение от Eovin

у меня условие -1<x<1, значит х не может равняться -1 и 1.

Ну так включите тогда концы интервала в условие, отсекающее неверно введённые значения x:

if ( x >= 1 || x <= -1 )

[Eovin]

Да, сама уже нашла))
Стыдно-то как - такую мелочь не заметила...)
Ну и в итоге все равно выводит не так как мне хочется

[Vago]

Сообщение от Eovin

Ну и в итоге все равно выводит не так как мне хочется

   int nDecimalDigits = 1 + (int) log10(1./eps);
   printf( "Znachenie funkcii = %.*g\n", nDecimalDigits, arctg );