用matlab编写程序求e-x^2（1-2x)的极值，不能使用函数

在MATLAB中，如果你想要寻找表达式 \( e^{-x^2(1-2x)} \) 的极值而不能直接使用内置函数如`max`、`min`或`fminbnd`等，你可以使用数值优化工具箱中的迭代方法，例如梯度下降法或牛顿法。这里我们可以尝试使用二阶导数测试来找到可能的极值点，然后手动计算函数值。

首先，你需要手动计算该函数的一阶导数和二阶导数：

一阶导数 \( f'(x) = -2xe^{-(1-2x)x^2} + 4x^2e^{-(1-2x)x^2} \)

二阶导数 \( f''(x) = 4(1 - 6x + 8x^3)e^{-(1-2x)x^2} - 4x^3e^{-(1-2x)x^2} \)

然后设置初始猜测点，比如 `x0 = 0`，并使用以下步骤：

1. 计算导数在给定点的值。
2. 如果一阶导数等于0，检查二阶导数判断是否是极大值还是极小值（如果二阶导数大于0，则是局部极小；小于0则可能是极大值）。
3. 如果一阶导数不为0，选择另一个方向调整 x 的值，直到一阶导数变为0，然后再按照上述步骤处理。

这个过程需要手动循环，并且可能会比较耗时，特别是对于复杂函数。下面是一个简单的示例代码，但这只是一个基础版的实现，实际应用可能需要更精细的算法控制：

syms x
% 定义函数及其导数
f = exp(-x.^2*(1 - 2*x));
df = diff(f, x);
ddf = diff(df, x);

% 初始猜测
x0 = 0;
tol = 1e-6; % 精度设定
dx = 1e-3; % 初始步长

while abs(dx) > tol
    df_x0 = df.subs(x, x0);
    
    if df_x0 == 0
        ddf_x0 = ddf.subs(x, x0);
        
        if ddf_x0 > 0
            % 极小值，向左移动
            x0 = x0 - dx;
        else
            % 极大值，向右移动
            x0 = x0 + dx;
        end
        
        dx = sign(ddf_x0)*dx/2; % 调整步长
    else
        % 导数不为零，继续搜索
        dx = -df_x0 / df.subs(x, x0);
    end
    
    fprintf('x = %f, f(x) = %f\n', x0, f.subs(x, x0));
end

% 输出结果
fprintf('\nPossible extremum at x = %.6f with value f(x) = %.6f.\n', x0, f.subs(x, x0));

%