Thursday, October 20, 2011

VPython Code

To get some practice in class with programing and specificallyVPython, Brian and I had to do exercise number 5 of the second handout.
So here is the code and the result:


from math import *


A=zeros(100,float)
B=zeros(100,float)
C=zeros(100,float)
Sum=0
D=float


for i in range(0,100):
    A[i]= sin(2*pi*i/100)
    B[i]= cos(2*pi*i/100)
    C[i]=A[i]*B[i]
for i in range(0,100):
    D=Sum+C[i]
    Sum=D
    
        
print A
print B
print C
print "sum of all values is", Sum 





[  0.00000000e+00   6.27905195e-02   1.25333234e-01   1.87381315e-01
   2.48689887e-01   3.09016994e-01   3.68124553e-01   4.25779292e-01
   4.81753674e-01   5.35826795e-01   5.87785252e-01   6.37423990e-01
   6.84547106e-01   7.28968627e-01   7.70513243e-01   8.09016994e-01
   8.44327926e-01   8.76306680e-01   9.04827052e-01   9.29776486e-01
   9.51056516e-01   9.68583161e-01   9.82287251e-01   9.92114701e-01
   9.98026728e-01   1.00000000e+00   9.98026728e-01   9.92114701e-01
   9.82287251e-01   9.68583161e-01   9.51056516e-01   9.29776486e-01
   9.04827052e-01   8.76306680e-01   8.44327926e-01   8.09016994e-01
   7.70513243e-01   7.28968627e-01   6.84547106e-01   6.37423990e-01
   5.87785252e-01   5.35826795e-01   4.81753674e-01   4.25779292e-01
   3.68124553e-01   3.09016994e-01   2.48689887e-01   1.87381315e-01
   1.25333234e-01   6.27905195e-02   1.22464680e-16  -6.27905195e-02
  -1.25333234e-01  -1.87381315e-01  -2.48689887e-01  -3.09016994e-01
  -3.68124553e-01  -4.25779292e-01  -4.81753674e-01  -5.35826795e-01
  -5.87785252e-01  -6.37423990e-01  -6.84547106e-01  -7.28968627e-01
  -7.70513243e-01  -8.09016994e-01  -8.44327926e-01  -8.76306680e-01
  -9.04827052e-01  -9.29776486e-01  -9.51056516e-01  -9.68583161e-01
  -9.82287251e-01  -9.92114701e-01  -9.98026728e-01  -1.00000000e+00
  -9.98026728e-01  -9.92114701e-01  -9.82287251e-01  -9.68583161e-01
  -9.51056516e-01  -9.29776486e-01  -9.04827052e-01  -8.76306680e-01
  -8.44327926e-01  -8.09016994e-01  -7.70513243e-01  -7.28968627e-01
  -6.84547106e-01  -6.37423990e-01  -5.87785252e-01  -5.35826795e-01
  -4.81753674e-01  -4.25779292e-01  -3.68124553e-01  -3.09016994e-01
  -2.48689887e-01  -1.87381315e-01  -1.25333234e-01  -6.27905195e-02]
