plot(y%3Dx%2Cy1%3Dx%5E2%2Cr%3Dcos(theta)%2Cr1%3Dsin(theta))
Want to compute something more complicated? Try a full Python/SymPy console at SymPy Live.