#include #include #include typedef double ftype; typedef ftype ** MATRIX; typedef ftype * VECTOR; typedef const ftype ** CONST_MATRIX; typedef const ftype * CONST_VECTOR; #define INIT_a 0.0 #define INIT_b 100.0 #define INIT_c 75.0 #define INIT_d 50.0 #define ERROR 0.000001 MATRIX matrix_alloc(int n){ MATRIX x = (MATRIX)malloc(sizeof(ftype *)*n+sizeof(ftype)*n*n); int i; x[0]=(ftype *)(x+n); for(i=1;i0)sum-=in[i-1]; if(i0)sum-=in[i-1][j]; if(i0)sum+=in[i-1][j]; if(i0)sum+=in[i-1]; if(i=1)h[i][i-1]=1.0; h[i][i]=-1.0; if(iN){ is_minus=-1; index-=N+1; } if(index==0){ T[i][j]=0; }else{ T[i][j]=sn[index-1]*mul*is_minus; } } for(i=0;i=ERROR){ failed|=1; } failed|=2; y[0][i]=0.0; }else{ y[0][i]=-R[i]/v; } for(j=1;j=1)sum-=x[N-1][i-1]; if(i