1 files changed, 144 insertions, 0 deletions
diff --git a/3845/CH33/EX33.3/Ex33_3.sce b/3845/CH33/EX33.3/Ex33_3.sce
new file mode 100644
index 000000000..42dfd9414
--- /dev/null
+++ b/3845/CH33/EX33.3/Ex33_3.sce
@@ -0,0 +1,144 @@
+//Example 33.3
+//See Table 33.2
+//Part (a)
+S_negative_Xi=-2;//Strangeness of negative Xi
+S_neutral_Lambda=-1;//Strangeness of neutral Lambda
+S_negative_Pi=0;//Strangeness of negative Pion 
+S_change=abs((S_negative_Pi+S_neutral_Lambda)-S_negative_Xi);//Change in strangeness
+printf('a.Strangeness changes by +%d',S_change)
+
+B_negative_Xi=+1;//Baryon number of negative Xi
+B_neutral_Lambda=+1;//Baryon number of neutral Lambda 
+B_negative_Pi=0;//Baryon number of negative Pion
+if  (B_negative_Xi==B_neutral_Lambda+B_negative_Pi)
+       printf('\n  Baryon number is conserved')
+else
+       printf('\n  Baryon number is not conserved')
+end
+
+Charge_negative_Xi=-1;//Charge of negative Xi
+Charge_neutral_Lambda=0;//Charge of neutral Lambda
+Charge_negative_Pi=-1;//Charge of negative Pion
+if Charge_negative_Xi==(Charge_neutral_Lambda+Charge_negative_Pi)
+    printf('\n  Charge is conserved')
+else
+    printf('\n  Charge is not conserved')
+end
+
+printf('\n  Lepton numbers:')
+Le_negative_Xi=0;//Electron family number of negative Xi
+Le_neutral_Lambda=0;//Electron family number of neutral Lambda 
+Le_negative_Pi=0;//Electron family number of negative Pion
+if  (Le_negative_Xi==Le_neutral_Lambda+Le_negative_Pi)
+       printf('\n  Electron family number is conserved')
+else
+       printf('\n  Electron family number is not conserved')
+end
+
+Ltau_negative_Xi=0;//Tau family number of negative Xi
+Ltau_neutral_Lambda=0;//Tau family number of neutral Lambda 
+Ltau_negative_Pi=0;//Tau family number of negative Pion
+if  (Ltau_negative_Xi==Ltau_neutral_Lambda+Ltau_negative_Pi)
+       printf('\n  Tau family number is conserved')
+else
+       printf('\n  Tau family number is not conserved')
+end
+
+Lmu_negative_Xi=0;//Muon family number of negative Xi
+Lmu_neutral_Lambda=0;//Muon family number of neutral Lambda 
+Lmu_negative_Pi=0;//Muon family number of negative Pion
+if  (Lmu_negative_Xi==Lmu_neutral_Lambda+Lmu_negative_Pi)
+       printf('\n  Muon family number is conserved')
+else
+       printf('\n  Muon family number is not conserved')
+end
+
+///////////////////////////////////////////////////////////////
+//Part (b)
+i=1;//An indicator; i is set to zero if any of the conditions fail
+
+//It is stated in the textbook that charge is conserved as s->d. The steps followed in this code analyze individual charges of the quarks composing the positive kaon and the leptons produced by its decay.
+//Positive Kaon is made of two quarks- up and strange antiquark, See Table 33.4
+charge_u=+2/3;//Charge of u quark, See Table 33.3
+charge_s_anti=+1/3//Charge of s antiquark, See Table 33.3
+charge_positive_Kaon=charge_u+charge_s_anti;//Charge of positive Kaon
+charge_positive_Muon=+1;//Charge of positive Muon
+charge_Mu_Neutrino=0;//Charge of Muon Neutrino
+if  (charge_positive_Kaon==charge_positive_Muon+charge_Mu_Neutrino)
+       printf('\n\nb.Charge is conserved')
+else
+       printf('\n\nb.Charge is not conserved')
+       i=0;
+end
+
+B_positive_Kaon=0;//Baryon number of positive Kaon
+B_positive_Muon=0;//Baryon number of positive Muon 
+B_Mu_Neutrino=0;//Baryon number of Muon Neutrino
+if  (B_positive_Kaon==B_positive_Muon+B_Mu_Neutrino)
+       printf('\n  Baryon number is conserved')
+else
+       printf('\n  Baryon number is not conserved')
+       i=0;
+end
+
+m_positive_Kaon=493.7;//Rest mass of positive Kaon (MeV/c^2)
+m_positive_Muon=105.7;//Rest mass of positive Muon (MeV/c^2)
+m_Mu_Neutrino=0;//Rest mass Muon Neutrino (MeV/c^2)
+//Decay can be spontaneous if positive Kaon has greater mass than the products of decay
+if (m_positive_Kaon>(m_positive_Muon+m_Mu_Neutrino))
+    printf('\n  Mass-energy is conserved')
+else
+    printf('\n  Mass-energy is not conserved')
+    i=0;
+end
+
+S_positive_Kaon=+1;//Strangeness of positive Kaon
+S_positive_Muon=0;//Strangeness of positive Muon
+S_Mu_Neutrino=0;//Strangeness of Muon Neutrino 
+if abs((S_positive_Muon+S_Mu_Neutrino)-S_positive_Kaon)<=1
+    printf('\n  Strangeness is conserved')
+else
+    printf('\n  Strangeness is not conserved')
+    i=0;
+end
+
+printf('\n  Lepton numbers:')
+Le_positive_Kaon=0;//Electron family number of positive Kaon
+Le_positive_Muon=0;//Electron family number of positive Muon 
+Le_Mu_Neutrino=0;//Electron family number of Muon Neutrino
+if  (Le_positive_Kaon==Le_positive_Muon+Le_Mu_Neutrino)
+       printf('\n  Electron family number is conserved')
+else
+       printf('\n  Electron family number is not conserved')
+       i=0;
+end
+
+Ltau_positive_Kaon=0;//Tau family number of positive Kaon
+Ltau_positive_Muon=0;//Tau family number of positive Muon 
+Ltau_Mu_Neutrino=0;//Tau family number of Muon Neutrino
+if  (Ltau_positive_Kaon==Ltau_positive_Muon+Ltau_Mu_Neutrino)
+       printf('\n  Tau family number is conserved')
+else
+       printf('\n  Tau family number is not conserved')
+       i=0;
+end
+
+Lmu_positive_Kaon=0;//Muon family number of positive Kaon
+Lmu_positive_Muon=-1;//Muon family number of positive Muon
+Lmu_Mu_Neutrino=+1;//Muon family number of Muon Neutrino
+if  (Lmu_positive_Kaon==Lmu_positive_Muon+Lmu_Mu_Neutrino)
+       printf('\n  Muon family number is conserved')
+else
+       printf('\n  Muon family number is not conserved')
+       i=0;
+end
+
+
+if i==1
+    printf('\n\n  The decay is allowed')
+else
+    printf('\n\n  The decay is not allowed')
+end
+//Some of the solution steps are not given in the textbook; the steps mentioned here may need to be check for theoretical correctness 
+//Openstax - College Physics
+//Download for free at http://cnx.org/content/col11406/latest