%% graficos de matrices

clearvars;
clear all;
close all;
clc;

%% incializacion de grilla
xmin = -3*pi;
dx = pi/10;
xmax = 3*pi;

ymin = -3*pi;
dy = pi/10;
ymax = 3*pi;

[x,y] = meshgrid(xmin:dx:xmax,ymin:dy:ymax);

%% funcion
z = 5*cos(x).*cos(y) + 0.05*(x.^2 + y.^2);

%% graficos
figure(1)
surf(x,y,z);
title('funcion z = 5cos(x)cos(y) + 0.05(x^2 + y^2)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('z (m)');
axis([ -3*pi,3*pi,-3*pi,3*pi, -1.5*max(max(abs(x))),1.5*max(max(abs(x))) ]);
daspect([1,1,1]);
grid on;
box on;
colormap jet;
colorbar;

figure(2)
mesh(x,y,z);
title('funcion z = 5cos(x)cos(y) + 0.05(x^2 + y^2)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('z (m)');
axis([ -3*pi,3*pi,-3*pi,3*pi, -1.5*max(max(abs(x))),1.5*max(max(abs(x))) ]);
daspect([1,1,1]);
grid on;
box on;
colormap jet;
colorbar;

figure(3)
contour(x,y,z,50);
title('funcion z = 5cos(x)cos(y) + 0.05(x^2 + y^2)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('z (m)');
axis([ -3*pi,3*pi,-3*pi,3*pi ]);
daspect([1,1,1]);
grid on;
box on;
colormap jet;
colorbar;

figure(4)
contour3(x,y,z,50);
title('funcion z = 5cos(x)cos(y) + 0.05(x^2 + y^2)');
xlabel('x (m)');
ylabel('y (m)');
zlabel('z (m)');
axis([ -3*pi,3*pi,-3*pi,3*pi ]);
daspect([1,1,1]);
grid on;
box on;
colormap jet;
colorbar;


for t = 0:0.01*pi:2*pi
    %% funcion

    zt = ( 5*cos(x).*cos(y) + 0.05*(x.^2 + y.^2) ).*cos(t);

    figure(5)
    mesh(x,y,zt);
    title('funcion zt = ( 5cos(x)cos(y) + 0.05(x^2 + y^2) )cos(t)');
    xlabel('x (m)');
    ylabel('y (m)');
    zlabel('z (m)');
    axis([ -3*pi,3*pi,-3*pi,3*pi, -1.5*max(max(abs(x))),1.5*max(max(abs(x))) ]);
    daspect([1,1,1]);
    grid on;
    box on;
    colormap jet;
    colorbar;

    pause(0.01);

end;

%% matrices especiales
N = 4;
M = 5;

M01 = rand(M,N);    %matriz de numeros aleatoreos entre 0 y 1
M02 = zeros(M,N);   %matriz de ceros
M03 = ones(M,N);    %matriz de unos
M04 = eye(M,M);     %matriz unos en la diagonal y cero en todo lo demas