integer q1,q4
      parameter(q1=1,q4=2)