Need staticstical analysis my model "Dark Energy" with DESI dataset

staticstical analysis my "Dark Energy" with DESI Dataset 1- Write a Python program that: reads the published DESI BAO mesurements, -Computes H(z) and D_A(z) - minimizes chi^2 -compares my model with lambdaCDM - plots the residuals and confidence regions Compute AIC and BIC to quantify whether the extra model complexity is justified.