function res= phibc0(phia,phib) res= [ phia(1) phib(1) ];