使用此标记来解决有关在科学中使用计算机的问题其中许多问题也适用于scicomp.stackexchange.com。
我正在使用 Jax 进行科学计算，具体来说，并行计算系统中所有粒子的成对相互作用力。这是为了模拟细菌的动力学......
问题描述： 我有与特定 conrer_ids 元素相对应的值。从这些值的多组中，我需要找到最大的值并生成包含 conrer_ids 的数组，
我正在尝试模拟基于菲克第二定律的简单扩散。 从 pylab 导入 * 将 numpy 导入为 np 网格点 = 128 定义轮廓（x）： 范围 = 2。 散乱=.1576 剂量=1 重新...
Solve_ivp 是来自 Scipy 的初始值问题求解器函数。 简单来说 scipy.integrate.solve_ivp（乐趣，t_span，y0，方法='RK45'，t_eval=无，dense_output=False，事件=无，矢量化=...
大家早上好，我有一个关于用 matplotlib 制作等高线图的问题。 如果我尝试用两个不同范围的变量绘制相同的函数（第二个包含另一个......
我正在尝试设置一些 YALMIP 约束，如下所示： 对于 i=1:q 约束 = [约束，trace(As(:, :, i)) == 1] 结尾 所以我本质上只是想断言所有（方形）ma...
显然，从 ODE 求解器获得非负解并非易事。在 Matlab 中，某些求解器可以使用 NonNegative 选项来获得非负解。有没有类似的选择...
我有一个三角信号（或更准确地说）是一个在一定区间（0, 20）内呈三角函数的函数（我不关心输出或输入信号在其他范围内的情况）。我想写...
我正在尝试编写一个单元测试来检查工程分析的输出。我有理论值，我想根据对一定数量的有效数字的分析来检查这些值。所以，...
我正在用厄米特（复）正半定矩阵进行一些计算。计算涉及取它们的指数。矩阵 A 的指数（以 2 为底），我将表示...
当 Const=False 和数字时，EXCEL LinEst 函数返回无意义 <1E-7
要重现我的问题，只需输入 A1: 9.1E-8 和 A2: 9.2E-8，然后标记这两个单元格并向下拖动（自动填充）到 A25。然后输入B1：“=1.1*A1”并将公式复制下来，...
我一直在阅读这篇文章《绘制正确分叉图的简单指南》，我想重现下图（图 10，第 2150011-7 页）： 我创建了这个程序，...
我现在正在研究两个多维数组 arr 和 cost。 arr 与尺寸 (width, height, m, n) 是密集的，而 cost 是与尺寸 (width, height, width, height) 稀疏的。 值：宽度和高度是...
V (cm3) T (°C) 1.032 10.0 1.094 29.5 1.156 50.0 1.215 69.5 1.273 90.0 编写一个 python 程序，通过假设求出系数 a 和 b 该液体体积随关系式 V = 1 + a...
我已经开始使用 valgrind，特别是用于代码分析的工具“callgrind”。 首先，我创建了一个用于数值积分的共享库（我已经实现了一个
我是 Rust 新手，我正在尝试使用 Rust 对多个矩阵执行顺序克罗内克积，即我想做 在Python中，我知道我能做到 将 numpy 导入为 np 从 functools 导入减少...
编码二维阵列中多个谐振子的四阶龙格-库塔方法。陷入单一振荡器
我之前在 Python 中使用四阶龙格库塔方法对一般（非线性）摆摆进行了数值求解。这是与此相关的问题，尝试求解第二个非线性摆......
我参考这个问题，已经有很好的答案了；但发现了不必要的操作（请参阅帖子中的讨论），我只是好奇我是否能成功消除这些...
尝试使用四阶龙格-库塔法求解非线性摆二阶微分方程，未得到预期结果
我是Python新手，所以我只将过程限制在matplotlib上，而不是进入NumPy。 我按照 Abhijit Kar Gupta 的书《Python 中的科学计算》，编写了 Python
我有一个多值函数，我想将其与 np.linspace 结合使用。每当函数中有 if/else 条件时，我都会收到值错误。 我期望应用一个函数...