FANDOM


マンデルブロ集合

前回の続き

#include <stdio.h>

#define ARS 1
#define ARE 4
#define AIS -1.5
#define AIE 1.5
#define XS 0.5
#define YS 0
#define AR_ORDER 0.002
#define AI_ORDER 0.002

int main(void)
	{
	FILE *fs;
	int count=1;
	double x0,xn,xn1,y0,yn,yn1,ar,ai;
	
	ar=ARS;
	ai=AIS;
	
	fs=fopen("plot_mandelbrot.dat","w+");
	
	while(ai <= AIE)
		{
		ar=ARS;
		
		while(ar <= ARE)
			{
			xn=XS;
			yn=YS;
			count=0;
			
			while(1)
				{
				xn1=ar*(xn-xn*xn+yn*yn)-ai*yn*(1-2*xn);
				yn1=ar*yn*(1-2*xn)+ai*(xn-xn*xn+yn*yn);
				
				if(xn1>1000||yn1>1000||xn1<-1000||yn1<-1000)
					{
					break;
					}
				else if(count>=1000)
					{
					fprintf(fs,"%lf %lf\n",ar,ai);
					
					break;
					}
				
				count++;
				
				xn=xn1;
				yn=yn1;
				}
			
			ar+=AR_ORDER;
			}
			
		ai+=AI_ORDER;
		}
	
	fclose(fs);
	
	return 0;
	}
Gnuplot mandelbrot-x0.5-y0.png

2008/07/24 12:29

プロフィール
  • 名前:Tommy6
  • ニックネーム:( ・ω・)
  • 性別:男
  • 生年月日:
  • 出身地:北
  • 血液型:B
2008年07月
1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30 31
最近の記事

Extension:DynamicPageList (DPL), version 2.3.0 : Error: No selection criteria found! You must use at least one of the following parameters: category, namespace, titlematch, linksto, uses, createdby, modifiedby, lastmodifiedby or their 'not' variants

アーカイブ

Extension:DynamicPageList (DPL), version 2.3.0 : Error: No selection criteria found! You must use at least one of the following parameters: category, namespace, titlematch, linksto, uses, createdby, modifiedby, lastmodifiedby or their 'not' variants

YouTube






編集

まだ今月分のポータルが作成されていません!

作成

まだ本日分のポータルが作成されていません!

作成

日記を書く

注意:トップページからのみ日記を作成してください。他のページから作ると記事空間がおかしくなります。

広告ブロッカーが検出されました。


広告収入で運営されている無料サイトWikiaでは、このたび広告ブロッカーをご利用の方向けの変更が加わりました。

広告ブロッカーが改変されている場合、Wikiaにアクセスしていただくことができなくなっています。カスタム広告ブロッカーを解除してご利用ください。

FANDOMでも見てみる

おまかせWiki