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(-x.^2).*cos(20*x)
g2=-402*exp(-x.^2).*cos(20.*x)+4*x.^2.*exp(-x.^2).*cos(20*x)+80*x.*exp(-x.^2).*sin(20*x)
G1=[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]
f1=[335.6814 897.9 1305.3 1475.3 1306.456]
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(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
%   1.0e+003 *

 % Columns 1 through 4 

  %  0.3028    0.8979    1.3053    1.4753

%  Column 5 

 %   1.4002