Assuming that the problem is continuous and not discrete (i.e. the producer pays half a dollar for 50 miles below the threshold), this question can be solved analytically (i.e. a formula based solution, like you said) by explicitly writing the distribution and integral.
The distribution of mileage is normal with known mean and variance, let's write $ f(x) = {1 \over \sqrt{2\pi}\sigma} exp(- {(x-\mu)^2 \over {2 \sigma^2}}) $
The cost function is 1 dollar per $m$ miles below some threshold $ x_0 $, so we can write $ c(x) = max( {{x_0-x} \over m} ,0) $.
The expected cost is the integral $ \int _{-\infty} ^{\infty} f(x)c(x)dx $. However, since the cost is 0 if $x$ is above $x_0$, we can simply write it as $ \int _{-\infty} ^{x_0} {{(x-x_0) \over m}} {1 \over \sqrt{2\pi}\sigma} exp(- {(x-\mu)^2 \over {2 \sigma^2}}) dx $.
From here, the solution is relatively simple.
Even if the problem is discrete (i.e. the producer pays only in discrete amount of dollars), you can verify your discrete result (e.g. by Viktor's answer) by computing the continuous result and see that they are not too far.
dnormdoes not return probabilities. It is a probability density function. – whuber Dec 04 '19 at 15:20