Даны действительные числа S, T, A0,.,A12. Получить P(X)= P(1) − P(T) + P^2 (S −T) − P^3(1) Где P(X) = A12X^12 + A11X^11 + . + A0 - C (СИ)

Формулировка задачи:

Пожалуйста посмотрите код, посоветуйте где что исправить, заранее спасибо) Даны действительные числа S, T, A0,..,A12. Получить P(X)= P(1) − P(T) + P^2 (S −T) − P^3(1) Где P(X) = A12X^12 + A11X^11 + ... + A0 (A0,..,A12 - это не степень, а порядок. X^12..X^1 - степень) Вот сам код:

#include<stdio.h>
#include<conio.h>
#include<math.h>
int a[13];
double p1 (double x, int i){
double rez=1; int j;
for (j=1;j<=i;j++) 
     rez*=x;
}       
double p (double x1)
       {
       double q=0; int j;
       for (j=1;j=12;j++)           
           {q+=a[j]*p1(x1,j);} 
       q+=a[0];
       }
main() 
{
int a[13],i; double rp,s,t;
for(int i=0;i<=12;i++)
         {
          printf("Vvedite elem. massiva A");
          scanf("%f",&a);
         }
printf("Vvedite s ");
scanf("%f",&s);
printf("Vvedite t ");
scanf("%f",&t);
rp=p*1-p*t+pow((p,2)*(s-t))-pow(p,3)*1;
printf("Otvet = ");
scanf("%f",&pr);
}

Уже несколько дней сижу над ней, голова совсем не варит, если увидете ошибки, скажите о них)

Решение задачи: «Даны действительные числа S, T, A0,.,A12. Получить P(X)= P(1) − P(T) + P^2 (S −T) − P^3(1) Где P(X) = A12X^12 + A11X^11 + . + A0»

textual
                            Листинг программы 
                        
                            rp=p*1-p*t+pow((p,2)*(s-t))-pow(p,3)*1;