矩阵元素的积分和矩阵图

Question

1。我用ODE求解器求解了非线性微分方程     需要与绘制答案的曲线拟合。我如何获得协调     EXCEL格式文件中的曲线上的点？  2.如何绘制积分     计算得出的答案？  3.在我的代码图中∫y（：，1）。我的代码是：

% landa_D=(epsilon_0*K*Te/n0+*e^2)^.5;
% kesi=x/landa_d;
%d/dt=d/d(kesi);
% eta=-e*phi/K*Te;
% u0=v0+/Cs0;
% Cs=Cs0/(1-alpha*(1-gama))^.5;
%syms z1 u eta gama kesi alpha
%constant
%y(1)=d(eta)/d(kesi);
%y(2)=d^2(eta)/d^2(kesi);
%y(3)=du/d(kesi);
function dydt = odefcn_1(eta,y,A,B,y0)
tspan = (0:.1: 15);
y0 = [5 1 2];
[eta,y] = ode23s(@ F , tspan, y0);
%z_1=quad(y(:,1));
plot(eta,y(:,1),'-o',eta,y(:,2),'-.',eta,y(:,3),'b');
end
function dydt_1=F(eta,y)
alpha=.08;
gama=80;
u0=2;
dydt_1 = zeros(3,1);
dydt_1(1) = y(1);
dydt_1(2) = u0/y(2)-((1-alpha)*exp(-eta)+alpha*exp(-gama*eta));
dydt_1(3)=y(1)/y(2);
end