读书笔记2-三体

弱小和无知不是生存的障碍,傲慢才是。

给岁月以文明,而不是给文明以岁月。

失去人性失去很多,失去兽性失去一切

毁灭你,与你何干?

黑暗森林中有一个叫人类的傻孩子,生了一堆火并在旁边高喊:“我在这儿,我在这儿!”——罗辑

宇宙很大,但生命更大。

新的文明在诞生,新的道德也在形成。

没有永恒的敌人或同志,只有永恒的责任。——章北海

你们是虫子。

把人类看做虫子的三体人似乎忘记了一个事实:虫子从来就没有被真正战胜过。

前进!前进!不择手段的前进!——维德

可怜虫们,准备去澳大利亚吧。

像坟墓一样简洁。

因为宇宙不是童话。

如果有墓碑,上面应该写:来了,爱了,给了她一颗星星,走了。——写云天明

孩子问,他们是烈士吗?
妈妈说,不是。
他们是敌人吗?
不是
那他们是什么?
他们是历史!

死亡是唯一一座永远亮着的灯塔,不管你向哪里航行,最终都得转向它指引的方向。一切都会逝去,只有死神永生。

唯一不可阻挡的是时间,它像一把利刃,无声的切开了坚硬和柔软的一切,恒定的向前推进着,没有任何东西能够使它的行进产生丝毫的颠簸,它确改变着一切。

而人类,已经无人能救了。

我需要一块二向箔,清理用。——歌者

从科学角度讲,毁灭一词并不准确,没有真正毁掉什么,更没有灭掉什么,物质总量一点不少都还在,角动量也还在,只是物质的组合方式变了变,像一副扑克牌,仅仅重洗而已......可生命是一手同花顺,一洗什么都没了。

我看到了我的爱恋
我飞到她的身边
我捧出给她的礼物
那是一小块凝固的时间
时间上有美丽的条纹
摸起来像浅海的泥一样柔软
她把时间涂满全身
然后拉起我飞向存在的边缘
这是灵态的飞行
我们眼中的星星像幽灵
星星眼中的我们也像幽灵
——歌者

 

转载于:https://www.cnblogs.com/zhulvbo/p/8746284.html

运动是指个物之间的相互作用和运动,可以使用Python进行模拟。其中,可以使用scipy库中的odeint函数求解微分方程组,来模拟运动。 以下是一个简单的运动模拟代码示例: ```python import numpy as np from scipy.integrate import odeint # 定义微分方程组 def three_body_equations(w, t, G, m1, m2, m3): x1, y1, vx1, vy1, x2, y2, vx2, vy2, x3, y3, vx3, vy3 = w r12 = np.sqrt((x2 - x1)**2 + (y2 - y1)**2) r13 = np.sqrt((x3 - x1)**2 + (y3 - y1)**2) r23 = np.sqrt((x3 - x2)**2 + (y3 - y2)**2) dx1dt = vx1 dy1dt = vy1 dvx1dt = G * m2 * (x2 - x1) / r12**3 + G * m3 * (x3 - x1) / r13**3 dvy1dt = G * m2 * (y2 - y1) / r12**3 + G * m3 * (y3 - y1) / r13**3 dx2dt = vx2 dy2dt = vy2 dvx2dt = G * m1 * (x1 - x2) / r12**3 + G * m3 * (x3 - x2) / r23**3 dvy2dt = G * m1 * (y1 - y2) / r12**3 + G * m3 * (y3 - y2) / r23**3 dx3dt = vx3 dy3dt = vy3 dvx3dt = G * m1 * (x1 - x3) / r13**3 + G * m2 * (x2 - x3) / r23**3 dvy3dt = G * m1 * (y1 - y3) / r13**3 + G * m2 * (y2 - y3) / r23**3 return dx1dt, dy1dt, dvx1dt, dvy1dt, dx2dt, dy2dt, dvx2dt, dvy2dt, dx3dt, dy3dt, dvx3dt, dvy3dt # 定义初始状态和参数 w0 = [1, 0, 0, 6, -1, 0, 0, -6, 0, 0, 0, 0] t = np.linspace(0, 10, 1000) G = 1 m1 = 1 m2 = 1 m3 = 1 # 求解微分方程组 wsol = odeint(three_body_equations, w0, t, args=(G, m1, m2, m3)) # 绘制轨迹图 import matplotlib.pyplot as plt plt.plot(wsol[:,0], wsol[:,1], label='Body 1') plt.plot(wsol[:,4], wsol[:,5], label='Body 2') plt.plot(wsol[:,8], wsol[:,9], label='Body 3') plt.legend() plt.show() ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值