a=[0:2.5:60]
a





x=[0
0.1
0.25
0.5
0.625
0.775
0.9
1
1.15
1.275
]
x=x*40 

y=[0.05775
6.53E-02
0.078125
0.0927291667
0.1066666667
0.080375
0.0443958333
0.0352708333
0.0164166667
0.0092041667
]


%coeff=polyfit(x,y,2);
%n=input('Valeur coeff poly')
%title('avec les donnees')
plot(x,y)



yi=spline(x,y,a)
plot(x,y,'o',a,yi)
r=input('valeur de x pour fonction')
yi=spline(x,y,r)
yi;
%yy=polyval(coeff,r)
%yy;
%title('avec fonction')
%yy=polyval(coeff,a)
%plot(a,yy)
%pid=fopen('test')
%fprintf(pid,yy)
%print -dpsc 'fig.ps'

