1
import "//support.s-script" : s
21
//x.spec["differential"] = true
22
//y.spec["differential"] = true
23
//z.spec["differential"] = true
24
//Vx.spec["differential"] = true
25
//Vy.spec["differential"] = true
26
//Vz.spec["differential"] = true
32
nxt = Knx * (Vc - sqrt(Vx ^ 2 + Vy ^ 2 + Vz ^ 2))
37
Vyz = s.limit(Vyz, -15.0, 25.0)
39
nyt = 1 + Kny * (Vyz - Vy)
41
nzt = Knz * (Vx * (Zz - z) - Vz * (Xz - x))
42
nzt = s.limit(nzt, -5.0, 5.0)
47
__rights__ = fn(self : __state, float : t) -> object
56
__d["Vx"] = (g/sqrt((Vx)^2+(Vy)^2+(Vz)^2))*Vx*(nxt-Vy*nyt/sqrt((Vx)^2+(Vz)^2))-g*Vz*nzt/sqrt((Vx)^2+(Vz)^2)
57
__d["Vy"] = (g/sqrt((Vx)^2+(Vy)^2+(Vz)^2))*Vy*nxt+(g/sqrt((Vx)^2+(Vy)^2+(Vz)^2))*nyt*sqrt((Vx)^2+(Vz)^2)-g
58
__d["Vz"] = (g/sqrt((Vx)^2+(Vy)^2+(Vz)^2))*Vz*(nxt-Vy*nyt/sqrt((Vx)^2+(Vz)^2))+g*Vx*nzt/sqrt((Vx)^2+(Vz)^2)