复数模板类、

/*  
02.* Copyright (c) 2013, 烟台大学计算机学院                      
03.* All rights reserved.                      
04.* 文件名称：test.cpp                      
05.* 作者：董万鹏                    
06.* 完成日期：2013 年 4 月 15 日                      
07.* 版本号：v1.0                    
08.*                      
09.* 输入描述：无                      
10.* 问题描述：                   
11.* 程序输出： 
12.* 问题分析：                     
13.* 算法设计：略                      
14.*/           

#include<iostream>
using namespace std;
template<class numtype>
class Complex
{
public:
    Complex( ){real=0;imag=0;}        
	Complex(numtype r,numtype i):real(r),imag(i){}
	Complex complex_add(Complex &);
	Complex complex_minus(Complex &);
	Complex complex_multiply(Complex &);
	Complex complex_divide(Complex &);
	void display();
private:
	numtype real;
	numtype imag;
};

template<class numtype> 
Complex<numtype> Complex<numtype>::complex_add(Complex<numtype> &c)
{
	Complex<numtype> c1;
	c1.real=real+c.real;
	c1.imag=imag+c.imag;
    return c1;
}

template<class numtype> 
void Complex<numtype>::display()
{
	cout<<real<<imag<<"i"<<endl;
}

template<class numtype>
Complex<numtype> Complex<numtype>::complex_minus(Complex<numtype> &c)
{
	Complex<numtype> c1;
	c1.real=real-c.real;
	c1.imag=imag-c.imag;
	return c1;
}

//复数相乘：(a+bi)(c+di)=(ac－bd)+(bc+ad)i
template<class numtype>
Complex<numtype> Complex<numtype>::complex_multiply(Complex<numtype> &c)
{
   Complex<numtype> c1;
   c1.real=real*c.real-imag*c.imag;
   c1.imag=imag*c.real+real*c.imag;
   return c1;
}

//复数相除：(a+bi)/(c+di)=(ac+bd)/(c^2+d^2) +(bc-ad)/(c^2+d^2)i  
template<class numtype>
Complex<numtype> Complex<numtype>::complex_divide(Complex<numtype> &c)
{
	  Complex<numtype> c1;
      numtype d=c.real*c.real+c.imag*c.imag; 
	  c1.real=(real*c.real+imag*c.imag)/d;      
      c1.imag=(imag*c.real-real*c.imag)/d;        
      return c1;   
}

int main( )
{
	Complex<int> c1(3,4),c2(5,-10),c3;  
	c3=c1.complex_add(c2);  
	cout<<"c1+c2="; 
	c3.display( );
	
	Complex<double> c4(3.1,4.4),c5(5.34,-10.21),c6;  
	c6=c4.complex_add(c5);  
	cout<<"c4+c5="; 
	c6.display( ); 
	
	c6=c4.complex_minus(c5);    
	cout<<"c4-c5=";  
	c6.display( );  
	
	c6=c4.complex_multiply(c5);    
	cout<<"c4*c5=";  
	c6.display( );  
	
	c6=c4.complex_divide(c5);    
	cout<<"c4/c5=";  
	c6.display( );  
	return 0;  
}