用向前欧拉法（古典显格式）求解一维非齐次热传导方程的matlab代码，非齐次项为x*exp(t)-6*x，初值条件为x.^3+x，左边界条件为0，右边界条件为1+exp（t）

下面是使用向前欧拉法求解一维非齐次热传导方程的 MATLAB 代码：

% 设置参数
N = 100;        % 空间离散化点数
T = 1;          % 时间总长度
M = 1000;       % 时间离散化点数
dx = 1/N;       % 空间步长
dt = T/M;       % 时间步长

% 初始化向量
x = linspace(0, 1, N+1);    % 空间网格
t = linspace(0, T, M+1);    % 时间网格
u = zeros(N+1, M+1);        % 数值解矩阵

% 设置初值条件
u(:, 1) = x.^3 + x;

% 设置边界条件
u(1, :) = 0;                    % 左边界条件
u(N+1, :) = 1 + exp(t);         % 右边界条件

% 使用向前欧拉法求解
for j = 1:M
    for i = 2:N
        u(i, j+1) = u(i, j) + dt/dx^2 * (u(i-1, j) - 2*u(i, j) + u(i+1, j)) + dt * (x(i)*exp(t(j)) - 6*x(i));
    end
end

% 绘制结果
[X, T] = meshgrid(x, t);
surf(X, T, u');
xlabel('x');
ylabel('t');
zlabel('u');
title('Solution of the Heat Equation');

请注意，这只是一个简单的示例代码，可能需要根据具体情况进行修改和调整。