Hi,
I'm using rk to solve a differential equation set. The script is
attached. There is a piecewise function I wrote in two ways: q(s) and
qq(s), where q(s) is simplified from pp(s) so actually they are
equivalent. But when I use them in the differential equation, one (q(s))
is solvable and the other (qq(s)) is not which results in a following error:
expt: undefined: 0 to a negative exponent.
Besides, the solved results are also not good (not smooth solution).
What causes these problems? Thanks.
Best regards,
Ning
--
Geotechnical Group
Department of Civil and Environmental Engineering
Hong Kong University of Science and Technology
Clear Water Bay, Kowloon, Hong Kong
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: stz2.wxm
URL: <http://www.math.utexas.edu/pipermail/maxima/attachments/20120413/fa2a9681/attachment.ksh>;