diff2m1.m

Micro Electronic Device Fundermental

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=2:h:2.1
y1=(1/24)*x.^6-(13/8)*x.^2
g2=(5/4)*x.^4-(13/4)

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=[194.0194 220.8087 231.1476 241.7920 230.7485]
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
%0.0601    0.0738    0.0882    0.1033

%  Columns 5 through 6 

%    0.1190    0.1355 %-10.4e%det3(ii)
    
%       16.7500    21.0601
%  Columns 1 through 4 

%   -5.6800   13.2300   15.4950   17.8500

%  Column 5 

%   -0.7351

%   0.0738    0.0882    0.1033    0.1190

%  Column 5 

%    0.1355