Технический форум
Вернуться   Технический форум > Программирование > Форум программистов > Delphi, Kylix and Pascal


Ответ
 
Опции темы Опции просмотра
Старый 27.03.2016, 02:18   #1 (permalink)
Алекс2
Member
 
Регистрация: 24.03.2016
Сообщений: 25
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию Интерполяционный полином Лагранжа. Паскаль

Выбрать интервал [a,b] і четыре узла интерполяцый в ньом х0, х1, х2, і х3, вычислить значение заданої функцыи у цих точках у0, у1, у2, у3. По этим даним сформировать интерполяционний полином Лагранжа L3;

Составить и отладить програму на паскаль, которая реализует введение входных даных; расчет значения функцыи, заданой таблично, при каком-либо значении n і хр; виведение результатов в удобной форме на екран и в файл;

помогите пожалуйста сделать таблицу значений из этого: y=ln(1+x)\(1-x)
а если можно, то и программу.
Заранее спасибо.
Алекс2 вне форума   Ответить с цитированием

Старый 27.03.2016, 02:18
Helpmaster
Member
 
Аватар для Helpmaster
 
Регистрация: 08.03.2016
Сообщений: 0

Несколько похожих на вашу тем уже создавались

Паскаль
Построить интерполяционный многочлен Лагранжа
Паскаль
Полином

Старый 27.03.2016, 09:21   #2 (permalink)
Vladimir_S
Специалист
 
Регистрация: 27.08.2008
Адрес: Санкт-Петербург
Сообщений: 27,807
Сказал(а) спасибо: 340
Поблагодарили 583 раз(а) в 208 сообщениях
Репутация: 113184
По умолчанию

Цитата:
Сообщение от Алекс2 Посмотреть сообщение
помогите пожалуйста сделать таблицу значений из этого: y=ln(1+x)\(1-x)
Уточните вид функции: знаменатель (1-х) находится под знаком логарифма, т.е.
y = ln[(1+x)/(1-x)]
или же так:
y = [ln(1+x)]/(1-x)?
Кстати, знак деления - это прямой слэш "/", а не обратный.
Ладно, будем считать, что правилен первый вариант. Тогда
Код:
Type
 Arr=Array[0..3] of Real;

Var
 Xi,Yi:Arr;
 a,b,D,X,Y,Step:Real;
 i,j:Byte;
 f:Text;

Function Fun(z:real):Real;
begin
 Fun:=Ln((1+z)/(1-z));
end;

Function Lagr(NL:Byte; XL,YL:Arr; ZL:Real):Real;
var
 S,T:Real;
 k,p:Byte;
begin
 S:=0;
 for k:=0 to NL do
  begin
   T:=1;
   for p:=0 to NL do
    if p<>k then T:=T*(ZL-XL[p])/(XL[k]-XL[p]);
   S:=S+YL[k]*T;
  end;
 Lagr:=S;
end;

Begin
 Writeln(' Attention! -1<a,b<+1; a<b');
 Repeat
  Write(' a = ');
  Readln(a);
  Write(' b = ');
  Readln(b);
 Until (a>-1) and (a<1) and (b>-1) and (b<1) and (a<b);
 D:=(b-a)/3;
 {Formation of Table}
 Writeln(' Number of point     Xi              Yi');
 for i:=0 to 3 do
  begin
   Xi[i]:=a+D*i;
   Yi[i]:=Fun(Xi[i]);
   Writeln(i:8,Xi[i]:16:3,Yi[i]:16:3);
  end;
 Writeln;
 {Calculation of Lagrange interpolation}
 Write(' Step (Step<(b-a))! = ');
 Readln(Step);
 Assign(f,'Result.txt');
 Rewrite(f);
 Writeln('     X       Lagrange interpolation      Accurate value');
 Writeln(f,'     X       Lagrange interpolation      Accurate value');
 X:=a;
 Repeat
  Writeln(X:8:3,Lagr(3,Xi,Yi,X):18:3,Fun(X):26:3);
  Writeln(f,X:8:3,Lagr(3,Xi,Yi,X):18:3,Fun(X):26:3);
  X:=X+Step;
 Until X>b;
 Close(f);
 Readln
End.
Миниатюры
aa01.jpg   lagr.jpg  
Vladimir_S вне форума   Ответить с цитированием
Старый 28.03.2016, 22:11   #3 (permalink)
Алекс2
Member
 
Регистрация: 24.03.2016
Сообщений: 25
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 10
По умолчанию

Большое спасибо за программу!!!
Алекс2 вне форума   Ответить с цитированием
Ads

Яндекс

Member
 
Регистрация: 31.10.2006
Сообщений: 40200
Записей в дневнике: 0
Сказал(а) спасибо: 0
Поблагодарили 0 раз(а) в 0 сообщениях
Репутация: 55070
Ответ

Опции темы
Опции просмотра

Ваши права в разделе
Вы не можете создавать новые темы
Вы не можете отвечать в темах
Вы не можете прикреплять вложения
Вы не можете редактировать свои сообщения

BB коды Вкл.
Смайлы Вкл.
[IMG] код Выкл.
HTML код Выкл.
Trackbacks are Вкл.
Pingbacks are Вкл.
Refbacks are Выкл.




Часовой пояс GMT +4, время: 21:57.

Powered by vBulletin® Version 6.2.5.
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.