In [1]:

    

%run src/ssy_discretized_test.py


    

In [2]:

    

import matplotlib.pyplot as plt


    

In [3]:

    

ssy = SSY()


    

In [4]:

    

D=3


    

In [5]:

    

test_val_spec_rad(ssy, K=D, I=D, J=D)


    

    
Out[5]:





0.999439176569547

In [6]:

    

D_vals = range(3, 14)
Lambda_vals = []

for D in D_vals:
    print(f"Calculating at D = {D}")
    Λ = test_val_spec_rad(ssy, K=D, I=D, J=D)
    Lambda_vals.append(Λ)


    

    




Calculating at D = 3
Calculating at D = 4
Calculating at D = 5
Calculating at D = 6
Calculating at D = 7
Calculating at D = 8
Calculating at D = 9
Calculating at D = 10
Calculating at D = 11
Calculating at D = 12
Calculating at D = 13

In [7]:

    

Lambda_vals


    

    
Out[7]:





[0.999439176569547,
 0.999383593701018,
 0.9993297989794658,
 0.999271704762543,
 0.9992039375594387,
 0.9991200758421778,
 0.9990109367116795,
 0.9988588974239533,
 0.9986168173417717,
 0.9982106640625895,
 0.9976187931764713]

In [ ]: