Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function cantr=cantor(l,a,x)
- if nargin <1
- tic
- l=0;
- N=1E7;
- cantr=ones(1,N);
- else
- cantr=a;
- l=l+1;
- endif
- if l==14
- return
- endif
- n=length(cantr);
- cantr(ceil(n/3)+1:2*ceil(n/3))=0;
- cantr(1:ceil(n/3))=cantor(l,cantr(1:ceil(n/3)));
- cantr(2*ceil(n/3)+1:end)=cantor(l,cantr(2*ceil(n/3)+1:end));
- if 1==0
- subplot (2,1,1)
- plot(cantr);
- imcan=zeros(1,length(cantr));
- for i=2:length(cantr)
- if cantr(i)
- imcan(i)=imcan(i-1)+dy;
- else
- imcan(i)=imcan(i-1);
- end
- end
- subplot (2,1,2)
- plot (imcan)
- toc
- endif
- endfunction
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement