简答3. 阅读以下函数文件代码 function project f=@(t,y)[y(2); -1-1/(y(1)+10)^2+y(2)^2]; %微
简答3. 阅读以下函数文件代码 function project f=@(t,y)[y(2); -1-1/(y(1)+10)^2+y(2)^2]; %微分方程 opts=odeset(’events’,@g); %定义events [t,y,tfinal] = ode45(f, [0,Inf],[1; 0] , opts); tfinal plot(t,y(:,1),’-’) %作图 axis([-0.1 tfinal+0.1 -2.1 1.1]) xlabel(’t’) ylabel(’y’) title(’Falling body’) %作图结束 function[gstop,isterminal,direction]=g(t,y) %events功能对应函数 gstop=y(1)+2; isterminal=1; direction=-1; end end 简述(1)f描述了一个怎样的微分方程?初值如何?(2)说明该代码中events功能的作用机制(解释其对应的函数)。(3)说明tfinal的意义。(4)大致描述作图区的内容 本题拍照上传。