求助CRC切线模量计算

凌云

想用切线模量法来考虑结构的初始应力缺陷，但不知道这个方程怎么解。希望高手能帮我使用matlab解一下，最好能给我解方程的文件。式中E=206000，fy=235.微分方程求得的是中间那段曲线。