Фильтр Баттерворта

void filter(std::vector<double>&input, std::vector<double>&output,std::string way_for_filter_koeff) {
    static std::vector<double> a, b;
    static bool first_loading = 1;
    int order = 10;
    
    if (first_loading) {
        
        std::ifstream input(way_for_filter_koeff);
        double koef;
        for (size_t i = 0; i <= order; i++) {
            input >> koef;
            b.push_back(koef);
        }
 
        for (size_t i = 0; i <= order; i++) {
            input >> koef;
            a.push_back(koef);
        }
        first_loading = 0;
    }
    
    order = a.size();
    static std::vector<double>z;
    
    for (int i  = 0;i<order;i++)
        z.push_back(0);      // Creates zeros if input z is omitted
    if (a.at(0) != 1) {
        for (int i = 0; i < order; i++) {
            b.at(i) /= a.at(0);  
            a.at(i) /= a.at(0);
        }
    }
    
    double new_sample = 0;
        new_sample = b.at(0) * *(input.end()-1) + z.at(0);
        
        for (int i = 1; i < order; i++) {
            z.at(i - 1) = b.at(i) * *(input.end() - 1) + z.at(i) - a.at(i) * new_sample;
            
        }
        z.erase(z.end() - 1);
 
        std::ofstream out("filtered.txt", std::ios::app);
        out << new_sample << std::endl;
        output.push_back(new_sample);
 
};