[求助] 将一个微分方程组的解应用到另一个方程

function dy=eq3(t,y)
dy=zeros(2,1);
dy(1)=[1.6+0.2*(1-y(1))^2+0.8*sin(200*pi*t)]/[1.6+0.4*(1-y(1))^2];
dy(2)=0.2*y(2)*[0.8+0.8*sin(200*pi*t)]*(y(1)-1)/[0.8+0.2*(1-y(1))^2]^2;

function f=fx(x,y);
f=(0.8+0.2*x.^2)^1.5*x.^2/y;
p=log(0.8+0.2*x.^2)-0.2857*log(0.0351+0.01348*quad(@(x)fx(x,y),x,1));%第二个方程，里面含有积分，积分范围【x 1】
H0=1;
[t,y]=ode45('eq3',[0 1],[0 H0]);
plot(x,p)