Метод простой итерации

[Artemm]

Помогите еще с таким заданием пожалуйста. Нужно решение в С++

Метод простой итерации найти все действительные корни уравнения
Оригинал: Метод простої ітерації знайти всі дійсні корені рівняння
( E=10^-5 ) ( x-tgx=0 ) xє [ 0, П/2]

[Artemm]

Ау люди , что никто не поможет ?

[Sazary]

Выкладывайте свой код. Иначе тема закроется.

[Artemm]

Вот такой код

#include<math.h>
#include<stdio.h>
float E,xk1,xk,x0;
float f(float x);
float modul(float x);
main()
{
	printf("vvedit E \n");
	scanf("%f",&E);
    printf("vvedit x0 \n");
	scanf("%f",&x0);


	xk=x0;
m1: xk1=f(xk);
	if(modul(xk1-xk)<E && modul(f(xk1)<E)) goto m2; else xk=xk1; goto m1;
	
m2: printf(" korin= %f", xk1);
	printf("\n f(xk1)= %f", x-sin(x)/cos(x));
	getchar();
	getchar();
	return 0;

}
float modul(float x)
{
	if (x<0) x=-x;
	return x;
}
float f(float x)
{
	return x-sin(x)/cos(x);
}

[Sazary]

Что это?
Зачем используете goto?
Зачем функция modul? Есть стандартная fabs.

Вот. Почитайте о методе итераций:
http://ru.wikipedia.org/wiki/Численн...B8.D1.82.D0.BC

В данном случае, как я понимаю, уравнение будет иметь вид: x(i+1) = arctg(x(i))

[Artemm]

нечего не понял , можете подробней обяснить, а лутше воплотить в код , или порядок действий написать , пожалуйста

[Sazary]

Там описан алгоритм.

Выбираете начальное приближение (например, pi/4).
x0 = M_PI/4;
Далее запускаете цикл. В теле цикла:
x1 = atan(x0);
ct = fabs(x1-x0);
x0 = x1;
Условие выхода - ct<E

[Artemm]

Незнаю , не тум тум , неполучяется вообще , О.

[Sazary]

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

Незнаю , не тум тум , неполучяется вообще , О.

Переводится как "напишите все за меня". Думать сами не хотите?
Все написал выше. Какая строчка непонятна?
http://alglib.sources.ru/equations/feqx.php

[Artemm]

Код в Pascal
помогите его переделать в С++ плиз

Program Iteracia;
Uses crt;
Label m1, m2;
Var xk, xk1, x0, eps, fk: real;
Function f1(x:real): real;
Begin
f1:= (x-sin(x)/cos(x));
End;
Function f2(x: real): real;
Begin
f2:=(2x-sin(x)/cos(x));
End;
Begin
Writeln(‘ Введіть значення параметрів x0, eps ’);
Readln( x0, eps);
xk:=x0;
m1: fk:=f2(xk);
xk1:=f2(xk);
If (abs(xk1-xk)<eps and abs(f1(xk1)<eps)) goto m2;
xk:=xk1;
goto m1;
m2: Writeln (‘x=’, xk1, ‘ f(x)=’, f1(xk1));
End.