x=[0.1 0.2 0.6 0.8 1];
y=[0.15 0.15 0.25 0.3 0.5];
a1=[1 x(1)
1 x(2)]
b1=[y(1)
y(2)]
c1=inv(a1)*b1
a2=[1 x(2)
1 x(3)];
b2=[y(2)
y(3)];
c2=inv(a2)*b2
a3=[1 x(3)
1 x(4)];
b3=[y(3)
y(4)];
c3=inv(a3)*b3;
a4=[1 x(4)
1 x(5)];
b4=[y(4)
y(5)];
c4=inv(a4)*b4
function[Ylin]=YL(t)
if t<=x(2) then
Ylin=c1(1)+c1(2)*t else
if t<=x(3) then
Ylin=c2(1)+c2(2)*t else
if t<=x(4) then
Ylin=c3(1)+c3(2)*t else
Ylin=c4(1)+c4(2)*t
end;
end;
end;
endfunction;
subplot(4,2,3);
plot2d(t,YL)
plot2d(x,y,-4)
xgrid();
xtitle('кусочно-линейная интерполяция','X','Y');
Yinterp=interpln(z,t);
subplot(4,2,4);
plot2d(t,Yinterp);
plot2d(x,y,-4)
xgrid();
xtitle('кус-лин интерполяция interpln(z,t)','X','Y')