在MATLAB中编程解决边界层问题通常涉及以下步骤：
选择合适的数学模型：
根据具体的边界层问题选择合适的数学模型，例如Blasius方程、Prandtl的混合长度理论等。
设置边界条件：
根据问题的具体要求设置边界条件，如速度、温度、压力等。
数值求解微分方程：
使用MATLAB提供的函数（如`ode45`、`fzero`、`fsolve`等）进行数值求解微分方程。这可能涉及到将高阶微分方程降阶、使用换元法简化问题等技巧。
处理边界条件：
在求解过程中，需要特别关注边界条件的应用，确保解的正确性。
结果可视化：
使用MATLAB的绘图功能（如`plot`、`surf`等）将求解结果可视化，帮助分析边界层的特性。
验证与调整：
通过与实验数据或其他数值方法的结果进行比较，验证所求解的正确性，并根据需要进行调整。
下面是一个简单的示例，展示如何使用MATLAB求解平板边界层问题：
```matlab
% 定义微分方程
function dy = rigid（t, y）
dy = zeros（3,1）；
dy（1） = y（2）；
dy（2） = y（3）；
dy（3） = -0.5*y（1）*y（3）；
end
% 设置求解的时间区间和初始条件
tspan = [0 5]；
y0 = [0 0 0]；
% 使用ode45函数求解微分方程
[T, Y] = ode45（@rigid, tspan, y0）；
% 绘制结果
plot（T, Y（:,1）,'-',T,Y（:,2）,'*',T,Y（:,3）,'+'）
```
在这个示例中，我们定义了一个简单的刚体运动方程作为微分方程，并使用`ode45`函数进行求解。通过绘制解随时间变化的曲线，可以观察到边界层的行为。
对于更复杂的边界层问题，可能需要结合使用多种数值方法和技巧，以及更高级的边界层理论。此外，还可以利用MATLAB的仿真工具箱（如Control System Toolbox、Fluid Flow Toolbox等）来辅助分析和可视化。