function t=FF(phi,E)

t=zeros(max([size(phi);size(E)]));
i1=find(E>1);
i2=find(E<=1);
E=E+eps;
invE=1./(E);
invsqrE=sqrt(invE);

if (length(phi)>1 & length(E)>1)
    t(i1)=F(phi(i1),invE(i1)).*invsqrE(i1);
    t(i2)=F(asin(sin(phi(i2)).*invsqrE(i2)),E(i2));
elseif (length(E)>1)
    t(i1)=F(phi,invE(i1)).*invsqrE(i1);
    t(i2)=F(asin(sin(phi).*invsqrE(i2)),E(i2));
else
    if(E>1)
    t=F(phi,invE).*invsqrE;
    else
    t=F(asin(sin(phi).*invsqrE),E);
    end
    
end    
