N := 6;

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

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

	h(5) := B;

        read(fn);
        eqs := ";

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

	ss := solve({eqs[3..13]},vars);
	Digits := 20;
        sB := fsolve(eqs[2],B);
	B := sB[n];
	h(5) := sB[n];

	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);