simulacion

clc
//define el paso a incrementar en x e y
Px=0.2;
Py=0.2;

//se define los intervalos y cada cuanto incrementa x e y a evaluar
x=-2:0.2:2;
y=-2:0.2:2;

//halla el tamañano de losvectores x e y
tamanox=size(x,"*")
tamanoy=size(y,"*")

//se define la malla a evaluar
[X,Y] = meshgrid(x,y);

//se define la funcion vectorial de 2 componentes
Fx=X.^1;
Fy=cos(-Y);

//tamaño ypocisión de la ventana de la figura del mapa de contornos
figure('Position',[100 200 400 400])

//se calcula la magnitud del campo F
MagF=sqrt(Fx.^2+Fy.^2)

//obtiene los valores maximo y minimode las matrices
MagFmin = min(MagF); MagFmax = max(MagF);

//mapa de colores para las graficas
xset("colormap",jetcolormap(64))

//se inserta la barra de colores del mapa de contornoscolorbar(MagFmin,MagFmax)

//se inserta el mapa de contornos
Sgrayplot(x,y,MagF)

//se inserta el titulo del mapa de contornos
xtitle("Mapa de contornos de la magnitud del campo vectorial F")//grafica el mapa de vectores
champ1(x, y, Fx, Fy)

//tamaño y pocisión de la ventana de la figura de los vectores del campo vectorial
figure('Position',[200 200 400 400])

//define el color delcampo vectorial
xset("colormap",jetcolormap(64))

//se inserta la barra de colores del campo vectorial
colorbar(MagFmin,MagFmax)

//se dibuja los vectores evaluados del campo vectorialchamp1(x,y,Fx,Fy,1)

//se inserta el titulo del mapa de contornos
xtitle("Campo vectorial de F")

//calculo la divergencia

//X(:,1)//fija a X

//Y(1,:)//fija a Y

//se organiza para garantizar quelos arreglos sean de las mismas dimensiones

//tamaño y pocisión de la ventana de la figura de los vectores del campo vectorial
figure('Position',[300 200 400 400])
DxFx=Fx;
DyFy=Fy;//realiza las derivadas parciales tanto en X como en Y
for j = 1:tamanoy-1
for i = 1:tamanox-1
DxFx(j,i)=Fx(j,i+1)-Fx(j,i)/Px;
DyFy(i,j)=Fy(i+1,j)-Fy(i,j)/Py;
end
end...