простой цикл

[3dg_fan]

привет всем, вот мучаюсь с такой простой задачей
вот код, но здесь ошибка a, b has non-integral type double
invalid types double[1][double] for array subscript

#include <iostream>
#include <cmath>
using namespace std;
	double factorial(double a)
	{
		if(a==1)
			return 1;
				else
					return a*factorial(a-1);
	}

	int main()
	{
		double n, ans,k;
		
		cin >> n;
		double a[n], b[n];
	

		for(k=1; k<n; k++)
		{
			a[k]=3*b[k-1] + 2*a[k-1];
			b[k]=2*a[k-1] + b[k-1];
			ans+=pow(2,k)/((1+a[k]*a[k] + b[k]*b[k])*factorial(k));
		
		}
	     cout<<ans;
	}

[_Alerter_]

double a[n], b[n];

n должна быть константой либо, если вы хотите, чтобы под массив выделялась память для заданного числа элементов, то

cin >> n;
double * a = new double[n]; //выделяем память под массив а
double * b = new double[n]; //выделяем память под массив b
... //код программы
//очищаем память
delete [] a; 
delete [] b;

Для a[0], b[0] и ans нужно задать значения, переменная k должна быть целочисленного типа, т.к. используется в качестве индекса массива

[3dg_fan]

спасибо, разобралась)