%Utt=Uxx+Uyy+sin(5x)cos(5y)sin(2t)+sin(3x)cos(2y)(cos(t)-g)
close all; clear all; clc

g = .0002

xVec = linspace(0,pi,55);
yVec = linspace(0,pi,55);
[xMat,yMat] = meshgrid(xVec,yVec);
zMat0 = 2*sin(3*xMat)*cos(2*yMat);
fMat0 = zeros(size(xMat));

subplot(1,3,1)
surfhandel1 = surf(xMat,yMat,zMat0);alpha(.5);
xlabel('x');ylabel('y');zlabel('Displacement');
titlehandle1 = title('')
axis([0,3.14159,0,3.14159,-.3,.3])

subplot(1,3,2)
surfhandel2 = surf(xMat,yMat,fMat0);alpha(.5); 
axis([0,3.14159,0,3.14159,-10,10])
xlabel('x');ylabel('y');zlabel('Force');title('Force')

subplot(1,3,3)
t = linspace(0,10000,150000);
Sum = zeros(size(zMat0));
for i = 1:length(t)
    zMat = (sin(0.2e1 * t(i)) / 0.46e2 - 0.1e1 / sqrt(0.50e2) * sin(sqrt(0.50e2) * t(i)) / 0.23e2) * sin(0.5e1 * xMat).*cos((5 * yMat)) + ((g / 0.13e2 - 0.1e1 / 0.12e2) * cos(sqrt(0.13e2) * t(i)) - g / 0.13e2 + cos(t(i)) / 0.12e2) * sin(0.3e1 * xMat).* cos((2 * yMat));
    fMat = -g*sin(xMat)+sin(5*xMat).*cos(5*yMat)*sin(2*t(i))+sin(3*xMat).*cos(2*yMat)*(cos(t(i))-g);
    set(surfhandel2, 'zdata', fMat);
    set(surfhandel1, 'zdata', zMat);
    title(sprintf('u(x,y,%d)', t(i)))
    Sum = Sum+zMat;
    imagesc(Sum);
    drawnow
    pause(.007)
end