Hallo Rafael,
your script is perfect, mine isn't. Where is the bug?
Regards,
Jens
clc(), clear(), mode(0),lines(0)
R=1.; //Radius Erde
a=50;//Abstand Mond - Erde
g0=1; //Schwerebeschleunigung des Mondes im Zentrum der Erde
VLev=[-0.02:0.002:0.02]; //Werte der Äquipotenziallinien-Kandidaten
n=400;//Anzahl der x- und y-Werte
xdel();
cmap=zeros(100,3); //Colormap
f = scf(); f.color_map = cmap;
x=linspace(a-R,a+R,n); y=linspace(-R,R,n); //Erdgebiet
[X,Y]=ndgrid(x,y);
pos=find((X-a).^2+Y.^2 > R^2);
X(pos)=%nan; Y(pos)=%nan;
V=-g0*(X+a^2../sqrt(X.^2+Y.^2)-2*a);
philev=contour2d(x,y,V, VLev)//Plot und dargestellte Äquipotenziallinien
ca=gca();
ca.isoview="on";
xset("fpf"," ") //Keine Levelanzeige, Reaktion erst nach dem 2. Lauf
-----------------------------------------------------------------------------------------------------------------------------
Am 29.04.2016 14:24, schrieb Rafael Guerra:
Hello,
Is the following example ok?
t=-%pi:0.02:%pi;
M=sin(t)'*cos(t);
n=length(t);
r=2;///radius/
fori=1:n
forj=1:n
ifsqrt(t(i)^2+t(j)^2)>rthenM(i,j)=%nan;end
end
end
_clf_();
_Sgrayplot_(t,t,M);
Otherwise, could you provide a sample of your pb script.
Regards,
Rafael
*From:*users [mailto:users-boun...@lists.scilab.org] *On Behalf Of
*Jens Simon Strom
*Sent:* Friday, April 29, 2016 1:00 PM
*To:* International users mailing list for Scilab.
<users@lists.scilab.org>
*Subject:* [Scilab-users] Circular contour plot
Hallo Scilab afficionados,
I would like to cut off all lines of a 2D contour plot which are
outside of a circle. I tried to set all external x and y to %nan.
This works however partially spoils the plot /inside/ the circle too.
Polar organisation of the x-y-points doesn't help either because x and
y have to be monotone. Could anyone give me a hint how to proceed here?
Kind Regards
Jens
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
