I am not sure how many people actually use this but l=1 case seems to be missing in Vcalc.f This is easy to fix by adding if(l.eq.1) F(2)=J*5.d0 right after F(0)=U Ilya