编写matlab代码实现：当k取值范围为0到pi时，Pin在0到20取值，Pout=Pin*cos(Pin)，找到Pout的最大值与最小值的差值定义为T，画出x轴为k、y轴为T的函数曲线。

在MATLAB中，你可以使用`for`循环、`cos`函数以及`max`和`min`函数来实现这个需求。首先，我们需要创建两个向量，一个是角度`k`，另一个是`Pout`。然后计算每个`Pin`对应的`Pout`值，并找出最大值和最小值之差。最后，我们可以用`plot`函数绘制函数曲线。下面是相应的MATLAB代码：

% 定义变量范围
Pin = 0:20; % Pin 的取值范围
k = linspace(0, pi, 180); % k 的取值从0到pi，180个等分点

% 初始化最大值和最小值
max_diff = NaN;
min_diff = Inf;

% 计算并存储Pout和差异
for i = 1:length(Pin)
    Pout(i) = Pin(i)*cos(Pin(i)); % 计算Pout
    diff = max(Pout(1:i)) - min(Pout(1:i)); % 计算差值
    if diff > max_diff
        max_diff = diff;
    elseif diff < min_diff
        min_diff = diff;
    end
end

% 找到最大差值对应的k值
[k_max_diff, idx] = max(diff);
[k_min_diff, ~] = min(diff);

% 绘制曲线
figure;
plot(k, diff, 'LineWidth', 2);
hold on;
plot([k_min_diff k_max_diff], [min_diff max_diff], 'r--', 'LineWidth', 1.5);
xlabel('k (弧度)');
ylabel('T (Pout的最大值减去最小值)');
title(['Max difference: T = ', num2str(max_diff), ', at k = ', num2str(k_max_diff)]);
legend('Function curve', ['Max difference (k = ', num2str(k_max_diff), ')']);
grid on;

运行这段代码后，你会得到一条函数曲线，x轴表示`k`，y轴表示`T`，并且会显示最大差值的位置。