X3 = 0;
 
for n=1:10
   X1 = X3;
   Y1 = Thermistor(X1);
   X2 = X1 + 0.01;
   Y2 = Thermistor(X2);
   X3 = X2 - (X2-X1)/(Y2-Y1)*Y2;
   disp([n, X1, Y1]);   
   X1 = X3;
   end