N := 8;

	# n := 4;  # n = 1..4

        fn := `gb.lp.4`;
        fn2 := `sol4_`.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),
		g(0), g(1), g(2), g(3), g(4), g(5), g(6), g(7) };

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

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

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