Метод прогонки

#include "stdafx.h"
#include"iostream"
#include<conio.h>
#include<stdlib.h>
#include<time.h>
#include <iomanip>
 
using namespace std;
int main()
{
    setlocale(LC_ALL,"");
    int *a,*b,*c,*d,*p,*q,*x,n,i,f;
    srand(time(NULL));
    cout<<"n= ";
    cin>>n;
    cout<<endl;
    cout<<endl<<endl;
    cout<<"Коэффициенты нижней диагонали a "<<endl;
    cout<<endl;
    a=new int[n];
    for(int i=0;i<n;i++)
    {
        cin>>a[i];
        cout<<" ";
        a[a[i]]++;
    }
    cout<<endl<<endl;
    cout<<"Коэффициенты основной диагонали b"<<endl;
    cout<<endl;
    
    b=new int[n];
    for(int i=0;i<n;i++)
    {
        cin>>b[i];cout<<" ";
        b[b[i]]++;
    }
    cout<<endl<<endl;
    cout<<"Коэффициенты верхней диагонали c";
 
    cout<<endl;
    c=new int[n];
    for(int i=0;i<n;i++)
    {
        cin>>c[i];cout<<" ";
        c[c[i]]++;
    }
    cout<<"Вектор свободных членов d";
    
    d=new int[n];
        for(int i=0;i<n;i++)
    {   cin>>d[i];
        cout<<" ";
        d[d[i]]++;
    }
    cout<<endl<<endl;
 
    p[2]=c[1]/b[1];
    q[2]=-d[1]/b[1];
    for (i=2;i=n;i++)
    {
        p[i+1]=c[i]/b[i]-a[i]*p[i];
        q[i+1]=a[i]*q[i]-d[i]/b[i]-a[i]*p[i];
    }
    for(f=n-1;f=1;f--)
    {
        x[n+1]=0;
        x[n]=q[n+1];
        cout<<x[n]<<" ";
    }
    getch();
    return 0;
}