clear;
f = input('请输入函数表达式:f(x) = ', 's');%注's',表明允许用户输入一个字符串
a = input('请输入区间左端值a:');
b = input('请输入区间右端值b:');
n = input('请输入区间等分值n:');
for i=1:n+1
x(i) = a + (b-a)/n*(i-1);
y(i) = eval(subs(f,'x(i)','x'));
end
n=n+1;
lamda(1)=1;%构造向量
miu(n)=1; %构造向量
h=diff(x); %若X为向量,Y = diff(X)= [X(2)-X(1),X(3)-X(2),...,X(n)-X(n-1)]
df=diff(y)/diff(x);
d(1)=6*(df(1)-1/2)/h(1);
d(n)=6*(0.5*81^-0.5-df(n-1))/h(n-1);
for j=2:n-1
lamda(j)=h(j)/(h(j-1)+h(j));
miu(j)=h(j-1)/(h(j-1)+h(j));
d(j)=6*(df(j)-df(j-1))/(h(j-1)+h(j));
end
miu=miu(2:end);
u=diag(miu,-1);r=