请问这怎么编程？本人不会

当轴心受压钢构件承受的外力增大到一定程度，构建会发生偏离其轴线发生而破坏，这种构件的破坏称为失稳

untitled.m运行入口脚本

Fi1 = Changxi_ZaiHeK(2,'guanjie');
Fi2 = Changxi_ZaiHeK(4,'hanjiekaikou');

Fall=[Fi1;Fi2];
%dlmwrite('Axialforce.txt',Fi1,'delimiter', '\t'); %只能存一次
dlmwrite('Axialforce.txt',Fall,'delimiter', '\t');%可以存多个

Changxi_ZaiHeK.m

function  K = Changxi_ZaiHeK(lam_Changxi, myType_GuanJie_0dot49_HanJieKaiKouJie_0dot76)

if(strcmp(myType_GuanJie_0dot49_HanJieKaiKouJie_0dot76,'guanjie')==1)
    alpha=0.49;
    lam0=0.40;
else if(strcmp(myType_GuanJie_0dot49_HanJieKaiKouJie_0dot76,'hanjiekaikou')==1)
    alpha=0.76;
    lam0=0.20;
    end
end

fi = 0.5*[1+alpha*(lam_Changxi - lam0)+lam_Changxi^2];

K = 1/( fi + sqrt(fi^2-lam_Changxi^2)   );

end