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

在MATLAB中，你可以使用循环结构和`max`和`min`函数来计算`Pout`数组的最大值和最小值差值。首先，我们需要创建一个二维数组`Pout`，然后遍历`k`和`Pin`的组合，绘制函数图形。以下是相应的代码：

% 定义变量范围
k = linspace(0, pi, 100); % k从0到pi均匀分布100个点
Pin = 0:20; % Pin从0到20步长为1

% 初始化结果矩阵
Pout = zeros(size(k));

% 计算并填充Pout数组
for i = 1:length(k)
    for j = 1:length(Pin)
        Pout(i,j) = Pin(j) * cos(Pin(j) * k(i));
    end
end

% 找到最大值和最小值
[T, idx_max] = max(Pout(:)); % 最大值及其索引
[T_min, idx_min] = min(Pout(:)); % 最小值及其索引

% 计算最大值和最小值的差值
T_diff = T - T_min;

% 绘制x轴为k，y轴为T的函数曲线
plot(k, Pout, 'LineWidth', 2);
hold on;
grid on;
xlabel('k');
ylabel('T (最大值和最小值差)');
title(sprintf('Pout = Pin * cos(Pin * k), T = %f', T_diff));

% 添加最大值和最小值点
plot([k(idx_max)], [T], 'ro', 'MarkerSize', 8, 'LineWidth', 2, 'DisplayName', 'Max');
plot([k(idx_min)], [T_min], 'bo', 'MarkerSize', 8, 'LineWidth', 2, 'DisplayName', 'Min');

% 显示图例
legend('Function', 'Max', 'Min', 'Location', 'Best');

% 清除之前的结果
hold off;

运行这段代码后，你应该会看到一个描绘了`Pout`随`k`变化的曲线，以及函数的最大值和最小值点。同时，`T_diff`将存储这两个值之间的差。