Q2=10 and Q2=1--PLB(Fig.4)

(个人学习笔记,仅供参考)

python 复制代码
import numpy as np
from scipy.special import kv, erfc
from scipy.integrate import dblquad
import matplotlib.pyplot as plt
import scipy.integrate as spi
python 复制代码
# Constants
w = 0.6198
g0_sq = 21.5989
rho = 0.782908
Q2 = 10.0
P2 = 0.01
z0 = 4.25
python 复制代码
for x in x_values:
    result, error = compute_integral(x)
    integral_results178.append(result)
    errors.append(error)

# Convert results to numpy arrays for further processing or plotting
integral_results178 = np.array(integral_results178)
errors = np.array(errors)
python 复制代码
integral_results178T = integral_results178 - integral_results178L

integral_results178R = integral_results178L/integral_results178T
python 复制代码
plt.figure(figsize=(5, 3))
plt.plot(x_values, integral_results1R, label='$Q^2=1$GeV$^2$', color='red')
plt.plot(x_values, integral_results178R, label='$Q^2=10$GeV$^2$', color='blue')
plt.xlabel('x')
plt.xscale('log')
plt.xlim(1e-5, 0.1)
plt.ylim(-0.05, 0)  # Adjust y-axis based on max result
plt.ylabel('$F_2^{\gamma}(x,Q^2)$')
plt.legend()
plt.show()
相关推荐
齐 飞3 分钟前
MongoDB笔记01-概念与安装
前端·数据库·笔记·后端·mongodb
云空4 分钟前
《Python 与 SQLite:强大的数据库组合》
数据库·python·sqlite
丫头,冲鸭!!!37 分钟前
B树(B-Tree)和B+树(B+ Tree)
笔记·算法
凤枭香1 小时前
Python OpenCV 傅里叶变换
开发语言·图像处理·python·opencv
测试杂货铺1 小时前
外包干了2年,快要废了。。
自动化测试·软件测试·python·功能测试·测试工具·面试·职场和发展
听忆.1 小时前
手机屏幕上进行OCR识别方案
笔记
艾派森1 小时前
大数据分析案例-基于随机森林算法的智能手机价格预测模型
人工智能·python·随机森林·机器学习·数据挖掘
小码的头发丝、1 小时前
Django中ListView 和 DetailView类的区别
数据库·python·django
Selina K2 小时前
shell脚本知识点记录
笔记·shell
2 小时前
开源竞争-数据驱动成长-11/05-大专生的思考
人工智能·笔记·学习·算法·机器学习