%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第七章例7－1－1程序
% 平面力系的簡化
% 西安電子科技大學出版社出版  陳懷琛編著 2007年6月
%
% 平面力系簡化
clear, N=input('輸入力的數目N= (書上取3) ')		%輸入力系的數據
for i=1:N
   i, F(i,:)=input('力F(i)的x,y兩個分量[Fx(i),Fy(i)]=(書上依次取[2,3],[-4,7],[3,-4]) ');
   r(i,:)=input('力F(i)的一個作用點的坐標r(i)=[rx,ry]=(書上依次取[-1,0],[1,-2],[1,2]) ');
end
ro=input('簡化中心ro的坐標ro=[xo,yo]=(書上取[-1,-1]) '); 		%輸入簡化中心的數據
F0=sum(F),			% 求主矢
for i=1:N								%計算各力對ro點的力矩
	m(i)=F(i,2)*(r(i,1)-ro(1))-F(i,1)*(r(i,2)-ro(2));
end
M0=sum(m)								% 相加求主矩
rt = M0/ [F0(2);-F0(1)] + ro				% 求合力作用點的坐標