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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Andrufka
0 / 0 / 0
Регистрация: 18.02.2010
Сообщений: 19
#1

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

10.09.2010, 14:17. Просмотров 1008. Ответов 9
Метки нет (Все метки)

Дамы и господа!!! Выручайте!!!
Надо найти площадь 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;
 
}
0
Миниатюры
Метод прямоугольника  
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
10.09.2010, 14:17
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод прямоугольника (C++):

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

Метод медиан из трех элементов VS улучшенный быстрый метод сортировки(метод Бентли-Макилроя) - C++
Здравствуйте! Дали весьма интересное задание. Сравнить два вышеуказанных метода сортировки для массива из 10000 элементов, результаты...

Мой код - метод бисекции, метод секущих (метод хорд) - C++
Всем привет!!! Изучаем в институте С++. Сделал код, и там, и там одна и та же проблема - при любых вбиваемых значениях программа делает...

Контур прямоугольника - C++
Доброй ночи. Подскажите пожалуйста, пытаюсь написать программку которая выводит на экран пустотелый прямоугольник. в чем ошибка? #include...

Построение прямоугольника - C++
Здравствуйте! Помогите пожалуйста построить прямоугольник шириной 30 и высотой 10! Заранее очень благодарна!

Площадь прямоугольника - C++
Возникла проблема с С++. Недавно начал изучать. Такая задача. Найти площадь прямоугольника, задав с клавиатуры значение длинны и ширины. ...

9
Евгений М.
1035 / 976 / 54
Регистрация: 28.02.2010
Сообщений: 2,829
Завершенные тесты: 2
11.09.2010, 06:57 #2
Цитата Сообщение от Andrufka Посмотреть сообщение
while (i<n) {
* * * * * * * * dy=f1(x)-exp(-x);
dy - это что? И причем здесь экспонента?

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

Добавлено через 54 секунды
Извиняюсь, плохо смотрел.
dy - разность между функциями на маленьком отрезке dx
0
Евгений М.
1035 / 976 / 54
Регистрация: 28.02.2010
Сообщений: 2,829
Завершенные тесты: 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().
0
Andrufka
0 / 0 / 0
Регистрация: 18.02.2010
Сообщений: 19
11.09.2010, 07:41  [ТС] #5
Спасибо огромное!!!
Это просто препод нам так обьяснял, либо я затупил )
Еще раз спасибо!
0
nikkka
Мат в 32 хода
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
11.09.2010, 21:37 #6
а я вообще "метод прямоугольника" в инете не нашёл...
0
Евгений М.
1035 / 976 / 54
Регистрация: 28.02.2010
Сообщений: 2,829
Завершенные тесты: 2
11.09.2010, 21:39 #7
Цитата Сообщение от nikkka Посмотреть сообщение
а я вообще "метод прямоугольника" в инете не нашёл...
? Видимо очень плохо искали.
0
silent_1991
Эксперт С++
4984 / 3041 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
11.09.2010, 21:40 #8
nikkka,
http://ru.wikipedia.org/wiki/Численное_интегрирование
Для начала. Первая ссылка в гугле...
0
nikkka
Мат в 32 хода
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
11.09.2010, 21:41 #9
Евгений М., линк пожалуйста. ;S

Добавлено через 48 секунд
silent_1991, аааа, численное интегрирование.. всеё ясно. я искал прямо "метод прямоугольника".
0
silent_1991
Эксперт С++
4984 / 3041 / 149
Регистрация: 11.11.2009
Сообщений: 7,027
Завершенные тесты: 1
11.09.2010, 21:42 #10
nikkka,
Я тоже. Только прямоугольников. Правда это одна из первых лаб по вычмату, так что мне, в принципе, и искать не надо было)))
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
11.09.2010, 21:42
Привет! Вот еще темы с ответами:

Площадь прямоугольника - C++
Написать программу вычисления площади прямоугольника.

движение прямоугольника - C++
как уменьшить скорость движения? #include &lt;stdio.h&gt; #include &lt;conio.h&gt; #include &lt;stdlib.h&gt; #include &lt;graphics.h&gt; #define BGIPATH...

Движение прямоугольника - C++
Здравствуйте! У меня проблема с этой работой. Программа идет. Но прямоугольник движется по диагонали,а мне нужно,чтобы он двигался по...

Исследовать итерационный метод- метод касательных для решения нелинейных уравнений - C++
прочитал много всего , но сам пример реализовать никак не могу , кто может помогите F(x) = x5+5x+1=0 с...


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

Или воспользуйтесь поиском по форуму:
10
Yandex
Объявления
11.09.2010, 21:42
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru