dãy Fibonacci

disp('tao day fibonacci ')
n=input(' nhap n ');
fb=[];
fb(1)=1;
fb(2)=2;
for a=3:(n+1)
    fb(a)=fb(a-1)+fb(a-2);
end
for a=1:(n+1)
    fb(a)
end

đường tròn nội tiếp tam giác

disp('nhap toa do diem a');
xa=input('xa=');
ya=input('ya=');
disp('nhap toa do diem b');
xb=input('xb=');
yb=input('yb=');
disp('nhap toa do diem c');
xc=input('xc=');
yc=input('yc=');
ab=sqrt((xa-xb)^2+(ya-yb)^2);
bc=sqrt((xb-xc)^2+(yb-yc)^2);
ca=sqrt((xc-xa)^2+(yc-ya)^2);
while ab+bc<=ca,ac+bc<=ab,bc+ca<=ac;
    disp('nhap lai toa do cac diem a,b,c');
    xa=input('xa=');
    ya=input('ya=');
    disp('nhap toa do diem b');
    xb=input('xb=');
    disp('nhap toa do diem c');
    xc=input('xc=');
    yc=input('yc=');
end
%ve tam giac
u=[xa xb xc xa];
v=[ya yb yc ya];
plot(u,v,'linewidth',2);
hold on
%ve duong tron noi tiep tam giac
disp('tam duong tron noi tiep tam giac,i=');
x1=(bc*xa+ca*xb+ab*xc)/(ab+bc+ca);
y1=(bc*ya+ca*yb+ab*yc)/(ab+bc+ca);
disp('tinh ban kinh truong tron noi tiep tam giac,r:');
p=(ab+bc+ca)/2;;
%nua chu vi tam giac
s=sqrt(p*(p-bc)*(p-ca)*(p-ab));
%dien tich tam giac
r=s/p;
t=0:0.01:360;
x=r*cos(t)+x1;
y=r*sin(t)+y1;
plot(x,y,'linewidth',2);
axis square equal

tam giác nội tiếp đường tròn

a=input('Nhap hoanh do tam duong tron ');
b=input('Nhap tung do tam duong tron ');
r=input('Nhap ban kinh duong tron ');
t=-2*pi:0.001:2*pi;
x=r*sin(t)+a;
y=r*cos(t)+b;
xa=input('Nhap hoanh do xa ');
while xa>a+r|xa<a-r
    xa=input('Nhap lai xa ')
end
ya=-sqrt(r^2-(xa-a)^2)+b;
xb=input('Nhap hoanh do xb ');
while(xb > a+r )| (xb < a-r)
    xa=input('Nhap lai xb ')
end
 
yb=sqrt(r^2 - (xb-a)^2)+b;
xc=input('Nhap hoanh do xc ');
while(xc > a+r) | (xc < a-r)
    xc=input('Nhap lai xc ')
end
yc=-sqrt(r^2 - (xc-a)^2)+b;
x1=xa:xb;
y1=((yb-ya)*x1+xb*ya-xa*yb)/(xb-xa);
x2=xb:xc;
y2=((yc-yb)*x2+xc*yb-xb*yc)/(xc-xb);
x3=xa:xc;
y3=((yc-ya)*x3+xc*ya-xa*yc)/(xc-xa);
ab=sqrt(((xa-xb)^2)+((ya-yb)^2))
bc=sqrt(((xc-xb)^2)+((yc-yb)^2))
ac=sqrt(((xa-xc)^2)+((ya-yc)^2))
plot(x,y,x1,y1,x2,y2,x3,y3)
grid on
text(a,b,'I')
text(xa,ya,'A')
text(xb,yb,'B')
text(xc,yc,'C')
axis square
axis equal

hệ pt tuyến tính

a1=input('nhap a1 ');
a2=input('nhap a2 ');
a3=input('nhap a3 ');
a4=input('nhap a4 ');
a5=input('nhap a5 ');
b1=input('nhap b1 ');
b2=input('nhap b2 ');
b3=input('nhap b3 ');
b4=input('nhap b4 ');
b5=input('nhap b5 ');
c1=input('nhap c1 ');
c2=input('nhap c2 ');
c3=input('nhap c3 ');
c4=input('nhap c4 ');
c5=input('nhap c5 ');
d1=input('nhap d1 ');
d2=input('nhap d2 ');
d3=input('nhap d3 ');
d4=input('nhap d4 ');
d5=input('nhap d5 ');

z=[a1 a2 a3 a4; b1 b2 b3 b4; c1 c2 c3 c4; d1 d2 d3 d4];
z1=[a5 a2 a3 a4; b5 b2 b3 b4; c5 c2 c3 c4; d5 d2 d3 d4];
z2=[a1 a5 a3 a4; b1 b5 b3 b4; c1 c5 c3 c4; d1 d5 d3 d4];
z3=[a1 a2 a5 a4; b1 b2 b5 b4; c1 c2 c5 c4; d1 d2 d5 d4];
z4=[a1 a2 a3 a5; b1 b2 b3 b5; c1 c2 c3 c5; d1 d2 d3 d5];

x1=(det(z1))./(det(z))
x2=(det(z2))./(det(z))
x3=(det(z3))./(det(z))
x4=(det(z4))./(det(z))

hệ pt bậc 1 2 ẩn