simulador polymath o matlab..en hp 50? es posible?

bueno gente la respuesta a ese problema en polymath es la siguiente:

con la ecuaciones antes mencionanda hace las siguiente ecuaciones:

dX/dt=F(X,t)=(k(1-t)*(Fbo*t-Nao*X)) / (vo+vo*t)

100 REM ... SEMIBATCH REACTOR PROGRAM ...
200 DIM CA(500) , CB(500) , NA(500), NB(500)
300 DIM NC(500) , RA(500) , X(500) , W(500) , T(500) , F(500)
400 REM ... INITIALIZE CONSTANTS ...
500 N1=2^-.5-.5
600 N2=1-2^-.5
700 N3=1+.5^.5
800 N4=.5^.5
900 CB0=.0025
950 FBO=.00125
1000 CAO=.05
1100 VO=5!
1200 UO=.05
1300 K=2.2
1400 DT=4
1500 NAO=CAO*VO
1600 NBI=0
1700 LPRINT '(ojo esto es una coma algo asi como elevada)INITIAL FLUID VOLUME IN THE VAT =VO = "VO' LITERS'
1800 LPRINT "VOLUME FLOW RATE = UO = "UO" LITERS/SEC"
1900 LPRINT "SPECIFIC REACTION RATE CONSTANT = K = "K" LITERS/SECOND/GRAM MOLE"
2000 LPRINT "MOLAR FLOW RATE OF B = FBO = "FBO" GRAM MOLES/SECOND"
2100 LPRINT "INITIAL MOLES IN VAT: NAI="NAO" NBI = "NBI" GRAM MOLES"
2200 LPRINT "THE STEPSIZE DT = "DT" SECONDS CAO= "CAO" CBO = "CB
2300 LPRINT " "
2400 LPRINT "I T X NA NC RA CB "
2500 LPRINT " SEC MOLES MOLES M/L/S M/L "
2600 I=1
2700 GOTO 3400
2800 TDUMMY=T(I)+DT
2900 I=I+1
3000 X(I)=XDUMMY
3100 T(I)=TDUMMY
3200 REM ....W= XT, THE PRODUCT OF CONVERSION AND TIME ....
3400 NA(I) =NAO*(1-X(I))
3500 NB(I) =NBI +FBO*T(I)-NAO*X(T)
3600 NC(I) =NAO*X(I)
3700 CA(I) =NA(I)/(VO+UO*T(I))
3800 CB(I) =NB(I)/(VO+UO*T(I))
3900 CA(I) =NC(I)/(VO+UO*T(I))
4000 RA(I) =K*CA(I)*CB(I)
4100 REM .... PRINT RESULTS ....
4200 LPRINT I T(I) X(I) NA(I) NC(I) RA(I) CB(I)
4300 REM .... USE RUNGE-KUTTA GILL TO SOLVE DIFFERENTIAL EQUATION ....
4400 REM .... CALCULATE K1 ....
4500 TDUMMY=T(I)
4600 XDUMMY=X(I)
4700 GOSUB 6900
4800 K1=F
4900 REM .... CALCULATE K2 ....
5000 TDUMMY=TDUMMY+.5*DT
5100 XDUMMY=XDUMMY+.5*DT*K1
5200 GOSUB 6900
5300 K2=F
5400 REM .... CALCULATE K3 ....
5500 XDUMMY=N1*DT*K1+N2*DT*K2+X(I)
5600 REM ....TDUMMY=TDUMMY ....
5700 GOSUB 6900
5800 K3=F
5900 REM ....CALCULATE K4 ....
6000 TDUMMY=T(I) +DT
6100 XDUMMY=X(I)-N4*DT*K2+N3*DT*K3
6200 GOSUB 6900
6300 K4=F
6400 REM .... CALCULATE X(I+1) ....
6500 XDUMMY=DT*(K1+2*N2*K2+2*N3*K3*+K4)/6+X(I)
6600 IF I=125 GOTO 6800
6700 GOTO 2800
6800 END
6900 REM .... CALCULATE RA(I)*V/NAO
7000 NA=NAO-NAO*XDUMMY
7100 NB=NBI+FBO*TDUMMY-NAO*XDUMMY
7200 V=VO+UO*TDUMMY
7300 F=K*NA*NB/V/NAO
7400 RETURN

BUENO EN FIN ESE ES LO QUE COLOCA EN EL PROGRAMA Y LUEGO RESUELVE. Y TIRA UNAS GRAFICAS..SON CURVAS DE CA VS t, Cb vs t y Cc vs t.

luego hay una parte que pone tambien en le programa que lo hace asi:

thh equations :
d(x)/d(t)=k*(1-x)*(fbo*t-nao*x)/(vi+vo*t)
nao=.25
vi=5
vo=.05
ca=nao*(1-x)/(vi+vo*t)
nbi=0
fbo=.00125
cb=(nbi+fbo*t-nao*x)/(vi+vo*t)
cc=nao*x/(vi+vo*t)
k=2.2
ra=-k*ca*cb
initial values: to=0.0, xo=0.0
final value: tf= 500.00

bue gente el que pueda hacer un programa similar a este..me salvaria un examen re dificil...

gracias adictos..amo este lugar....