#!/usr/bin/env python
# coding: utf-8

# In[ ]:


import numpy as np

def Glycoscillator():
    IC_inside = [ 0.5325, 0.4767, 0.0618, 0.3173, 0.2075, 1.5311, 0.0572]
    IC_limit = [1.21440084, 0.27349879, 0.04031784, 0.05089061, 0.14671042, 1.88564678, 0.00505542]
    IC_outside = [ 0.4820, 1.9893, 0.0644, 0.3065, 0.1984, 2.6602, 0.0539]
    IC_1 = [ 0.875, 1.175, 0.12, 0.225, 0.19, 1.405, 0.075] 
    
    IC_test = [3.5, 2, 1, 1.5, 2.3, 1.8, 1.4]

    IC = np.transpose([IC_inside, IC_limit, IC_outside, IC_1, IC_test])
    return IC

def Repressilator():
    IC_inside = [3.79156715, 5.65839438, 4.45129228, 3.87443063, 5.07860964, 4.89017854]
    IC_limit = [ 0.48708115, 4.65643912, 60.66488088, 2.75453539,  1.88284155, 71.56193699]
    IC_outside = [200, 0, -50, 120, 0, -10]
    IC_1 = [40, 60, 70, 50, 55, 65]

    IC_test = [80, 40, 20, 25, 50, 60]

    IC = np.transpose([IC_inside, IC_limit, IC_outside, IC_1, IC_test])
    return IC

