帮我用MWORKS的语言编写以下程序： 将波浪频散方程，采用迭代法进行求解，将其求解波浪频散方程编写为一个函数dispersion，并计算T＝2：10，水深h＝20：-1：1m的波长、波速。

在MWorks（MathWorks Wave Optics Simulation Environment）环境中，你可以使用MATLAB语言编写这个程序，因为MWorks通常用于光波模拟，但它支持MATLAB脚本。以下是一个简化的示例，假设我们使用了数值积分方法如欧拉方法进行迭代求解：

function [lambda, velocity] = dispersion(Ts, depths)
% 定义波浪频散方程的函数
dispersion_func = @(k) -1i * k * sqrt(1 + (kh)^2); % 例子，这是一个简单的线性频散方程，kh代表深度/波数

% 波数范围
kmin = 2*pi/Ts; % 波长最小值对应的波数
kmax = 2*pi/(Ts*1e-3); % 波长最大值对应的波数，这里假设T=10ms

% 水深范围
kh_range = depths ./ linspace(Ts, Ts*1e-3, length(depths)); % kh = depth / wave number

% 初始化结果矩阵
lambda = zeros(size(depths));
velocity = real(dispersion_func(kmax));

for i = 1:length(depths)
    % 迭代求解每个深度下的波长
    % 这里假设使用欧拉方法，实际应用中可能需要选择更精确的数值方法
    kh = kh_range(i);
    lambda(i) = 2*pi ./ kh; % 波长 = 2pi/k
    % 进行一次迭代更新波速（简化）
    velocity(i) = velocity * exp(-dispersion_func(kmin) * (kh_range(i) - kh));
end
end

% 调用函数并计算指定范围的波长和波速
Ts_values = 2:10; % 时间间隔范围
depths_values = 20:-1:1; % 深度范围
[lambda_results, velocity_results] = dispersion(Ts_values, depths_values);

% 输出结果
disp("波长:")
disp(lambda_results)
disp("波速:")
disp(velocity_results)

%