-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathranges_excel.py
51 lines (36 loc) · 1.08 KB
/
ranges_excel.py
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
# -*- coding: utf-8 -*-
"""
Created on Wed Nov 22 09:47:05 2023
@author: rnv
test file for range_excel function
"""
import time
import ray
import pandas as pd
import numpy as np
from scipy.stats import skewnorm
from scipy.stats import norm, kurtosis,skew, tmean,tvar
from srim import Layer, Target
from ranges import ranges_width
from wait_for_task import wait_for_tasks
from auto_srim import auto_srim
path=r"C:\Users\Documents\SRIMsim\"
ranges1,width1 = ranges_width(path)
import numpy as np
from scipy.stats import skewnorm
import matplotlib.pyplot as plt
kurt2=kurtosis(test[0:57,1]/np.sum(test[0:57,1]))
skewness2=skewnorm.skew(test[0:57,1])
mean=skewnorm.mean(test[0:57,1])
var=tvar(test[0:57,1])
pdffunc = extras.pdf_mvsk([1323, 904.7, -0.2128, 2.6706])
rangeh = np.arange(0, 5700, 0.1)
plt.plot(rangeh, pdffunc(rangeh))
plt.show()
fig, ax = plt.subplots(1, 1)
ax.hist(test[0:57,1],bins=58)
rng = np.random.default_rng(19680801)
N_points = 100000
n_bins = 20
# Generate two normal distributions
dist1 = rng.standard_normal(N_points)