clear all close all clc S=solve( 'r1/M0/(r1/k11)*(Ri/kt0)^(1/2)*24+C=0.13' , ... '(r1*0.21^2+2*0.21*o.79+r2*0.79^2)/M0/(r1*0.21/k11+r2*0.79/k22)*(Ri/kt0)^(1/2)*24+C=0.09' , ... '(r1*0.2^2+2*0.2*o.8+r2*0.8^2)/M0/(r1*0.2/k11+r2*0.8/k22)*(Ri/kt0)^(1/2)*24+C=0.08' , ... '(r1*0.09^2+2*0.09*o.91+r2*0.91^2)/M0/(r1*0.09/k11+r2*0.91/k22)*(Ri/kt0)^(1/2)*24+C=0.44' , ... '(r1*0.05^2+2*0.05*o.95+r2*0.95^2)/M0/(r1*0.05/k11+r2*0.95/k22)*(Ri/kt0)^(1/2)*31+C=0.4' , ... '(r1*0.08^2+2*0.08*o.92+r2*0.92^2)/M0/(r1*0.08/k11+r2*0.92/k22)*(Ri/kt0)^(1/2)*24+C=0.34' , ... '(r1*0.02^2+2*0.02*o.98+r2*0.98^2)/M0/(r1*0.02/k11+r2*0.98/k22)*(Ri/kt0)^(1/2)*24+C=0.43' , ... '(r1*0.15^2+2*0.15*o.85+r2*0.85^2)/M0/(r1*0.15/k11+r2*0.85/k22)*(Ri/kt0)^(1/2)*28+C=0.44' ) S=[S.r1 S.MO S.k11 S.Ri S.kt0 S.C S.r2 S.k22