Изпечатай

[koleto]

Здравейте. Искам да намеря синус на даден ъгъл на Java, без да ползвам Math.sin. За целта реших да използвам Ред на Тейлър.



Ето какво написах от тази сума:

Код

GeSHi (Java):
	/*
	 * Calculating sin x
	 * 
	 * @param x - size of the angle in degrees
	 * @param p - accuracy of calculation
	 * @return sin of x
	 */
	public static double sin(double x, int p){
		double s = 0;
		for(int i=0; i<=p; i++){
			s += Math.pow(-1, i)*(Math.pow(x, 2*i+1)/fac(2*i+1));
		}
		return s ;		 
	}
 
	/*
	 * Calculating factorial
	 * 
	 * @param x - number to factorial
	 * @return x factorial
	 */
	public static double fac(double x){
	if(x<=1)     
            return 1;
        else
            return x * fac(x-1);
	}
 

Проблема е че получавам грешен резултат. Нещо от рода на:

Math.sin(x)=0.5000000000000004
sin(x)=6.19795175229794E19

Ще се радвам ако някой ми помогне да си открия грешката. Мерси предварително.

[shoshon]

Аз не съм кой знае колко добър с математиката, но в ред на Тейлор нямаше ли остатък?

Интересното е, че формулата ти е изключително неоптимизирана и отнема безкрайно много време да се изчисли. Което ме навежда на мислата - ти колко слагаш на p? със сисгурност трябва да е по-голямо от 200-300  за да получиш нещо точно.

От друга страна гледам, че разликата идва още в първите цифри, което говори за грешка от началото на итерацията.

Ама много ме мързи да го пиша...мда...

[clovenhoof]

Нямаш грешка в коденето. Сега го пробах на С++ правилно си смята.

double r = sin(0.52359, 10); // = 0.49999 което е правилно

[clovenhoof]

бе ти да не въвеждаш градуси?

[b2l]

Код

GeSHi (Java):
...
s += (Math.pow(-1, i)/fac(2*i+1))*Math.pow(x, 2*i+1);
...

[shoshon]

Код

GeSHi (Java):
...
s += (Math.pow(-1, i)/fac(2*i+1))*Math.pow(x, 2*i+1);
...

Да те хванат, че да те цапнат...

[b2l]

Да те хванат, че да те цапнат...

Пии ти до 4 па да те видим  тогава...

[bop_bop_mara]

Айде малко офтопик - кой ще го напише на Scheme/Lisp? 

[b2l]

Айде малко офтопик - кой ще го напише на Scheme/Lisp? 

Айде още малко офтопик: Защо да се пише на Scheme като няма смисъл от това?

[lkr]

Айде малко офтопик - кой ще го напише на Scheme/Lisp? 

Код

GeSHi (Scheme):
 

[b2l]

Радвам се, че се намират ентусиасти и че има хора без работа...

[chen_dzen]

а за имането на хора без работа не би трябвало 

[cdiem]

@koleto - ако сложиш повече скоби, мисля, че ще проработи
@mara, @backtolife - трия кода; няма смисъл от него (ще си го запазя за мен си, щото ми е първата работеща програма на този език ;b)

[b2l]

Много ме кефите като почнете да си ги мерите - особено на Scheme. Някой на някой друг такъв безсмислен език няма ли да се включи? Нещо на руски или на китайски ала-бала.

[bop_bop_mara]

Много ме кефите като почнете да си ги мерите - особено на Scheme. Някой на някой друг такъв безсмислен език няма ли да се включи? Нещо на руски или на китайски ала-бала.

Айде дай да ти видим Ц-острото като искаш нещо безсмислено

lkr, cdiem