code finaly worked, and give results consistent with Ingrid's.