sigma_from_rphi.c

最大似然估计算法

#include "timeseries.h"
                                                                                                                                     
void sigma_from_rphi(int ndim, double *phi, double *unit) {

int j, k;

for (j = 0; j < ndim; j++) {
	unit[j] = 1.0;
	for (k = 0; k < j; k++) unit[j] *= sin(phi[k]*M_PI/180.0);
	if (j < ndim-1) unit[j] *= cos(phi[j]*M_PI/180.0);
}
	
		
}