pist:=proc(fu,n)
	local h,ma,j;
	h:=evalf(1/n);
	ma:=Matrix(n+1,2);
	for j to n+1 do
		ma[j,1]:=(j-1)*h;
		ma[j,2]:=evalf(fu((j-1)*h));
	end;
	ma;
end proc;

linspl:=proc(fu,n)
	local ma ;
	ma:=pist(fu,n);
	Spline(ma,x,degree=1);
end proc;

kuspl:=proc(fu,n)
	local ma ;
	ma:=pist(fu,n);
	Spline(ma,x);
end proc;


pist1:=proc(fu,xv)
	local n,ma,j;
	n:=Dimension(xv);
	ma:=Matrix(n,2);
	for j to n do
		ma[j,1]:=xv[j];
		ma[j,2]:=evalf(fu(xv[j]));
	end;
	ma;
end proc;

linspl1:=proc(fu,xv)
	local ma ;
	ma:=pist1(fu,xv);
	Spline(ma,x,degree=1);
end proc;

kuspl1:=proc(fu,xv)
	local ma ;
	ma:=pist1(fu,xv);
	Spline(ma,x);
end proc;