function pdm("real[0..s]" x, "real" qe＝0) "//initial running error is zero"