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);
J41 =0.00145.*p-0.001764542.*p.^3.*(0.577215664901533+log(0.06042.*p)).*(-0.1666666666675-0.0124999999575.*(0.12084*p).^2-0.0003720241492.*(0.12084*p).^4-0.0000060272126.*(0.12084*p).^6-0.0000000627385.*(0.12084*p).^8)-0.001764542*p.^3.*(0.1388888888910+0.0181249998963.*(0.12084*p).^2+0.0006731867710.*(0.12084*p).^4+0.0000124726259.*(0.12084*p).^6+0.0000001428666.*(0.12084*p).^8)+0.001701781.*p.^3.*(0.577215664901533+log(0.059695.*p)).*(-0.1666666666675-0.0124999999575.*(0.11939*p).^2-0.0003720241492.*(0.11939*p).^4-0.0000060272126.*(0.11939*p).^6-0.0000000627385.*(0.11939*p).^8)+0.001701781.*p.^3.*(0.1388888888910+0.0181249998963.*(0.11939*p).^2+0.0006731867710.*(0.11939*p).^4+0.0000124726259.*(0.11939*p).^6+0.0000001428666.*(0.11939*p).^8);
J041 =0.00145.*p-0.001764542.*p.^3.*(0.577215664901533+log(0.06042.*p)).*(-0.1666666666675-0.0124999999575.*(0.12084*p).^2-0.0003720241492.*(0.12084*p).^4-0.0000060272126.*(0.12084*p).^6-0.0000000627385.*(0.12084*p).^8)-0.001764542*p.^3.*(0.1388888888910+0.0181249998963.*(0.12084*p).^2+0.0006731867710.*(0.12084*p).^4+0.0000124726259.*(0.12084*p).^6+0.0000001428666.*(0.12084*p).^8)+0.001701781.*p.^3.*(0.577215664901533+log(0.059695.*p)).*(-0.1666666666675-0.0124999999575.*(0.11939*p).^2-0.0003720241492.*(0.11939*p).^4-0.0000060272126.*(0.11939*p).^6-0.0000000627385.*(0.11939*p).^8)+0.001701781.*p.^3.*(0.1388888888910+0.0181249998963.*(0.11939*p).^2+0.0006731867710.*(0.11939*p).^4+0.0000124726259.*(0.11939*p).^6+0.0000001428666.*(0.11939*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)).^2.*cos(0.*p);
 f=b3p.*f1.*J41.*J041./p.^6;
  pp=spline(p,f);
 int_pp=fnint(pp);
 Ssp=ppval(int_pp,[0.001,9.302])*[-1;1];
 si=0.00145*0.160;
 M4041=8*4e-7*pi*Ssp/si^2;
 M4041;
 
 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);
 J42 =0.000725.*p-exp(-0.12084.*p).*(0.12084)^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.12084*p)-0.6950910682./(0.12084*p).^2+1.1666191735./(0.12084*p).^3-2.9952418887./(0.12084*p).^4+9.4332266362./(0.12084*p).^5-30.6134088398./(0.12084*p).^6+90.4117679352./(0.12084*p).^7-226.0598597708./(0.12084*p).^8+460.4696268423./(0.12084*p).^9-747.4466697296./(0.12084*p).^10+951.1027599993./(0.12084*p).^11-932.4652438657./(0.12084*p).^12+688.1998022368./(0.12084*p).^13-369.0779601792./(0.12084*p).^14+135.5755240815./(0.12084*p).^15-30.4795722146./(0.12084*p).^16+3.1605958754./(0.12084*p).^17)+exp(-0.11939.*p).*0.11939^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11939*p)-0.6950910682./(0.11939*p).^2+1.1666191735./(0.11939*p).^3-2.9952418887./(0.11939*p).^4+9.4332266362./(0.11939*p).^5-30.6134088398./(0.11939*p).^6+90.4117679352./(0.11939*p).^7-226.0598597708./(0.11939*p).^8+460.4696268423./(0.11939*p).^9-747.4466697296./(0.11939*p).^10+951.1027599993./(0.11939*p).^11-932.4652438657./(0.11939*p).^12+688.1998022368./(0.11939*p).^13-369.0779601792./(0.11939*p).^14+135.5755240815./(0.11939*p).^15-30.4795722146./(0.11939*p).^16+3.1605958754./(0.11939*p).^17);
 J042 =0.000725.*p-exp(-0.12084.*p).*(0.12084)^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.12084*p)-0.6950910682./(0.12084*p).^2+1.1666191735./(0.12084*p).^3-2.9952418887./(0.12084*p).^4+9.4332266362./(0.12084*p).^5-30.6134088398./(0.12084*p).^6+90.4117679352./(0.12084*p).^7-226.0598597708./(0.12084*p).^8+460.4696268423./(0.12084*p).^9-747.4466697296./(0.12084*p).^10+951.1027599993./(0.12084*p).^11-932.4652438657./(0.12084*p).^12+688.1998022368./(0.12084*p).^13-369.0779601792./(0.12084*p).^14+135.5755240815./(0.12084*p).^15-30.4795722146./(0.12084*p).^16+3.1605958754./(0.12084*p).^17)+exp(-0.11939.*p).*0.11939^(1/2).*p.^(1/2).*(1.2533141373+1.0966488851./(0.11939*p)-0.6950910682./(0.11939*p).^2+1.1666191735./(0.11939*p).^3-2.9952418887./(0.11939*p).^4+9.4332266362./(0.11939*p).^5-30.6134088398./(0.11939*p).^6+90.4117679352./(0.11939*p).^7-226.0598597708./(0.11939*p).^8+460.4696268423./(0.11939*p).^9-747.4466697296./(0.11939*p).^10+951.1027599993./(0.11939*p).^11-932.4652438657./(0.11939*p).^12+688.1998022368./(0.11939*p).^13-369.0779601792./(0.11939*p).^14+135.5755240815./(0.11939*p).^15-30.4795722146./(0.11939*p).^16+3.1605958754./(0.11939*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)).^2.*cos(0.*p);
  f=b3p.*f1.*J42.*J042./p.^6;
  pp=spline(p,f);
 int_pp=fnint(pp);
 Ssp=ppval(int_pp,[9.302,93.0233])*[-1;1];
  si=0.00145*0.160;
 M4042=8*4e-7*pi*Ssp/si^2;
 M4042;
 M404=M4041+M4042;