这是一种快速且非迭代的椭圆拟合。 用法： A = EllipseDirectFit(XY) 输入：XY(n,2)是n个点的坐标数组x(i)=XY(i,1), y(i)=XY(i,2) 输出：A = [abcdef]' 是系数向量最佳拟合椭圆的方程： ax^2 + bxy + cy^2 + dx + ey + f = 0， 要将此向量转换为几何参数（半轴、中心等），请使用标准公式，例如 Wolfram Mathworld 中的 (19) - (24)： http://mathworld.wolfram.com/Ellipse.html 这种椭圆拟合是在文章中提出的AW Fitzgibbon, M. Pilu, RB Fisher “椭圆的直接最小二乘拟合” IEEE 翻译帕米，卷。 21，第 476-480 页（1999 年） 作者将其称为“直接椭圆拟合”。 我的代码基于数

在Dyson-Schwinger方程的框架内，并通过多次反射扩展，我们研究了有限体积对球形手性相变的影响，并特别讨论了其对临界终点（CEP）可能位置的影响 ）。 根据我们的计算，当我们使用球体而不是立方体时，有限体积对相变的影响并不像先前计算的那么重要。 例如，随着球形体积的半径从无限减小到2 fm，临界温度 在零化学势和有限温度下，仅轻微下降。 在有限的化学势和有限的温度下，CEP的位置朝着较小的温度和较高的化学势移动，但变化幅度不超过20％。 结果，我们发现不仅体积的大小，而且体积的形状对相变都有很大的影响。

在本文中，我们找到了具有广义多态状态方程（GPEoS）的Einstein–Maxwell方程的精确解。 为此，我们考虑具有带电各向异性物质分布的球对称物体。 我们通过Durgapal（Phys Rev D 27：328，1983）引入的变换将场方程重写为简单形式，然后解析求解这些方程。 对于这些解决方案的物理可接受性，我们绘制了物理量，例如能量密度，各向异性，声速，切向和径向压力。 我们发现所有解决方案均满足所需的物理条件。 结论是，我们所有的结果都简化为带有线性，二次态和多态状态方程的各向异性带电物质分布的情况。

AdS 4中较高自旋电流相互作用的形式是从Weyl 0形式扇形中的完全非线性较高自旋方程导出的。 明确确定了由标量和自旋轴建立的自旋一电流前的耦合常数以及Yukawa耦合。 所有其他高自旋电流相互作用的耦合均被隐式确定。 显示所有耦合都独立于非线性高自旋理论的相位参数。 顶点对较高自旋相位参数的适当全息依赖性显示为整体场上的边界条件所致。

代码主要为四阶龙格库塔求解四自由度动力学模型，可无缝衔接时变刚度的导入以及后续振动加速度、位移的提取，可出相图等非线性结果。稍加修改方程即可完成简单的六自由度动力学模型求解。主要适用于刚学习，齿轮动力学同学。

热传导 差分方程

拉普拉斯方程数学代码跌落模拟器 该存储库包含与本文相关的代码： Liimatainen等。 “在生物和合成憎水表面上绘制微尺度的润湿变化”，《自然通讯》 （2017年），1798年。 doi：10.1038 / s41467-017-01510-7 请参考该论文引用代码。 它是什么？ 该存储库包含用于计算两个表面之间的液滴的粘附力的代码。 解决的问题是轴对称（2D）。 液滴满足Young-Laplace方程。 降落的体积和边界条件也应该是已知的。 该力是拉普拉斯项和毛细管项的总和。 使用射击方法解决了边值问题。 有关所有详细信息，请参见本文的补充说明。 该代码使用Matlab R2016b编写并经过测试，并使用Simulink和Parallel计算工具箱。 如果需要，可以通过用for循环替换所有parfor循环，非常轻松地消除对并行计算工具箱的依赖，但代价是代码的运行速度严重降低。 Simulink是solveIVP.m代码不可或缺的一部分，但应该可以仅使用ode45来实现相同的行为。 如何从纸上重新创建数字 从下载代码 运行脚本（代码需要对/output目录的

自洽-肖丁格-泊松 二维薛定谔-泊松方程的自洽解

吉尔数值算法，求解刚性微分方程组，稳定好，收敛好

运用能量估计的方法,在临界Sobolev空间H1/2(R3)中,研究了三维不可压磁微极流体方程小初值整体强解的渐进性质.设(u,ω,b)是三维不可压磁微极流体方程在临界Sobolev空间H1/2(R3)中小初值(u0,ω0,b0)∈H1/2(R3)对应的整体强解,那么解的H1/2(R3)范数‖u,ω,b‖H1/2关于时间t是非增函数,且当t→+∞时,极限为0;并且使得整体强解(u,ω,b)存在的小初值(u0,ω0,b0)构成的集合是空间H1/2(R3)中的开集.