Theory_Transition_Saturation_Curves.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Wed Aug 25 10:32:58 2021

@author: oem
"""


import numpy as np
import matplotlib.pyplot as plt

def SaturationCurve(S, gamma, det):
    return 0.5*S*(1/(1+S+4*((det**2)/(gamma**2))))

gamma = 2*np.pi*1e6*22
det = 2*np.pi*1e6*10

Svec = np.arange(0, 5, 0.01)

SatCurve = np.array([SaturationCurve(s, gamma, det) for s in Svec])

plt.plot(Svec, SatCurve)