MatLab / hermite_input

function [xii, f, df]=hermite_input(xi,funcc,proizv)
xii=zeros(1,(xi(2)-xi(1))/(pi/6));
xii(1)=xi(1);
for l=2:size(xii,2)+1
xii(l)=xii(l-1)+pi/6;
end
% xii=linspace(xi(1),xi(2),xi(2)-xi(1))
f=subs(funcc,xii);
df=zeros(proizv,size(f,2));
for i=1:proizv
for j=1:size(f,2)
df(i,j)=subs(diff(funcc,'z',i),xii(j));
end
end
% f=eval(f);