subroutine pda_deflt(alg, iv, liv, lv, v)
  
    ***  supply ***sol (version 2.3) default values to iv and v  ***
  
    ***  alg = 1 means regression constants.
    ***  alg = 2 means general unconstrained optimization constants.