% == Input Data ==== y = [0 1 2 2 1 0]; % == Algorithm ==== x = 1:length(y); s=length(y)-1; % степень интерполяционного полинома k = polyfit(x,y,s); % коэффициенты интерполяционного полинома Y = @(t) polyval(k, t); % анонимная функция от переменной t (полином с коэффициентами k) dY = @(t) polyval(polyder(k), t); % производная integrand_L = @(t) sqrt(1 + dY(t).^2); % подынтегральная функция для расчета длины кривой L = integral(integrand_L, min(x), max(x)) % длина кривой S = integral(Y, min(x), max(x)) % площадь под графиком функции