Недавно начал изучать эту штуку. Мне нужно вычислить коэффициенты ряда Фурье. Пошарил под доке нашел функцию fourier, которая вроде как должна возвращать список коэффициенты.
Вот функция которую нужно разложить в ряд Фурье
T:40; A:0.32; flsq(x) := (v:sin(2*%pi*x/T), if (v >= 0) then 0 else A);
Тут пытаюсь вычислить foursimp(fourier(flsq(x), x, T/2)); Однако на выходе он только выполняет подстановку функции flsq, но не вычисляет конкретных значений.
Ещё пытался сделать вручну, т.е. так: fourcosn(n) := ((2/T) * quad_qags(flsq(x) * cos(n*2*%pi*x/T), x, 0, T))[1];
Однако тут для n != 1, выдаётся сообщение об ошибке 2. На сколько я понял это связано с потерей точности.
Заранее спасибо
