func6p4.m

数值分析最常用的四十种算法

[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