The formulas above will be slightly different to charts as charts are projected onto a spheriod, normally the WGS84 ellipsopid.
The formula results above will agree directly with unprojected data, or an unprojected map.
The WGS84 ellipsoid is flatter at the poles than a perfect sphere, hence the spherical geometry answers you received above will have a small error as a result of the projection.