I have to evaluate an integral over a triangular area:
int_0^1 int_0^y [f(x,y)] dx,dy
Initially I though I could use dblquad, but is seems the integration limits have to be numbers, not variables. I have backed of to using quadgk to do the inner integral, then summing over all y. Does anyone know a better way?