function y=fel(x);
% file name: fel.m
% residual function to be used for fitting of 
% the Arrhenius rate function to measured data
global rate temp;
a=rate-x(1)*exp(x(2)-x(2)./temp); 
y=a'*a;


