func6p4.m

来自「数值分析最常用的四十种算法」· M 代码 · 共 7 行

[t,y]=ode23('sub6p4',[0 exp(1)],1);   
f=linspace(1,2,14)
for i=1:14
    f(i)=(1-log(y(i)))*y(i)
end
plot(t,y,'-r',y,f,'.-')
%Please also the figure,v(e)=6.8798