C++総和の質問です

質問

Microsoft Visual Basic 6.0の質問です．
色々試してもうまくいかない箇所があるので質問させていただきます．
この式の中で，am=1となっているものをam=1からam=ｘ（xはscanfで入力）
とし，最後のRRR上で1からxまでのRRRの総和を計算させることは可能でしょうか？
基本的な総和と違い，RRRより前に計算しているKmやKn上にそれぞれのamを代入し，最終的にRRRの総和を計算するという形にするやりかたがわからない状態です．

プログラムが無駄に長くて申し訳ありませんが，わかる方いらっしゃいましたらよろしくお願いいたします．

以下プログラムです

//hello2.cpp
#define _COMPLEX_DEFINED
#include <iostream>
#include <string>
#include "complex.h"
#include "bessel.h"
#include <math.h>
#include <cmath>
#include <stdio.h>
#define RADIAN(ARC) ((ARC) * 3.14159 / 180)
using namespace std;

int main(void)
{
const double pi=3.141592654;
const double I0=2829;
const double r0=0;
const double z=0;
const double B1=6.07*pow(10,7);
int k;
double omega,freq,x;
double R,Rk;
double lambdag,ag;
double lambda1,a1,l1;
double lambda2,a2,l2;
double p,p2,pg;
double a11,a21,ag1,Rk1,l11,l21,p11,p21;
double Km,Kn,Kg,am;
complex r1,g,S2;
complex sigma1,sigma2,sigmag;
complex G,P,B,H;
complex RR,NI,NI2,NI3,BE,RRR;
double r,r00,r001,r01;
double imag(),real();

printf("接触熱抵抗");
scanf("%lf",&R);
printf("加熱光半径(10^-6)×");
scanf("%lf",&Rk1);
Rk=Rk1*pow(10,-6);
printf("1層目熱伝導率");
scanf("%lf",&lambda1);
printf("1層目温度伝導率(10^-6)×");
scanf("%lf",&a11);
a1=a11*pow(10,-6);
printf("1層目厚み(10^-6)×");
scanf("%lf",&l11);
l1=l11*pow(10,-6);
printf("2層目熱伝導率");
scanf("%lf",&lambda2);
printf("2層目温度伝導率(10^-6)×");
scanf("%lf",&a21);
a2=a21*pow(10,-6);
printf("2層目厚み(10^-6)×");
scanf("%lf",&l21);
l2=l21*pow(10,-6);
printf("バッキング熱伝導率");
scanf("%lf",&lambdag);
printf("バッキング温度伝導率(10^-6)×");
scanf("%lf",&ag1);
ag=ag1*pow(10,-6);
printf("1層目半径境界(10^-6)×");
scanf("%lf",&p11);
p=p11*pow(10,-6);
printf("2層目半径境界(10^-6)×");
scanf("%lf",&p21);
p2=p21*pow(10,-6);
printf("バッキング半径境界");
scanf("%lf",&pg);

x=2;
freq=pow(10,x);
omega=2*pi*freq;

am=1;
Km=am/p;
Kn=am/p2;
Kg=am/pg;
sigma1=cpow(pow(Km,2)+complex(0,omega/a1),0.5);
sigma2=cpow(pow(Kn,2)+complex(0,omega/a2),0.5);
sigmag=cpow(pow(Kg,2)+complex(0,omega/ag),0.5);

r1=B1/sigma1;
g=(lambdag*sigmag)/(lambda1*sigma1);
S2=(lambda2*sigma2)/(lambda1*sigma1);

G = (1/(2*(S2+g)))*((S2-g)*(-1+R*lambda2*sigma2+S2)*cexp(-2*sigma2*l2)-(S2+g)*(1+R*lambda2*sigma2+S2));
P = (1/(2*(S2+g)))*((S2-g)*(1-R*lambda2*sigma2+S2)*cexp(-2*sigma2*l2)-(S2+g)*(-1-R*lambda2*sigma2+S2));
B = (1/(2*(S2+g)))*((S2-g)*((1-R*lambda2*sigma2)*r1+S2)*cexp(-2*sigma2*l2)-(S2+g)*(-1*(1+R*lambda2*sigma2)*r1+S2));
H = ((g+r1)*cexp(sigma1*l1)*G+(1-g)*cexp(-1*B1*l1)*B)/((1-g)*cexp(-1*sigma1*l1)*P+(1+g)*cexp(sigma1*l1)*G);
r=0;
r00=0;
r001=1;
r01=0.00005;
BE=1;
NI= exp(-1*(2*pow(r01,2)/pow(Rk,2)))*j0(Km*r)*r01;
NI2=2*pi*r001*j0(Km*r00);
NI3=cpow(j1(am),-2);
RR= (B1*I0/(2*lambda1*(cpow(sigma1,2)-pow(B1,2))))*2*pow(p,-2)*cpow(BE,-2)*NI*(1/(1-g))*(-2*H+r1+1)*NI2*cexp(sigmag*z);
RRR=-atan(((RR.imag())/(RR.real()))*(180/pi));

printf("%e\n",R*Rk);
printf("%e\n",Rk);
printf("%e\n",lambda1);
printf("%e\n",B1);
printf("%e\n",sigma1);
printf("%e\n",sigmag);
printf("%e\n",g);
printf("%e\n",r1);
printf("%e\n",j0(am));
printf("%e\n",G);
printf("%e\n",P);
printf("%e\n",B);
printf("%e\n",H);
printf("%e\n",NI3);
printf("%e",RR.real());
printf("+");
printf("%e",RR.imag());
printf("i""\n");
printf("%e\n",RRR);

scanf("%d",&k);//

return 0;
}

回答No.2

f272 · ベストアンサー

いくつか訂正。

よく見たらamもdoubleにしてるのね。
それなら
for (ami=1; am<=x; i++) {
この辺もこのままではまずいね。（ついでに言えば明らかな間違いもあるし...）
ループの変数はint型にするようにして
for (i=1; i<=x; i++) {
am=(double)i;
のようにするんだね。

お礼

ありがとうございます．
ただ，書き間違いがありました！申し訳ありません．
１からｘの総和のｘはプログラム内で定義してあるｘとは別物です…

プログラムにあるのとは別に（仮に）ｘ１として，ｘ１をscanfで入力し，１からｘ１の総和を計算する形になります．

回答No.1

f272

am=1からam=xまでの総和を計算したいということはxは1以上の整数だと想定しているのかな？
だとすればxはint型で定義しておくこと。
その場合
freq=pow(10,x);
の部分もそれに合わせて書きなおすこと。

/* ループに入る前に和を0にしておく */
sumRRR=0.0;
/* ここからループに入る */
for (ami=1; am<=x; i++) {

/* ここでKm, Kn, ...やRRRを計算する */

/* 最後に和を計算する */
sumRRR+=RRR;
}
/* ループを抜けた時点でRRRの総和がsumRRRに入っている */