untitled3.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=2:h:2.12
y1=(1/24)*x.^6-(13/8)*x.^2
g1=1/4*x.^5-13/4.*x
m1=(1/(2*h))*(((1/24)*(x+h).^6-(13/8)*(x+h).^2)-((1/24)*(x-h).^6-(13/8)*(x-h).^2))
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,x,m1),grid