[ create a new paste ] login | about

Link: http://codepad.org/CCTtwU2C    [ raw code | output | fork ]

C++, pasted on Oct 19:
#include <stdio.h>
//#include <conio.h>
#include <math.h>
 
double f(double x, double eps)
{
    double v, sum = x;
    long n = 1;
    double fn = pow(x,3)/6.0;
    while(eps < fabs(fn))
    {
        sum += fn;
        n   += 1;
        fn *= (2*n - 1)*(2*n - 1)*x*x;
        fn /= 2*n*(2*n + 1);
    }
    return sum;
}
 
int main()
{
    char sign[] = " ";
    double x, fx, dx = 0, xn = 0, xk = 0, eps = 0;
    printf("Enter xn : ");/*scanf("%lf",&xn)*/ xn = -0.9;
    printf("Enter xk : ");/*scanf("%lf",&xk)*/ xk = 0.9;
    printf("Enter dx : ");/*scanf("%lf",&dx)*/ dx = 0.1;
    printf("Enter eps: ");/*scanf("%lf",&eps)*/eps= 1E-3;
    printf("|-----------------------------------|\n");
    printf("|    x   |  f(x)  |  asin |(sin-f)|\n");
    printf("|-----------------------------------|\n");
    for(x = xn; x <= xk; x += dx)
    {
        fx = f(x,eps);
        if(x < 0)
            sign[0] = '\0';
        else
            sign[0] = ' ';
        printf("|%s%.4f |%s%.4f |%s%.4f|%s%.4f  |\n",
            sign,x,sign,fx,sign,asin(fx),sign, fabs(sin(x) - fx));
    }
    printf("|-----------------------------------|\n");
    printf("Press any key to continue\n");
    getchar();
    return 0;
}


Output:
Enter xn : Enter xk : Enter dx : Enter eps: |-----------------------------------|
|    x   |  f(x)  |  asin |(sin-f)|
|-----------------------------------|
|-0.9000 |-1.1165 |nan|0.3331  |
|-0.8000 |-0.9253 |-1.1817|0.2079  |
|-0.7000 |-0.7747 |-0.8862|0.1305  |
|-0.6000 |-0.6431 |-0.6985|0.0784  |
|-0.5000 |-0.5232 |-0.5506|0.0438  |
|-0.4000 |-0.4107 |-0.4232|0.0212  |
|-0.3000 |-0.3045 |-0.3094|0.0090  |
|-0.2000 |-0.2013 |-0.2027|0.0027  |
|-0.1000 |-0.1000 |-0.1002|0.0002  |
|-0.0000 |-0.0000 |-0.0000|0.0000  |
| 0.1000 | 0.1000 | 0.1002| 0.0002  |
| 0.2000 | 0.2013 | 0.2027| 0.0027  |
| 0.3000 | 0.3045 | 0.3094| 0.0090  |
| 0.4000 | 0.4107 | 0.4232| 0.0212  |
| 0.5000 | 0.5232 | 0.5506| 0.0438  |
| 0.6000 | 0.6431 | 0.6985| 0.0784  |
| 0.7000 | 0.7747 | 0.8862| 0.1305  |
| 0.8000 | 0.9253 | 1.1817| 0.2079  |
| 0.9000 | 1.1165 | nan| 0.3331  |
|-----------------------------------|
Press any key to continue


Create a new paste based on this one


Comments: