将上题代码汇合在一起发给我

当然，下面是将你之前提供的两段代码结合在一起的版本：

f=@(x) 0.4*x(2)+x(1)^2+x(2)^2-x(1)*x(2)+1/30*x(1)^3;
x0=[0.5;0.5];
A=[-1,-0.5;-0.5,-1];
b=[-0.4;-0.5];
lb=[0;0];
option=optimset('Display','off');
[xmin,fmin]=fmincon(f,x0,A,b,[],[],lb,[],[],option);

[t,x]=ode45(@fun,[0,20],[0;1]);
subplot(2,2,1);plot(t,x(:,1));
subplot(2,2,2);plot(x(:,1),x(:,2));

这段代码首先定义了一个函数 f，然后设置了最小化问题的初始点 x0，线性不等式约束 A 和 b，以及下界限制 lb。接着使用 fmincon 函数求解最小化问题，并将最优解和最优值存储在 xmin 和 fmin 中。

然后，使用 ode45 函数求解常微分方程，并将结果存储在变量 t 和 x 中。最后，使用 subplot 和 plot 函数绘制两个子图，分别显示时间 t 和 y(1) 的关系，以及 y(1) 和 y(2) 的关系。

这段代码将求解最小化问题和求解常微分方程的过程整合在一起，并绘制了相关的曲线图。