diff --git a/Schrick-Noah_Project-1.R b/Schrick-Noah_Project-1.R
index e1f7fb2..a21d15d 100644
--- a/Schrick-Noah_Project-1.R
+++ b/Schrick-Noah_Project-1.R
@@ -7,8 +7,24 @@
 setwd(dirname(rstudioapi::getActiveDocumentContext()$path))
 
 ## Part A: compute binding energy
+computeBindingEnergy <- function(A, Z){
+  if (A %% 2 == 1) {
+    spin <- 0
+  }else if (A %% 2 == 0 && Z %% 2 == 0){
+    spin <- 12
+  }
+  else {
+    spin <- -12
+  }
+  
+  B <- (15.67*A) - (17.23*(A^(2/3))) - ((0.75)*((Z^2)/(A^(1/3)))) - ((93.2)*(((A-(2*Z))^2)/A)) + (spin/(A^(1/2)))
+  
+  return(B)
+}
 
 ## Part B: compute Z=28 and A=58
+computeBindingEnergy(58,28)
+(15.67*58) - (17.23*(58^(2/3))) - ((0.75)*((28^2)/(58^(1/3)))) - ((93.2)*(((58-(2*28))^2)/58)) + (12/(58^(1/2)))
 
 ## Part C: table with two columns: nucleons and binding energy per nucleon