• Sign Up
• Login
• API
• FAQ
• Tools
• Archive
SHARE
TWEET

# ZMIENIENIONE LAB3

a guest Sep 11th, 2019 110 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. #include <iostream>
2. #include <math.h>
3. #include <stdlib.h>
4. #include <cstdio>
5.
6. using namespace std;
7.
8. double funct_1(double x)
9. {
10.     return sin(x);
11. }
12.
13. double funct_2(double x)
14. {
15.     return cos(x);
16. }
17.
18. double funct_3(double x)
19. {
20.     return -sin(x);
21. }
22.
23.
24. int main(void)
25. {
26.     double a=0.0, b=18.0;
27.     int N=800;
28.     double h=(b-a)/(double)N;
29.     double pochodna_pierwsza,pochodna_druga;
30.     double x;
31.     double err1,err2;
32.     FILE *file1=fopen("pochodna_pierwsza.txt","w");
33.     FILE *file2=fopen("pochodna_druga.txt","w");
34.     fprintf(file1,"X: \t f = Cos(x): \t\t  Pochodna 1: \t\t Blad: \n");
35.     fprintf(file2,"X: \t f = Cos(x): \t\t  Pochodna 2: \t\t Blad: \n");
36.
37.     for(int i=0;i<N;i++)
38.     {
39.         x=a+(double)i*h;
40.         pochodna_pierwsza=(-funct_1(x+2.0*h)+4.0*funct_1(x+h)-3.0*funct_1(x))/(2.0*h);
41.         err1=fabs(funct_2(x)-pochodna_pierwsza);
42.         fprintf(file1,"%.5lf\t%.15lf\t%.15lf\t%.15lf\n",x,funct_2(x),pochodna_pierwsza,err1);
43.         pochodna_druga=(funct_1(x+h)-2.0*funct_1(x)+funct_1(x-h))/(h*h);
44.         err2=fabs(funct_3(x)-pochodna_druga);
45.         fprintf(file2,"%.5lf\t%.15lf\t%.15lf\t%.15lf\n",x,funct_2(x),pochodna_druga,err2);
46.     }
47.     fclose(file1);
48.     fclose(file2);
49.     return 0;
50. }
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top