function [retval] = s103021110 (input1, input2)
format long;
d = [0.3 exp(-.6); .4 exp(-.8); .5 exp(-1); .6 exp(-1.2)];
dd = [0.3 .5488; .4 .4493; .5 .3679; .6 .3012];
l = zeros(4,4);

%%%%%%%%%%%%%%%%%%%%%%%%%
for i = 1:4
  l(i,1) = d(i,2);
end

for i = 1:3
  l(i,2) = ( l(i+1,1)-l(i,1) )/( d(i+1,1)-d(i,1) );
end

for i = 1:2
  l(i,3) = ( l(i+1,2)-l(i,2) )/( d(i+2,1)-d(i,1) );
end

for i = 1:1
  l(i,4) = ( l(i+1,3)-l(i,3) )/( d(i+3,1)-d(i,1) );
end
%disp(l);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

f = @(x) l(1,1) + l(1,2)*(x-d(1,1)) + l(1,3)*(x-d(1,1))*(x-d(2,1)) + l(1,4)*(x-d(1,1))*(x-d(2,1))*(x-d(3,1))-x;

%%%%%%%%%%%
yy = f(d(1,1));
y = f(1.3);
%disp(yy);

for i = 04:10^-6:2
  if f(i)<0
    disp(i);
    break
  end
end

r = f(i-10^-6);
disp(r);
  
endfunction