将这段代码改写为matlab语言

将这段Python代码转换成MATLAB语言，我们需要注意语法差异以及MATLAB特定的库和功能。以下是代码的MATLAB版本：

% 导入所需工具箱
import numpy as np;
import math;
syms x;
que_f = 1 ./ (1 + 25 * x^2);

% MATLAB中无需定义符号变量，直接使用
f = sym('1 ./ (1 + 25 * x^2)');

% 定义拉格朗日插值函数
function L_N = langrange(N, x_in)
    L_N = 0;
    for i = 1:N
        l_i = 1;
        for j = 1:N
            if j ~= i
                l_i = l_i * (x - x_in(j)) ./ (x_in(i) - x_in(j));
            end
        end
        L_N = L_N + evalin(symphonyEngine,'subs',que_f,x,x_in(i)) .* l_i;
    end
end

% 定义插值点数组和插值次数列表
N_num = [5, 8, 10, 15];
% 设置字体和显示格式
set(0, 'DefaultTextFontName', 'Times New Roman');
set(0, 'DefaultAxesFontName', 'Times New Roman');
set(groot, 'DefaultFigurePaperType', 'a4');

% 变量初始化
xx = linspace(-1, 1, 1000);
LG_E_N = zeros(size(N_num));
LGR_E_N = zeros(size(N_num));
LGL_E_N = zeros(size(N_num));

% 循环计算并存储各方法误差
for k = 1:length(N_num)
    N = N_num(k);
    x_h = [-1 + (0:N)/N; 1]; % 注意MATLAB中需要指定端点
    L = langrange(N+1, x_h);
    
    % 计算误差并更新列表
    for j = 1:numel(xx)
        LG_E_N(k) = max(LG_E_N(k), abs(que_fsubs(x, xx(j)) - subs(L, x, xx(j))));
        LGR_E_N(k) = max(LGR_E_N(k), abs(que_fsubs(x, xx(j)) - subs(L, x, xx(j))));
        LGL_E_N(k) = max(LGL_E_N(k), abs(que_fsubs(x, xx(j)) - subs(L, x, xx(j))));
    end
end

% 创建误差曲线图
figure;
plot(N_num, log10([LG_E_N(1:end-1); 1e-14]), 'DisplayName', 'LG');
hold on;
plot(N_num, log10([LGL_E_N(1:end-1); 1e-14]), 'DisplayName', 'LGL');
plot(N_num, log10([LGR_E_N(1:end-1); 1e-14]), 'DisplayName', 'LGR');
hold off;
xlabel('N');
ylabel('log10(误差)');
title('三种插值节点的误差曲线');
legend;

% 更改字体设置
set(gca, 'FontName', 'Times New Roman');

这个MATLAB代码实现了与Python代码类似的功能，包括拉格朗日插值、其他插值方法和误差评估。请注意MATLAB的语法和函数调用方式与Python有所不同。