| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 		
			
			 | 
		#4 | 
| 
			
			 普通会员 
			
			
			
			注册日期: 2009-02-22 
				
				
				
					帖子: 91
				 
				
				
				声望力: 18 ![]()  | 
	
	
	
		
		
			
			 
			
			不太理解你的意思,看下程序可否: 
		
		
		
		
		
		
		
	syms x format long e ss=zeros(10,10); for k=1:3 for l=1:3 ak=(k-1/2)*pi; al=(l-1/2)*pi; aa=sinh(ak*x)*sinh(al*x); bb=sin(al*x)*sinh(ak*x); pp=ak*al*(aa+bb); ss(k,l)=quadl(inline(pp),0,1); end end 不过,上界1应小点,不然会出现warning!  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 |