From 7a73162607544204032aa66cce755daf21edebda Mon Sep 17 00:00:00 2001
From: Matt Strapp <matt@mattstrapp.net>
Date: Tue, 24 May 2022 11:18:46 -0500
Subject: Graduate

Signed-off-by: Matt Strapp <matt@mattstrapp.net>
---
 ee1301/wk6/lab5/vectorArray.cpp | 43 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)
 create mode 100644 ee1301/wk6/lab5/vectorArray.cpp

(limited to 'ee1301/wk6/lab5/vectorArray.cpp')

diff --git a/ee1301/wk6/lab5/vectorArray.cpp b/ee1301/wk6/lab5/vectorArray.cpp
new file mode 100644
index 0000000..92d23bd
--- /dev/null
+++ b/ee1301/wk6/lab5/vectorArray.cpp
@@ -0,0 +1,43 @@
+#include <iostream>
+
+double randVec();
+double momentum (double velocity);
+int main() {
+  srand(time(NULL));
+  double sum[3]={0,0,0};
+  double* randVel;  double* momArray;
+  randVel = new double[1000];
+  momArray = new double[1000];
+  for (int i=0; i<1000; i++) {
+    randVel[i]=randVec();
+    momArray[i]=momentum(randVel[i]);
+    int j=i%3;
+    sum[j]+=momArray[i];
+  }
+  std::cout << '<' << sum[0]/1000.0 << ',' <<  sum[1]/1000.0 << ',' << sum[2]/1000.0 << '>' << std::endl;
+  delete[] randVel;
+  delete[] momArray;
+}
+
+
+double randVec() {
+  double vector[3];
+  for (int i=0; i<3; i++) {
+    if (rand()%2+1==1) {
+      vector[i]=rand()%100+1.0;
+    } else {
+      vector[i]=rand()%100*-1.0+1.0;
+    }
+  }
+  return *vector;
+}
+
+double momentum(double velocity) {
+  double mass;
+  if (rand()%2+1==0) {
+    mass=rand()%10+1.0;
+  } else {
+    mass=rand()%10*-1.0+1.0;
+  }
+  return mass*velocity;
+}
\ No newline at end of file
-- 
cgit v1.2.3