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