1 #include<bits/stdc++.h>

 2 using namespace std;

 3 int t,n,m,dx[4]={1,0,-1,0},dy[4]={0,1,0,-1};

 4 double f[2005][2005];

 5 int id(int x,int y){

 6     if ((x<1)||(y<1)||(x>n)||(y>m))return 0;

 7     return (x-1)*m+y;

 8 }

 9 double guess(){

10     for(int ii=n;ii;ii--){

11         for(int i=id(ii,1);i<=id(ii,m);i++){

12             double s=abs(f[i][i]);

13             for(int j=i+1;j<=id(ii,m);j++)s=max(s,abs(f[j][i]));

14             for(int j=i;j<=id(ii,m);j++)

15                 if (s==abs(f[j][i])){

16                     s=f[j][i];

17                     for(int k=id(ii-1,1);k<=id(ii,m);k++)swap(f[j][k],f[i][k]);

18                     swap(f[j][id(n,m)+1],f[i][id(n,m)+1]);

19                     break;

20                 }

21             for(int j=id(ii-1,1);j<=id(ii,m);j++)f[i][j]/=s;

22             f[i][id(n,m)+1]/=s;

23             for(int j=id(ii,1);j<=id(ii,m);j++){

24                 if (j==i)continue;

25                 s=f[j][i];

26                 for(int k=id(ii-1,1);k<=id(ii,m);k++)f[j][k]-=f[i][k]*s;

27                 f[j][id(n,m)+1]-=f[i][id(n,m)+1]*s;

28             }

29         }

30         if (ii==1)return f[1][id(n,m)+1];

31         for(int i=id(ii-1,1);i<=id(ii-1,m);i++)

32             for(int j=id(ii,1);j<=id(ii,m);j++){

33                 for(int k=id(ii-1,1);k<=id(ii-1,m);k++)f[i][k]-=f[i][j]*f[j][k];

34                 f[i][id(n,m)+1]-=f[i][j]*f[j][id(n,m)+1];

35             }

36     }

37 }

38 int main(){

39     while (scanf("%d%d",&n,&m)!=EOF){

40         if ((!n)&&(!m))return 0;

41         memset(f,0,sizeof(f));

42         for(int i=1;i<=n;i++)

43             for(int j=1;j<=m;j++)

44                 f[id(i,j)][id(i,j)]=f[id(i,j)][id(n,m)+1]=-1;

45         for(int k=0;k<4;k++)

46             for(int i=1;i<=n;i++)

47                 for(int j=1;j<=m;j++)

48                     scanf("%lf",&f[id(i,j)][id(i+dx[k],j+dy[k])]);

49         f[id(n,m)][id(n,m)+1]=0;

50         printf("%.8f\n",guess());

51     }

52 }