solar computation for hj andrews in matlab