a=[0:2.5:60]
a





x=[0.15
0.2 
0.275 
0.35	
0.4	
0.475 
0.55	
0.6	
0.7	
0.75	
0.825	
0.875
0.95	
1.05	
1.2	
1.325]
x=x*40 

y=[	1350
	1375
	1400
	1450
	1575
	1950
	2175
        2075
	1800
	1575
	1325
	1100
	850
	575
        400	
        350]


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'

