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)
g2=-2./(x.^3).*exp(-1./x)+1./x.^4.*exp(-1./x)
G2=[10 1 0 0 0;1 10 1 0 0;0 1 10 1 0;0 0 1 10 1;0 0 0 1 10]
%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]
f2=[-4.0465 -4.3254 -4.2352  -4.1444 -3.7307]
m2=G2^(-1)
m3=m2*f2'


%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(g2(ii+1)-m3(ii));
    fprintf('%-10.4f       %-10.4f      %-10.4f       %-10.4e\n',x(ii),g2(ii+1),m3(ii),det1(ii));
end
%plot(x,y1,x,m2),grid
 %   -4.4144   -4.3254   -4.2352   -4.1444

%  Columns 5 through 7 

 %  -4.0534   -3.9626