function b_spline()

fid=fopen('B_spline.out','r');

x=zeros(101,1);
y=zeros(101,1);

a=zeros(101,2);

[a]=fscanf(fid,'%lf %lf');


for i=1:101
x(i,1)=a(2*i-1,1);
y(i,1)=a(2*i,1);
end


plot(x,y)
hold on;

plot (1550,.031,'+')
hold on;
plot (1100,.044,'+')
hold on;
plot (830,.0498,'+')
hold on;
plot (500,.04,'+')