请用matlab写出改程序

以下是用MATLAB编写的求解非线性方程的程序：

function [f] = solve_nonlinear_equation(c, N, epsilon)
    h = 1/N;
    
    % 初始猜测解
    f = ones(N, 1);
    
    % 牛顿迭代
    while true
        F = zeros(N, 1);
        J = zeros(N, N);
        
        % 计算非线性方程和雅可比矩阵
        for i = 2:N-1
            F(i) = c^2 * (f(i+1) - 4*f(i) + 6*f(i-1) - 4*f(i-2) + f(i-2))/h^4 - (f(i+1) - 2*f(i) + f(i-1))/h^2 - exp(f(i));
            J(i, i-2:i+1) = [1/h^4, -4/h^4 + 2/h^2, 6/h^4 - 2/h^2 - exp(f(i)), -4/h^4];
        end
        
        % 边界条件
        F(1) = f(1) - 10;
        F(N) = f(N) - 20;
        J(1, 1) = 1;
        J(N, N) = 1;
        
        % 牛顿迭代更新
        delta_f = J \ (-F);
        f = f + delta_f;
        
        % 判断是否收敛
        if norm(delta_f) < epsilon
            break;
        end
    end
end

你可以调用该函数来求解非线性方程。例如，通过以下方式调用函数来求解c=1，N=100和epsilon=1e-6的情况：

c = 1;
N = 100;
epsilon = 1e-6;
f = solve_nonlinear_equation(c, N, epsilon);

请注意，这只是一个示例程序，你可以根据自己的需求进行修改和调整。