вот уравнение: tg(0,47x+0,2)=x*x;
eps=0.0001;
вот текст хорды:
{$n+}
Uses Crt;
Var
x, Eps, a, b, c: Double;
n: Integer;
Function f(x: Double): Double;
Begin
f := (sin(0.47*x+0.2))/(cos(0.47*x+0.2))-x*x
End;
Function Chord(a, b: Double): Double;
Begin
Chord := (f(b)*a - f(a)*b) / (f(b) - f(a))
End;
begin
ClrScr;
Writeln('vvedu a и b');
Read(a, b);
WriteLn('vvedu Eps');
Read(Eps);
ClrScr;
n := 0;
Repeat
C := Chord(a, b);
If f(a)*f(c) > 0 Then a := c
Else b := c;
Inc(n)
Until Abs(Chord(a, b) - C) < Eps;
x := c;
WriteLn(' x=', x:10:7);
WriteLn('kol-vo iteraciy = ',n);
ReadKey
end.
вот итерация:
program ITERAT;
uses crt;
const max_iter=100;
var
i :integer;
x,x0,eps,M :real;
function F(x:real):real;
begin
F:=(sin(0.47*x+0.2))/(cos(0.47*x+0.2))-x*x;
end;
begin
clrscr;
write('vvedu x='); readln(x);
write('vvedu eps='); readln(eps);
i:=0;
repeat
M:=-(F(x+eps)-F(x-eps))/(2*eps); {коэффициент для улучшения сходимости}
x0:=x;
x:=x0+F(x0)/M; {сердце метода итераций}
inc(i);
writeln('---iteraciya',i:3,' x=',x);
writeln('F(x)=',F(x),' tochnost=',abs(x-x0));
until (abs(x-x0)<=eps)or(i>max_iter);
if (abs(x-x0)<=eps) then writeln('Ответ: X=',x)
else writeln('otvet ne naiden! za ',max_iter:0,' wagov iteraciya ne sowlas.');
end.
Препод говорит что ответ должен быть один и тот же..максимум он должен отличаться на 0.0001((
но ответы у меня разные(((
помогите студенту плз(((надо до понедельника...
