<div dir="ltr"><div>Hi all,</div><div> </div><div>I try to imetate the code for Thomae function to graph like Thomae function F(p/q)=sin(1/sqrt(q)).</div><div>What I got does not look right and I'm not sure what did I do wrong.</div>
<div> </div><div>Here is my try</div><div> </div><div>\documentclass{article}<br>\usepackage{pstricks-add}<br>\usepackage{multido}</div><div>\makeatletter<br>\def\psThomae{\pst@object{psThomae}}<br>\def\psThomae@i(#1,#2)#3{%<br>
\addbefore@par{dotsize=1pt}<br> \begin@ClosedObj<br> \addto@pscode{<br> 1 1 #3 {<br> dup<br> /ipSave ED % save loop value<br> /ip ED % dito<br> 1 1 #3 { <br> dup <br> /iqSave ED % sabve loop value<br>
/iq ED % dito<br> {<br> iq 0 le { exit } if<br> ip iq mod<br> /ip iq def<br> /iq ED<br> } loop<br> ip 1 eq {<br> /xVal ipSave iqSave div def<br> xVal #1 ge { xVal #2 le {<br>
\psk@dotsize<br> \@nameuse{psds@\psk@dotstyle}<br> \pst@usecolor\pslinecolor xVal 1 iqSave sqrt div RadtoDeg sin \tx@ScreenCoor </div><div> 2 copy moveto Dot } if } if<br>
} if<br> } for<br> } for<br> }%<br> \end@ClosedObj%<br>}</div><div>\makeatother</div><div>\begin{document}<br>\begin{center}<br> \psset{unit=4cm,dotsize=1pt 0}<br> \begin{pspicture}(2.5,1.5)<br> \psaxes{->}(0,0)(2.5,1.5)<br>
\psThomae(0,2){300}<br> \end{pspicture}<br>\end{center}</div><div>\end{document}</div><div> </div><div> </div><div>Thanks a lot </div><div> </div><div> </div><div>-- <br></div><div dir="ltr"><div>Regards,</div>
<div> </div>
<div>Hamed Al-Sulami</div></div><br>
</div>