return C;
}
array plnpow(double X,int mpow){//stvarennie masiva stupiniej X da 'mpow' uklucxna
array tmp(mpow+1);
double pw=1.0;
for(unsigned int i=0;i<tmp.n;i++){
tmp.V[i]=pw;
pw=pw*X;
}
return tmp;
}
double plnigl01(array A){//integravannie palinoma u miezxach ad 0 da 1
double cum=0.0;
for(unsigned long int i=0;i<A.n;i++)cum+=A.V[i]/((double)i+1.0);
return cum;
}
matrix plndif(matrix A){//diferencavannie niekalkich palinomaux zapisanych slupkami u matrycu
matrix B;
B.m=A.m -1;B.n=A.n;
B.M=new array[B.m];
for(unsigned long int i=0;i<B.m;i++)B.M[i](B.n);
for(i=0;i<B.m;i++)
for(unsigned long int j=0;j<B.n;j++)
B.M[i].V[j]=((double)i+1.0)*A.M[i+1].V[j];
return B;
}
matrix svior(matrix A){//vylicxenni "sviortki" matrycy (integrala palinaminalnaha pieramnazxennia slupkoux)
matrix B;
B.m=A.n;B.n=A.n;
B.M=new array[B.m];
for(unsigned long int i=0;i<B.m;i++)B.M[i](B.n);
for(i=0;i<B.m;i++)
for(unsigned long int j=0;j<B.n;j++)
B.M[i].V[j]=plnigl01(plnprd(vector(A,i),vector(A,j)));
return B;
}
matrix svior(matrix A, unsigned long int lCl, unsigned long int rCl){//vylicxenni "sviortki" matrycy (integrala palinaminalnaha pieramnazxennia slupkoux), ale nie usich slupkoux, a ad 'lCl' da 'rCl'
matrix B;
B.m=rCl-lCl+1;B.n=B.m;
B.M=new array[B.m];
for(unsigned long int i=0;i<B.m;i++)B.M[i](B.n);
for(i=0;i<B.m;i++)
for(unsigned long int j=0;j<B.n;j++)
B.M[i].V[j]=plnigl01(plnprd(vector(A,lCl+i),vector(A,lCl+j)));
return B;
}
Распечатка сладов презентации
| | |
| | |
| | |
| | |
| | |
| | |
| | |