format short
%syms x
%plot(x,y1),grid
%plot(x,y1,x,g1),grid
%h=((48*0.5e-4)/135)^(1/4)
h=0.02
x=1:h:1.12
y1=exp(-1./x)
g1=(1./x.^2).*exp(-1./x)
G1=[4 1 0 0 0;1 4 1 0 0;0 1 4 1 0;0 0 1 4 1;0 0 0 1 4]
%f1=[(3/h)*(((1/24)*(x+h).^6-(13/8)*(x+h).^2)-((1/24)*(x-h).^6-(13/8)*(x-h).^2))-g20 (3/h)*(((1/24)*(x+2h).^6-(13/8)*(x+2h).^2)-((1/24)*(x).^6-(13/8)*(x).^2)) (3/h)*(((1/24)*(x+3h).^6-(13/8)*(x+3h).^2)-((1/24)*(x+h).^6-(13/8)*(x+h).^2)) (3/h)*(((1/24)*(x+4h).^6-(13/8)*(x+4h).^2)-((1/24)*(x+2h).^6-(13/8)*(x+2h).^2)) (3/h)*(((1/24)*(x+5h).^6-(13/8)*(x+5h).^2)-((1/24)*(x+3h).^6-(13/8)*(x+3h).^2))-g25]
%f1=[(3/h*0.0738-16.7500) 3/h*0.0882 3/h*0.1033 3/h*0.1190 (3/h*0.1355-21.0601)]
%f1=[-5.6800 13.2300 15.4950 17.8500 -0.7351]
f1=[1.7958 2.1209 2.0790 2.0380 1.6716]
m2=G1^(-1)
m3=m2*f1'


%m2=(1/(h^2))*(((1/24)*(x+h).^6-(13/8)*(x+h).^2)-2*((1/24)*x.^6-(13/8)*x.^2)+((1/24)*(x-h).^6-(13/8)*(x-h).^2))
%det1=abs(g1-m3)
j=length(x);
fprintf('------------------------------------------------------------\n')
fprintf('  x               g1               m1                det1\n');
fprintf('------------------------------------------------------------\n')
for ii=1:5,
    det1(ii)=abs(g1(ii+1)-m3(ii));
    fprintf('%-10.4f       %-10.4f      %-10.4f       %-10.4e\n',x(ii),g1(ii+1),m3(ii),det1(ii));
end
%plot(x,y1,x,m2),grid
%   2.1637    2.1209    2.0790    2.0380

 % Columns 5 through 7 

 %   1.9980    1.9588    1.9205