参数微分方程组--绘制3d曲面图
 

现有两参数D和P，微分方程组有6个，自变量时间
D范围0到1
P范围1到100

求问如何绘制出y（3）关于D和P的曲面图
即求D和P的值 使y（3）在t趋近于无穷时 最小

请问这种情况如何用ode45解？

方程如下：
function [dy] = ode(t, y)

k1 = 31;
k2 = 0.02;
k3 = 0.0001875;
k4 = 0.45;
k5 = D*k4;
k6 = 200;
k7 = 30;

dy = zeros (6,1);
dy(1) = k1 - k2*y(1) - k3*y(3)*y(1) - k3*y(6)*y(1);
dy(2) = k3*y(3)*y(1) - k4*y(2);
dy(3) = k6*k4*y(2) - k7*y(3) + D*k6*k5*y(5);
dy(4) = k3*y(6)*y(1) - k2*y(4) - k3*y(3)*y(4);
dy(5) = k3*y(3)*y(4) - k5*y(5);
dy(6) = P*P*D*k6*k5*y(5) - k7*y(6);

end