%Rysowanie funkcji DropWave
x1min=-2; x1max=2;
x2min=-2; x2max=2;
R=1500;
x1=x1min:(x1max-x1min)/R:x1max;
x2=x2min:(x2max-x2min)/R:x2max;

for j=1:length(x1)

    for i=1:length(x2)
        frac1 = 1+cos(12*sqrt(x1(j).^2+x2(i).^2));
        frac2 = 0.5*(x1(j).^2+x2(i).^2)+2;
        f(i)=-frac1/frac2;
    end

    f_tot(j,:)=f;

end

figure(1)
meshc(x1,x2,f_tot);colorbar;set(gca,'FontSize',10);
zlabel('f(x)');
title('DropWave','FontSize',20,'FontWeight','bold');