50737

Формули Нютона через кінцеві різниці

Лабораторная работа

Информатика, кибернетика и программирование

Мета. Навчитися обчислити значення функції при даному значенні аргумента, використовуючи формули Н’ютона через кінцеві різниці. Обладнання. Лист формату А4, ручка, олівець, програмне забезпечення С++.

Украинкский

2014-01-29

108.5 KB

1 чел.

М Чуприна С. 1ПМ-08

Лабороторна робота №6

Тема. Формули Н’ютона через кінцеві різниці

Мета. Навчитися обчислити значення функції при даному значенні аргумента, використовуючи формули Н’ютона через кінцеві різниці.

Обладнання. Лист формату А4, ручка, олівець, програмне забезпечення С++.

Хід роботи

     3. Індивідуальна робота

Використовуючи першу і другу інтерполяційну формули Н’ютона, обчислити значення функції при даних значеннях аргумента. При складанні таблиць різниць контролювати обчислення

Таблиця8

x

y

0,01

0,991824

0,06

0,951935

0,11

0,913650

0,16

0,876905

0,21

0,841638

0,26

0,807789

0,31

0,775301

0,36

0,744120

0,41

0,714193

0,46

0,685470

0,51

0,657902

0,56

0,631442

№ варіанта

х1

х2

18

0,1243

0,492

#include<iostream.h>

#include<math.h>

double N(int n,double* y){

int i,k,L=n;

double*z=new double[n];

for(i=0;i<n;i++)z[i]=y[i+1]-y[i];

if(n>1)

for(i=1;i<n;i++){

for(k=0;k<L;k++)z[k]=z[k+1]-z[k];

L=L-1;}

return z[0];

delete []z;}

double M(int n,int I,double* y){

int i,k,L=n-I;

double*z=new double[I+1];

for(i=n-1;i>=0;i--)z[i]=y[i+1]-y[i];

if(I>1)

for(i=1;i<n;i++){

for(k=0;k<L;k++)z[k]=z[k+1]-z[k];

L=L-1;}

return z[n-I];

delete []z;}

unsigned long p(int n)

{if(n<=1)

return 1;

return n*p(n - 1);}

void main(){

int n,i,k,l=1;

double P,Y,X,h;

cout<<"Введите количество входных данных:";cin>>n;

double*x=new double[n];

double*y=new double[n];

cout<<"Введите узлы интерполЯции:\n";

for(i=0;i<n;i++){cout<<" X["<<i+1<<"]: ";cin>>x[i];}

cout<<"Введите соответствующие значениЯ:\n";

for(i=0;i<n;i++){cout<<" Y["<<i+1<<"]: ";cin>>y[i];}

cout<<"Введите значение Х:";cin>>X;

if(X<x[1]&&X>x[0]){

P=y[0];

h=x[1]-x[0];

for(i=1;i<n;i++){

Y=N(i,y);

for(k=0;k<i;k++)Y*=(X-x[k]);

P+=Y/(p(i)*pow(h,i));}

cout<<"По первой интерполЯционной\n формуле f("<<X<<") = "<<P;}

Висновок: я навчилася обчислюваити значення функції при даному значенні аргумента, використовуючи формули Н’ютона через кінцеві різниці.

if(X>x[n-1]&&X<x[n]){

P=y[n];

h=x[1]-x[0];

for(i=1;i<n;i++){

Y=M(n,i,y);

for(k=n;k>=n-i;k--)Y*=(X-x[k]);

P+=Y/(p(i)*pow(h,i));}

cout<<"По второй интерполЯционной\n формуле f("<<X<<")="<<P;}

delete []x;

delete []y;}


 

А также другие работы, которые могут Вас заинтересовать

72144. Подведомственность и подсудность 22.5 KB
  Между ними дела распределяются следующим образом. К ведению общих судов относятся независимо от того из каких правоотношений возник спор все дела одной из сторон в которых выступает гражданин. Арбитражные суды рассматривают преимущественно дела связанные с разрешением хозяйственных споров...
72145. Стадии гражданского процесса 23 KB
  Стадия гражданского процесса – совокупность процессуальных действий, направленных к одной ближайшей цели: принятие заявления, подготовка дела к судебному разбирательству, судебное разбирательство и т.д.
72146. Понятие и принципы уголовно-процессуального права 36 KB
  В наибольшей степени регламентирована правовыми нормами уголовно-процессуальная деятельность. Это обусловлено, с одной стороны, тем, что из всех правонарушений именно преступления представляют наибольшую опасность для общества и особенно важно их раскрытие и наказание виновных...
72147. Стадии уголовного процесса 26.5 KB
  Процессуальные стадии это сменяющие друг друга части ступени развития юридического процесса. Так в уголовном процессе обычно выделяют следующие стадии: возбуждения уголовного дела; предварительного расследования; подготовки к судебному разбирательству; судебного разбирательства; кассационного рассмотрения...
72148. Володимир Володимирович Маяковський 79 KB
  Маяковському, щоб описати краще сучасників, необхідно було вчитися майстерності. І він вирішує покинути ряди партії, щоб знаходитися на легальному становищі. Першим ділом береться за живопис, навчається у Жуковського. Через рік починає освоювати рукоділля у Келіна.
72149. Основные виды специальных ПР-мероприятий 42 KB
  Товар на выставке можно показать в действии рассказать о нем подробно и именно так как хотите вы. Показ нового только что созданного вами изделия на коммерческой выставке простой и дешевый способ запустить его в производство найти покупателей идеи или инвесторов.
72151. Архивное хранение документов и дел 29.39 KB
  Система хранения информации стала формироваться одновременно с системой делового письма. Первая опись датирована 1288 годом. (Ипатьевская летопись) В период приказов места хранения документов не были стабильными.
72152. Движение документов в организации 19.21 KB
  Порядок движения документов в организации выделяются несколько этапов: Экспедиционная обработка документов поступающих в организацию. Предварительное рассмотрение документов службой ДОУ. Организация рационального движения документов внутри организации.