非线弹性材料的应力于弯曲关系，matlab求法。

梁（矩形梁，纯弯曲）的材料非线性应力应变曲线如图2所示（数据见右表），现对该梁两端施加M=1000N.mm的弯矩，梁的尺寸为高度20mm，长度50mm，厚度为5mm。
要求：根据材料的应力应变曲线，写出截面应力与弯矩之间弯矩的解析表达式和离散模型表达式，采用数值方法计算该梁截面应力分布，并绘图，并带上单位，调试出结果，配程序源代码。
提示：该材料没有明确的弹性模量，需要根据应变差值得到应力，即应力与应变的比值不是材料力学中学习到的弹性模量的概念。