Matlab - Ecuaciones derivadas parciales Urgente

   
Vista:

Ecuaciones derivadas parciales Urgente

Publicado por Cesar (2 intervenciones) el 20/01/2014 01:29:27
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
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
-1
Responder
Imágen de perfil de JOSE JEREMIAS CABALLERO

Ecuaciones derivadas parciales Urgente

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3441 intervenciones) el 20/01/2014 04:12:14
prueba con valor pequeños de iteracion
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar