How to graphically represent this function?

[scidam]

You need to compute values of your function for all values in S. Just replace line #28 in your code with the following:

bc_P = [CON_call_price(120,100, 1, 0.25, 0.2, 0.03, 1, s) for s in S]