function T = INV_init(x_min,del_x,x_max,T_min,f_init,T_max)
%

x = [x_min:del_x:x_max];
num_pts = length(x);
T = zeros(num_pts,1);
T(1) = T_min;
T(num_pts) = T_max;
for p = 2:num_pts-1
    T(p) = feval(f_init,x(p));
end
