The wording of this question confused me a little. In particular that we want the fraction of unfolded at different values of external force. I interpret this just as the probability of being in the unfolded state because you are either in the unfolded state or the folded state. Thus the probability of the unfolded state gives you the fraction in the folded state compared to the number in either state so the fraction in the folded state. This is the first graph below. That being said I went to office hours and asked about this and was told that you wanted the ratio of the probability of the unfolded state over the probability of the folded state. Thus I have provided this as the second graph.
/shared-libs/python3.7/py-core/lib/python3.7/site-packages/ipykernel_launcher.py:6: RuntimeWarning: divide by zero encountered in log
below K=1 we start to get negative entropy.
b) We get a non physical result because we let q_tran be less than 1. This doesn't make any sense because this says that there are no accessible states for the system to be in however it has to be in some state so q_tran should never be lower than 1. One reason this error came up is because we assumed that the temperature was high enough so we could approximate the summation in the partition function by an integral. Here however we let the temperature get arbitrarily low and thus this approximation makes no sense.
/shared-libs/python3.7/py-core/lib/python3.7/site-packages/ipykernel_launcher.py:7: RuntimeWarning: divide by zero encountered in true_divide import sys