N := 10;

	# n := 1;  # n = 1..8

        fn := `gb.lp.2`;
        fn2 := `sol2_`.n.`.m`;
        lprint(fn);
        lprint(fn2);

        read(fn):
        eqs := ":

  	vars := { h(0), h(1), h(2), h(3), h(4), h(5), h(6), h(7), h(8), h(9),
		g(0), g(1), g(2), g(3), g(4), g(5), g(6), g(7), g(8), g(9) };

	ss := solve({eqs[2..21]},vars):
	Digits := 40;
        sA := fsolve(eqs[1],A);
	A  := sA[n];

	interface(screenwidth=600);
	interface(quiet=true);
	writeto(fn2);

        lprint(`function [h,g] = sol2_`.n);
        
        lprint(`h = [`);
        k := 'k':
        for k from 0 to N-1 do
                lprint(subs(ss,h(k)));
        od;
        lprint(`]'*sqrt(2);`);

        lprint(`g = [`);
        k := 'k':
        for k from 0 to N-1 do
                lprint(subs(ss,g(k)));
        od;
        lprint(`]'*sqrt(2);`);


        writeto(`terminal`);
	interface(quiet=false);
	interface(screenwidth=80);