Example 1.4-1 using Maple functions with the Matlab program

 

Note that units are not inserted in the realtions.

> restart;
> x:=T->kapoe*T;A simple Maple function of one variable
> mu:=(T,omegmu)->.1*2.6693e-5*sqrt(M*T)/(sigma^2*omegmu);
A more complex Maple function of two variables
> kapoe:=1./190; sigma:=3.996; M:=44.01;
The CO2 Parameters
> x(200.);x(800.);
Getting the arguments for omegacalc in Matlab for 
T = 200K and 800K.  Omegamu was found to be 1.5471 
at 200K and 0.9595 at 800K
> mu(200,1.5471);mu(800,0.9595);
Executing our mu function for T = 800
[Maple Math] [Maple Math]
> 

BS&L found 1.013e-4 g/(cm*s) at 200 K and 3.269e-4 g/(cm*s) at 800 K