Schoonschip, 68000 version of June 27, 1991. Public version. Date: Wed May 12 1993 18:05:55. Memory: start 00020008, length 476860. Command line: lf1 lf1.out C SMALL CHARGED SPHERE: Lorentz Force F(mu,nu)*U(nu) File: lf1 Input: tau1.o, tau2.o, tau3.o from dtau.z RdotA from RdotA.z denom2, denom3 from denoms.z Output: f1 in file lf1.z Starting date: November 8, 1987 Last revision: December 2, 1989 Seems correct on 2 Dec 89. Enter dtau.z ! Load tau1.o, tau2.o, tau3.o Enter denoms.z Enter RdotA.z Names tau1.o,tau2.o,tau3.o,RdotA,denom2,denom3 A dtau=4,eps=4 A c I mu,nu V u,a,j,n F D1y,D2y,D3y,D1u,D2u,D1a A norm Common tau1.o,tau2.o,tau3.o Common RdotA,denom2,denom3 X R(mu) = - D1y(mu,n)*eps - D2y(mu,n,n)*eps^2/2 - D3y(mu,n,n,n)*eps^3/6 - u(mu)*dtau - D1u(mu,n)*eps*dtau - D2u(mu,n,n)*eps^2*dtau/2 - a(mu)*dtau^2/2 - D1a(mu,n)*eps*dtau^2/2 - j(mu)*dtau^3/6 X U(mu) = u(mu) + D1u(mu,n)*eps + D2u(mu,n,n)*eps^2/2 + a(mu)*dtau + D1a(mu,n)*eps*dtau + j(mu)*dtau^2/2 X A(mu) = a(mu) + D1a(mu,n)*eps + j(mu)*dtau *fix Begin. Time 0 sec. A eps=3 ! eps^3 is leading nonsingular order in this section ! it won't be done right B eps X F(mu,nu) = { R(mu)*U(nu) - R(nu)*U(mu) } * d3 * ( c^2 - RdotA ) + { R(mu)*A(nu) - R(nu)*A(mu) } * d2 Z f = F(mu,u) ! Lorentz force density L 4 Id, uDu = c^2 L 5 Id, uDa = 0 L 6 Id, uDj = -aDa L 7 Id, D1u(u,l~) = 0 L 8 Id, D1a(u,l~) = -D1u(a,l) C Id, D2u(u,l~,m~) = -D1u(nu,l)*D1u(nu,m) L 9 Id, dtau = tau1.o*eps + tau2.o*eps^2 + tau3.o*eps^3 *yep L 1 Id, Commu, D1y,D2y,D3y,D1u,D2u,D1a *yep L 1 Id, d2 = denom2 *yep L 1 Id, Commu, D1y,D2y,D3y,D1u,D2u,D1a *yep L 1 Id, d3 = denom3 *yep L 1 Id, Commu, D1y,D2y,D3y,D1u,D2u,D1a L 2 Id, eps = 0 ! superfluous order from RdotA *yep L 1 Id, Addfa, eps ! ready to kill nonsingular order *yep L 1 Id, eps = 0 ! kill nonsingular order *yep L 1 Id, Addfa, eps^-1 ! restore eps powers L 2 Id, D1y(mu,n) = eta(mu,n) + D1y(u,n)*u(mu)*c^-2 *yep L 1 Id, uDu = c^2 Keep f *next f = + a(mu)*eps^-1 * ( - 1/2*c^-1*norm^-1 ) + D1y(Naa,n)*eta(mu,n)*D1y(u,n)*D1u(Naa,n)*eps^-1 * ( 3*c^-1*norm^-5 ) + D1y(Naa,n)*eta(mu,n)*D2y(Naa,n,n)*eps^-1 * ( - 3/2*c*norm^-5 ) + D1y(u,n)*D1y(u,n)*a(mu)*eps^-1 * ( - 1/2*c^-3*norm^-3 ) + D1y(u,n)*D1u(mu,n)*eps^-1 * ( c^-1*norm^-3 ) + D2y(mu,n,n)*eps^-1 * ( - 1/2*c*norm^-3 ) + D2y(u,n,n)*u(mu)*eps^-1 * ( 1/2*c^-1*norm^-3 ) + eta(mu,n)*eps^-2 * ( - c*norm^-3 ) + eta(mu,n)*D1y(u,n)*D1y(u,n)*D1y(a,n)*eps^-1 * ( - 3/2*c^-3*norm^-5 ) + eta(mu,n)*D1y(u,n)*D2y(u,n,n)*eps^-1 * ( 3/2*c^-1*norm^-5 ) + eta(mu,n)*D1y(a,n)*eps^-1 * ( 1/2*c^-1*norm^-3 ) + 0. Common f1 Z f1(mu) = f L 2 Id, eps^-2 = 0 B eps Keep f *next f1(mu) = + a(mu)*eps^-1 * ( - 1/2*c^-1*norm^-1 ) + D1y(Naa,n)*eta(mu,n)*D1y(u,n)*D1u(Naa,n)*eps^-1 * ( 3*c^-1*norm^-5 ) + D1y(Naa,n)*eta(mu,n)*D2y(Naa,n,n)*eps^-1 * ( - 3/2*c*norm^-5 ) + D1y(u,n)*D1y(u,n)*a(mu)*eps^-1 * ( - 1/2*c^-3*norm^-3 ) + D1y(u,n)*D1u(mu,n)*eps^-1 * ( c^-1*norm^-3 ) + D2y(mu,n,n)*eps^-1 * ( - 1/2*c*norm^-3 ) + D2y(u,n,n)*u(mu)*eps^-1 * ( 1/2*c^-1*norm^-3 ) + eta(mu,n)*D1y(u,n)*D1y(u,n)*D1y(a,n)*eps^-1 * ( - 3/2*c^-3*norm^-5 ) + eta(mu,n)*D1y(u,n)*D2y(u,n,n)*eps^-1 * ( 3/2*c^-1*norm^-5 ) + eta(mu,n)*D1y(a,n)*eps^-1 * ( 1/2*c^-1*norm^-3 ) + 0. Z f2(mu) = f L 2 Id, eps^-1 = 0 B eps *next f2(mu) = + eta(mu,n)*eps^-2 * ( - c*norm^-3 ) + 0. Delete tau1.o,tau2.o,tau3.o,denom2,denom3,RdotA Write lf1.z *end End run. Time 1 sec.