【正文】
break。 end if l~=0 break。 for s=l1:l if s~=k a=a*(zx0(s))/(x0(k)x0(s))。 end y(i)=sum。y(%d)=%f\nx1=%.3f y1=%.5f,x2=%.3f y2=%.5f\n\n39。%輸出插值結(jié)果和所需節(jié)點(diǎn)endendend分段二次插值%分段二次插值function y=piece_square(x0,y0,x)% x0,y0為已知點(diǎn),x為待求點(diǎn)n=length(x0)。m=length(x)。Error! Please input again!\n39。% x0和y0長度不等時(shí),報(bào)錯(cuò)elsefor i=1:m z=x(i)。 l=0。Error!x(%d) is out of range!\n39。 break。 if zp l=j。 end%一旦l有非零值,則終止循環(huán),選出合適的l end if l==0 l=n1。 for s=l1:l+1 if s~=k a=a*(zx0(s))/(x0(k)x0(s))。 end y(i)=sum。y(%d)=%f\nx1=%.3f y1=%.5f\nx2=%.3f y2=%.5f\nx3=%.3f y3=%.5f\n\n39。%輸出插值結(jié)果與所需節(jié)點(diǎn)endendend拉格朗日全區(qū)間插值%拉格朗日全區(qū)間插值function y=lagrange(x0,y0,x)% x0,y0為已知點(diǎn),x為待求點(diǎn)n=length(x0)。m=length(x)。Error! Please input again!\n39。%x0和y0長度不等時(shí),報(bào)錯(cuò)elsefor i=1:m z=x(i)。 if zx