Sunday, July 31, 2011

Graduate school could be very boring at times

Matlab code for Batman's equation:



clear; clc;
[x,y] = meshgrid(3.01:0.01:10,-3*sqrt(33)/7+0.01:0.01:10);
z=x.*x/49.0 .* sqrt( abs( abs(x)-3)./( abs(x)-3 ))+ y.*y/9.0.*sqrt( abs(y+3*sqrt(33)/7)./(y+3*sqrt(33)/7)) -1;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-10:0.01:-3.01,-3*sqrt(33)/7+0.01:0.01:10);
z=x.*x/49.0 .* sqrt( abs( abs(x)-3)./( abs(x)-3 ))+ y.*y/9.0.*sqrt( abs(y+3*sqrt(33)/7)./(y+3*sqrt(33)/7)) -1;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-4:0.01:4,-5:0.01:5);
z=abs(x)/2-(3*sqrt(33)-7)/112*x.*x-3+sqrt( 1-( abs(abs(x)-2)-1 ).^2 )-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(0.76:0.01:.99,-5:0.01:5);
z=9*sqrt(abs((1-abs(x)).*(abs(x)-0.75))./((1-abs(x)).*(abs(x)-0.75)))-8*abs(x)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-0.99:0.01:-.76,-5:0.01:5);
z=9*sqrt(abs((1-abs(x)).*(abs(x)-0.75))./((1-abs(x)).*(abs(x)-0.75)))-8*abs(x)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(0.51:0.01:.74,-5:0.01:5);
z=0.75*sqrt(abs((0.75-abs(x)).*(abs(x)-0.5))./((0.75-abs(x)).*(abs(x)-0.5)))+3*abs(x)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-0.76:0.01:-0.49,-5:0.01:5);
z=0.75*sqrt(abs((0.75-abs(x)).*(abs(x)-0.5))./((0.75-abs(x)).*(abs(x)-0.5)))+3*abs(x)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-0.51:0.01:0.49,-5:0.01:5);
z=2.25*sqrt(abs((0.5-x).*(x+0.5))./((0.5-x).*(x+0.5)))-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(1.01:0.01:3.000,-5:0.01:5);
z=6*sqrt(10)/7+(1.5-0.5*abs(x)).*sqrt(abs(abs(x)-1)./(abs(x)-1))-6*sqrt(10)/14*sqrt(4-(abs(x)-1).^2)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
[x,y] = meshgrid(-3:0.01:-0.99,-5:0.01:5);
z=6*sqrt(10)/7+(1.5-0.5*abs(x)).*sqrt(abs(abs(x)-1)./(abs(x)-1))-6*sqrt(10)/14*sqrt(4-(abs(x)-1).^2)-y;
contour(x,y,z,[0 0],'b-'); hold on;clear x y z;
axis equal; xlim([-7 7]); ylim ([-3 3]);