function outputs(t,E,A,l,uz,w,cab_num,x,y,z,r,s0,Etype,bianhao,cab_n,filename1)


for i=1:1:cab_num
    if Etype(i)==0
        s(i,1)=t(i)*s0(i)/(E(i)*A(i))+s0(i);
        derts(i,1)=s(i)-s0(i);
    elseif Etype(i)==1
        derts(i,1)=t(i)*l(i)*(l(i)+uz(i)^2/l(i)+w(i)^2*r(i)^2/(12*t(i)^2*l(i)))/(E(i)*A(i)*r(i));
        s(i,1)=fs(l(i),r(i),uz(i),s0(i),E(i),A(i),w(i),t(i))+s0(i)+derts(i,1);
    end
    %     derts(i)=h(i)*(l(i)+uz(i)^2/l(i))/(E(i)*A(i));
end
X=[bianhao,x,y,z]';

smodelfile=strrep(filename1,'.txt','變形后坐標.txt');
fid=fopen(smodelfile,'wt');     %自應力模態文件名smodel
if fid<0
    errordlg({'文件打開出錯'});
else
    fprintf(fid,'           單元號    單元初始長度     變形后長度       變化長度       桿件徑向力\n');
    fprintf(fid,'%15d%15.5f%15.5f%15.5f%15.5f\n',[cab_n,s0,s,derts,t]');
    fprintf(fid,'\n\n\n          節點號碼      X坐標值       Y坐標值         Z坐標值\n');
    fprintf(fid,'%15d%15.4f%15.4f%15.4f\n',X);
end
fclose(fid);