For j from 1 by 1 to 10
    do
        plot({[f(x,j),g(x,j), x=x_min..x_max]}, options)
    end do;