<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>