可以这样做
q0=solve('(f*(28+0.42*q)+1132*c)/(0.9*q-13)-(f*(171+0.3*q)+1624*c*c)/(152-0.95*q)','q')；
simple（q0)
-(0.0014947683109118086696562032884903*(537700.0*c + 56380.0*f + 14.142135623730950488016887242097*(2670343200.0*c^4 + 3929511600.0*c^3 + 482644680.0*c^2*f + 1445606450.0*c^2 + 878709340.0*c*f + 37565777.0*f^2)^(1/2) + 730800.0*c^2))/f
-(0.0014947683109118086696562032884903*(537700.0*c + 56380.0*f - 14.142135623730950488016887242097*(2670343200.0*c^4 + 3929511600.0*c^3 + 482644680.0*c^2*f + 1445606450.0*c^2 + 878709340.0*c*f + 37565777.0*f^2)^(1/2) + 730800.0*c^2))/f
得到了q关于f与c的表达式
把随机数f=rand(100);c=rand(100);带入q的表达式中去
q=[-(0.0014947683109118086696562032884903*(537700.0*c + 56380.0*f + 14.142135623730950488016887242097*(2670343200.0*c^4 + 3929511600.0*c^3 + 482644680.0*c^2*f + 1445606450.0*c^2 + 878709340.0*c*f + 37565777.0*f^2)^(1/2) + 730800.0*c^2))/f
-(0.0014947683109118086696562032884903*(537700.0*c + 56380.0*f - 14.142135623730950488016887242097*(2670343200.0*c^4 + 3929511600.0*c^3 + 482644680.0*c^2*f + 1445606450.0*c^2 + 878709340.0*c*f + 37565777.0*f^2)^(1/2) + 730800.0*c^2))/f];
