用matlab編寫程序計算正六面體各面的輻射傳熱量。下面是我自己編的，運行不出來，哪位高人能給我改出來

a=5.67e+08;%黑體輻射常數

x=0.2;%角系數

A=1;

b=input('please input;b=');%b為灰度

T1=input('please input;T1=');%表面1的溫度

Eb1=b*a*T1^4;

T2=input('please input;T2='); %表面2的溫度

Eb2=b*a*T2^4;

T3=input('please input;T3=');% 表面3的溫度

Eb3=b*a*T3^4;

T4=input('please input;T4='); %表面4的溫度

Eb4=b*a*T4^4;

T5=input('please input;T5='); %表面5的溫度

Eb5=b*a*T5^4 ;

T6=input('please input;T6='); %表面6的溫度

Eb6=b*a*T6^4;

[Eb1,Eb2,Eb3,Eb4,Eb5,Eb6];

syms J1 J2 J3 J4 J5 J6;

[J1,J2,J3,J4,J5,J6]=solve('(Eb1-J1)/(1/(A*b)-1/A)+(J2-J1)/(1/A*x)+(J3-J1)/(1/A*x)+(J4-J1)/(1/A*x)+(J5-J1)/(1/A*x) +(J6-J1)/(1/A*x)=0',...

'(Eb2-J2)/(1/(A*b)-1/A)+(J1-J2)/(1/A*x)+(J3-J2)/(1/A*x)+(J4-J2)/(1/A*x)+(J5-J2)/(1/A*x) +(J6-J2)/(1/A*x)=0',...

'(Eb3-J3)/(1/(A*b)-1/A)+(J2-J3)/(1/A*x)+(J1-J3)/(1/A*x)+(J4-J3)/(1/A*x)+(J5-J3)/(1/A*x) +(J6-J3)/(1/A*x)=0',...

'(Eb4-J4)/(1/(A*b)-1/A)+(J2-J4)/(1/A*x)+(J3-J4)/(1/A*x)+(J1-J4)/(1/A*x)+(J5-J4)/(1/A*x) +(J6-J4)/(1/A*x)=0',...

'(Eb5-J5)/(1/(A*b)-1/A)+(J2-J5)/(1/A*x)+(J3-J5)/(1/A*x)+(J4-J5)/(1/A*x)+(J1-J5)/(1/A*x) +(J6-J5)/(1/A*x)=0',...

'(Eb6-J6)/(1/(A*b)-1/A)+(J2-J6)/(1/A*x)+(J3-J6)/(1/A*x)+(J4-J6)/(1/A*x)+(J5-J6)/(1/A*x) +(J1-J6)/(1/A*x)=0',...

' J1','J2','J3','J4','J5','J6');

[J1,J2,J3,J4,J5,J6] ;

Q1=(Eb1-J1)/(1/(A*b)-1/A);

Q2=(Eb2-J2)/(1/(A*b)-1/A);

Q3=(Eb3-J3)/(1/(A*b)-1/A);

Q4=(Eb4-J4)/(1/(A*b)-1/A);

Q5=(Eb5-J5)/(1/(A*b)-1/A);

Q6=(Eb6-J6)/(1/(A*b)-1/A);

Q1=eval(Q1)

Q2=eval(Q2)

Q3=eval(Q3)

Q4=eval(Q4)

Q5=eval(Q5)

Q6=eval(Q6)

請使用上面的程序段，問題自己看看對比壹下就行了！

祝妳學習愉快！