!    include "lapack.f90"
    use lapack95
    real*8 a(5,10),b(5,5)
    do 1 i=1,5
      b(i,1)=1.D0
1     a(i,i)=1.D0
    call dgels(a,b,'N',info)
    print *,"info=",info
101 format(I3,5F14.7)
    do 2 i=1,5
2        print 101,i,a(i,1:5)
    do 3 I=1,5
3        print 101,i,b(i,1:5)
    end