运动学控制方程（17-1）为二阶常微分方程，可通过ode45函数求解。编制自定义函数function k=DD_fun（～,y,～,m1 ,m2,g,11,12,Jc），并存为文件DD_fun.m，以方便GUI主程序调用。代码如下：
