I have a critical problem that I wish to solve it with maple
suppose we have a list like following: y_obs=(2,4,8,7,9,52,35,478,52) and corresponding variance σy=(.2,.3,.5,.87,.1.2,.22,.78,.99,1.5)
we know y as the function of x described such as y_theoric=x+p and minimizing X is
X=Sigma [(y_theoric-y_obs)^2]/σy which includes the sum of nine numbers...
the question is:
How we can find p from likelihood function and plot general behavior of y versus of x through two above series?
for example this solution used in article under the names Hubble parameter data constraints on dark energy by Yun Chen and Bhatra Ratra (Physics Letters B)