Решение алгебраических (линейных) уравнений

@АКиреев · Регистрация: 31.01.2012

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

Здесь вышла ошибка "вещественное деление на 0 (Program3.pas, строка 17)". Кто сможет поправить? Заранее большое спасибо!

Pascal

program lab_3;      {метод Лина}
const n=4;
var a4,a3,a2,a1,a0,x,p,q,b4,b3,b2,b1,b0:real;
 
begin
 a4:=1;   {начальные значения коэффициентов}
 a3:=1;
 a2:=0;
 a1:=2;
 a0:=1;
 repeat
  b4:=a4;
  b3:=a3-p*b4;
  b2:=a2-p*b3-q*b4;
  b1:=a1-p*b2-q*b3;
  b0:=a0-p*b1-q*b2;
  p:=(a1-q*b3)/b2;
  q:=a0/b2;
 until ((ABS(a4-b4)<n) AND (ABS(a3-p*b4)-(b3)<n) AND (ABS(a2-p*b3-q*b4)-(b2)<n) AND (ABS(a1-p*b2-q*b3)-(b1)<n) AND (ABS(a0-p*b1-q*b2)-(b0)<n));
 x:=a4*(x*x*x*x)+a3*(x*x*x)+a2*(x*x)+a1*x+a0*1;
   writeln('x=',x:0:5);
end.

@Erea · 13.05.2013, 00:45

При первой попытке выполнения цикла repeat твои переменные принимают значения:
b4 = a4 = 1
b3 = a3 - p*b4 = 1 - 0*1 = 1
b2 = a2 - p*b3 - q*b3 = 0 - 0*1 - 0*1 = 0 (! вот тут собственно и ошибка, потому что потом идет p = (a1 - q*b3)/b2 = (2 - 0*1)/0 - вот и оно, деление на ноль)

Решение - рассчитать или узнать еще каким-либо способом начальные значения p и q и сообщить программе перед началом выполнения цикла.

@barbudo59 · 13.05.2013, 00:47

Сейчас уже не могу вспомнить суть метода Лина (давно это было!), но вот что пишет Его Величество Интернет: "Метод Лина может не привести к нахождению делителя"
Похоже, с оглядкой на строчки 17 и 18 надо сделать оговорку - если b2=0, то решения нет

@Erea 130 / 54 / 27 Регистрация: 13.03.2013 Сообщений: 162
	13.05.2013, 00:45	2
	При первой попытке выполнения цикла repeat твои переменные принимают значения: b4 = a4 = 1 b3 = a3 - pb4 = 1 - 01 = 1 b2 = a2 - pb3 - qb3 = 0 - 01 - 01 = 0 (! вот тут собственно и ошибка, потому что потом идет p = (a1 - qb3)/b2 = (2 - 01)/0 - вот и оно, деление на ноль) Решение - рассчитать или узнать еще каким-либо способом начальные значения p и q и сообщить программе перед началом выполнения цикла. 0

@barbudo59 Я не экстрасенс 382 / 339 / 34 Регистрация: 22.01.2013 Сообщений: 1,126
	13.05.2013, 00:47	3
	Сейчас уже не могу вспомнить суть метода Лина (давно это было!), но вот что пишет Его Величество Интернет: "Метод Лина может не привести к нахождению делителя" Похоже, с оглядкой на строчки 17 и 18 надо сделать оговорку - если b2=0, то решения нет 0

Опции темы