Некорректное выполнение программы

#include<iostream>
#include <iomanip>
#include<stdio.h>
#include<math.h>
using namespace std;
double func(int i,double c, double d);
int main(void)
{
setlocale(LC_ALL,"Russian");
double h=0.1;//постоянный шаг
double x[11];//массив значений x
double y[11];//массив значений y
double k[4];
double f[11];
//задаем разности
double df[11];
double d2f[11];
double d3f[11];
x[0]=0;//левая граница отрезка
for(int i=0;i<11;i++)
{x[i+1]=x[i]+h;}//прибавляем шаг
y[0]=0;
k[0]=func(0,x[0],y[0]);
for(int i=1;i<3;i++)
{k[i]=func(0,x[0]+h/2,y[0]+k[i-1]/2);}
k[3]=func(0,x[0]+h,y[0]+k[2]);
y[1]=y[0]+h*(k[0]+2*k[1]+2*k[2]+k[3])/6;//вычисление первого отрезка методом Рунги-Кутта
for(int i=0;i<11;i++)
{f[i]=func(i,x[i],y[i]);}
for(int i=1;i<11;i++)
//вычисляем разности
{df[i]=f[i]-f[i-1];}//первая разность
for(int i=2;i<11;i++)
{d2f[i]=f[i]-2*f[i-1]+f[i-2];}//вторая разность
for(int i=3;i<11;i++)
{d3f[i]=f[i]-3*f[i-1]+3*f[i-2]-f[i-3];}//третья разность
for(int i=1;i<11;i++)
{
y[i+1]=y[i]+h*f[i]+h*h*df[i+1]/2+5*h*h*h*d2f[i+2]/12+3*h*h*h*h*d3f[i+3]/8;//разностная схема метода Адамса
}
cout<<"Таблица вычислений:"<<endl;
for(int i=0;i<11;i++)
{cout<<setw(10)<<x[i]<<setw(10)<<y[i]<<endl;}
getchar();
getchar();
return 0;
}
//пользовательская функция
double func(int i,double a, double b)
{
double f;
f=(1-b*b)*cos(a)+0.6*b;//исходное дифференциальное уравнение
return f;
}