N := 8;

	# n := 6;  # n = 1..6

        fn := `gb.lp.1`;
        fn2 := `sol1_`.n.`.m`;
        lprint(fn);
        lprint(fn2);

	h(7) := B;

        read(fn):
        eqs := ":

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

	ss := solve({eqs[3..17]},vars):
	Digits := 40;
        sB := fsolve(eqs[2],B);
	B  := sB[n];
	h(7) := B;

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

        lprint(`function [h,g] = sol1_`.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);`);


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