You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
69 lines
2.0 KiB
C++
69 lines
2.0 KiB
C++
|
|
|
|
/*-------------- DIVIDED DIFFERENCES GENERATION FOR INTERPOLATION -----*/
|
|
|
|
/* THE PROGRAM GENERATES A DIVIDED DIFFERENCES TABLE FROM GIVEN
|
|
|
|
DATA. THE VALUES OF x AND CORRESPONDING y = f(x) ARE TO BE
|
|
|
|
ENTERED IN THE ARRAY FORMAT.
|
|
|
|
INPUTS : 1) Number of entries of the data.
|
|
|
|
2) Values of 'x' & corresponding y = f(x).
|
|
|
|
OUTPUTS : divided difference table array. */
|
|
|
|
/*------------------------------ PROGRAM --------------------------*/
|
|
|
|
#include<stdio.h>
|
|
#include<math.h>
|
|
#include<stdlib.h>
|
|
#include<conio.h>
|
|
|
|
void main()
|
|
{
|
|
double y[20][20],x[20];
|
|
/* ARRAY OF y[n][n] ELEMENTS FOR DIVIDED DIFFERENCE TABLE */
|
|
int i,j,k,n;
|
|
clrscr();
|
|
printf("\n\tDIVIDED DIFFERENCES GENERATION FOR INTERPOLATION");
|
|
printf("\n\n The form of equation is y = f(x)\n");
|
|
printf("\n\nEnter the number of entries (max 20) = ");
|
|
/* ENTER THE NUMBER OF ENTRIES IN THE TABLE */
|
|
scanf("%d",&n);
|
|
for(i = 0; i < n; i++)
|
|
{
|
|
/* LOOP TO GET x AND y = f(x) IN THE TABLE */
|
|
printf("x%d = ",i);
|
|
scanf("%lf",&x[i]);
|
|
printf(" y%d = ",i);
|
|
scanf("%lf",&y[i][0]);
|
|
}
|
|
k = n;
|
|
for(j = 1; j < n; j++)
|
|
{
|
|
/* LOOP TO CALCULATE DIVIDED DIFFERENCES IN THE TABLE*/
|
|
k = k - 1;
|
|
for(i = 0; i < k; i++)
|
|
{
|
|
y[i][j] = (y[i+1][j-1] - y[i][j-1])/(x[i+j]-x[i]);
|
|
}
|
|
}
|
|
k = n;
|
|
clrscr();
|
|
printf("\n The divided difference table is as follows . . .\n");
|
|
printf("\n\tx\ty\tDy\tD2y\tD3y\tD4y\tD5y\tD6y\n");
|
|
for(i = 0; i < n; i++)
|
|
{
|
|
/* LOOP TO PRINT DIVIDED DIFFERENCES IN THE TABLE */
|
|
printf("\nx%d = %4.2lf",i,x[i]);
|
|
for(j = 0; j < k-i; j++)
|
|
{
|
|
printf("\t%4.2lf ",y[i][j]);
|
|
}
|
|
printf("\n");
|
|
}
|
|
}
|
|
/*-------------------------------End of program -----------------------*/
|