```
#include <iostream>
using namespace std;
double polynome3(double a0,double a1,double a2, double a3)
`
double y[21];
double z[21];
for (int k=0;k<21;k++)
{
y[0]=1;
y[1]=1;
y[2]=1;
/* Bernoulli's algorithm */
y[k+3]=-(((a1*y[k+2])+(a2*y[k+1])+(a3*y[k]))/a0);
// cout<<y[10]/y[9]<<endl;
}
double alpha1=y[20]/y[19];
cout<< "The first solution is alpha1 = " << alpha1<<endl;
double b0=a0;
double b1=a1+alpha1*b0;
double b2=-(a3/alpha1);
cout<<b0<<" " <<b1<<" "<<b2<<endl;
for(int j=0;j<21;j++)
{
z[0]=0;
z[1]=1;
z[j+2]=-((b1*z[j+1]+b2*z[j])/b0);
}
double alpha2=(z[20]/z[19]);
cout<<" The second solution is alpha2= " <<alpha2<<endl;
double c0=b0;
double c1=-(b2/alpha2);
double alpha3=-(c1/c0);
cout<<c0<<" " <<c1<<endl;
cout<<" The third solution will be alpha3 = "<< alpha3<< "\n"<<endl;
return 0;
}
int main()
{
cout<<"Enter your coefficient : "<<endl;
double x,y,z,t;
cin>>x;
cin>>y;
cin>>z;
cin>>t;
polynome3(x,y,z,t);
/* methode de Bernoulli */
return 0;
}
# Heading Here #
```

## Recommended Answers

Jump to PostAnd?

Jump to PostSome of the things you will see are that you have looped over the end of your arrays. e.g. line 18 and 36. You have written:

`double y[21]; // stuff... for(int k=0;k<21;k++) { y[k+3]=-((a1*y[k+2])+(a2*y[k+1])+(a3*y[k]))/a0); }`

Note that k+3 means that k goes to 24. That will cause memory …

## All 7 Replies

NathanOliver
429
Veteran Poster
Featured Poster

rubberman
1,355
Nearly a Posting Virtuoso
Featured Poster

Luisa_1
0
Newbie Poster

NathanOliver
429
Veteran Poster
Featured Poster

Luisa_1
0
Newbie Poster

rubberman
1,355
Nearly a Posting Virtuoso
Featured Poster

StuXYZ
731
Practically a Master Poster

ddanbe
commented:
Deep insights.
+15

Be a part of the DaniWeb community

We're a friendly, industry-focused community of 1.19 million developers, IT pros, digital marketers, and technology enthusiasts learning and sharing knowledge.