Skip to content
Snippets Groups Projects
Commit b04e9503 authored by Andres Felipe Pabon Orozco's avatar Andres Felipe Pabon Orozco
Browse files

Upload New File

parent 44bfb45e
Branches main
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 469 additions and 0 deletions
Untitled.ipynb 0 → 100644
+ 469
− 0
View file @ b04e9503
%% Cell type:markdown id:6a48eeee-a178-4cbb-9aa2-1653b9026e21 tags:
grafica de primos hasta el 1000
%% Cell type:code id:3f5faedb-41fa-4d68-ba05-9414218927e2 tags:
``` python
import numpy as np
import matplotlib.pyplot as plt
```
%% Cell type:code id:7655ceec-69a8-4ed6-9768-08bf522c5d83 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:6c485dc9-ac8a-4a64-8705-c1d8e1b4ac7a tags:
``` python
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:84c78578-d238-4092-a1b1-d397a030c47f tags:
``` python
len(num_prim)
```
%% Output
35
%% Cell type:code id:18fd005b-9e5d-43f5-b16d-0fca7332193c tags:
``` python
def num_prim_func(N):
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:
break
if cond !=0:
#print(num[i], "es primo")
#print(num[j], "no 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:9f9e616d-c830-4942-a429-a5fbea2d1242 tags:
``` python
```
%% Cell type:code id:e6b54d00-23ab-442b-9b75-cc874c8d0809 tags:
``` python
counting_primes(1e4)
```
%% Output
1229
%% Cell type:code id:0d8521cc-193f-4440-a77d-068cb661c021 tags:
``` python
nume=np.arange(2,1001)
y=[]
for n in np.arange (0,len(nume)):
y.append(counting_primes(num[n]))
```
%% Cell type:code id:65ad1293-5d8d-4ec9-afd2-105e50564242 tags:
``` python
```
%% Cell type:code id:0a8820d0-104e-4258-83ce-9d5336c8d149 tags:
``` python
plt.figure()
plt.plot (nume,y,".")
plt.title("cantidad de primos")
plt.plot(nume,np.rint(nume/np.log(nume)))
```
%% Output
[<matplotlib.lines.Line2D at 0x223a4f2bb50>]
%% Cell type:code id:6d8e19f2-f0bb-47e1-92db-fa4a22c9f37f tags:
``` python
f=np.rint(nume/np.log(nume))
error_p=np.abs(y-f)/y*100
```
%% Cell type:code id:b8210c4d-b6f6-40d9-b666-b0dbf28d6e42 tags:
``` python
plt.figure()
plt.plot (nume[10:],error_p[10:])
```
%% Output
[<matplotlib.lines.Line2D at 0x223a4ecf670>]
%% Cell type:code id:55867cc5-f03d-4fef-b4cc-c345b7768b09 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