#include<iostream.h>
#include<conio.h>
void main(){
clrscr();
int nstuds[]={31,42,51,35,31};
float p=(45-40)/10;
float dy[5][5];
float sum=0;
for(int i=0;i<5;i++){
	sum=sum+nstuds[i];
	dy[0][i]=sum;
}
float yp=0;
for(i=1;i<5;i++){
	for(int j=0;j<5-i;j++){
		dy[i][j]=dy[i-1][j+1]-dy[i-1][j];
	}
}
//float fp=1,fi=1;
for(i=0;i<5;i++){
	float fp=1,fi=1;
	for(int j=0;j<1;j++){
		fp=fp*(p-(float)j);
		cout<<fp<<",";
		fi=fi*((float)(i-j));
		cout<<fi;
	}
	yp=yp+((fp/fi)*dy[i][0]);
}
cout<<"Number of students: "<<yp;
getch();
}