[  1.00000000e+00   9.98026728e-01   9.92114701e-01   9.82287251e-01
   9.68583161e-01   9.51056516e-01   9.29776486e-01   9.04827052e-01
   8.76306680e-01   8.44327926e-01   8.09016994e-01   7.70513243e-01
   7.28968627e-01   6.84547106e-01   6.37423990e-01   5.87785252e-01
   5.35826795e-01   4.81753674e-01   4.25779292e-01   3.68124553e-01
   3.09016994e-01   2.48689887e-01   1.87381315e-01   1.25333234e-01
   6.27905195e-02   6.12323400e-17  -6.27905195e-02  -1.25333234e-01
  -1.87381315e-01  -2.48689887e-01  -3.09016994e-01  -3.68124553e-01
  -4.25779292e-01  -4.81753674e-01  -5.35826795e-01  -5.87785252e-01
  -6.37423990e-01  -6.84547106e-01  -7.28968627e-01  -7.70513243e-01
  -8.09016994e-01  -8.44327926e-01  -8.76306680e-01  -9.04827052e-01
  -9.29776486e-01  -9.51056516e-01  -9.68583161e-01  -9.82287251e-01
  -9.92114701e-01  -9.98026728e-01  -1.00000000e+00  -9.98026728e-01
  -9.92114701e-01  -9.82287251e-01  -9.68583161e-01  -9.51056516e-01
  -9.29776486e-01  -9.04827052e-01  -8.76306680e-01  -8.44327926e-01
  -8.09016994e-01  -7.70513243e-01  -7.28968627e-01  -6.84547106e-01
  -6.37423990e-01  -5.87785252e-01  -5.35826795e-01  -4.81753674e-01
  -4.25779292e-01  -3.68124553e-01  -3.09016994e-01  -2.48689887e-01
  -1.87381315e-01  -1.25333234e-01  -6.27905195e-02  -1.83697020e-16
   6.27905195e-02   1.25333234e-01   1.87381315e-01   2.48689887e-01
   3.09016994e-01   3.68124553e-01   4.25779292e-01   4.81753674e-01
   5.35826795e-01   5.87785252e-01   6.37423990e-01   6.84547106e-01
   7.28968627e-01   7.70513243e-01   8.09016994e-01   8.44327926e-01
   8.76306680e-01   9.04827052e-01   9.29776486e-01   9.51056516e-01
   9.68583161e-01   9.82287251e-01   9.92114701e-01   9.98026728e-01]
[  0.00000000e+00   6.26666168e-02   1.24344944e-01   1.84062276e-01
   2.40876837e-01   2.93892626e-01   3.42273553e-01   3.85256621e-01
   4.22163963e-01   4.52413526e-01   4.75528258e-01   4.91143625e-01
   4.99013364e-01   4.99013364e-01   4.91143625e-01   4.75528258e-01
   4.52413526e-01   4.22163963e-01   3.85256621e-01   3.42273553e-01
   2.93892626e-01   2.40876837e-01   1.84062276e-01   1.24344944e-01
   6.26666168e-02   6.12323400e-17  -6.26666168e-02  -1.24344944e-01
  -1.84062276e-01  -2.40876837e-01  -2.93892626e-01  -3.42273553e-01
  -3.85256621e-01  -4.22163963e-01  -4.52413526e-01  -4.75528258e-01
  -4.91143625e-01  -4.99013364e-01  -4.99013364e-01  -4.91143625e-01
  -4.75528258e-01  -4.52413526e-01  -4.22163963e-01  -3.85256621e-01
  -3.42273553e-01  -2.93892626e-01  -2.40876837e-01  -1.84062276e-01
  -1.24344944e-01  -6.26666168e-02  -1.22464680e-16   6.26666168e-02
   1.24344944e-01   1.84062276e-01   2.40876837e-01   2.93892626e-01
   3.42273553e-01   3.85256621e-01   4.22163963e-01   4.52413526e-01
   4.75528258e-01   4.91143625e-01   4.99013364e-01   4.99013364e-01
   4.91143625e-01   4.75528258e-01   4.52413526e-01   4.22163963e-01
   3.85256621e-01   3.42273553e-01   2.93892626e-01   2.40876837e-01
   1.84062276e-01   1.24344944e-01   6.26666168e-02   1.83697020e-16
  -6.26666168e-02  -1.24344944e-01  -1.84062276e-01  -2.40876837e-01
  -2.93892626e-01  -3.42273553e-01  -3.85256621e-01  -4.22163963e-01
  -4.52413526e-01  -4.75528258e-01  -4.91143625e-01  -4.99013364e-01
  -4.99013364e-01  -4.91143625e-01  -4.75528258e-01  -4.52413526e-01
  -4.22163963e-01  -3.85256621e-01  -3.42273553e-01  -2.93892626e-01
  -2.40876837e-01  -1.84062276e-01  -1.24344944e-01  -6.26666168e-02]
sum of all values is -4.56579218877e-15


the exercise asked to define two set of arrays and see if they are orthogonal to each other by getting the dot product of them. Since the final value is -4.56579218877e-15 which is very close to zero we can assume these two vectors are perpendicular to each other.

    

   

No comments:

Post a Comment