<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.6002.18071" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Hello</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I am fairly new to using Scilab (and Matlab) but
need to convert some existing files. I have attached an example which I can't
get to work and the image of the correct plot from Matlab.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Original Matlab file:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>rho_m=3300; % in
kg/m^3<BR>rho_infill=2700; % in
kg/m^3<BR>E=1e11; % in kg/m/s^2<BR>v=0.25;<BR>g=9.81; % in
m/s^2<BR>l1=1e3:2e3:50e3; % in
m<BR>l2=75e3:25e3:1000e3; % in m<BR>l=[l1 l2]; % in
m<BR>k=1./l; % in 1/m</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>for Te=0:10e3:100e3 % in
m<BR> D=E*Te^3/(12*(1-v^2));<BR>
phi=1./((D.*k.^4./((rho_m-rho_infill)*g))+1);<BR>
semilogx(k,phi)<BR> hold on<BR> title('\Phi
_e(k) for various T_e')<BR> xlabel('Wavenumber
k')<BR> ylabel('\Phi
_e(k)')<BR>end<BR>text(1.7e-6,0.5,'T_e=100km')<BR>text(3.6e-5,0.5,'T_e=10km')<BR>text(8e-5,0.95,'T_e=0km')</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Output from the conversion below:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>// Display mode<BR>mode(0);</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>// Display warning for floating point
exception<BR>ieee(1);</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>rho_m = 3300;<BR>// in kg/m^3<BR>rho_infill =
2700;<BR>// in kg/m^3<BR>E = 100000000000;<BR>// in kg/m/s^2<BR>v = 0.25;<BR>g =
9.81;<BR>// in m/s^2<BR>l1 = 1000:2000:50000;<BR>// in m<BR>l2 =
75000:25000:1000000;<BR>// in m<BR>l = [l1,l2];<BR>// in m<BR>k = 1 ./l;<BR>//
in 1/m</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>for Te = 0:10000:100000<BR> // in m<BR>
D = (E*(Te^3))/(12*(1-v^2));<BR> phi = 1 ./((D .*(k .^4))
./((rho_m-rho_infill)*g)+1);<BR> // !! L.14: Matlab function semilogx not
yet converted, original calling sequence used.<BR>
semilogx(k,phi)<BR> set(gca(),"auto_clear","off")<BR> title("\Phi
_e(k) for various T_e")<BR> xlabel("Wavenumber k")<BR> ylabel("\Phi
_e(k)")<BR>end;<BR>// !! L.20: Matlab function text not yet converted, original
calling sequence used.<BR>text(0.0000017,0.5,"T_e=100km")<BR>// !! L.21: Matlab
function text not yet converted, original calling sequence
used.<BR>text(0.000036,0.5,"T_e=10km")<BR>// !! L.22: Matlab function text not
yet converted, original calling sequence
used.<BR>text(0.00008,0.95,"T_e=0km")</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Hopefully someone can see what is needed to get
things working. Thanks in advance all.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Lester Anderson</FONT></DIV></BODY></HTML>