function c=paral(a,b) c=a.*b./(a+b);