... diem tua pt=mt_1*g; size(pt); %ma tran tham so t(0->1) voi buoc step %n la so diem noi suy trong 1 phan doan t=0:step:1; n=round(1/step)+1; %ma tran he so M(const) trong (1.8) m=[2 -2 ... mt_1=inv(mt); %ma tran ve phai cua 1.19 for i=2:(r-1) g(i,:)=3*(pp(i+1,:)-pp(i-1,:)); 101end g(1,:)=3*(pp(2,:)-pp(1,:)); g(r,:)=3*(pp(r,:)-pp(r-1,:)); %ma tran cac vector tiep tuyen ... [r,c]=size(pp); %ma tran he so cua 1.19 v=linspace(1,1,r-1); triu=diag(v,1); tril=diag(v,-1); trim=4*eye(r); trim(1,1)=2; trim(r,r)=2; mt=triu+tril+trim; %ma tran nghich dao cua m4t...