import numpy linalg import numpy random from matplotlib import pyplot

  • Anonymous
  • Python
  • 22 Jan 2020 
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
import numpy.linalg 

import numpy.random

from matplotlib import pyplot as plt

import numpy as np





def calculate_b_c(p):

    A = [[p[0][0]*p[0][1], p[0][1]*p[0][1], p[0][0], p[0][1], 1],

         [p[1][0]*p[1][1], p[1][1]*p[1][1], p[1][0], p[1][1], 1],

         [p[2][0]*p[2][1], p[2][1]*p[2][1], p[2][0], p[2][1], 1],

         [p[3][0]*p[3][1], p[3][1]*p[3][1], p[3][0], p[3][1], 1],

         [p[4][0]*p[4][1], p[4][1]*p[4][1], p[4][0], p[4][1], 1]]

    

    B = [[p[0][0]*p[0][0], p[0][1]*p[0][1], p[0][0], p[0][1], 1],

         [p[1][0]*p[1][0], p[1][1]*p[1][1], p[1][0], p[1][1], 1],

         [p[2][0]*p[2][0], p[2][1]*p[2][1], p[2][0], p[2][1], 1],

         [p[3][0]*p[3][0], p[3][1]*p[3][1], p[3][0], p[3][1], 1],

         [p[4][0]*p[4][0], p[4][1]*p[4][1], p[4][0], p[4][1], 1]]

    

    C = [[p[0][0]*p[0][0], p[0][0]*p[0][1], p[0][0], p[0][1], 1],

         [p[1][0]*p[1][0], p[1][0]*p[1][1], p[1][0], p[1][1], 1],

         [p[2][0]*p[2][0], p[2][0]*p[2][1], p[2][0], p[2][1], 1],

         [p[3][0]*p[3][0], p[3][0]*p[3][1], p[3][0], p[3][1], 1],

         [p[4][0]*p[4][0], p[4][0]*p[4][1], p[4][0], p[4][1], 1]]

    

    D = [[p[0][0]*p[0][0], p[0][0]*p[0][1], p[0][1]*p[0][1], p[0][1], 1],

         [p[1][0]*p[1][0], p[1][0]*p[1][1], p[1][1]*p[1][1], p[1][1], 1],

         [p[2][0]*p[2][0], p[2][0]*p[2][1], p[2][1]*p[2][1], p[2][1], 1],

         [p[3][0]*p[3][0], p[3][0]*p[3][1], p[3][1]*p[3][1], p[3][1], 1],

         [p[4][0]*p[4][0], p[4][0]*p[4][1], p[4][1]*p[4][1], p[4][1], 1]]

    

    E = [[p[0][0]*p[0][0], p[0][0]*p[0][1], p[0][1]*p[0][1], p[0][0], 1],

         [p[1][0]*p[1][0], p[1][0]*p[1][1], p[1][1]*p[1][1], p[1][0], 1],

         [p[2][0]*p[2][0], p[2][0]*p[2][1], p[2][1]*p[2][1], p[2][0], 1],

         [p[3][0]*p[3][0], p[3][0]*p[3][1], p[3][1]*p[3][1], p[3][0], 1],

         [p[4][0]*p[4][0], p[4][0]*p[4][1], p[4][1]*p[4][1], p[4][0], 1]]

    

    F = [[p[0][0]*p[0][0], p[0][0]*p[0][1], p[0][1]*p[0][1], p[0][0], p[0][1]],

         [p[1][0]*p[1][0], p[1][0]*p[1][1], p[1][1]*p[1][1], p[1][0], p[1][1]],

         [p[2][0]*p[2][0], p[2][0]*p[2][1], p[2][1]*p[2][1], p[2][0], p[2][1]],

         [p[3][0]*p[3][0], p[3][0]*p[3][1], p[3][1]*p[3][1], p[3][0], p[3][1]],

         [p[4][0]*p[4][0], p[4][0]*p[4][1], p[4][1]*p[4][1], p[4][0], p[4][1]]]

    

    a = numpy.linalg.det(A)

    b = numpy.linalg.det(B)

    c = numpy.linalg.det(C)

    d = numpy.linalg.det(D)

    e = numpy.linalg.det(E)

    f = numpy.linalg.det(F)

    

    

    

    return [a,b,c,d,e,f]

    

points = np.array([[numpy.random.randint(-10,10),numpy.random.randint(-10,10)] for _ in range(5)])

print(points)

print (calculate_b_c(points))



plt.scatter(points[:,0], points[:,1])
← newer dump
older dump →