In [1]:
;type "lin_interp.jl"
In [2]:
include("lin_interp.jl")
Out[2]:
In [3]:
grid = [0, 2, 4, 6, 8, 10]
vals = [1, 4, 5, 8, 9, 11]
f = my_lin_interp(grid,vals)
Out[3]:
In [4]:
println(f([3.5,10,-5,20]))
In [5]:
# 齊藤一織です。
In [6]:
g(x)=2 .*cos(6x) .+sin(14x) .+ 2.5
Out[6]:
In [7]:
c_grid = 0:.2:1
Out[7]:
In [8]:
Af = my_lin_interp(c_grid,g(c_grid))
Out[8]:
In [9]:
residual(x)=g(x)-Af(x)
Out[9]:
In [10]:
using PyPlot
In [11]:
plot(f([0,1,2,3,4,5,6,7,8,9,10,11]),"b-")
Out[11]:
In [17]:
plot(residual([-3,-2,-1,0,1,2,3]),"b-")
Out[17]:
In [ ]:
In [ ]: