VPF::Интеграл с использованием комплексных переменных

Vasya2904

Дата 14.6.2016, 22:43 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 8
Регистрация: 20.8.2011

Репутация: нет
Всего: нет

Здравствуйте! Имеется функция f(z)=x*intgrl(z)
Решать надо методом средних прямоугольников. Не доходит до меня как правильно(
Мой код

Код


#include <iostream>
#include <cmath>
#include <complex>
#include <fstream>
long double pi = 3.1415;
using namespace std;
int main() {
    long double l,h, k, Z1, q,V, Vk, z1=0, alpha = 0,beta, n=10.00, dz,zz,h2, ed=1, pii = pi/100, m;
    complex <long double>sh(1,pi);
    complex <long double>hh;
    complex <long double>da;
    complex <long double>dj;
    complex <long double> i (0,1);
    complex <long double>p1;
    complex <long double>p2;
    complex <long double>p3;
    complex <long double>p4;
    complex <long double>p5;
    complex <long double>p6;
    complex <long double>z;
    cout << "ENter shirinu: ";
    cin >> l;
    cout <<endl << "ENter visotu ";
    cin >> h;
    cout << endl << "ENter chislo kavitacii ";
    cin >> k;
    cout << endl << "ENter skorost ";
    cin >> V;
    //z=(2*(alpha*(alpha*alpha+1))/pi)*(q/v);
    alpha = l/h;
    //beta=0.9;
    beta = sqrt((pow(alpha,2))*((k+1+(1/pow(alpha,2)))/(k+1+pow(alpha,2))));
cout << beta << endl;
q=h*V; //Расход жидкости в половине струи
zz=q/Vk;
z1 = 2*(alpha*(pow(alpha, 2)+1)/pi);  // Первая чкасть формулы
Vk=sqrt((k+1))*V; //Скорость кавитации
ofstream fout("file.txt");        //Открытие файлов
ofstream fout1("file1.txt");
hh=sh/n;
h2=h/2;


for ( dz=0; dz <= pi; dz = dz+pii)   // циклиз менения dz от 0 до pi
    {
        for (m = 10; m>=1; m--)
        {
    if (beta== dz)
    cout <<  "Error! "<< endl;
    else{
        da = exp(i*dz);   //r*exp(i,w)
        dj=da+h2;
p1 = pow(beta,2)*pow(dj,2)-ed;  
p2 = pow(dj,2)-pow(beta,2);      
p3 = p1/p2;                                
p4 = sqrt(p3);                            
p5 = z1*p4*(pow(dj,2)-ed)/(pow(alpha,2)-pow(dj,2)*(pow(alpha,2)*pow(dj,2))-ed);
//cout<< sqrt((pow(beta,2)*pow(dz,2)-1)) << endl;
//0cout<< sqrt(pow(dz,2)-pow(beta,2));
p6+=p5; //    сумма значений от 0 до pi
z=p6*hh;
cout<<z<< endl;
fout << z.real()<< endl; 
fout1 << z.imag()<< endl; 
}
}//cout<<"\n"<<z;    
    }
    
    
    
    cout<< alpha<<endl<<Vk<<endl<<z1<<endl<< hh;
fout.close();
fout1.close();
//cout << "\n"<<z1;
    return 0;
}

Главное то, что посередине. Программа работает, но графики не такие как надо. Если нетрудно объясните просто какой цикл нужен и я буду очень благодарен:)

frostlatte

Дата 16.6.2016, 14:57 (ссылка)

(голосов: 0)

Загрузка ...

Unregistered

Модератор: Сообщение скрыто.

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Для новичков \| Следующая тема »