Здравствуйте, помогите перевести функцию с матлаба на си... Код (Text): function f=n1yink(u,n,dt,m) % f=n1yink(u,n,dt,m) - Вычисление определенного интеграла % по квадратурным формулам Ньютона-Котеса % u - массив значений интегрируемой функции % n - число точек в массиве % dt - шаг интегрирования % m - число точек на интервале (порядок точности+1) c=[ 1 0 0 0 2 1 4 0 0 6 1 3 0 0 8 7 32 12 0 90 19 75 50 0 288 41 216 27 272 840 751 3577 1323 2989 17280]; f=0;f1=0; if any([n<=1 m<1]) break, end % невозможно интегрировать m=min([m n 8]); m1=m-1; % порядок точности if m1==0 % специальный случай - формула прямоугольников f=sum(u(1:n-1))*dt; break end % Коэффициенты квадратурной формулы d=[c(m1,1:round(m/2)) c(m1,fix(m/2):-1:1)]; % строка коэф-тов dd=c(m1,5); % Вычисление интеграла s=zeros(m,1); % столбец значений частичных сумм for i=1:m ix=(0:m1:n-m)+i; s(i)=sum(u(ix)); end f=m1*dt/dd*(d*s); n1=n-max(ix); % число неиспользованных точек if n1==0 break, end % нет неиспользованных точек f1=n1yink(u(n-n1:n),n1+1,dt,n1+1); f=f+f1; Заранее очень благодарен!
