/*myComplex.h*/
void complexADD(double *result, double *a, double *b)
{
	*result = *a + *b ;
	*(result+1) = *(a+1) + *(b+1); 
}
void complexDEC(double *result, double *a, double *b)
{
	*result = *a - *b ;
	*(result+1) = *(a+1) - *(b+1); 
}
void complexMUL(double *result, double *a, double *b)
{
	*result = *a*(*b) - *(a+1)*(*(b+1)) ;
	*(result+1) = *(a+1)*(*b) + *a*(*(b+1)) ; 
}
void complexDIV(double *result, double *a, double *b)
{
	*result = ( *a*(*b) + *(a+1)*(*(b+1)) ) / ( *b*(*b) + *(b+1)*(*(b+1)) ) ;
	*(result+1) = ( *(a+1)*(*b) - *a*(*(b+1)) ) / (*b*(*b) + *(b+1)*(*(b+1))) ; 
}