function [y,y_n]=fusion(p1,p1_n,p2,p2_n)  %p1_attack,p1_normal,p2_attack,p2_normal
y=(p1+p2-p1_n*p2-p1*p2_n-p1*p2)/(1-p1*p2_n-p1_n*p2);
y_n=(p1_n+p2_n-p1_n*p2-p1*p2_n-p1_n*p2_n)/(1-p1*p2_n-p1_n*p2);