fork download
copy 
  1. #include <iostream>
  2. #include "math.h"
  3.  
  4. using namespace std;
  5.  
  6. int main()
  7. {
  8.    double h=0.01, A=-2.5/(h*h),B=-2.5/(h*h), F, P=3.14, C;
  9.    int N=1/h;
  10.  
  11.    double *xi=new double [N+1];
  12.    double *vi=new double [N+1];
  13.  
  14.    xi[0]=1;
  15.    vi[0]=0;
  16.  
  17.  
  18.    for(int i=1; i<=N; i++)
  19.    {
  20.        double x=i*h;
  21.        C=-5/(h*h)-sin(3*x);
  22.        F=sin(P*x)/2;
  23.        xi[i]=B/(C-A*xi[i-1]);
  24.        vi[i]=(A*vi[i-1]-F)/(C-A*vi[i-1]);
  25.    }
  26.    double *zi=new double [N+1];
  27.  
  28.    zi[N]=2;
  29.    zi[N-1]=xi[N]*zi[N]+vi[N];
  30.  
  31.    for(int i=N-2; 0<i; i--)
  32.        zi[i]=xi[i]*zi[i+1]+vi[i];
  33.  
  34.    for(int i=0; i<N; i++)
  35.    {
  36.        cout <<"z["<<i<<"]="<<zi[i]<<endl;
  37.    }
  38.  
  39.    return 0;
  40. }
  41.  
Success #stdin #stdout 0s 3472KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
z[0]=0
z[1]=1.77167
z[2]=1.77168
z[3]=1.77168
z[4]=1.77169
z[5]=1.77171
z[6]=1.77174
z[7]=1.77178
z[8]=1.77184
z[9]=1.77191
z[10]=1.772
z[11]=1.77211
z[12]=1.77224
z[13]=1.7724
z[14]=1.77258
z[15]=1.7728
z[16]=1.77304
z[17]=1.77331
z[18]=1.77362
z[19]=1.77397
z[20]=1.77435
z[21]=1.77477
z[22]=1.77523
z[23]=1.77574
z[24]=1.77629
z[25]=1.77689
z[26]=1.77753
z[27]=1.77823
z[28]=1.77897
z[29]=1.77977
z[30]=1.78062
z[31]=1.78153
z[32]=1.78249
z[33]=1.78352
z[34]=1.7846
z[35]=1.78574
z[36]=1.78694
z[37]=1.78821
z[38]=1.78954
z[39]=1.79093
z[40]=1.7924
z[41]=1.79392
z[42]=1.79552
z[43]=1.79718
z[44]=1.79891
z[45]=1.80071
z[46]=1.80259
z[47]=1.80453
z[48]=1.80654
z[49]=1.80863
z[50]=1.81078
z[51]=1.81301
z[52]=1.81531
z[53]=1.81769
z[54]=1.82013
z[55]=1.82265
z[56]=1.82524
z[57]=1.82791
z[58]=1.83065
z[59]=1.83346
z[60]=1.83634
z[61]=1.83929
z[62]=1.84232
z[63]=1.84541
z[64]=1.84858
z[65]=1.85181
z[66]=1.85512
z[67]=1.85849
z[68]=1.86193
z[69]=1.86544
z[70]=1.86901
z[71]=1.87265
z[72]=1.87635
z[73]=1.88011
z[74]=1.88394
z[75]=1.88782
z[76]=1.89177
z[77]=1.89577
z[78]=1.89983
z[79]=1.90394
z[80]=1.90811
z[81]=1.91233
z[82]=1.9166
z[83]=1.92092
z[84]=1.92529
z[85]=1.9297
z[86]=1.93416
z[87]=1.93865
z[88]=1.94319
z[89]=1.94776
z[90]=1.95238
z[91]=1.95702
z[92]=1.9617
z[93]=1.9664
z[94]=1.97114
z[95]=1.9759
z[96]=1.98068
z[97]=1.98548
z[98]=1.99031
z[99]=1.99515
https://ideone.com/kfddVK
language:
C++14 (gcc 8.3)
created:
9 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!