數學史上你認為最醜陋的公式是什麼？

12-27

姊妹篇：數學史上你認為最美的公式是什麼？

[1]就瘮人來說的話,我認為是馬丁迭代公式

該公式生成的圖像看得我起雞皮疙瘩.尤其是馬丁迭代是一個生長的過程,它會一圈一圈地變粗.

[2]就邪惡來說的話,我認為是"骷髏五星"(這是我起的名字,如果有人知道它真正的名字請告知),它的公式我不會描述,只能用語言來說下它的生成方式:

生成正五邊形的5個頂點，和一個當前點設置為原點；

a.隨機選擇五邊形的某一個頂點，計算出它與當前點的中點位置；
b.將計算出的中點做為當前點，重新執行操作a
由此迭代處理上千萬次後，即能生成骷髏五星的像素數據。

[3]就混亂來說的話,我認為是"N體",每一個星體的運動軌跡都是一團亂麻.
下圖為一個運動物體在六個固定物體的引力作用下的運動軌跡:

該圖像的生成公式是我寫的,也是一團亂麻:

p0x=[static]r*sin(0.0)
p1x=[static]r*sin(PI/3)
p2x=[static]r*sin(PI*2/3)
p3x=[static]r*sin(PI*3/3)
p4x=[static]r*sin(PI*4/3)
p5x=[static]r*sin(PI*5/3)
p0y=[static]r*cos(0.0)
p1y=[static]r*cos(PI/3)
p2y=[static]r*cos(PI*2/3)
p3y=[static]r*cos(PI*3/3)
p4y=[static]r*cos(PI*4/3)
p5y=[static]r*cos(PI*5/3)
m0=[static]1.0
m1=[static]1.0
m2=[static]1.0
m3=[static]1.0
m4=[static]1.0
m5=[static]1.0
k=[static]0.5*w*(u*u + v*v)
d0x=[static]p0x - x
d0y=[static]p0y - y
q0=[static]-g*m0*w/sqrt(d0x*d0x + d0y*d0y)
d1x=[static]p1x - x
d1y=[static]p1y - y
q1=[static]-g*m1*w/sqrt(d1x*d1x + d1y*d1y)
d2x=[static]p2x - x
d2y=[static]p2y - y
q2=[static]-g*m2*w/sqrt(d2x*d2x + d2y*d2y)
d3x=[static]p3x - x
d3y=[static]p3y - y
q3=[static]-g*m3*w/sqrt(d3x*d3x + d3y*d3y)
d4x=[static]p4x - x
d4y=[static]p4y - y
q4=[static]-g*m4*w/sqrt(d4x*d4x + d4y*d4y)
d5x=[static]p5x - x
d5y=[static]p5y - y
q5=[static]-g*m5*w/sqrt(d5x*d5x + d5y*d5y)
e=[static]k + q0 + q1 + q2 + q3 + q4 + q5

f0=(p0x - x)*(p0x - x) + (p0y - y)*(p0y - y)
r0=g*m0/f0
f1=(p1x - x)*(p1x - x) + (p1y - y)*(p1y - y)
r1=g*m1/f1
f2=(p2x - x)*(p2x - x) + (p2y - y)*(p2y - y)
r2=g*m2/f2
f3=(p3x - x)*(p3x - x) + (p3y - y)*(p3y - y)
r3=g*m3/f3
f4=(p4x - x)*(p4x - x) + (p4y - y)*(p4y - y)
r4=g*m4/f4
f5=(p5x - x)*(p5x - x) + (p5y - y)*(p5y - y)
r5=g*m5/f5
i=r0*(p0x - x)/sqrt(f0)
j=r0*(p0y - y)/sqrt(f0)
i=i + r1*(p1x - x)/sqrt(f1)
j=j + r1*(p1y - y)/sqrt(f1)
i=i + r2*(p2x - x)/sqrt(f2)
j=j + r2*(p2y - y)/sqrt(f2)
i=i + r3*(p3x - x)/sqrt(f3)
j=j + r3*(p3y - y)/sqrt(f3)
i=i + r4*(p4x - x)/sqrt(f4)
j=j + r4*(p4y - y)/sqrt(f4)
i=i + r5*(p5x - x)/sqrt(f5)
j=j + r5*(p5y - y)/sqrt(f5)
x=x + u*t + 0.5*i*t*t
y=y + v*t + 0.5*j*t*t
u=u + i*t
v=v + j*t
d0x=p0x - x
d0y=p0y - y
q0=-g*m0*w/sqrt(d0x*d0x + d0y*d0y)
d1x=p1x - x
d1y=p1y - y
q1=-g*m1*w/sqrt(d1x*d1x + d1y*d1y)
d2x=p2x - x
d2y=p2y - y
q2=-g*m2*w/sqrt(d2x*d2x + d2y*d2y)
d3x=p3x - x
d3y=p3y - y
q3=-g*m3*w/sqrt(d3x*d3x + d3y*d3y)
d4x=p4x - x
d4y=p4y - y
q4=-g*m4*w/sqrt(d4x*d4x + d4y*d4y)
d5x=p5x - x
d5y=p5y - y
q5=-g*m5*w/sqrt(d5x*d5x + d5y*d5y)
k=e - q0 - q1 - q2 - q3 - q4 - q5
k=if(k &< 0.0, 0, k)
s=sqrt(2*k/w)/sqrt(u*u + v*v)
u=u*s
v=v*s