Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <iostream>
- #include <cfloat>
- #include <cmath>
- #define N 10
- using namespace std;
- double dx = 0.1;
- double Y0 = -1;
- double yp = 2;
- class DifferentialEquation{
- double precision;
- double init;
- double derivative;
- public:
- DifferentialEquation(double _p, double _i, double _d){
- precision = _p;
- init = _i;
- derivative = _d;
- }
- double * calculatingY(){
- double values[N] = {0};
- values[0] = init;
- for(int i=1; i<N; i++){
- values[i] = values[i-1] + precision*pow(values[i-1], derivative);
- }
- double *result = values;
- return result;
- }
- };
- int main(){
- DifferentialEquation R(dx, Y0, yp);
- double *primaryFunc;
- primaryFunc = R.calculatingY();
- //for(int i=0; i<N; i++){
- // cout << primaryFunc[i] << endl;
- //}
- // Nie mam pojęcia dlaczego powyższy zakomentowany sposób wypisania danych daje inne wartości niż poniższy
- cout << primaryFunc[0] << endl << primaryFunc[1] << endl << primaryFunc[2] << endl << primaryFunc[3] << endl << primaryFunc[4] << endl << primaryFunc[5] << endl << primaryFunc[6] << endl << primaryFunc[7] << endl << primaryFunc[8] << endl << primaryFunc[9];
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement