p=0.001:0.0001:9.302;
 I01p =1.0000000027+24.9999990983.*(0.01075*p).^2+156.2500502179.*(0.01075*p).^4+434.0266800860.*(0.01075*p).^6+678.1808514731.*(0.01075*p).^8+678.0845424100.*(0.01075*p).^10+471.3119673418.*(0.01075*p).^12+239.2391892805.*(0.01075*p).^14+95.9019523737.*(0.01075*p).^16+26.2421075769.*(0.01075*p).^18+9.6664625536.*(0.01075*p).^20+0.1439945401.*(0.01075*p).^22+0.6688315150.*(0.01075*p).^24;
 I11p =0.1075.*p.*(0.50000000009+6.24999996950.*(0.01075*p).^2+26.04166836566.*(0.01075*p).^4+54.25343507056.*(0.01075*p).^6+67.81726179804.*(0.01075*p).^8+56.51119244161.*(0.01075*p).^10+33.65156773642.*(0.01075*p).^12+14.98220505757.*(0.01075*p).^14+5.28380368332.*(0.01075*p).^16+1.35566931981.*(0.01075*p).^18+0.41187059355.*(0.01075*p).^20+0.01604459097.*(0.01075*p).^22+0.02411174292.*(0.01075*p).^24);
 K01p =(0.577215664901533+log(0.1075/2.*p)).*(-1.00000000013-0.24999999338.*(0.1075*p).^2-0.01562505298.*(0.1075*p).^4-0.00043387943.*(0.1075*p).^6-0.00000695123.*(0.1075*p).^8)+0.01155625.*p.^2.*(0.250000000002+0.023437499887.*(0.1075*p).^2+0.004795718494.*(0.1075*p).^4+0.000014125984.*(0.1075*p).^6+0.000000157733.*(0.1075*p).^8);
J31 =0.00145.*p-0.001672868.*p.^3.*(0.577215664901533+log(0.059355.*p)).*(-0.1666666666675-0.0124999999575.*(0.11871*p).^2-0.0003720241492.*(0.11871*p).^4-0.0000060272126.*(0.11871*p).^6-0.0000000627385.*(0.11871*p).^8)-0.001672868*p.^3.*(0.1388888888910+0.0181249998963.*(0.11871*p).^2+0.0006731867710.*(0.11871*p).^4+0.0000124726259.*(0.11871*p).^6+0.0000001428666.*(0.11871*p).^8)+0.001612314.*p.^3.*(0.577215664901533+log(0.05863.*p)).*(-0.1666666666675-0.0124999999575.*(0.11726*p).^2-0.0003720241492.*(0.11726*p).^4-0.0000060272126.*(0.11726*p).^6-0.0000000627385.*(0.11726*p).^8)+0.001612314.*p.^3.*(0.1388888888910+0.0181249998963.*(0.11726*p).^2+0.0006731867710.*(0.11726*p).^4+0.0000124726259.*(0.11726*p).^6+0.0000001428666.*(0.11726*p).^8);
J121 =0.00145.*p-0.001584426.*p.^3.*(0.577215664901533+log(0.05829.*p)).*(-0.1666666666675-0.0124999999575.*(0.11658*p).^2-0.0003720241492.*(0.11658*p).^4-0.0000060272126.*(0.11658*p).^6-0.0000000627385.*(0.11658*p).^8)-0.001584426*p.^3.*(0.1388888888910+0.0181249998963.*(0.11658*p).^2+0.0006731867710.*(0.11658*p).^4+0.0000124726259.*(0.11658*p).^6+0.0000001428666.*(0.11658*p).^8)+0.001526038.*p.^3.*(0.577215664901533+log(0.057565.*p)).*(-0.1666666666675-0.0124999999575.*(0.11513*p).^2-0.0003720241492.*(0.11513*p).^4-0.0000060272126.*(0.11513*p).^6-0.0000000627385.*(0.11513*p).^8)+0.001526038.*p.^3.*(0.1388888888910+0.0181249998963.*(0.11513*p).^2+0.0006731867710.*(0.11513*p).^4+0.0000124726259.*(0.11513*p).^6+0.0000001428666.*(0.11513*p).^8);
 b3p=59/60.*I01p.*I11p.*0.1075.*p./(K01p.*I11p.*p.*0.1075*59/60+1/60);
  f1=sin(0.160/2.*p).*sin(0.165/2.*p).*cos(0.4075.*p);
 f=b3p.*f1.*J31.*J121./p.^6;
  pp=spline(p,f);
 int_pp=fnint(pp);
 Ssp=ppval(int_pp,[0.001,9.302])*[-1;1];
 s1=0.00145*0.160;
 s12=0.00145*0.165;
 M3121=8*4e-7*pi*Ssp/(s1*s12);
 M3121;
 
    p=9.302:0.0001:93.0233;
 I01p =1.0000000027+24.9999990983.*(0.01075*p).^2+156.2500502179.*(0.01075*p).^4+434.0266800860.*(0.01075*p).^6+678.1808514731.*(0.01075*p).^8+678.0845424100.*(0.01075*p).^10+471.3119673418.*(0.01075*p).^12+239.2391892805.*(0.01075*p).^14+95.9019523737.*(0.01075*p).^16+26.2421075769.*(0.01075*p).^18+9.6664625536.*(0.01075*p).^20+0.1439945401.*(0.01075*p).^22+0.6688315150.*(0.01075*p).^24;
 I11p =0.1075.*p.*(0.50000000009+6.24999996950.*(0.01075*p).^2+26.04166836566.*(0.01075*p).^4+54.25343507056.*(0.01075*p).^6+67.81726179804.*(0.01075*p).^8+56.51119244161.*(0.01075*p).^10+33.65156773642.*(0.01075*p).^12+14.98220505757.*(0.01075*p).^14+5.28380368332.*(0.01075*p).^16+1.35566931981.*(0.01075*p).^18+0.41187059355.*(0.01075*p).^20+0.01604459097.*(0.01075*p).^22+0.02411174292.*(0.01075*p).^24);
 K01p=exp(-0.1075.*p)./((0.1075)^(1/2).*p.^(1/2)).*(1.2533141373-0.1566640293./(0.1075*p)+0.0880947319./(0.1075*p).^2-0.0909977576./(0.1075*p).^3+0.1303595356./(0.1075*p).^4-0.2065010623./(0.1075*p).^5+0.2970331052./(0.1075*p).^6-0.3390507775./(0.1075*p).^7+0.2809180905./(0.1075*p).^8-0.1557795499./(0.1075*p).^9+0.0512812921./(0.1075*p).^10-0.0075446371./(0.1075*p).^11);
 J32 =0.000725.*p-exp(-0.11871.*p).*(0.11871)^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11871*p)-0.6950910682./(0.11871*p).^2+1.1666191735./(0.11871*p).^3-2.9952418887./(0.11871*p).^4+9.4332266362./(0.11871*p).^5-30.6134088398./(0.11871*p).^6+90.4117679352./(0.11871*p).^7-226.0598597708./(0.11871*p).^8+460.4696268423./(0.11871*p).^9-747.4466697296./(0.11871*p).^10+951.1027599993./(0.11871*p).^11-932.4652438657./(0.11871*p).^12+688.1998022368./(0.11871*p).^13-369.0779601792./(0.11871*p).^14+135.5755240815./(0.11871*p).^15-30.4795722146./(0.11871*p).^16+3.1605958754./(0.11871*p).^17)+exp(-0.11726.*p).*0.11726^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11726*p)-0.6950910682./(0.11726*p).^2+1.1666191735./(0.11726*p).^3-2.9952418887./(0.11726*p).^4+9.4332266362./(0.11726*p).^5-30.6134088398./(0.11726*p).^6+90.4117679352./(0.11726*p).^7-226.0598597708./(0.11726*p).^8+460.4696268423./(0.11726*p).^9-747.4466697296./(0.11726*p).^10+951.1027599993./(0.11726*p).^11-932.4652438657./(0.11726*p).^12+688.1998022368./(0.11726*p).^13-369.0779601792./(0.11726*p).^14+135.5755240815./(0.11726*p).^15-30.4795722146./(0.11726*p).^16+3.1605958754./(0.11726*p).^17);
 J122 =0.000725.*p-exp(-0.11658.*p).*(0.11658)^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11658*p)-0.6950910682./(0.11658*p).^2+1.1666191735./(0.11658*p).^3-2.9952418887./(0.11658*p).^4+9.4332266362./(0.11658*p).^5-30.6134088398./(0.11658*p).^6+90.4117679352./(0.11658*p).^7-226.0598597708./(0.11658*p).^8+460.4696268423./(0.11658*p).^9-747.4466697296./(0.11658*p).^10+951.1027599993./(0.11658*p).^11-932.4652438657./(0.11658*p).^12+688.1998022368./(0.11658*p).^13-369.0779601792./(0.11658*p).^14+135.5755240815./(0.11658*p).^15-30.4795722146./(0.11658*p).^16+3.1605958754./(0.11658*p).^17)+exp(-0.11513.*p).*0.11513^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11513*p)-0.6950910682./(0.11513*p).^2+1.1666191735./(0.11513*p).^3-2.9952418887./(0.11513*p).^4+9.4332266362./(0.11513*p).^5-30.6134088398./(0.11513*p).^6+90.4117679352./(0.11513*p).^7-226.0598597708./(0.11513*p).^8+460.4696268423./(0.11513*p).^9-747.4466697296./(0.11513*p).^10+951.1027599993./(0.11513*p).^11-932.4652438657./(0.11513*p).^12+688.1998022368./(0.11513*p).^13-369.0779601792./(0.11513*p).^14+135.5755240815./(0.11513*p).^15-30.4795722146./(0.11513*p).^16+3.1605958754./(0.11513*p).^17);
 b3p=59/60.*I01p.*I11p.*0.1075.*p./(K01p.*I11p.*p.*0.1075*59/60+1/60);
 f1=sin(0.160/2.*p).*sin(0.165/2.*p).*cos(0.4075.*p);
  f=b3p.*f1.*J32.*J122./p.^6;
  pp=spline(p,f);
 int_pp=fnint(pp);
 Ssp=ppval(int_pp,[9.302,93.0233])*[-1;1];
  s1=0.00145*0.160;
 s12=0.00145*0.165;
 M3122=8*4e-7*pi*Ssp/(s1*s12);
 M3122;

 M312=M3121+M3122;
 M312;