для уравнения 15*x^4-4*x^3-6*x^2-4*x-1=0
написать решение методом касательных и хорд с точность 10^-5
производная первого порядка 60*x^3-12*x^2-12*x-4=0
Полная версия: задачи на численное решение нелинейных уравнений
Что, очень трудно посмотреть FAQ ПРЕЖДЕ чем задать вопрос? Обязательно плодить одни и те же темы еще и еще раз? 
FAQ: Численное решение и визуализация
FAQ: Численное решение и визуализация
эту тему в FAQ я видел, но я не понимаю куда вставить свое уравнение. а если взять решение из Faq то прога выдает ощибку 8087 в
извините, но в паскале не разбираюсь
Код
Begin
F := Sqr(Arg) - 2
End;
F := Sqr(Arg) - 2
End;
извините, но в паскале не разбираюсь
напиши первой строчкой в программе
подключение математического сопроцессора для работы с вещественными числами.
{$n+}
подключение математического сопроцессора для работы с вещественными числами.
исправте если неправильно пожайлуста. метод хорд для уравнения выше.
Исходный код
program xord;
uses Crt;
const
eps = 0.00001;
var
lB, rB, x, C, C2 : real;
function Fun(x : real) : real;
begin
Fun := (15*x*x*x*x) - (4*x*x*x)-(6*x*x)-(4*x) - 1;
end;
function GetC(a,b : real) : real;
begin
GetC := (Fun(B)*a - Fun(a)*B) / (Fun(B) - Fun(a));
end;
begin
lB := -1;
rB := 0;
ClrScr;
repeat
C := GetC(lB,rB);
if (Fun(lB)*Fun© > 0)
then lB := C
else rB := C;
WriteLn('X = ',C:10:10);
C2 := GetC(lB,rB);
until (abs(C2 - C) < eps);
ReadLn;
end.
uses Crt;
const
eps = 0.00001;
var
lB, rB, x, C, C2 : real;
function Fun(x : real) : real;
begin
Fun := (15*x*x*x*x) - (4*x*x*x)-(6*x*x)-(4*x) - 1;
end;
function GetC(a,b : real) : real;
begin
GetC := (Fun(B)*a - Fun(a)*B) / (Fun(B) - Fun(a));
end;
begin
lB := -1;
rB := 0;
ClrScr;
repeat
C := GetC(lB,rB);
if (Fun(lB)*Fun© > 0)
then lB := C
else rB := C;
WriteLn('X = ',C:10:10);
C2 := GetC(lB,rB);
until (abs(C2 - C) < eps);
ReadLn;
end.
Что может быть не правильно ? Все что тебе надо было сделть, это вставить свое уравнение и граници поиска. 
Если ты не изменял ничего кроме этого, то все правильно.
Если ты не изменял ничего кроме этого, то все правильно.
" x " получается очень много, уравнение надо простое брать или производную от уравнения?. и еще вопрос :"что такое приближение в методе Ньютона"
Цитата(mafia1232 @ 29.05.05 11:34)
" x " получается очень много, уравнение надо простое боать или производную от уравнения?
Если бы исходный текст не был изменен, таких вопросов не было бы... А в измененном тексте разбирайтесь самостоятельно... Сумели изменить - ...
справился с хордми по вашему примеру, но так и не понял что такое начальное приближение для метода касательных и как его вычислить. Начальное приближение случайно не начальная точка отсчета до корня или начало отрезка в который входит корень?
Цитата(mafia1232 @ 29.05.05 11:52)
начало отрезка в который входит корень?
:yes:
Слушай, ну что с тобой ? там же дальше черным по белому написано
Метод Ньютона(касательных)
прошу прощения, не допонял, ты не не нашел, а не понял в чем там соль да ? ;)
Метод Ньютона(касательных)
прошу прощения, не допонял, ты не не нашел, а не понял в чем там соль да ? ;)
тогда последний вопрос: что дописать к методу Ньютона (основа из вашего макета) что бы он показывал корень, а не переходил к окну программы команда
а вот сама прога
Код
readln;
не помогает. а вот сама прога
Исходный код
{$n+}
Program metodnewton;
Function newton(start, Eps: Extended): Extended;
Var
X, prev: Extended;
{ function }
Function F(Arg: Extended): Extended;
Begin
F :=(15*x*x*x*x)-(4*x*x*x)-(6*x*x)-(4*x)-1
End;
{ proizvodnaya }
Function Deriv(Arg: Extended): Extended;
Begin
Deriv := (60*x*x*x)-(12*x*x)-(12*x)-4
End;
Begin
X := start;
Repeat
prev := X;
X := prev - F(prev) / Deriv(prev);
Until Abs(X - prev) <= Eps;
newton := X
End;
Var a, Eps: Extended;
begin
WriteLn('vvedite priblegenie a');
Read(a);
WriteLn('vvedite to4nost Eps');
Read(Eps);
WriteLn('koren raven x= ', newton(a, Eps)
7);
readln;
end.
Program metodnewton;
Function newton(start, Eps: Extended): Extended;
Var
X, prev: Extended;
{ function }
Function F(Arg: Extended): Extended;
Begin
F :=(15*x*x*x*x)-(4*x*x*x)-(6*x*x)-(4*x)-1
End;
{ proizvodnaya }
Function Deriv(Arg: Extended): Extended;
Begin
Deriv := (60*x*x*x)-(12*x*x)-(12*x)-4
End;
Begin
X := start;
Repeat
prev := X;
X := prev - F(prev) / Deriv(prev);
Until Abs(X - prev) <= Eps;
newton := X
End;
Var a, Eps: Extended;
begin
WriteLn('vvedite priblegenie a');
Read(a);
WriteLn('vvedite to4nost Eps');
Read(Eps);
WriteLn('koren raven x= ', newton(a, Eps)
7);readln;
end.
Все работает, скорее всего ты вводишь "плохие данные", и где-нибудь получается делиние на 0 или корень из отрицательного ...
а по поводу привильно ли счетает, это опять же, смотря что ты исправлял.
а по поводу привильно ли счетает, это опять же, смотря что ты исправлял.
ты не понял он выходит не посредством ошибки, а просто корень не показывает, а показывает только если ты прогу второй раз запускаешь. А считает она правильно.
Щас в Борланде запущу.
поставь
хмм странно, не получислось, тогда
поставь
uses crt;
второй строчкойхмм странно, не получислось, тогда
uses crt;
в начале, а вместо
readln;
-
repeat until keypressed;
вот прога pascal
Спасибо klem4 ! наконец работает как должна. а куда вставить clrscr?
Или поменяй Read на ReadLn... Это тоже решит проблему...
Цитата(mafia1232 @ 29.05.05 12:35)
Спасибо klem4 ! наконец работает как должна. а куда вставить clrscr?
После основного
Begin
Цитата
Или поменяй Read на ReadLn... Это тоже решит проблему
там как раз и был readln
:no: Ты не понял... Здесь:
WriteLn('vvedite priblegenie a');
Read(a); { <-- Поставь ReadLn }
WriteLn('vvedite to4nost Eps');
Read(Eps); { <-- Поставь ReadLn }
Не мучайте модуль crt
Проще нажать Alt+F5 и все увидите! :D
Проще нажать Alt+F5 и все увидите! :D
Цитата(hiv @ 31.05.05 15:17)
Не мучайте модуль crt
Проще нажать Alt+F5 и все увидите! :D
Проще нажать Alt+F5 и все увидите! :D
тут не согласен, его желательно всегда подключать :yes: , моного лажи вылезает, когда его нету.
Он всетаки стандартный модуль для клавиатуры и дисплея.
Случай выше - пример, программа работает не так как написно.
У меня нормально работает и без модуля CRT
Это текстовая версия — только основной контент. Для просмотра полной версии этой страницы, пожалуйста, нажмите сюда.