compute lat/long given lat/long plus distance and bearing

Has anyone got a formula that:
given a lat/long starting point plus a distance and bearing
will output the new lat/long position?

I found this but it doesnt seem to work right

double curLat=40.12345;
double curLong=-122.12345;
double bearing=90; // east
double dist=1609.344; // 1 mile in meters

d_lat = asin(sin(curLat)* cos(dist)+ cos(curLat)* sin(dist)*cos(bearing));
d_lon = atan2(sin(bearing)*sin(dist)*cos(curLat),
cos(dist)-sin(altitude)*sin(curLat));
d_long = fmod( curLong - d_lon + pi, 2*pi )- pi;

printf("new coordinates: %lg x %lg\n", d_lat, d_long);

my results look like this:
new coordinates: -0.384715 x -6.1364

Thanks
Eric