Normal Derivative values at boundary

Dear All

I am trying to get the normal derivative value at a specified boundary.

for this I have written code as:

real Nu  = int1d(Th, hot)( dx(T)*N.x + dy(T)*N.y); 

for( int n1 = 0; n1<=20 ; n1++)
{  court<< Nu(20/n1, 0) << endl;
}

Am I right for Nu values at a boundary hot.

varf vNeuman(used,v) = int1d(Th)(v*(dx(Uex)*N.x+dy(Uex)*N.y)/lenEdge);
Vh vn; // Finite element function to store the Normal derivative
vn[]= vNeuman(0,Vh);

plot(vn,wait=1);