diff --git a/dipdip-dirty/lrhankel_recspace_dirty.c b/dipdip-dirty/lrhankel_recspace_dirty.c
index 05b013e..5648e4b 100644
--- a/dipdip-dirty/lrhankel_recspace_dirty.c
+++ b/dipdip-dirty/lrhankel_recspace_dirty.c
@@ -14,20 +14,21 @@ typedef complex double (*lrhankelspec)(double, double, double,
 		const complex double *,
 		const complex double *,
 		const complex double *,
+		const complex double *,
 		const complex double *);
 // complex double fun(double c, double k0, double k, ccd *a, ccd *b, ccd *d, ccd *e)
 
 
 complex double fk5q1n0l(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (FF*e[0]-5*e[1]+10*e[2]-10*e[3]+5*e[4]-e[5])/k0;
 }
 complex double fk5q1n1l(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (-FF*d[0]+5*d[1]-10*d[2]+10*d[3]-5*d[4]+d[5])/(k0*k);
 }
 complex double fk5q1n2l(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	double t = 2/(k*k);
 	return (     (FF*e[0] - t*a[0] + FF*t*d[0]*a[0])
 		-5 * (e[1] - t*a[1] + t*d[1]*a[1])
@@ -37,12 +38,19 @@ complex double fk5q1n2l(double c, double k0, double k,
 		-    (e[5] - t*a[5] + t*d[5]*a[5])
 	       )/k0;
 }
-complex double fk5q2n0l(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return 0; // FIXME
+complex double fk5q2n0(double c, double k0, double k,
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
+	return (
+		 -     ash[0]
+		 + 5 * ash[1]
+		 -10 * ash[2]
+		 +10 * ash[3]
+		 - 5 * ash[4]
+		 +     ash[5]
+	       ) / (k0*k0);
 }
 complex double fk5q2n1l(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (   FF *b[0]*a[0]
 		  - 5 *b[1]*a[1]
 		  +10 *b[2]*a[2]
@@ -52,7 +60,7 @@ complex double fk5q2n1l(double c, double k0, double k,
 	       )/(k*k0*k0);
 }
 complex double fk5q2n2l(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return  (       b[0]*a[0]*a[0]
 		  + 5 * b[1]*a[1]*a[1]
 		  -10 * b[2]*a[2]*a[2]
@@ -62,16 +70,29 @@ complex double fk5q2n2l(double c, double k0, double k,
 		) / (k*k*k0*k0);
 }
 
+complex double fk5q3n0l(double c, double k0, double k, 
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) { // FIXME	
+	return  ( /* FIXME */
+		  -     k*b[0] + a[0] * ash[0]
+		  + 5 * k*b[1] + a[1] * ash[1]
+		  -10 * k*b[2] + a[2] * ash[2]
+		  +10 * k*b[3] + a[3] * ash[3]
+		  - 5 * k*b[4] + a[4] * ash[4]
+		  +     k*b[5] + a[5] * ash[5]
+		)/(k0*k0*k0);
+}
+
+
 complex double fk5q1n0s(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (e[0]-5*e[1]+10*e[2]-10*e[3]+5*e[4]-e[5])/k0;
 }
 complex double fk5q1n1s(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (-d[0]+5*d[1]-10*d[2]+10*d[3]-5*d[4]+d[5])/(k0*k);
 }
 complex double fk5q1n2s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	double t = 2/(k*k);
 	return (     (e[0] - t*a[0] + t*d[0]*a[0])
 		-5 * (e[1] - t*a[1] + t*d[1]*a[1])
@@ -81,12 +102,11 @@ complex double fk5q1n2s(double c, double k0, double k,
 		-    (e[5] - t*a[5] + t*d[5]*a[5])
 	       )/k0;
 }
-complex double fk5q2n0s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return 0; // FIXME
-}
+
+lrhankelspec fk5q2n0s = fk5q2n0, fk5q2n0l = fk5q2n0;
+
 complex double fk5q2n1s(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return (   FF *b[0]*a[0]
 		  - 5 *b[1]*a[1]
 		  +10 *b[2]*a[2]
@@ -96,7 +116,7 @@ complex double fk5q2n1s(double c, double k0, double k,
 	       )/(k*k0*k0);
 }
 complex double fk5q2n2s(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
+		const complex double *a, const complex double *b, const complex double *d, const complex double *e, const complex double *ash) {
 	return  (  FF * b[0]*a[0]*a[0]
 		  + 5 * b[1]*a[1]*a[1]
 		  -10 * b[2]*a[2]*a[2]
@@ -106,49 +126,13 @@ complex double fk5q2n2s(double c, double k0, double k,
 		) / (k*k*k0*k0);
 }
 
-#if 0
-complex double fk5q1n0s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return (p[0]-5*p[1]+10*p[2]-10*p[3]+5*p[4]-p[5])/k0;
-}
-complex double fk5q1n1s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return ;
-}
-complex double fk5q1n2s(double c, double k0, double k, 
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return ;
-}
-complex double fk5q2n0s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return ;
-}
-complex double fk5q2n1s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return ;
-}
-complex double fk5q2n2s(double c, double k0, double k,
-		const complex double *a, const complex double *b, const complex double *d, const complex double *e) {
-	return ;
-}
-
-static lrhankelspec transfuns_n[MAXKAPPA+1][MAXQM+1][MAXN+1] = {
-	{NULL,NULL,NULL},{NULL,NULL,NULL},{NULL,NULL,NULL},
-	{NULL,NULL,NULL},{NULL,NULL,NULL},{NULL,NULL,NULL},
-	{NULL,NULL,NULL},{NULL,NULL,NULL},{NULL,NULL,NULL},
-	{NULL,NULL,NULL},{NULL,NULL,NULL},{NULL,NULL,NULL},
-	{NULL,NULL,NULL},{NULL,NULL,NULL},{NULL,NULL,NULL},
-	{TODO,TODO,TODO},{TODO,TODO,TODO},{TODO,TODO,TODO}
-};
-#endif
-
 static lrhankelspec transfuns_f[MAXKAPPA+1][MAXQM+1][MAXN+1] = {
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
-	{{fk5q1n0l,fk5q1n1l,fk5q1n2l},{fk5q2n0l/*FIXME*/,fk5q2n1l,fk5q2n2l}}
+	{{fk5q1n0l,fk5q1n1l,fk5q1n2l},{fk5q2n0,fk5q2n1l,fk5q2n2l}}
 };
 
 static lrhankelspec transfuns_n[MAXKAPPA+1][MAXQM+1][MAXN+1] = {
@@ -157,7 +141,7 @@ static lrhankelspec transfuns_n[MAXKAPPA+1][MAXQM+1][MAXN+1] = {
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
 	{{NULL,NULL,NULL},{NULL,NULL,NULL}},
-	{{fk5q1n0s,fk5q1n1s,fk5q1n2s},{fk5q2n0s/*FIXME*/,fk5q2n1s,fk5q2n2s}}
+	{{fk5q1n0s,fk5q1n1s,fk5q1n2s},{fk5q2n0,fk5q2n1s,fk5q2n2s}}
 };
 
 void lrhankel_recpart_fill(complex  double *target,
@@ -186,19 +170,21 @@ int main() {
 
 	for (double k = kmin; k <= kmax; k += kstep) {
 		printf("%f ", k);
-		complex double a[kappa+1], b[kappa+1], d[kappa+1], e[kappa+1];
+		complex double a[kappa+1], b[kappa+1], d[kappa+1], e[kappa+1], ash[kappa+1];
 		for (size_t sigma = 0; sigma <= kappa; ++sigma) {
 			a[sigma] = (sigma * c - I * k0);
 			b[sigma] = csqrt(1+k*k/(a[sigma]*a[sigma]));
 			d[sigma] = 1/b[sigma];
 			e[sigma] = d[sigma] / a[sigma];
+			ash[sigma] = casinh(a[sigma]/k);
 		}
 		for (size_t qm = 0; qm <= MAXQM; ++qm)
 			for (size_t n = 0; n <= MAXN; ++n) 
-			if (!((qm==1)&&(n==0))){ // skip q==2, n=0 for now
+			if (/*!*/((qm==1)&&(n==0))){ //  not skip q==2, n=0 for now
 // complex double fun(double c, double k0, double k, ccd *a, ccd *b, ccd *d, ccd *e)
 				complex double result = 
-					transfuns_f[kappa][qm][n](c,k0,k,a,b,d,e);
+					//transfuns_f[kappa][qm][n](c,k0,k,a,b,d,e,ash);
+				 	fk5q2n0s(c,k0,k,a,b,d,e,ash);
 				printf("%.16e %.16e ", creal(result), cimag(result));
 			}
 		printf("\n");