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

Квадратуре Симпсон - C++

Восстановить пароль Регистрация
 
_abba_
0 / 0 / 0
Регистрация: 31.05.2013
Сообщений: 5
31.05.2013, 12:06     Квадратуре Симпсон #1
Помогите пожалуйста
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
#include <iostream>
#include <conio.h> 
#include <math.h>
#include <iomanip>
#include <fstream>
 
using namespace std; 
const float a=-1,b=1,c=0,d=2, hx=0.1, hy=0.1;
 
main()
{     
      ofstream out("Otvet9.txt");
      int nx=(b-a)/hx;
      int ny=(d-c)/hy;
      float f[300][300],y[300],x[300],g1[300][300],g2[300][300];
      float k[300][300],u[300][300];
      x[0]=a;y[0]=c;
      for (int i=0;i<=nx;i++)
      x[i]=x[0]+i*hx;
      for (int j=0;j<=ny;j++)
      y[j]=y[0]+j*hy;
      for (int i=0;i<=nx;i++)
        { 
         k[i][0]=0;k[i][ny]=0;
        }
      for (int j=0;j<=ny;j++)
      {
       k[0][j]=0;k[nx][j]=0;      
      }
      for (int i=0;i<=nx;i++)
      for (int j=0;j<=ny;j++)
      f[i][j]=y[j]*cos(M_PI*x[i]/2)*sin(M_PI*y[j]/2);
      for (int i=0;i<=nx;i++)
      for (int j=0;j<=ny;j++)
      u[i][j]=cos(M_PI*x[i]/2)*sin(M_PI*y[j]/2);
      for (int i=1;i<nx;i++)
      k[i][1]=(f[i][2]-f[i][0])/(4*hy)-(f[i+1][0]-2*f[i][0]+f[i-1][0]+f[i+1][1]-2*f[i][1]+f[i-1][1])*y[1]/(6*hx*hx);
      for (int i=1;i<nx;i++)
      for (int j=2;j<ny;j++)
      { float Sum1 = 0; 
      for (int j=1;j<ny;j++)
        
        {g1[i][2*j-1]=4*(f[i+1][j]-2*f[i][j]+f[i-1][j]);
             if(g1[i][2*j-1]!=g1[nx][ny-1])
             {
                Sum1 += g1[i][2*j-1];
             }
        }
           float Sum2 = 0;
      for (int j=1;j<ny;j++)
         
           {
            g2[i][2*j]=2*(f[i+1][j]-2*f[i][j]+f[i-1][j]);
            if(g2[i][2*j]!=g2[nx][ny])
            {
                Sum2 += g2[i][2*j];
            }
           }
      for (int j=2;j<ny;j++)
      k[i][j]=(f[i][j+1]-f[i][j-1])/(4*hy)-(f[i+1][0]-2*f[i][0]+f[i-1][0]+Sum1+Sum2+f[i+1][j]-2*f[i][j]+f[i-1][j])*y[j]/(6*j*hx*hx);
      
      }
      
      
      for (int i=1;i<nx;i++)
      for (int j=1;j<ny;j++)
      
      out<<setprecision(10)<<"u["<<i<<","<<j<<"]"<<"   "<<u[i][j]<<"        "<<k[i][j]<<"        "<<fabs(u[i][j]-k[i][j])<<endl;
      out.close();
      getch();
  }
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
31.05.2013, 12:06     Квадратуре Симпсон
Посмотрите здесь:

Pascal Симпсон + Графика
Построение графика. Перевод кода с С++. Симпсон MathCAD

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Tulosba
:)
Эксперт С++
4378 / 3221 / 297
Регистрация: 19.02.2013
Сообщений: 9,044
31.05.2013, 14:04     Квадратуре Симпсон #2
Цитата Сообщение от _abba_ Посмотреть сообщение
Помогите пожалуйста
Помочь Вам сформулировать вопрос, или что?
Yandex
Объявления
31.05.2013, 14:04     Квадратуре Симпсон
Ответ Создать тему
Опции темы

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