FUNCTION Errorcalc, ast_spec, star_spec, ast_err, star_err
sigmat = FltArr(2,335)
spart = (star_err[1,*] * star_err[1,*]) / (star_spec[1,*] * star_spec[1,*])
apart = (ast_err[1,*] * ast_err[1,*]) / (ast_spec[1,*] * ast_spec[1,*])
divpart = (ast_spec[1,*] / star_spec[1,*]) * (ast_spec[1,*] / star_spec[1,*])
fincalc = sqrt((apart + spart) * divpart)
sigmat[0,*] = star_spec[0,*]
sigmat[1,*] = fincalc
RETURN, sigmat
END