help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

still problems, version? (was: matlab code: isovalues and isosurface)


From: Uwe Brauer
Subject: still problems, version? (was: matlab code: isovalues and isosurface)
Date: Sat, 16 Jul 2011 20:21:02 +0200
User-agent: Gnus/5.110011 (No Gnus v0.11) XEmacs/21.5-b29 (linux)

   > Am Samstag, 16. Juli 2011, 17:25:25 schrieb Uwe Brauer:
   > You simply type (octave 3.4)
   > graphics_toolkit fltk
   > at the beginning of your script or octave session. In 3.2 it was "backend 
fltk" 
   > as far as I remember.


Thanks but the suggested code does not work for me.
I have attached the code with the your and Liam's
modification below.


In Linux backend fltk works, but when running the code again
I obtain 


,----
| error: number of columns must match (0 != 6)
| error: called from:
| error:   /usr/share/octave/3.2.3/m/plot/__patch__.m at line 46, column 12
| error:   /usr/share/octave/3.2.3/m/plot/patch.m at line 55, column 5
| error:   
/home/oub/ALLES/tex/vorlesungen/11-12/Calculo-Comput_11_12/Cuantica/newcode.m 
at line 15, column 13
| 
`----

In windows where I also have 3.2.4 installed the same errors
appears. It seems that 3.4. Is not available as precompiled
binaries, neither for Linux nor for Windows. Am I correct?

Uwe 

,----
| [x,y,z]=meshgrid(-20:.5:20);
| function Psi=FdOn320(x,y,z)
|  r=sqrt(x.^2+y.^2+z.^2);
|  Psi=exp(-r./3).*(3*z.^2-r.^2).*2*sqrt(15/(120*pi))/81;
| endfunction 
|  Psi=FdOn320(x,y,z);
|  m=min(Psi(:));
|  M=max(Psi(:));
|  color=[1 1 0;1 0 1; 0 1 1;0 1 0;1 0 0; 0 0 1; 0 0 0];
|  n=5;
|  for i=1:n-1
|      isovalue=m+i*(M-m)/n;
|      fv = isosurface(x,y,z,Psi,isovalue);
|      hpatch = patch(fv);
|      isonormals(x,y,z,Psi,hpatch);
| %     Alphalevel=0.2+0.8*i/n;
|      set(hpatch,'FaceColor',color(i,:),'EdgeColor','black')
| 
| %       
set(hpatch,'FaceAlpha',Alphalevel,'FaceColor',color(i,:),'EdgeColor','none')
| %       set(hpatch,'FaceColor',color(i,:),'EdgeColor','none')
|       daspect([1,1,1])
|       view(3); 
|       axis tight
|  %     camlight left; 
|  end
| 
| 
`----



reply via email to

[Prev in Thread] Current Thread [Next in Thread]