Oblicz wartość funkcji w punktach - problem z sinusem

Dzień dobry, piszę program obliczający wartość funkcji w punktach x=a+i*h, gdzie i=0,..., M;
Funckja określona wzorem

f=1/3*sin(x)+1/2

a,h - typ rzeczywisty
M -całkowite

Na chwilę obecną mój program wygląda tak, ale ewidentnie gdzieś jest błąd ponieważ wyrzuca mi wynik 0. Niestety nie wiem co robię źle i serdecznie proszę o jakies naprowadzenie.

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <cmath>
#include <ctime>

using namespace std;

int main()
{   float a,h;
    float funkcja;
    int M;
    cout << "Podaj wartosc a: " <<endl;
    cin >> a;
    cout << "Podaj wartosc h: " <<endl;
    cin >> h;
    cout << "Podaj wartosc M: " <<endl;
    cin >> M;
    for (int i=0; i<=M; i++)
    {
           funkcja = 1/3*sin(a+i*h)+1/2;
       }
    cout << "Funkcja f wynosi " << funkcja <<  endl;
    return 0;
}

Dziękuję bardzo! Taka drobnostka, a teraz działa.