使用matlab编写程序：有一半径为4a的大圆，里面有一个半径为a的小圆，现在，大圆固定而小圆在大圆内相切滚动，设起点M的坐标为(4a.0)，确定小圆上一点M的轨迹曲线，并动态模拟内切小圆的滚动轨迹。

好的，这是一个很有趣的问题。我会尽力回答你的问题并提供一个MATLAB程序来模拟内切小圆的滚动轨迹。

首先，让我们考虑小圆在大圆内滚动的情况。我们可以将小圆的圆心坐标表示为(x,y)，其中x表示小圆圆心的横坐标，y表示小圆圆心的纵坐标。我们可以用下面的公式来计算小圆圆心的运动轨迹：

x = 4a + a*cos(theta)
y = a*sin(theta)

其中theta是小圆圆心相对于起点M的运动角度。我们可以通过改变theta的值来模拟小圆的运动轨迹。

接下来，我们需要绘制大圆和小圆的图形，并模拟小圆在大圆内的滚动。下面是一个MATLAB程序来实现这个任务：

% 定义大圆半径和小圆半径
R = 4;
r = 1;

% 定义起点M的坐标
x0 = R + r;
y0 = 0;

% 定义theta的取值范围
theta = linspace(0, 2*pi, 200);

% 循环绘制小圆的运动轨迹
for i = 1:length(theta)
    % 计算小圆圆心坐标
    x = x0 + r*cos(theta(i));
    y = y0 + r*sin(theta(i));
    
    % 绘制大圆
    t = linspace(0, 2*pi, 100);
    plot(R*cos(t), R*sin(t), 'k');
    hold on;
    
    % 绘制小圆
    plot(x + r*cos(t), y + r*sin(t), 'r');
    
    % 绘制小圆圆心
    plot(x, y, 'ro');
    
    % 设置坐标轴范围
    axis([-5 9 -5 5]);
    axis equal;
    
    % 暂停一段时间，以便观察运动轨迹
    pause(0.05);
    
    % 清除图形
    clf;
end

在这个程序中，我们使用了一个for循环来计算小圆圆心的位置，并绘制大圆、小圆和小圆圆心的图形。我们还使用了pause函数来暂停程序的执行，以便观察运动轨迹。最后，我们使用clf函数清除图形，以便绘制下一帧。

使用上述MATLAB程序可以模拟内切小圆的滚动轨迹。你可以在MATLAB中运行代码并观察运动轨迹。