Question

error C2297: '*' : illegal, right operand has type 'double *'

#include <cstdlib>
#include <iostream>
#include <cmath>
/*
利用Levinson递推演算法求解n阶Toeplitz型的联立方程式
*/
using namespace std;
int tlvs(double *t,int n,double *b,double *x)
//int n;
//double t[],b[],x[];
{ int i,j,k;
double a,beta,q,c,h,*y,*s;
s=(double *)malloc(n*sizeof(double));
y=(double *)malloc(n*sizeof(double));
a=t[0];
if (fabs(a)+1.0==1.0)
{ free(s); free(y); printf("fail\n"); return(-1);}
y[0]=1.0; x[0]=b[0]/a;
for (k=1; k<=n-1; k++)
{ beta=0.0; q=0.0;
for (j=0; j<=k-1; j++)
{ beta=beta+y[j]*t[j+1];
q=q+x[j]*t[k-j];
}
if (fabs(a)+1.0==1.0)
{ free(s); free(y); printf("fail\n"); return(-1);}
c=-beta/a; s[0]=c*y[k-1]; y[k]=y[k-1];
if (k!=1)
for (i=1; i<=k-1; i++)
*s=y[i-1]+c*y[k-i-1];
a=a+c*beta;
if (fabs(a)+1.0==1.0)
{ free(s); free(y); printf("fail\n"); return(-1);}
h=(b[k]-q)/a;
for (i=0; i<=k-1; i++)
{ x=x+h*s; y=s;}
x[k]=h*y[k];
}
free(s); free(y);
return(1);
}
int main(int argc, char *argv[])
{
int i;
static double x[6];
static double t[6]={6.0,5.0,4.0,3.0,2.0,1.0};
static double b[6]={11.0,9.0,9.0,9.0,13.0,17.0};
if (tlvs(t,6,b,x)>0)
for (i=0; i<=5; i++)
cout<<"x("<<i<<")="<<x<<"\n";
system("PAUSE");
return EXIT_SUCCESS;
}

Answer 1

{ x=x+h*s; y=s;}
s是个指针，参与*运算是什么意思？

追问

怎么改？