ipar(1) = 0
    ipar(2) = 1
    ipar(3) = 1
    ipar(4) = n * 7
    ipar(5) = 10
    ipar(6) = 10000
    fpar(1) = 1.0E-6
    fpar(2) = 1.0E-8
    fpar(11) = 0.0

    allocate(wk(n * 7))

10    call bcg(n, right, result, ipar, fpar, wk)

    print*, ipar(1)
    if (ipar(1).eq.1) then
        call amux(n, wk(ipar(8)), wk(ipar(9)), Ax, Ai, Ap)
        goto 10
    else if (ipar(1).eq.3 .or. ipar(1).eq.5) then
         call lusol(n, wk(ipar(8)), wk(ipar(9)), alu, jlu, ju)
         goto 10
    end if