Необработанное исключение по адресу 0x4050D302 в ConsoleApplication17.exe: 0xC0000005: Access violation (параметры: 0x00000008)

#include <iostream>
 
#include <stdio.h>
 
#include <math.h>
using namespace std;
 
long int t=0,a=0;
 
long double fx(long double x1,long double x2);
 
long double gradfx(long double x1,long double x2,short int r);
 
int main(){
 
cout<<"gradientnii metod s drobleniem shaga."<<endl;
double x1[562];
 
double x2[562];
 
x1[0]=-1;
 
x2[0]=0;
 
long double alfa=0.39,e=0.0001,delta=0.1;
 
long double minx1,minx2,minf;
 
int k=0;
 
long double *grad = new long double [5];
 
long double D,y,z1;
 
grad[1]=gradfx(x1[0],x2[0],1);
 
grad[2]=gradfx(x1[0],x2[0],2);
 
D=pow(grad[1],2)+pow(grad[2],2);
 
metka2:
 
y=fx(x1[k],x2[k]);
 
metka:
 
z1=fx(x1[k]-alfa*grad[1],x2[k]-alfa*grad[2]);
 
if((z1-y)>-(alfa*delta*D)) {alfa=alfa/2; goto metka; }
 
else{
 
x1[k+1]=x1[k]-alfa*grad[1];
 
x2[k+1]=x2[k]-alfa*grad[2];
 
grad[1]=gradfx(x1[k+1],x2[k+1],1);
 
grad[2]=gradfx(x1[k+1],x2[k+1],2);
 
D=pow(grad[1],2)+pow(grad[2],2);}
 
if(sqrt(D)>e) { k++; goto metka2; }
 
else{
 
minx1=x1[k+1];
 
minx2=x2[k+1];
 
minf=fx(minx1,minx2);}
 
cout<<"minx1="<<minx1<<" minx2="<<minx2<<" minf="<<minf<<endl;
 
cout<<"N="<<t+a<<endl;
 
cout<<"k="<<k<<endl;
 
return(0);
 
}
 
long double fx(long double x1,long double x2){
 
a++; long double f;
 
f=x1-1.4*x2+exp(0.01*x1+0.11*x2);
 
return(f);}
 
long double gradfx(long double x1,long double x2,short int r){
 
::t++;
 
long double gradf;
 
if (r==1)
 
    gradf=1+0.01*exp(0.01*x1+0.11*x2);
 
if (r==2)
 
    gradf=-1.4+0.11*exp(0.01*x1+0.11*x2);
 
return(gradf);
 
}