Помощь - Поиск - Пользователи - Календарь
Полная версия: хорда и итерация
Форум «Всё о Паскале» > Pascal, Object Pascal > Задачи
candLe
Вопрос вот в чем.написано 2 проги..одна методом хорд другая методом итераций(взяты коды с вашего сайта)

вот уравнение: 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((
но ответы у меня разные(((

помогите студенту плз(((надо до понедельника... wacko.gif
volvo
И где ты нашел разные ответы?

Метод итераций:
Цитата
vvedu x=1
vvedu eps=0.0001
---iteraciya 1 x= 8.31784301934810E-001
F(x)=-2.09492134166151E-002 tochnost= 1.6821569806519032E-0001
---iteraciya 2 x= 8.10451272488642E-001
F(x)=-3.58424238759207E-004 tochnost= 2.1333029446167662E-0002
---iteraciya 3 x= 8.10073328559432E-001
F(x)=-1.13206256414607E-007 tochnost= 3.7794392921042252E-0004
---iteraciya 4 x= 8.10073209112537E-001
F(x)=-1.11717818156165E-014 tochnost= 1.1944689515264173E-0007
Ответ: X= 0.8100732

Метод хорд:
Цитата
vvedu a и b
0 1
vvedu Eps
0.0001
x= 0.8100538
kol-vo iteraciy = 6


Все, что я сделал - заставил в методе итераций печатать ответ в том же формате, в котором он печатается в первой программе:
 if (abs(x-x0)<=eps) then writeln('Ответ: X=',x:10:7)
 else writeln('otvet ne naiden! za ',max_iter:0,' wagov iteraciya ne sowlas.');

candLe
спасибо большое за помощь !zdarov.gif 1.gif
Это текстовая версия — только основной контент. Для просмотра полной версии этой страницы, пожалуйста, нажмите сюда.