I have the following formula for KummerU:
U = gamma(1-b)*KummerMz(a,b,z)/gamma(a-b+1) + gamma(b-1)*z.^(1-b).*KummerMz(a-b+1,2-b,z)/gamma(a)
from [https://pmc.polytechnique.fr/pagesperso/dg/confluent/KummerU.m][1]
I am generating it in Excel using a Gamma and KummerM function. But I do not get the desired score of .953742 for a=.25, b=.25, z=.5.
Any input to where the calculation is going wrong?