新手求助！麻烦编个程序。

[legendmq]

用范德华方程计算，333K，100000KPa下，氮的压缩因子。要求用roots命令，分析并找出最后对的根。