Number fields

See nf_in_gp

Quaternion algebra

F=nfinit(y^2-2)
A = alginit(F, [-1,-1]);