h = 0.000008;
f = @(x) exp(x);
f_00 = f(0)
f_01 = f(h)
f_02 = f(2*h)
f_03 = f(3*h)               
f_10 = (f_01 - f_00)/h   %<---f'(0)
f_11 = (f_02 - f_01)/h   %<---f'(h)
f_12 = (f_03 - f_02)/h   %<---f'(2h)
f_20 = (f_11 - f_10)/h   %<------f''(0)   
f_21 = (f_12 - f_11)/h   %<------f''(h)
f_30 = (f_21 - f_20)/h   %<----------f'''(0)
fprintf('%f %f %f %f\n',f_00, f_10, f_20, f_30);
%first_index = 0; second_index = 0:3        
%first_index = 1; second_index = 0:2
%first_index = 2; second_index = 0:1
%first_index = 3; second_index = 0:0
for second_index  = 0:3
    F(0,second_index) = f(second_index * h);
end

for first_index = 1:3
    for second_index = 3 - first_index
        F(first_index, second_index) =  %....
            
    end
end
        
