曲柄滑块运行程序通常包括以下步骤：
输入已知数据
清除之前的计算结果。
定义已知参数，如连杆长度（l1, l2）、偏心距（e）、转动惯量（hd）、角速度（omega1）、角加速度（alpha1）等。
主程序
使用for循环遍历每一时刻（n1从1到721），计算每一时刻的曲柄角度（theta1）。
调用子程序`slider_crank`，输入当前时刻的曲柄角度、角速度、角加速度以及连杆长度，计算并返回曲柄滑块机构的位移（s3）、角速度（omega2）、线速度（v3）、角加速度（alpha2）和线加速度（a3）。
输出计算结果，并绘制曲线。
子程序
`slider_crank`函数接收曲柄角度（theta1）、角速度（omega1）、角加速度（alpha1）、连杆长度（l1, l2）和偏心距（e）作为输入参数。
根据曲柄滑块机构的运动学方程，计算并返回所需的位移、速度、加速度等参数。
```matlab
% 输入已知数据
clear；
l1 = 0.1；
l2 = 0.3；
e = 0；
hd = pi/180；
du = 180/pi；
omega1 = 10；
alpha1 = 0；
% 主程序
for n1 = 1:721
theta1（n1） = （n1-1）*hd； % 计算每一时刻的曲柄角度
[theta2（n1）, s3（n1）, omega2（n1）, v3（n1）, alpha2（n1）, a3（n1）] = slider_crank（theta1（n1）, omega1, alpha1, l1, l2, e）；
end
% 绘制结果
plot（theta1, s3, theta1, v3, theta1, a3）；
xlabel（'Time （degrees）'）；
ylabel（'Displacement （m）'）；
title（'slider_crank Motion'）；
grid on；
% 子程序
function [theta2, s3, omega2, v3, alpha2, a3] = slider_crank（theta1, omega1, alpha1, l1, l2, e）
% 根据曲柄滑块机构的运动学方程进行计算
% 这里省略具体计算过程，实际代码应根据运动学方程进行实现
end
```
请注意，上述代码仅为示例，实际应用中需要根据具体的物理模型和计算需求进行详细实现。