pro funct_decline_regress,x,a,f,pder
;
cval=0.308446            ; need to update to ensure a smooth join
z=x-1900
f=cval+a(0)*z+a(1)*z*z
;
if n_params() ge 4 then begin
  pder=[[z],[z*z]]
endif
;
end