Здравствуйте, Erop, Вы писали:
E>1) Сообщу на всякий случай, что double *x и double x[3] в описании ПАРАМЕТРОВ функции ОБОЗНАЧАЮТ ОДНО И ТО ЖЕ
E>Так что лучше писать ПОНЯТНЕЕ
E>E>double dot_product(const double x[], const double y[])
E>{
E> return x[0]*y[0] + x[1]*y[1] + x[2]*y[2];
E>}
E>
Вот она —
красота
/* Subroutine int dgesvd_(char *jobu, char *jobvt, integer *m, integer *n,
doublereal *a, integer *lda, doublereal *s, doublereal *u, integer *
ldu, doublereal *vt, integer *ldvt, doublereal *work, integer *lwork,
integer *info)
*/
dgesvd_( &JOBU, &JOBVT, &M, &N, a, &LDA, s, uu,
&LDU, vt, &LDVT, wk, &LWORK, &INFO);