added lecture 18

lecture_17/lecture_17.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 65,
+   "execution_count": 4,
    "metadata": {
     "collapsed": true
    },
@@ -13,7 +13,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 66,
+   "execution_count": 5,
    "metadata": {
     "collapsed": true
    },
@@ -2218,7 +2218,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 93,
+   "execution_count": 6,
    "metadata": {
     "collapsed": false
    },
@@ -2385,7 +2385,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 94,
+   "execution_count": 7,
    "metadata": {
     "collapsed": false
    },
@@ -2409,7 +2409,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 95,
+   "execution_count": 8,
    "metadata": {
     "collapsed": false
    },

lecture_18/.ipynb_checkpoints/lecture_18-checkpoint.ipynb

@@ -0,0 +1,6 @@
+{
+ "cells": [],
+ "metadata": {},
+ "nbformat": 4,
+ "nbformat_minor": 2
+}

lecture_18/Newtint.m

@@ -0,0 +1,33 @@
+function yint = Newtint_bak(x,y,xx)
+% Newtint: Newton interpolating polynomial
+% yint = Newtint(x,y,xx): Uses an (n - 1)-order Newton
+%   interpolating polynomial based on n data points (x, y)
+%   to determine a value of the dependent variable (yint)
+%   at a given value of the independent variable, xx.
+% input:
+%   x = independent variable
+%   y = dependent variable
+%   xx = value of independent variable at which
+%        interpolation is calculated
+% output:
+%   yint = interpolated value of dependent variable
+
+% compute the finite divided differences in the form of a
+% difference table
+n = length(x);
+if length(y)~=n, error('x and y must be same length'); end
+b = zeros(n,n);
+% assign dependent variables to the first column of b.
+b(:,1) = y(:); % the (:) ensures that y is a column vector.
+for j = 2:n
+  for i = 1:n-j+1
+    b(i,j) = (b(i+1,j-1)-b(i,j-1))/(x(i+j-1)-x(i));
+  end
+end
+% use the finite divided differences to interpolate
+xt = 1;
+yint = b(1,1);
+for j = 1:n-1
+  xt = xt*(xx-x(j));
+  yint = yint+b(1,j+1)*xt;
+end

lecture_18/challenger_oring.csv

@@ -0,0 +1,24 @@
+Flight#,Temp,O-Ring Problem
+1,53,1
+2,57,1
+3,58,1
+4,63,1
+5,66,0
+6,66.8,0
+7,67,0
+8,67.2,0
+9,68,0
+10,69,0
+11,69.8,1
+12,69.8,0
+13,70.2,1
+14,70.2,0
+15,72,0
+16,73,0
+17,75,0
+18,75,1
+19,75.8,0
+20,76.2,0
+21,78,0
+22,79,0
+23,81,0