zg.cpp

计算方法中的上机题

#include<iostream.h>
#include<stdlib.h>
#include<math.h>
#define n 4
static float a[n]={0,0,0,0};
static float b[n]={0,0,0,0};
static float c[n]={0,0,0,0};
static float d[n]={0,0,0,0};
static float x[n]={0,0,0,0};
void main()
{
	cout<<"please give a:";
	int i;
	for(i=1;i<n;i++)
	{
		cin>>a[i];
	}
	cout<<"please give b:";
	for(i=1;i<=n;i++)
	{
		cin>>b[i];
	}
	cout<<"please give c:";
	for(i=1;i<n;i++)
	{
		cin>>c[i];
	}
	cout<<"please give d:";
	for(i=1;i<=n;i++)
	{
		cin>>d[i];
	}
	for(i=2;i<=n;i++)
	{
		a[i]=a[i]/b[i-1];
		b[i]=b[i]-a[i]*c[i-1];
	}
	for(i=2;i<=n;i++)
	{
		d[i]=d[i]-a[i]*d[i-1];
	}
	x[n]=d[n]/b[n];
	for(i=n-1;i>=1;i--)
	{
		x[i]=(d[i]-c[i]*x[i+1])/b[i];
	}
	for(i=1;i<=n;i++)
	{
		cout<<x[i]<<"\n";
	}
}