#/usr/bin/env python
# encoding: utf-8

from __furure__ import division, print_function
import numpy as np
from scipy.integrate import quad
from scipy.optimize import fixed_point
from math import exp, sin, cos

def F(m, r, D, beta, J):
    a = (1+D)/2
    b = (1-D)/2
    X = np.vstack([[1, J],[J, 1]])
    M = (a*X).dot(m)
    R = (b*X).dot(r)


def Fb():
    pass

def Za():
pass

def Zb():
    pass

def ma():
    pass

def mb():
    pass

def m2a():
    pass

def m2b():
    pass

def ra():
    pass

def rb():
    pass

def mr():
    pass

def
