function [I] = f3_s_2y(q,scaling,F_rad,S_phi,S_Z1,S_K1,S_Z2,S_K2)
Q=q*2*F_rad;
S=CalTYSk(S_Z1,S_Z2,S_K1,S_K2,S_phi,Q);
