diff1c2.m

来自「Micro Electronic Device Fundermental」· M 代码 · 共 19 行

format short
%syms x
%plot(x,y1),grid
%plot(x,y1,x,g1),grid
%h=((3*0.5e-4)/45)^(1/3)
h=0.02
x=1:h:1.12
y1=exp(-1./x)
g1=(1./x.^2).*exp(-1./x)
m1=(1/(2*h))*(exp(-1./(x+h))-exp(-1./(x-h)))
det1=abs(g1-m1)
j=length(x);
fprintf('------------------------------------------------------------\n')
fprintf('  x               g1               m1                det1\n');
fprintf('------------------------------------------------------------\n')
for ii=1:j,
    fprintf('%-10.4f       %-10.4f      %-10.4f      %-10.4e\n',x(ii),g1(ii),m1(ii),det1(ii));
end
plot(x,y1),grid