Skip to content
Snippets Groups Projects
Commit 15386646 authored by Oscar Danilo Lopez Jaime's avatar Oscar Danilo Lopez Jaime
Browse files

Subir nuevo archivo

parent 29c663cf
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Primos_clase_2junio.ipynb 0 → 100644
+ 399
− 0
View file @ 15386646
%% Cell type:code id:628f22b5-0cec-45ac-aaf9-facef0fbfe12 tags:
``` python
import numpy as np
import matplotlib. pyplot as plt
```
%% Cell type:code id:39a4ff49-3acd-409b-ad76-4e387de2e1ee tags:
``` python
N= 150
num =np.arange(2, N+1)
num_prim = [num[0]]
for i in np.arange(1,len(num)):
for j in np.arange(0,i):
cond = num[i]%num[j]
if cond== 0:
print(num[i], "no es primo")
break
if cond != 0:
print(num[i], "es primo")
num_prim.append(num[i])
```
%% Output
3 es primo
4 no es primo
5 es primo
6 no es primo
7 es primo
8 no es primo
9 no es primo
10 no es primo
11 es primo
12 no es primo
13 es primo
14 no es primo
15 no es primo
16 no es primo
17 es primo
18 no es primo
19 es primo
20 no es primo
21 no es primo
22 no es primo
23 es primo
24 no es primo
25 no es primo
26 no es primo
27 no es primo
28 no es primo
29 es primo
30 no es primo
31 es primo
32 no es primo
33 no es primo
34 no es primo
35 no es primo
36 no es primo
37 es primo
38 no es primo
39 no es primo
40 no es primo
41 es primo
42 no es primo
43 es primo
44 no es primo
45 no es primo
46 no es primo
47 es primo
48 no es primo
49 no es primo
50 no es primo
51 no es primo
52 no es primo
53 es primo
54 no es primo
55 no es primo
56 no es primo
57 no es primo
58 no es primo
59 es primo
60 no es primo
61 es primo
62 no es primo
63 no es primo
64 no es primo
65 no es primo
66 no es primo
67 es primo
68 no es primo
69 no es primo
70 no es primo
71 es primo
72 no es primo
73 es primo
74 no es primo
75 no es primo
76 no es primo
77 no es primo
78 no es primo
79 es primo
80 no es primo
81 no es primo
82 no es primo
83 es primo
84 no es primo
85 no es primo
86 no es primo
87 no es primo
88 no es primo
89 es primo
90 no es primo
91 no es primo
92 no es primo
93 no es primo
94 no es primo
95 no es primo
96 no es primo
97 es primo
98 no es primo
99 no es primo
100 no es primo
101 es primo
102 no es primo
103 es primo
104 no es primo
105 no es primo
106 no es primo
107 es primo
108 no es primo
109 es primo
110 no es primo
111 no es primo
112 no es primo
113 es primo
114 no es primo
115 no es primo
116 no es primo
117 no es primo
118 no es primo
119 no es primo
120 no es primo
121 no es primo
122 no es primo
123 no es primo
124 no es primo
125 no es primo
126 no es primo
127 es primo
128 no es primo
129 no es primo
130 no es primo
131 es primo
132 no es primo
133 no es primo
134 no es primo
135 no es primo
136 no es primo
137 es primo
138 no es primo
139 es primo
140 no es primo
141 no es primo
142 no es primo
143 no es primo
144 no es primo
145 no es primo
146 no es primo
147 no es primo
148 no es primo
149 es primo
150 no es primo
%% Cell type:code id:023fa8d7-dce7-40db-b3c8-5799ffd2b8c9 tags:
``` python
print(num_prim)
```
%% Output
[2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149]
%% Cell type:code id:365d17da-08d7-498a-b76d-cdb27ce21741 tags:
``` python
len(num_prim)
```
%% Output
35
%% Cell type:code id:ad9ff5e9-c7f2-4893-947f-5bbf57fdafd9 tags:
``` python
def num_prim_func(N):
"""retorna los numeros primos"""
num =np.arange(2, N+1)
num_prim = [num[0]]
for i in np.arange(1,len(num)):
for j in np.arange(0,i):
cond = num[i]%num[j]
if cond== 0:
#print(num[i], "no es primo")
break
if cond != 0:
#print(num[i], "es primo")
num_prim.append(num[i])
return num_prim
def counting_primes(N):
num_prim = num_prim_func(N)
return len(num_prim)
```
%% Cell type:code id:f7192b33-7448-4e2b-a98a-94ba6363e174 tags:
``` python
counting_primes(10000)
```
%% Output
1229
%% Cell type:code id:7dfbcb81-064e-4fad-b544-d65d09c0c79f tags:
``` python
nums=np.arange(2,1001)
y=[]
for i in np.arange(0,len(nums)):
y.append(counting_primes(nums[i]))
y = np.array(y)
```
%% Cell type:code id:05b80860-aa12-435d-b593-cae20c413a74 tags:
``` python
plt.figure()
plt.plot(nums,y,".")
plt.plot(nums,np.rint(nums/np.log(nums)))
```
%% Output
[<matplotlib.lines.Line2D at 0x1e8bcae22c0>]
%% Cell type:code id:c1ae88a4-4a65-4971-8646-cb0b3968eb39 tags:
``` python
f=np.rint(nums/np.log(nums))
error_p= np.abs(y - f)/y*100
```
%% Cell type:code id:7503f174-54fe-40be-a0c0-36a2155a1190 tags:
``` python
plt.figure
plt.plot(nums[10:],error_p[10:])
```
%% Output
[<matplotlib.lines.Line2D at 0x1e8bc2eba30>]
%% Cell type:code id:01f20788-ff7e-4d32-a606-4ff5202cff71 tags:
``` python
```
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment