friendship

friendship

Thursday, December 8, 2011

Membuat Tabung posisi mendatar.


Tabung Mendatar
clear;
clc;
disp('-=- Matlab Tabung Mendatar -=-');
t=input('Panjang Tabung = ');
r=input('Diameter Tabung = ');
x=input('Koordinat Sumbu X = ');
y=input('Koordinat Sumbu Y = ');
grid on;
hold on;
axis square;
xlabel ('sumbu X'), ylabel ('sumbu Y');
axis([0 99 0 99]);
title('Matlab Tabung Mendatar');
a=r;
b=r/2;
elips(y,x,a,b);
garis(x+a,y,x+a,y-t);
garis(x-a,y,x-a,y-t);
elips(x,y-t,a,b);

function elips(x,y,a,b)
yi=y-b:0.1:y+b;
xi=(a/b)*(sqrt(b^2-(yi-y).^2))+x;
plot(yi,xi);
plot(yi,2*x-xi);

function garis(xa,ya,xb,yb)
hold on;
axis([0 100 0 100]);
axis square;
if ya==yb
   if xa<xb
      for xi=xa:0.01:xb;
         yi=xi+ya-xi;
         plot(yi,xi);
      end
   else
      for xi=xa:-0.01:xb;
         yi=xi+ya-xi;
         plot(yi,xi);
      end
   end
else
   m=(xb-xa)/(yb-ya);
   if ya<yb
      for yi=ya:.0.01:yb;
         xi=m*(yi-ya)+xa;
         plot(yi,xi);
      end
   else
      for yi=ya:-0.01:yb;
         xi=m*(yi-ya)+xa;
         plot (yi,xi);
       end
   end
end