VPF::Найти ошибку - Форум программистов

norman

Дата 17.4.2006, 23:20 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 11
Регистрация: 17.4.2006

Репутация: нет
Всего: нет

Метод Зейделя для системы нелинейных уравнений, Помогите найти ошибку...

При запуске выдаёт ошибку "Floating point operation"...хотя, я никаких немыслимых вычисленей здесь не проводил...
Если кто-нить работал с этим методом или чё-нить про него знает, помогите плиз...

Код

program metod_zaidelya;
uses crt;
var
x,y,x0,y0,e:real;

function f1 (x,y:real):real;
begin
f1:=(sin(x+y)/cos(x+y))-sqr(x)+0.1;
end;

function f2 (x,y:real):real;
begin
f2:=0.6*sqr(x)+2*sqr(y)-1;
end;

procedure solve (x0,y0,e:real); {ОШИБКА, КАК МНЕ КАЖЕТСЯ, ВОЗНИКАЕТ ВОТ В ЭТОЙ ПРОЦЕДУРЕ}
var x1,y1:real;
const k1=0.001; k2=-0.001;
begin
repeat
x1:=x0+k1*f1(x0,y0);
y1:=y0+k2*f2(x1,y0);
x0:=x1;
y0:=y1;
until abs(f1(x0,y0))+abs(f2(x0,y0))<e; {е=0.0000001}
x:=x0;
y:=y0;
end;

begin
clrscr;
write('Введите начальное приближение корня x0:');
read(x0);
write('Введите начальное приближение корня y0:');
read(y0);
write('Задайте точность вычислений e:');
read(e);
writeln;
solve (x0,y0,e);
writeln('|     x    |     y    |    e     |    f1(x,y)     |     f2(x,y)    |');
writeln('|',x:2:10,'|',y:2:10,'|',e:2:10,'|',f1(x,y):2:10,'|',f2(x,y):2:10,'|');
end.

Fixin

Дата 18.4.2006, 18:07 (ссылка)

(нет голосов)

Загрузка ...

Ёжик

Профиль
Группа: Комодератор
Сообщений: 1357
Регистрация: 6.1.2004

Репутация: 5
Всего: 18

Цитата(norman @ 17.4.2006, 23:20

)

Floating point operation

Операция с палавающей точкой. Это разве error?
Я думаю, компилятор выдает ошибку не из-за применения

Цитата(norman @ 17.4.2006, 23:20

)

Метод Зейделя для системы нелинейных уравнений

--------------------

Пользуйтесь ссылкой <<решенный\нерешенный вопрос>>
IRC-канал Винграда

nostromo

Дата 19.4.2006, 09:13 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 194
Регистрация: 23.3.2006

Репутация: 4
Всего: 10

А что в отладчике?
Сколько раз repeat/until прокручивается?

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »