Solving newton interpolation with c

dabillionaire

Hi,

i have no idea how to start up my program on newton interpolation. i have solved it manually but im clueless on how 2 start my program. below represnts the input which i want 2 solve. every 4 lines represnts one problem. the 1st line- n,m where n is no of data points and m is no of evaluation pts. 2nd and 3rd line represent x and y of n pts. 4th line is m @ which it is to be evaluated. the function is sinx. can anyone help me out??


3 2
0.1 0.4 0.5
3.1 3.2 4.1
1.5 1.6
3 5
1.0 2.0 3.0
1.7 1.8 1.1
1.1 1.2 1.4 1.3 1.8
15 5
0.1 0.4 0.5 1.2 2.0 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.95
3.1 3.2 4.1 2.7 1.9 1.7 1.8 1.1 1.7 1.8 1.1 1.7 1.8 1.1 2.3
1.0 1.2 1.6 2.0 2.5

ma output is supposed to be in d form

e1 fN(e1)
e2 fN(e2)
.
.
em fN(em)

dabillionaire

e1, e2 and em represent the fourth line of each problem

mgb_phys

Numerical recipes ( http://www.nr.com/ ) is the classic book on doing this sort of thing.
The code is available online for free. The book is sometimes available depending on what they feel like that day.

You can usually pick up the older editions of the Fortran/Pascal or C versions cheaply.

KTC

Is it?? They have the most restrictive terms for its code one can find for a book. Is the code for the old version now free (in whatever sense)? I can't find it on their website, only links to buy its code cd-rom.

LMAO. How true, how true.

mgb_phys

It's not there anymore, sometimes you can read the whole book online, sometimes you can download the code - they are an amazing bunch.

The licence for the code used to say something like -
"you can't copy this code from the disk, you can type it out of the book but you can't change anything. We own the copyright on the algorithms - you can't use the algorithms and change the code."

There was a bunch of lecturers who boycotted the book and lobbyed CUP to stop publishing it. It would be a bit embarrasing when CUP would have to sue the Cambridge University Library.

KTC

Yeah. It was noted in an ACCU review of the 2nd edition C version that one can't even typed the code into one's own computer if you didn't own the book. i.e. you got it out of the library / workplace etc.

The 2nd edition C version said:
So you have to own the book, you're not allowed to have another person type the code in for you. It's only available for personal use, and only on 1 computer each time you type the code in, specifically the one you typed the code in on. If you want to use the code (for personal use) on more than one computer, you have to type it in each time.

mgb_phys

And worst of all they have the C arrays start from 1 (or they do in my 1988 'C' edition)

jim mcnamara

There were some weird C compilers back in early DOS days that used something called an option-base setting, defined in a header file: 0 = normal C; 1= arrays start at 1 like BASIC. I dunno why.

Since everything in Numerical Recipes is available thru the GNU Scientific library, just pitch the damn thing. GSL is free.

mgb_phys

I think my 'C' edition was a too close copy of the Fortran one.

The explanations of the algorithms are still quite good. Even if you are going to use NAG or GSL it's a good first place to look for which algorithm you should use and what the problems might be.