update to provide standard errors

the routine now returns standard errors estimated from a linearized
model near the solution