Некорректно работает прогрмма - C (СИ)

#include <stdio.h>
#include <math.h>
 
#define a1=21.333;
 
float Pvalue(float a, float b, float c)
      {
 
      if(c<b)
      b=c;
      return a+b;
 
      }
 
int main()
{
 
float a, b, c, h, proiz, deltak, p, x, ak, xo;
int m, n, k, i;
 
printf("vvedite a,b,c, x0, m, h, n\n");
scanf("%f%f%f%f%d%f%d", &a, &b, &c, &xo, &m, &h, &n);
 
p = Pvalue(a,b,c);
 
printf("%f\n", p);
for (i=1; i<=m; i++)
{
    x=xo+(i+1)*h;
     proiz = 1;
    printf("\nPosledovatelnost #%d\n", i);
    for (k=2; k<=n; k++)
        {
         deltak=(p*p*sin(x+3)*sin(x+3)*(k+1))/(k*k+2*k);
         ak=ak*deltak;
         if (k%2==0)
         {
         proiz=proiz+ak;
         }
         printf("a%d = %f\n", k,ak);
         }
         printf("\nProizvedenie chetnih chlenov = %f", proiz);
    }
 
}