plot[ exp( (.09-.03)*(c-1)*15 + (1 - c^2)*.22^2 * 15/2 ), c, 0, 5]