/* mul2_fix.mod - binary model with fixed relative position */

observation b1, b2, b3, b4, b5;
parameter   x1, y1, par, xdot, ydot, A1, x2, y2, A2;
data        obs, it, t, fx, fy, fp;
variable    mas = 4.8481369E-9, M1 = 0.71, M2 = 0.2485, p1, p2;
main()
{
    while ( import() ) {
        export(b1 - (A1+A2));
    }
    while ( import() ) {
        p1 = (fx*(x1+t*xdot) + fy*(y1+t*ydot) + fp*par)*mas;
        p2 = (fx*(x2+t*xdot) + fy*(y2+t*ydot) + fp*par)*mas;
        export(b2 - M1*(A1*cos(p1)+A2*cos(p2)));
        export(b3 - M1*(-A1*sin(p1)-A2*sin(p2)));
        export(b4 - M2*(A1*cos(2*p1)+A2*cos(2*p2)));
        export(b5 - M2*(-A1*sin(2*p1)-A2*sin(2*p2)));
    }
}