N := 6;

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

        fn := `gb.lp.3`;
        fn2 := `sol3_`.n.`.m`;
        lprint(fn);
        lprint(fn2);

        read(fn);
        eqs := ";

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

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

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

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