function plot_deformed_2DMESH_matlab(nodeCoordinates,elementNodes,...
    displacement,magnification) 
clf
[numelements nodes]=size(elementNodes);
node=zeros(nodes,1);
coordx=zeros(nodes,1);
coordy=zeros(nodes,1);

for ielem=1:numelements
	for inodo=1:4
		node(inodo)=elementNodes(ielem,inodo);
		coordx(inodo)=nodeCoordinates(node(inodo),1);
		coordy(inodo)=nodeCoordinates(node(inodo),2);
		dispx(inodo)=displacement(node(inodo)*2-1);
		dispy(inodo)=displacement(node(inodo)*2);
	end
	xini=[coordx(1),coordx(2),coordx(3),coordx(4)];
	yini=[coordy(1),coordy(2),coordy(3),coordy(4)];
	xfin=xini+magnification*dispx;
	yfin=yini+magnification*dispy;
    xini_yini=[xini' yini']
    xfin_yfin=[xfin' yfin']
    patch([xini],[yini],1,'EdgeColor','black','Facecolor','none')
	patch([xfin],[yfin],1,'EdgeColor','white','Facecolor','blue')
end
return
end