Matlab - Ecuaciones derivadas parciales Urgente

Hola, estoy montando un programa de matlab para resolver una ecuacion de derivadas parciales parabolica. del tipo dC/dt= d^2C/dx + d^2C/dy Esto es lo que llevo. Si alguien tiene un codigo general para este tipo de ecuaciones me ayudaria muchisimo....El error que me arroja es Attempted to access C(2001,2000,1); index out of bounds because
size(C)=[2000,10000,2]. en la linea 2...Aun no consigo mi error cualquier ayuda o codigo general me sirve GRACIASS




h=0.0005;
f=1/h;
c=5/h;
%C=zeros(f,c);
dx=h;
dy=dx;
r=0.5;
g=1;
dt=r.*dx.^2./g;
C(f,c)=1;
for n=1:10000
for i=2:f
for j=2:c
if i==1
C(i-1,j,n)=C(i+1,j,n)-2.*dx.*0.1.*C(i,j,n);
end
if j==1
C(i,j-1,n)=C(i,j+1,n);
end
if i==f
C(i+1,j,n)=C(i-1,j,n);
end
if j==c
C(i,j+1,n)=C(i,j-1,n)
end
C(i,j,n+1)=C(i,j,n)+r.*(C(i-1,j,n)+C(i,j-1,n)-4.*C(i,j,n)+C(j+1,j,n)+C(i,j+1,n));
end
end
t(n+1)=t(n)+dt;
if C(1,f)==0.01
break
disp(t')
end
end

prueba con valor pequeños de iteracion