הנה שתיים
תוכנית ראשונה - מציירת את עקומת הילברט:
clear s=input('Enter stage (1-9): '); f=[]; for i=1:s f=[pi/2-f pi/2 f 0 f 3*pi/2 3*pi/2-f]; end x=zeros(1,4^s); y=x; for i=1:4^s-1 x(i+1)=x(i)+cos(f(i)); y(i+1)=y(i)+sin(f(i)); end clf plot(x,y) axis square axis off
שימו לב שכאשר אתם מריצים את התוכנה היא מצפה לקלט - מספר שלם בין 1 ל - 9. תוכנית שנייה - נחשו לבד מה היא מציירת
:
clear clf a=[.849 .197 -.15 0]; b=[.037 -.226 .283 0]; c=[-.037 .226 .26 0]; d=[.849 .197 .237 .16]; e=[.075 .399 .572 .497]; f=[.183 .049 -.084 0]; x=zeros(1,1000000); y=x; x(1)=e(4); for i=2:1000000 k=rand; r=1+(k>.73)+(k>.86)+(k>.97); x(i)=a(r)*x(i-1)+b(r)*y(i-1)+e(r); y(i)=c(r)*x(i-1)+d(r)*y(i-1)+f(r); end plot(y,x,'k.','MarkerSize',1) axis([0 1.1 .3 .8]) axis off
בשתי התוכניות תעשו זום על הציור המתקבל.