Форум программистов, компьютерный форум, киберфорум
Наши страницы

курсач с интегралом.разными методами - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Массив: Поменять местами последний столбец со столбцом, содержащим наибольшее по абсолютной величине элемент матрицы. http://www.cyberforum.ru/cpp-beginners/thread891410.html
В матриці з дійсних чисел поміняти місцями останній стовпець зі стовпцем, що містить найбільший за абсолютною величиною елемент матриці. В матрице из действительных чисел поменять местами...
C++ Блок-схемы Формумчане, у меня нестандартный вопрос: кто-нибудь из вас может мне помочь в составлении блок-схем к программам? Необходимо это для зачета. Программы есть, а вот с составлением схем у меня проблемы.... http://www.cyberforum.ru/cpp-beginners/thread891400.html
C++ как можно сделать покрасивее?
/Задание. Разделение массива на несколько Дан массив i состоящий из n чисел типа int (случайное равномерное заполнение выполняется программой), n=10. Произвести сохранение всех элементов массива...
Как найти локальные максимумы? C++
Подскажите пожалуйста как найти локальные максимумы??
C++ максимальное количество подряд идущих положительных элементов http://www.cyberforum.ru/cpp-beginners/thread891370.html
Помогите с переводом с pascal на С++ Program prog525; Uses Crt; {Массив А - вводимый массив Массив В - текущий массив Массив С - массив для последовательности с наиб. кол. полож. эл. ...
C++ При вычисление выдает nan Здравствуйте, дорогие форумчане. Небольшая программка для вычисления полупериметра треугольника и т.д., реализованная через классы. Полупериметр считает нормально, а вот при подсчете других формул... подробнее

Показать сообщение отдельно
Настенка69
0 / 0 / 0
Регистрация: 15.05.2013
Сообщений: 26

курсач с интегралом.разными методами - C++

04.06.2013, 23:17. Просмотров 265. Ответов 0
Метки (Все метки)

код есть но не совсем рабочий просьба посмотреть
может что подскажите дельногобуду очень благодарна интеграл вида sin ^3 х dx от 0 до п/2
C++
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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
#include "stdafx.h"
#include "stdio.h"
#include "math.h"
 
void func_5(float a, float b, int n, float eps)
{
    int i;
    float h,x,I0,I1;
    h=(b-a)/n;
    x=a+h/2;
    I0=x/(1+(pow(x,2)));
    for (i=1;i<=n-1;i++) {
        x=x+h;
        I0=I0+pow(sin(x),3);
    }
    I0=I0*h;
    I1=0;
    while (fabs(I0-I1)>=eps) {
        I1=I0;
        n=2*n;
        h=(b-a)/n;
        x=a+h/2;
        I0=pow(sin(x),3);
        for (i=1;i<=n-1;i++) {
            x=x+h;
            I0=I0+pow(sin(x),3);
        }
        I0=I0*h;
    }
    printf("integral_5 = %f\n",I0);
}
 
void func_4(float a, float b, int n, float eps)
{
    int i;
    float h,x,I0,I1;
    h=(b-a)/n;
    x=a;
    I0=0;
    for(i=0;i<=n-1;i++) {
        I0=I0+pow(sin(x),3);
        x=x+h;
    }
    //I0=h*(pow(sin(a),3)+pow(sin(b),3))/2+I0);
    I1=0;
    while(fabs(I0-I1)>eps)
    {
        I1=I0;
        n=n*2;
        h=(b-a)/n;
        x=a;
        I0=0;
        for(i=0;i<=n-1;i++)    {
            I0=I0+pow(sin(x),3);
            x=x+h;
        }
        //I0=h*(pow(sin(a),3)+pow(sin(b),3))/2+I0);
    }
    //printf("I=%f\n",I0);
    printf("integral_3 = %f\n",I0);
}
 
void func_3(float a, float b, int n, float eps)
{
    int i;
    float h,x,I0,I1,f2i_1,f2i;
    //если n - нечетное, то умножаем на 2
    if(n%2==1)
        n=n*2;
    printf("eps=");
    scanf("%f",&eps);
    f2i_1=0; //для "нечетных" точек
    f2i=0; //для "четных" точек
    h=(b-a)/n; x=a+h;
    for(i=1;i<=(n-1)/2;i++) {
        f2i_1=f2i_1+pow(sin(x),3); x=x+h;
        f2i=f2i+pow(sin(x),3); x=x+h;
    } 
    f2i_1=f2i_1+pow(sin(x),3);
    //I0=h*pow(sin(a),3)+pow(sin(b),3)+4*f2i_1+2*f2i)/3;
    I1=0;
    while(fabs(I0-I1)>=eps)
    {    
        I1=I0;
        n=n*2;
        f2i_1=0; //для "нечетных" точек
        f2i=0; //для "четных" точек
        h=(b-a)/n; x=a+h;
        for(i=1;i<=(n-1)/2;i++) {
            f2i_1=f2i_1+pow(sin(x),3);  x=x+h;
            f2i=f2i+pow(sin(x),3);  x=x+h;
        }
        f2i_1=f2i_1+sin(3*x);
        //I0=h*pow(sin(a),3)+pow(sin(b),3)+4*f2i_1+2*f2i)/3;
    }
    //printf("I=%f\n",I0);
    printf("integral_3 = %f\n",I0);
}
 
void func_2(float a, float b, int n, float eps)
{
    int i;
    float h,x,I0,I1;
    h=(b-a)/n;
    x=a+h;
    I0=pow(sin(a),3);
    for (i=2;i<=n;i++) {
        x=x+h;
        I0=I0+pow(sin(x),3);
    }
    I0=I0*h;
    I1=0;
    while (fabs(I0-I1)>=eps) {
        I1=I0;
        n=2*n;
        h=(b-a)/n;
        x=a+h;
        I0=pow(sin(x),3);
        for (i=2;i<=n;i++) {
            x=x+h;
            I0=I0+pow(sin(x),3);
        }
        I0=I0*h;
    }
    printf("integral_2 = %f \n",I0);
}
 
void func_1(float a, float b, int n, float eps)
{
    int i;
    float h,x,I0,I1;
    h=(b-a)/n;
    x=a;
    I0=pow(sin(x),3);
    for (i=1;i<=n-1;i++) {
        x=x+h;
        I0=I0+pow(sin(x),3);
    }
    I0=I0*h;
    I1=0;
    while (fabs(I0-I1)>=eps) {
        I1=I0;
        n=2*n;
        h=(b-a)/n;
        x=a;
        I0=pow(sin(x),3);
        for (i=0;i<=n-1;i++)        {
            x=x+h;
            I0=I0+pow(sin(x),3);
        }
        I0=I0*h;
    }
    printf("integral_1 = %f\n",I0);
}
 
int main ( )
{
    int i,n;
    float a,b,h,x,I0,I1,eps;
    printf("a=");
    scanf_s("%f",&a);
    printf("b=");
    scanf_s("%f",&b);
    printf("n=");
    scanf_s("%i",&n);
    
    printf("eps=");
    scanf_s("%f",&eps);
    func_1(a, b, n, eps);
    func_2(a, b, n, eps);
    func_3(a, b, n, eps);
    func_4(a, b, n, eps);
    func_5(a, b, n, eps);
}
Добавлено через 16 минут
никто значит не знает да??
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru