h1(2) := A;
vars := {h0(0),h0(1),h1(0),h1(1),h1(3)};


read(`gb.lp.fact.2`);
eqs := ";
ss := solve({eqs[1..3],eqs[5..6]},vars);
sA := solve(eqs[4]);
assign(ss);

# ---------- print out filter coefficients ----------

writeto(`coeff2.m`);

lprint(`function [h0,h1] = coeff2(k)`);
lprint(`if k == 0`);
lprint(`   A = `,sA[1],`;`);
lprint(`elseif k == 1`);
lprint(`   A = `,sA[2],`;`);
lprint(`end`);

lprint(`h0 = [`);
        lprint(h0(0));
        lprint(h0(1));
        lprint(h0(0));
lprint(`]/sqrt(2);`);
lprint(`h0 = [h0; 0; 0; 0; 0; 0];`);

lprint(`h1 = [`);
        lprint(h1(0));
        lprint(h1(1));
        lprint(h1(2));
        lprint(h1(3));
        lprint(h1(2));
        lprint(h1(1));
        lprint(h1(0));
lprint(`]/sqrt(2);`);
lprint(`h1 = [h1; 0];`);

writeto(`terminal`);