clear all, close all

K=0.9;     % try K from 0.7 to 1

c(1,:)=[-K -K -K];
c(2,:)=[-K +K -K];
c(3,:)=[-K -K +K];
c(4,:)=[-K +K +K];
c(5,:)=[+K -K -K];
c(6,:)=[+K +K -K];
c(7,:)=[+K -K +K];
c(8,:)=[+K +K +K];

x(1)=rand()        % random starting point
y(1)=rand()
z(1)=rand()

for j=1:40000  
   
    cc=c(ceil(8*rand()),:);
   
    x(j+1) = cc(1) + (x(j)-cc(1))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
    y(j+1) = cc(2) + (y(j)-cc(2))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
    z(j+1) = cc(3) + (z(j)-cc(3))/((sqrt((x(j)-cc(1))^2+(y(j)-cc(2))^2+(z(j)-cc(3))^2))^2);
end


plot3(x,y,z,'.','MarkerSize',4) , hold on

楼: Originally posted by jackzhang19 at 2012-03-07 21:52:31:
应该会有偏微分方程吧。你先把方程给出来吧，不是每个人都对这个方程熟悉的。