Chebyshev coefficients: Difference between revisions

{{draft task|Mathematics}}

=={{header|C}}==

C99.

#include <string.h>

#include <math.h>

 

return 0;

{{trans|C++}}

using System.Collections.Generic;

using System.Linq;

}

}

{{out}}

<pre>Coefficients:

0.900 0.62160996827066 0.62160996827111 4.444223E-013

0.950 0.58168308946388 0.58168308946379 -8.992806E-014</pre>

 

=={{header|D}}==

This imperative code retains some of the style of the original C version.

 

/// Map x from range [min, max] to [min_to, max_to].

writefln("%1.3f % 10.10f % 10.10f % 4.2e", x, f, approx, approx - f);

}

{{out}}

<pre>Coefficients:

0.900 0.62160996827066445648 0.62160996827066445674 2.71e-19

0.950 0.58168308946388349416 0.58168308946388349403 -1.63e-19</pre>

 

=={{header|Go}}==

 

Two variances here from the WP presentation and most mathematical presentations follow other examples on this page and so keep output directly comparable. One variance is that the Kronecker delta factor is dropped, which has the effect of doubling the first coefficient. This simplifies both coefficient generation and polynomial evaluation. A further variance is that there is no scaling for the range of function values. The result is that coefficients are not necessarily bounded by 1 (2 for the first coefficient) but by the maximum function value over the argument range from min to max (or twice that for the first coefficient.)

 

import (

}

return x1*t - s + .5*c.c[0]

{{out}}

<pre>

1.0 0.54030231 0.54030231 -4.476e-13

</pre>

 

=={{header|J}}==

From 'J for C Programmers: Calculating Chebyshev Coefficients [[http://www.jsoftware.com/learning/a_first_look_at_j_programs.htm#_Toc191734318]]

chebft =: adverb define

:

(2 % x) * +/ f * 2 o. o. (0.5 + i. x) *"0 1 (i. x) % x

)

Calculate coefficients:

10 (2&o.) chebft 0 1

1.64717 _0.232299 _0.0537151 0.00245824 0.000282119 _7.72223e_6 _5.89856e_7 1.15214e_8 6.59629e_10 _1.00227e_11

=={{header|Java}}==

Partial translation of [[Chebyshev_coefficients#C|C]] via [[Chebyshev_coefficients#D|D]]

{{works with|Java|8}}

import java.util.function.Function;

 

System.out.println(d);

}

 

<pre>Coefficients:

6.59630183807991E-10

-1.0021913854352249E-11</pre>

 

=={{header|Julia}}==

{{works with|Julia|0.6}}

{{trans|Go}}

 

c::Vector{Float64}

min::Float64

approx = evaluate(c, x)

@printf("%.1f %12.8f %12.8f % .3e\n", x, computed, approx, computed - approx)

 

{{out}}

0.9 0.62160997 0.62160997 -4.449e-13

1.0 0.54030231 0.54030231 -4.476e-13</pre>

=={{header|Kotlin}}==

{{trans|C}}

 

typealias DFunc = (Double) -> Double

System.out.printf("%1.3f %1.8f %1.8f % 4.1e\n", x, f, approx, approx - f)

}

 

{{out}}

1.000 0.54030231 0.54030231 4.5e-13

</pre>

 

=={{header|Microsoft Small Basic}}==

{{trans|Perl}}

pi=Math.pi

a=0

EndIf

TextWindow.WriteLine(i+" : "+t+cheby[i])

{{out}}

<pre>

9 : -0,0000000000100189955816952521

</pre>

 

 

 

 

 

 

 

{{out}}

-0.232299371615172

 

{{trans|C}}

 

 

 

}

 

{{out}}

-2.3229937e-01

-5.3715115e-02

+1.1521427e-08

+6.5962992e-10

=={{header|Phix}}==

{{trans|Go}}

{{out}}

<pre>

1.0 0.54030231 0.54030231 -4.477e-13

</pre>

=={{header|Python}}==

{{trans|C++}}

 

def test_func(x):

return None

 

{{out}}

<pre>Coefficients:

0.900 0.6216099683 0.6216099683 4.46e-13

0.950 0.5816830895 0.5816830895 -8.99e-14</pre>

=={{header|Racket}}==

 

{{trans|C}}

 

(: chebft (Real Real Nonnegative-Integer (Real -> Real) -> (Vectorof Real)))

(define (chebft a b n func)

(module+ test

(chebft 0 1 10 cos))

 

{{out}}

6.596299173544651e-010

-1.0022016549982027e-011)</pre>

 

=={{header|REXX}}==

 

The numeric precision is dependent on the number of decimal digits specified in the value of '''pi'''.

numeric digits length( pi() ) - length(.) /*DIGITS default is nine, but use 71. */

parse arg a b N . /*obtain optional arguments from the CL*/

if b=='' | b=="," then b= 1 /*B " " " " " */

if N=='' | N=="," then N= 10 /*N " " " " " */

Dma= (b-a) / 2 /*calculate one─half of the difference.*/

Dpa= (b+a) / 2 /* " " " " sum. */

if a=pi then return -1; if a=pi*.5 | a=pi*2 then return 0; pit= pi/3; z= 1

if a=pit then return .5; if a=pit*2 then return -.5; q= x*x; _= 1

/*──────────────────────────────────────────────────────────────────────────────────────*/

pi: pi=3.1415926535897932384626433832795028841971693993751058209749445923078164;return pi

{{out|output|text=&nbsp; when using the default inputs:}}

<pre>

</pre>

{{out|output|text=&nbsp; when using the following input of: &nbsp; &nbsp; <tt> , &nbsp; , &nbsp; 20 </tt>}}

0 Chebyshev coefficient is: 1.647169475390313686961473816799

1 Chebyshev coefficient is: -0.232299371615171942121038341150

18 Chebyshev coefficient is: -3.976201538410589537318561880598E-27

19 Chebyshev coefficient is: 2.859065292763079576513213370136E-29

 

=={{header|Scala}}==

{{Out}}Best seen running in your browser either by [https://scalafiddle.io/sf/DqRNe2A/0 ScalaFiddle (ES aka JavaScript, non JVM)] or [https://scastie.scala-lang.org/M5Ye6h8ZRkmTCNzexUh3uw Scastie (remote JVM)].

 

object ChebyshevCoefficients extends App {

c.foreach(d => println(f"$d%23.16e"))

 

=={{header|Sidef}}==

 

 

 

}

 

 

{{out}}

</pre>

 

=={{header|VBScript}}==

{{trans|Microsoft Small Basic}}

To run in console mode with cscript.

Dim coef(10),cheby(10)

pi=4*Atn(1)

If cheby(i)<=0 Then t="" Else t=" "

WScript.StdOut.WriteLine i&" : "&t&cheby(i)

{{out}}

<pre>

8 : 6,59629917354465E-10

9 : -1,0022016549982E-11

</pre>

 

=={{header|zkl}}==

{{trans|C}}{{trans|Perl}}

fcn chebft(a,b,n,func){

bma,bpa,fac := 0.5*(b - a), 0.5*(b + a), 2.0/n;

})

}

{{out}}

<pre>