Форум программистов, компьютерный форум CyberForum.ru

Метод прямоугольника - C++

Восстановить пароль Регистрация
 
Andrufka
0 / 0 / 0
Регистрация: 18.02.2010
Сообщений: 19
10.09.2010, 14:17     Метод прямоугольника #1
Дамы и господа!!! Выручайте!!!
Надо найти площадь Q способом Монте Карло и способом
прямоугольника!
Метод Монте карло работает вроде правильно, а прямоугольника нет
Помогите разобраться!
Можно и на паскале!
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
#include "stdafx.h"
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <math.h>
 
#define n_MAX 100000L
 
double Karlo();
double MRect();
 
//////////////////////////////////////////////////////////////////////
double f1(double x){return 1+sqrt(1-x*x);}
double f2(double x){return 1-sqrt(1-x*x);}
 
//////////////////////////////////////////////////////////////////////
double main(int argc, char* argv[])
{
//  printf("Hello World!\n");
//  printf("\n%x",RAND_MAX);
//  getch();
 
    double z,s;
    z=Karlo();
    printf("\nz=%f",z);
    s=MRect();
    printf("\ns=%f",s);
    printf("\n");
    return 0;
}
 
//////////////////////////////////////////////////////////////////////
 
double Karlo(){
    long k,n;
    long t;
    double x,y,s;
    double a,b,c,d;
    double r=1; //radiys
    //------------------
    a=-1; b=1; c=0; d=2;
    srand(time(&t)%37);
    n=k=0;
    while (n<n_MAX){
        x=(double)rand()/RAND_MAX;
        y=(double)rand()/RAND_MAX;
        x=(b-a)*x+a;
        y=(d-c)*y+c;
        if (((x*x+pow(y-1,2)<=1) && (x>0) && (y<exp(-x))) ||
            ((x*x+y*y<=1) && (x<0) && (y>exp(-x)))) k++;
        n++;
    }
    s=(double)k/n*(b-a)*(d-c);
    return s;
}
/////////////////////////////////////////////////////////////////////////
double MRect(){
    double a,a2,b,b2,x,dx,dy,s1,s2;
    int i,n;
    n=10000;
    a=-1; b=0;
    dx=(b-a)/n;
    s1=0; x=a; i=0;
    while (i<n) {
        dy=f1(x)-exp(-x);
        if (dy>0) {s1=s1+dy*x;}
        x=x+dx; i++;
    }
    a2=0; b2=1;
    dx=(b-a)/n;
    s2=0; x=a2; i=0;
    while (i<n) {
        dy=f2(x)-exp(-x);
        if (dy>0) {s2=s2+dy*x;}
        x=x+dx; i++;
    }
    return s1+s2;
 
}
Миниатюры
Метод прямоугольника  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Евгений М.
1033 / 974 / 53
Регистрация: 28.02.2010
Сообщений: 2,817
Завершенные тесты: 2
11.09.2010, 06:57     Метод прямоугольника #2
Цитата Сообщение от Andrufka Посмотреть сообщение
while (i<n) {
* * * * * * * * dy=f1(x)-exp(-x);
dy - это что? И причем здесь экспонента?

Добавлено через 54 секунды
Извиняюсь, плохо смотрел.
Andrufka
0 / 0 / 0
Регистрация: 18.02.2010
Сообщений: 19
11.09.2010, 07:00  [ТС]     Метод прямоугольника #3
Цитата Сообщение от Евгений М. Посмотреть сообщение
dy - это что? И причем здесь экспонента?

Добавлено через 54 секунды
Извиняюсь, плохо смотрел.
dy - разность между функциями на маленьком отрезке dx
Евгений М.
1033 / 974 / 53
Регистрация: 28.02.2010
Сообщений: 2,817
Завершенные тесты: 2
11.09.2010, 07:10     Метод прямоугольника #4
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
double MRect(){
        double a,a2,b,b2,x,dx,dy,s1,s2;
        int i,n;
        n=1000000;
        a=-1; b=0;
        dx=(b-a)/n;
        s1=0; x=a; i=0;
        while (i<n) {
                dy=f1(x)-exp(-x);
        if (dy>0) {s1=s1+dy;}
                x=x+dx; i++;
        }
        s1*=dx;
        a2=0; b2=1;
        dx=(b-a)/n;
        s2=0; x=a2; i=0;
        while (i<n) {
                dy=f2(x)-exp(-x);
                if (dy>0) {s2=s2+dy;}
                x=x+dx; i++;
        }
        s2*=dx;
        return s1+s2;
 
}
Исправил следующее:
C
1
s1=s1+dy*x;
На
C
1
s1=s1+dy;
И добавил
C
1
s1*=dx
Аналогично с s2.
Результат более менее близок к результату Karlo().
Andrufka
0 / 0 / 0
Регистрация: 18.02.2010
Сообщений: 19
11.09.2010, 07:41  [ТС]     Метод прямоугольника #5
Спасибо огромное!!!
Это просто препод нам так обьяснял, либо я затупил )
Еще раз спасибо!
nikkka
Мат в 32 хода
 Аватар для nikkka
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
11.09.2010, 21:37     Метод прямоугольника #6
а я вообще "метод прямоугольника" в инете не нашёл...
Евгений М.
1033 / 974 / 53
Регистрация: 28.02.2010
Сообщений: 2,817
Завершенные тесты: 2
11.09.2010, 21:39     Метод прямоугольника #7
Цитата Сообщение от nikkka Посмотреть сообщение
а я вообще "метод прямоугольника" в инете не нашёл...
? Видимо очень плохо искали.
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
11.09.2010, 21:40     Метод прямоугольника #8
nikkka,
http://ru.wikipedia.org/wiki/Численное_интегрирование
Для начала. Первая ссылка в гугле...
nikkka
Мат в 32 хода
 Аватар для nikkka
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
11.09.2010, 21:41     Метод прямоугольника #9
Евгений М., линк пожалуйста. ;S

Добавлено через 48 секунд
silent_1991, аааа, численное интегрирование.. всеё ясно. я искал прямо "метод прямоугольника".
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
11.09.2010, 21:42     Метод прямоугольника
Еще ссылки по теме:

C++ Контур прямоугольника
Площадь прямоугольника C++
C++ СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
silent_1991
Эксперт C++
4938 / 3014 / 149
Регистрация: 11.11.2009
Сообщений: 7,024
Завершенные тесты: 1
11.09.2010, 21:42     Метод прямоугольника #10
nikkka,
Я тоже. Только прямоугольников. Правда это одна из первых лаб по вычмату, так что мне, в принципе, и искать не надо было)))
Yandex
Объявления
11.09.2010, 21:42     Метод прямоугольника
Ответ Создать тему
Опции темы

Текущее время: 16:54. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru