http://goo.gl/ug9HQI

\documentclass{article}

\usepackage{luamplib}

\begin{document}

\begin{mplibcode}

beginfig(1);

def compute_curve(suffix f)(expr xmin,xmax,xinc) = ( (xmin,f(xmin))

    for x=xmin+xinc step xinc until xmax: ..(x,f(x)) endfor )

enddef;   

vardef f(expr x)=sind(x*180/3.14) enddef;

path g;

g:=compute_curve(f,-1.57,9.44,0.1);

draw g xscaled 8pt yscaled 0.9cm withpen pencircle scaled 2bp;

draw (-1cm,0)--(3cm,0);

draw (0,1cm)--(0,-1cm);

draw fullcircle scaled 1.8cm;

label(btex$y$etex, (0,1.2cm));

label(btex$x$etex, (3.2cm,0));

endfig;

\end{mplibcode}

\end{document}