www.pudn.com > fractal.rar > juliaTP.m

clear all; close all; clc;  col=40; m=200; cx=0; cy=0; l=1.5; x=linspace(cx-l,cx+l,m); y=linspace(cy-l,cy+l,m); [X,Y]=meshgrid(x,y); c= -.745429; Z=X+i*Y; for k=1:col; Z=Z.^2+c; W=-abs(Z); A = ( W <-.1); B = 1-A; V = abs(C.*C); D = B.*sin(-0.1*V); end E = D-1.5; colormap copper(256); mesh(0.05*A+E); shading flat; axis([ 0 140 0 140 -2 -0.8]) axis('square','equal','off');