[pstricks] \pst-func

Patrice Mégret Patrice.Megret at fpms.ac.be
Wed Dec 10 14:00:09 CET 2008



I will generate the points with a mathematical software and then use\dataplot, \fileplot or \listplot.

Best regards,


-----Message d'origine-----
De : pstricks-bounces at tug.org [mailto:pstricks-bounces at tug.org] De la part de Herbert Voss
Envoyé : mercredi 10 décembre 2008 13:37
À : Graphics with PSTricks
Objet : Re: [pstricks] \pst-func

Patrice Mégret schrieb:

> I am using pst-func to plot functions with \psIntegral.
> In the documentation, (page 21) <x t function> and <x function> should be written <t x function> and <t function> to be compatible with the relations 16 to 18 and with the stack order.
> My functions are:
> F(x,t)=60*\int_0^pi (cos(2pixcos(t))-cos(2pix))^2/sin(t) dt
> and: 
> D = 2*(cos^2(2pix))/F
> I have no problem to draw F with \psIntegral:
> \documentclass[11pt,twoside]{article}
> \usepackage{pstricks}
> \usepackage{pst-func}
> \usepackage{pstricks-add}
> \begin{document}
> \pstVerb{/F {exch dup 3 -1 roll 2 mul Pi mul dup 3 -1 roll COS mul COS exch COS sub 2 exp exch SIN div 60 mul} bind def}
> % ce code définit la fonction [cos(2pix cos(t))-cos(2pix)]^2 / sin(t) avec x=h/lambda
> \psset{llx=-1.5cm,lly=-0.5cm,urx=0.2cm,ury=0.2cm,xAxisLabel={$h/\lambda$},xAxisLabelPos={0.5,-45},yAxisLabel={$R_r$ en ohms},yAxisLabelPos={-0.1,150}}
> \begin{psgraph}[Dy=50,Dx=0.1,xticksize=300 0,yticksize=1 0, comma=true,axesstyle=frame](0,0)(1,300){10cm}{5cm}
> \psIntegral[Simpson=100,plotpoints=200,linecolor=orange,linewidth=1.5pt]{0}{1}(0.01,3.13){F}
> \end{psgraph}
> \end{document}
> but I don't see how I can draw D? Maybe I should use the Simpson integration routine inside pst-func.pro but it is not clear for me how to use it. 

there seems no quick solution for that problem
and I am limited with my time until end of this year ...
so I cannot help until then.

> By the way, there is a small trouble with the labeling on the x-axis: if Dx=0.1, the last label 1.0 is not printed, if Dx=0.2, the last label 1.0 is printed. To solve the problem, in my final version I have put 1.001 as the max x-axis and Dx=0.1.

that is the old problem with TeX and rounding errors. The labels
are set on TeX side. But it is no problem to fix it in your way.

PSTricks mailing list
PSTricks at tug.org

More information about the PSTricks mailing list