|
|
Регистрация Восстановить пароль |
Регистрация | Задать вопрос |
Заплачу за решение |
Новые сообщения |
Сообщения за день |
Расширенный поиск |
Правила |
Всё прочитано |
|
Опции темы | Поиск в этой теме |
14.12.2013, 16:16 | #1 |
Пользователь
Регистрация: 24.11.2013
Сообщений: 13
|
Интерполяционный полином ньютона
Доброго времени суток! Мне нужно задать разности(у1,у2,у3) циклом перед основной формулой(N3),но не получается...помогите пожалуйста. заранее спасибо.
uses crt; var xi,yi:array [1..5] of real; y1: array [1..4] of real; y2: array [1..3] of real; y3: array [1..2] of real; i,k: integer; N3,x,q,w: real; Begin clrscr; writeln ('Заданные значения узлов:'); Writeln ('№ | x | y '); writeln ('---------------'); xi[1]:=1.6; yi[1]:=1.6416; xi[2]:=1.7; yi[2]:=2.3961; xi[3]:=1.8; yi[3]:=3.3536; xi[4]:=1.9; yi[4]:=4.5441; xi[5]:=2.0; yi[5]:=6.0000; for i:=1 to 5 do writeln (i-1,' | ', xi[i]:2:2,' | ', yi[i]:2:3); writeln ('Разности первого порядка:'); for i:=1 to 4 do begin y1[i]:=yi[i+1]-yi[i]; write (' y1=', y1[i]:4:4); end; writeln; writeln ('Разности второго порядка:'); for i:=1 to 3 do begin y2[i]:=y1[i+1]-y1[i]; write(' y2=',y2[i]:4:4); end; writeln; Writeln ('Разности третьего порядка:'); for i:=1 to 2 do begin y3[i]:=y2[i+1]-y2[i]; write(' y3=',y3[i]:4:4); end; writeln; repeat writeln ('Введите значение переменной х:'); readln (x); if (x>xi[1]) or (x<xi[5]) then writeln('Ошибка.Заданное значение не попадает в интервал. '); until (x>xi[1]) and (x<xi[5]); for i:=2 to 4 do if (x>xi[i-1]) and (x<xi[i]) then begin k:=i-1; w:=xi[i-1]; end; xi[1]:=xi[k]; for i:=2 to 5 do xi[i]:=xi[1]+(i-1)*0.1; writeln('Новые значения X: '); for i:=1 to 5 do write(xi[i]:5:2); writeln; q:= (x - w)/0.1; writeln('q= ',q:2:3); for i:= 1 to 5 do N3:=yi[1]+q*y1[1]+y2[1]*((q*(q-1))/2)+y3[1]*((q*(q-1)*(q-2))/6); writeln ('Значение интерполируемой переменной = ', N3:4:4); end. |
Похожие темы | ||||
Тема | Автор | Раздел | Ответов | Последнее сообщение |
интерполяционный полином Ньютона построить в делфи | Klik_1602 | Помощь студентам | 2 | 16.06.2011 06:33 |
Интерполяционный полином ньютона-где то ошибка в коде | DiGris | Общие вопросы C/C++ | 1 | 04.05.2010 20:45 |
Интерполяционный полином Ньютона - Visual Basic | gred | Помощь студентам | 1 | 26.02.2009 20:12 |
Интерполяционный полином Ньютона | Skay | Общие вопросы Delphi | 2 | 19.10.2007 22:18 |