我正在尝试在平方区域上积分独立的二元正态分布。数值积分与蒙特卡罗模拟不匹配。这是怎么回事？ import numpy as np from scipy import integrate from s...

我正在尝试使用函数 quadde 在 Julia 中使用 DoubleExponentialFormulas.jl 计算积分。但是，在将结果与公式中的 Mathematica 进行比较时，我得到的相对...

我正在为我的数值方法科目做一项任务，他们要求使用 Runge-Kutta 和 Taylor 方法来估计微分系统（Chua 电路）的解，其场由下式确定： F(x,y,z)=\[2\*(y-0.5\*x...

对于一个学校项目，我正在 Python 中实现一种有限差分方法，以数值求解以下反应扩散偏微分方程组： 到目前为止，我已经实现了以下目标： import numpy as np import mat...

我尝试了一种非常幼稚的方法来逼近 Lisp 中函数的一阶导数，并得出了这样的东西： (defparameter *delta-x* 0.00001) (defun diff (f x) (/ (...

我正在阅读我的前辈同事编写的用于计算大圆距离的 matlab 函数。 地球表面两点之间的距离应使用以下公式计算： d = r * arccos[(sin(lat1) * sin(lat2)) + c...

我对浮点数的表示方式略知一二，但恐怕还不够。 一般问题是： 对于给定的精度（就我而言，以 10 为基数的精确小数位数），16 位、32 位和 64 位 IEEE-754 系统可以表示多大范围的数字...

我有这个代码 rm(list=ls()) N = 20000 xvar <- runif(N, -10, 10) e <- rnorm(N, mean=0, sd=1) yvar <- 1 + 2...

假设您想运行一个平分算法/二进制搜索，精确到机器精度，由于浮点数范围的指数减半，该算法/二进制搜索将始终在几百步内终止。为此，应满足以下条件：让浮点数 A<B 和 M = （A+B）/2。然后 当且...

有人告诉我，1e100 不能完全表示为 64 位浮点数。使用 64 位浮点数，我们可以在有效数中有 16 位有效数字，并操纵指数（最多约 10^308）。因此，由于 1e100 只有一个有效数字，为什...