Решение уравнений метод Зейделя

01.12.2008, 13:53. **Показов** 4903. **Ответов** 8

Студворк — интернет-сервис помощи студентам

Мне надо методом Зейделя решить такое уравнение
x1^2=x2^3
x2+3*sinx1=0.
Вроде запускается, но после нескольких запусков выдает ошибку в строчке x11:=sqrt(x2*x2*x2);

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
var
E,x1,x2,F1,F2,x11,x21:real;
k:integer;
label M,M1;
 
begin
clrscr;
writeln(' METOD Zeidelya ');
writeln;
write('Vvedite kolichestvo iteracii k= ');
readln(k);
write('Vvedite tochnost E= ');
readln(E);
x1:=0;
x2:=0;
x1:=x1+k;
x2:=x2+k;
M1: x11:=sqrt(x2*x2*x2);
x21:=-3*sin(x1);
F1:=x11*x11-x21*x21*x21;
F2:=x21+3*sin(x1);
if(abs(F1)<=E) and (abs(F2)<=E) then goto M;
x1:=x11;
x2:=x21;
k:=k+1;
goto M1;
M:
writeln;
writeln('F1= ',F1);
readln;
writeln('F2= ',F2);
readln;
end.

01.12.2008, 14:59

Если кто-то знает напишите плиззззз

@mik-a-el · 01.12.2008, 15:08

Какую ошибку выдает?

01.12.2008, 15:16

"Invalid floating point operation"

@Puporev · 01.12.2008, 16:12

Тут у Вас что-то с математикой.
В определенный момент x21:=-3*sin(x1); становится отрицательным (при положительном sin(x1).
x2:=x21; тоже отрицательно.
x11:=sqrt(x2*x2*x2); не может быть выполнено, поэтому здесь ошибка.

01.12.2008, 16:26

А не подскажете как этого избежать?

@Puporev · 01.12.2008, 16:52

1. Нужно установитьить область определиня sin(x1)<=0, это если не ошибаюсь
[(2n-1)*pi,2n*pi] и писать если то, то это, а если не то, тогда другое.
2. Непонятно зачем Вы к переменной прибавляете К-количество итераций? По-моему надо прибавлять что-то типа шага, я просто не знаю этот метод, но суть все-равно не меняется, зачем прибавлять 1,2,3... Итерации считают, а не прибавляют к переменной, они у Вас даже разного типа. Короче разбирайтесь.

03.12.2008, 09:55

Переделала по теории. Но теперь у меня не меняется вообще значение F. Что я сделала не правильно?
Можете помочь?

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
var
E,x1,x2,F1,F2,x11,x21:real;
k:integer;
label M,M1;
 
begin
write('Vvedit kilkist iteracii k= ');
readln(k);
write('Vvedit tochnist E= ');
readln(E);
k:=0;
 
M1: x1:=0;
x2:=0;
x11:=sqrt(x2*x2*x2);
x21:=-3*sin(x1);
F1:=x1*x11-x21*x21*x21;
F2:=x21+3*sin(x11);
 
if(abs(F1)<=E) and (abs(F2)<=E) then
begin
goto M;
end
else
x1:=x11;
x2:=x21;
k:=k+1;
goto M1;
 
M:
writeln('F1= ',F1);
readln;
writeln('F2= ',F2);
readln;
end.

03.12.2008, 14:00

Подскажите, пожалуйста, от чего должны зависить х1 и х2

Новые блоги и статьи Все статьи Все блоги /
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .	Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .
Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .

KatyaKKK
	01.12.2008, 14:59
	Если кто-то знает напишите плиззззз

@mik-a-el Администратор 87822 / 53143 / 249 Регистрация: 10.04.2006 Сообщений: 13,761
	01.12.2008, 15:08
	Какую ошибку выдает? 0

@Puporev Почетный модератор 64314 / 47610 / 32743 Регистрация: 18.05.2008 Сообщений: 115,168
	01.12.2008, 16:12
	Тут у Вас что-то с математикой. В определенный момент x21:=-3sin(x1); становится отрицательным (при положительном sin(x1). x2:=x21; тоже отрицательно. x11:=sqrt(x2x2*x2); не может быть выполнено, поэтому здесь ошибка. 0

KatyaKKK
	01.12.2008, 16:26
	А не подскажете как этого избежать?

@Puporev Почетный модератор 64314 / 47610 / 32743 Регистрация: 18.05.2008 Сообщений: 115,168
	01.12.2008, 16:52
	1. Нужно установитьить область определиня sin(x1)<=0, это если не ошибаюсь [(2n-1)pi,2npi] и писать если то, то это, а если не то, тогда другое. 2. Непонятно зачем Вы к переменной прибавляете К-количество итераций? По-моему надо прибавлять что-то типа шага, я просто не знаю этот метод, но суть все-равно не меняется, зачем прибавлять 1,2,3... Итерации считают, а не прибавляют к переменной, они у Вас даже разного типа. Короче разбирайтесь. 0

KatyaKKK
	03.12.2008, 14:00
	Подскажите, пожалуйста, от чего должны зависить х1 и х2