如何用Matlab修改或者写泊肃叶方程的代码

Question

式1中的方程对应的代码如下：
global M g p R r

M=0.5;%千克

g=9.8;

p=1000;%水的密度

R=0.03;

r=0.005;

t1=[0,75];

H1=[0;0;0];

[t,H] = ode45(@dH,t1,H1);

plot(t,H);

function dHdt=dH(t,H)

global M g p R r

dHdt=zeros(3,1);

dHdt(1)=H(2);

dHdt(2)=M*g-p*g*pi*(H(1)-H(3))*(R^2-r^2);

dHdt(3)=sqrt(-4*g*r^4/R^4*(H(3)-H(1)));

end

现在把1式修改为3式 对应的代码应该怎么改？
添加的常量值为0.004

图2中的两个方程应该怎么写？