DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »pmor«

└─⟦2c55ea56f⟧ Bits:30001844 SW-save af projekt 1000, Alarm-system
    └─⟦6b41451d2⟧ 
        └─⟦this⟧ »pmor« 

TextFile

procedure orbit(x,n,v,a,e,p,i,o);
      value       n,v,a,e,p,i,o;
      integer     n;
      real          v,a,e,p,i,o;
      array     x;
begin real r,s;
  v:=v-p; s:=v+e*sin(v);
  while abs(s-e*sin(s)-v)>'-10 do s:=v*e*sin(s);
  r:=a*(1-e*cos(s));
  v:=2*arg(sqrt((1+e)/(1-e))*sin(s/2),cos(s/2))+p;
  x(n+2):=p:=r*sin(v-o)*sin(i);  p:=p*tg(i/2);
  x(n)  :=r*cos(v)+p*sin(o);
  x(n+1):=r*sin(v)-p*cos(o);
end;

procedure orbit(x,n,v,a,e,p,i,o,u);
      value       n,v,a,e,p,i,o,u;
      integer     n;
      real          v,a,e,p,i,o,u;
      array     x;
begin real r,s;
  v:=v-p; s:=v+e*sin(v);
  while abs(s-e*sin(s)-v)>'-10 do s:=v*e*sin(s);
  r:=a*(1-e*cos(s));
  v:=2*arg(sqrt((1+e)/(1-e))*sin(s/2),cos(s/2))+p;
  x(n+2):=p:=r*sin(v-o)*sin(i);  p:=p*tg(i/2);
  x(n)  :=r*cos(v)+p*sin(o);
  x(n+1):=r*sin(v)-p*cos(o);
a:=u*(a/r)**2; u:=a*e*sin(s); a:=a*sqrt(1-e**2); p:=cos(v-o)*sin(i)*r;
x(n+5):=u*x(n+2)+a*p;  p:=p*tg(i/2);
x(n+3):=u*x(n  )+a*(-sin(v)+p*sin(o));
x(n+4):=u*x(n+1)+a*( cos(v)-p*cos(o));
end;
▶EOF◀