matlab中p p1 p2,多元函数f=f(t,u,P1,P2,P3)的求解和数值提取和作图的问题

syms t u P1  P2  P3;

a12=21.98*10^7;

a111=2.57*10^9;

a112=6.95*10^9;

a123=13.13*10^9;

S11=11*10^-12;

S12=-4.1*10^-12;

S44=27.5*10^-12;

Q11=0.084;

Q12=-0.025;

Q44=0.035;

>> A13=a12+Q12*(Q11+Q12)/(S11+S12);

>> A12=a12-[(Q11^2+Q12^2)*S12-2*Q11*Q12*S11]/(S11^2-S12^2)+Q44^2/2/S44;

>>  eq1=sym('A1=(3.255*t-1960)*10^5-u*(Q11+Q12)/(S11+S12)');

eq2=sym('A3=(3.255*t-1960)*10^5-2*u*Q11/(S11+S12)');

eq3=sym('A11=(-0.3775*t-53.92)*10^5-0.5*[(Q11^2+Q12^2)*S11-2*Q11*Q12*S12]/(S11^2-S12^2)');

eq4=sym('A33=(-0.3775*t-53.92)*10^5+Q12^2/(S11+S12)');

>> eq5=sym('0=2*A1*(P1+P2)+4*A11*(P1^3+P2^3)+2*A13*(P1*P3^2++P2*P3^2+P3*P1^2+P3*P2^2)+2*A3*P3+4*A33*P3^3+2*A12*(P1*P2^2+P1^2*P2)+6*a111*(P1^5+P2^5+P3^5)+2*a112*(2*P1^3*P2^2+2*P1^3*P3^2+P1^4*P2+P3^4*P2+P1^4*P3+2*P3^3*(P1^2+P2^2)+P2^4*P3)+2*a123*(P1*P2^2*P3^2+P1^2*P2*P3^2+P1^2*P2^2*P3);');

>> [t u P1  P2  P3]=solve(eq1,eq2,eq3,eq4,eq5,t,u,P1,P2,P3);

Warning: Explicit solution could not be found.

> In solve at 160

主要想做出来f=f(t,u,P1,P2,P3)=0时，的取值   谢谢啊