---

awk 'NR==1 {print; next} {if ($1!="H") count++; $2=sin($2); $3=sin($3); $4=sin($4); print} END {print "Число атомов != H:", count}' atoms.xyz