### C C++ CODE : Shooting method for solving boundary value problem

Working C C++  Source code program for Shooting method for solving boundary value problem
``````#include<iostream.h>
#include<conio.h>
float f1(float ,float,float);
float f2(float, float ,float);
int main()
{
int i;    float x1,x2,y1,y2,h,z,sy1,sy2,sz1,sz2,temp1,temp2,c[2],d[2];
clrscr();
cout<<"Enter x1,y1,x2,y2 ? " ;cin>>x1>>y1>>x2>>y2;
cout<<"Enter h ? ";cin>>h;
temp1=x1;
temp2=y1;

for(i=0;i<2;i++)
{
cout<<"enter z ? ";cin>>z;
d[i]=z;
while(x1<=x2-h)
{
sy1=f1(x1,y1,z);
sz1=f2(x1,y1,z);
sy2=f1(x1+h,y1+h*sy1,z+h*sz1);
sz2=f2(x1+h,y1+h*sy1,z+h*sz1);
y1+=h*(sy1+sy2)/2;
z+=h*(sz1+sz2)/2;
x1+=h;
cout<<endl<<x1<<"   "<<y1<<"   "<<z<<"  ";
}
c[i]=y1;
x1=temp1;
y1=temp2;
}

//for(i=0;i<2;i++)
//cout<<endl<<d[i]<<"   "<<c[i];
z=(y2-c[0])*(d[1]-d[0])/(c[1]-c[0])+d[0];
cout<<endl<<z;

while(x1<=x2-h)
{

cout<<endl<<x1<<"   "<<y1<<"   "<<z<<"  ";
sy1=f1(x1,y1,z);
sz1=f2(x1,y1,z);
sy2=f1(x1+h,y1+h*sy1,z+h*sz1);
sz2=f2(x1+h,y1+h*sy1,z+h*sz1);
y1+=h*(sy1+sy2)/2;
z+=h*(sz1+sz2)/2;
x1+=h;
}
getch();
return 0;
}

float f1(float x,float y,float z)
{
return z;
}

float f2(float x,float y,float z)
{
return 3*x+4*y;
}

``````