Fix special functions on large intervals

tidy linsolves in specialfunctions (#632)

it's better to change the relative tolerance than to scale and rescale (because this avoids dividing by 0 and then multiplying by 0 again)