[watcher]

如果按
rk4('dadiab',7,theta,dtheta,y);
这种格式调用dadiab，那么dadiab好像应该是函数，即以function开头

lz的做了很多工作，但是遗憾得告诉你，matlab中ode系列函数可以完善的解决你的问题。也就是说，除非你是故意想锻炼一下你的matlab编程技巧，完全没必要自己编写runge-kutta求解函数。
建议先看一下matlab中关于ode求解的帮助文档（通常我会从函数ode45开始）。