# The DE is S'=rS-k, with r=0.04 and k=(150)*(12), with S(0)=9000: Eqn1:=diff(S(t),t)=0.04*S(t)-(150)*12;

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJTRzYiNiNJInRHRilGKywmRickIiIlISIjISUrPSIiIg==

AmtMoney:=rhs(dsolve( { Eqn1, S(0)=9000},S(t)));

LCYiJitdJSIiIi1JJGV4cEc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjLCRJInRHRiojRiQiI0QhJitnJA==

# Paid in full? Set the equation to zero- Since we're in Maple, we can also solve it. solve(AmtMoney=0,t);

LCQtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjIyIiJiIiJSIjRA==

# The DE will be P'=kP with P(0)=50. The solution is easy enough to write directly, but in Maple: Eqn2A:=diff(P(t),t)=k*P(t);

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJQRzYiNiNJInRHRilGKyomSSJrR0YpIiIiRidGLg==

P1:=rhs(dsolve({Eqn2A,P(0)=50},P(t)));

LCQtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2IyomSSJrR0YoIiIiSSJ0R0YoRiwiI10=

# Now solve for k- We know that t=2 gives 100 fish: Eqn2B:=subs(t=2,P1)=100;

LywkLUkkZXhwRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMsJEkia0dGKSIiIyIjXSIkKyI=

kk:=solve(Eqn2B,k);

LCQtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjIiIjIyIiIkYq

# Final population of fish so far: P2:=subs(k=kk,P1);

LCQtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2IywkKiYtSSNsbkdGJTYjIiIjIiIiSSJ0R0YoRjAjRjBGLyIjXQ==

# If we harvest removing k fish per year, then: Eqn2C:=diff(P(t),t)=kk*P(t)-k;

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJQRzYiNiNJInRHRilGKywmKiYtSSNsbkc2JEYlSShfc3lzbGliR0YpNiMiIiMiIiJGJ0Y0I0Y0RjNJImtHRikhIiI=

# Since we're using Maple, we can also go ahead and solve (this is linear and separable): dsolve({ Eqn2C,P(0)=50 },P(t));

Ly1JIlBHNiI2I0kidEdGJSwmKiYtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYjIiIjISIiSSJrR0YlIiIiRjAqJilGMCwkRicjRjNGMEYzLCYiI11GM0YpISIjRjNGMw==

# Interesting that Maple chose to move into base 2 instead of base e! #Cannonball: mv'=mg-gam v, or v'=g-(gam/m) v m:=10; g:=9.8; gam:=5;

IiM1

JCIjKSohIiI=

IiIm

Eqn3:=diff(v(t),t)=g-(gam/m)*v(t);

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJ2RzYiNiNJInRHRilGKywmJCIjKSohIiIiIiJGJyNGLyIiIw==

# Note that initial velocity is -10 (minus for "up") v1:=rhs(dsolve({Eqn3,v(0)=-10},v(t)));

LCYjIiMpKiIiJiIiIi1JJGV4cEc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjLCRJInRHRiwjISIiIiIjIyEkWyJGJQ==

# Just for fun, we can plot the velocity to see if it makes sense (and it seems to!). plot(v1,t=0..10);

LSUlUExPVEc2KC0lJ0NVUlZFU0c2JDdkdzckJCIiISEiIiQhJCsiISIiNyQkIjFHYzdEXWdkXyEjPCQhMTpiPDdKLEsjKiEjOjckJCIxS2hBWCE+QSQpKiEjPCQhMUhmWjE5LSFlKSEjOjckJCIxKFF4YTQmbyhcIiEjOyQhMTpzIkgpZVFreSEjOjckJCIyMTE3Q1tlYiwjISM8JCExMyVmIylRZ0I7KCEjOjckJCIyKClvUHZdcTRgIyEjPCQhMjo2KW8iPSpbIlsnISM7NyQkIjEvMztLTSMpM0khIzskITF5ZldyKT1kJ2UhIzo3JCQiMkNWJ0dkQ2guTiEjPCQhMW5vPDssV1ZfISM6NyQkIjFzVihbKFJLOlMhIzskITEoZid5NyFwZWglISM6NyQkIjIkeWM4RldSRFghIzwkITFNdF0lKXkzMVMhIzo3JCQiMXRYIkhlayswJiEjOyQhMnVvUWhrcFtSJCEjOzckJCIxViYzPE0lPjdiISM7JCEyajxoQmJJJ3BHISM7NyQkIjF2WyhcKkhXS2chIzskITJkP0s8SnhFSCMhIzs3JCQiMClmPlJ5I1tiJyEjOiQhMTFpdUguRUc8ISM6NyQkIjEwMztLQ0NlcSEjOyQhMlpQWiRbRDYpPiIhIzs3JCQiMWtGYjU2Ujp2ISM7JCExQ3g4Jz5JNkcoISM7NyQkIjFRd18wXikqZSEpISM7JCEydHBhXSoqPi8kPSEjPDckJCIxbEloQTBbPiYpISM7JCIxWy5NTjhhc0UhIzs3JCQiMVoiSGU7WV4wKiEjOyQiMU0nKmYjKSplPHkoISM7NyQkIjFHYzdENUlIJiohIzskIjI7K2ouJjQ6PjchIzs3JCQiMiNRd18wQyZcKyIhIzskIjInUWpiMU40InAiISM7NyQkIjEwNT9TKipbYTUhIzokIjJrSExwakUjSEAhIzs3JCQiMi4tLzM7eGg1IiEjOyQiMjhRTGI4XVxkIyEjOzckJCIyIzM7S2tAa2A2ISM7JCIyRi07Rk1TVShIISM7NyQkIjJHZDlIZVJbPyIhIzskIjEvYFZQclYlUiQhIzo3JCQiMjtJZz9KPiFlNyEjOyQiMk0pKlxTWWknPlEhIzs3JCQiMmAxOEVVN1ZJIiEjOyQiMXhHNCczRzI9JSEjOjckJCIybEReLUQ1Vk4iISM7JCIxVDdcSlxUaFghIzo3JCQiMkNZI1wpKkcnZlMiISM7JCIxaW8jPSVRJFslXCEjOjckJCIyQVcpb1BcXGM5ISM7JCIxLCMpUiQpSFo1YCEjOjckJCIxKWY+UkcoUTA6ISM6JCIwYDVnKj5jYmMhIzk3JCQiMkJZI1wpKlFuZjohIzskIjFhaioqeSVwKkdnISM6NyQkIjJYITQ9T0dYMzshIzskIjFHLXgjKWYmZk4nISM6NyQkIjJkND5Rd00wbSIhIzskIjE0JypleDxSJ3AnISM6NyQkIjF6ZDpKKEd4cSIhIzokIjE6UDBCT0ooKnAhIzo3JCQiMmIyOkk/QyRmPCEjOyQiMWMkUiVSbEY9dCEjOjckJCIyd18wNl9yeSE9ISM7JCIxO0g+WSg0R2goISM6NyQkIjIpPlJ5Y3VoZT0hIzskIjE0OTltVDg4eiEjOjckJCIyMDZBVykzQjM+ISM7JCIxN2B4aDBbKj4pISM6NyQkIjIxOEVfYXIsJz4hIzskIjEzI3oqKUhXPFwpISM6NyQkIjFzVihbZCc+NT8hIzokIjE3UTUlcFZodykhIzo3JCQiMlcjWydIUmE4MSMhIzskIjEkZjhbdlooUiEqISM6NyQkIjJGWSNcKWYpMzdAISM7JCIxJD4yJypmaFVJKiEjOjckJCIyJDM7S2shMyhlQCEjOyQiMU0iNG4weDlhKiEjOjckJCIyJUhlO0wjUkBAIyEjOyQiMTosYkonUW0hKSohIzo3JCQiMiUpb1B2cUcqZkEhIzskIjJfKFxBd3V5LjUhIzo3JCQiMiVwUXhhQSkzSiMhIzskIjJBKW9NIyk0JXktIiEjOjckJCIyYC0wNUleJ2ZCISM7JCIyTWEhPlRqSF01ISM6NyQkIjJEWydIZkkhUVQjISM7JCIxczpVLnBmdTUhIzk3JCQiMmIyOkldJXBnQyEjOyQiMkBFZ11OOV40IiEjOjckJCIybEZiNXI7Ul4jISM7JCIwKG9nIVFFeTYiISM4NyQkIjEoUnplKEhVaUQhIzokIjJDLDAqSGcrUTYhIzo3JCQiMk1tS2w/Ils6RSEjOyQiMm03eGwzRSZmNiEjOjckJCIxaEBWJ1toN20jISM6JCIyTyZwKSlINGt4NiEjOjckJCIyYC4yOSkpKlI4RiEjOyQiMWNPTV5HeCg+IiEjOTckJCIyWiNcKXAqSHdqRiEjOyQiMmtqLmZFRm5AIiEjOjckJCIyIyp5ZDpAJDQ5RyEjOyQiMillISkqXFQpPk43ISM6NyQkIjFmPE5xI1FVJ0chIzokIjJPcTcmKjRYSkQiISM6NyQkIjJWJG9PdD5UN0ghIzskIjJrOVJaZm4qcDchIzo3JCQiMWAwNkEhKltrSCEjOiQiMi1LLFw2LnhHIiEjOjckJCIxVCJHY0FBVCwkISM6JCIyJltNRz4+PS84ISM6NyQkIjIoemU8TillajEkISM7JCIxRyFbQS0qM0A4ISM5NyQkIjFWJ0dkYVhPNiQhIzokIjFwI1tCJXksTzghIzk3JCQiMUZgMThpIWY7JCEjOiQiMkxxMz9fNkBOIiEjOjckJCIyMjlHYyNHJ2ZAJCEjOyQiMkA+KSlSN1ByTyIhIzo3JCQiMXZcKiopZigqZUUkISM6JCIwRzNDNWM8USIhIzg3JCQiMk5lO0xtYSE9TCEjOyQiMkAyblAjNGsnUiIhIzo3JCQiMjA0PU9BJDRtTCEjOyQiMiUqeiFRKT02K1QiISM6NyQkIjIvLC0vRyJIOk0hIzskIjJkQytDTHZMVSIhIzo3JCQiMUJYITR5NidwTSEjOiQiMiczVTFXU3ZQOSEjOjckJCIxPU9zV1h6PU4hIzokIjJRIylHKXAsV105ISM6NyQkIjImSGU7TEo0cE4hIzskIjIiKT0hUVBhNGo5ISM6NyQkIjIxMDU/XVgtaSQhIzskIjJpXiYqUkBWY1oiISM6NyQkIjF0WCJIW2tzbSQhIzokIjJNa2JwUSgqb1siISM6NyQkIjIoXClwUj4mUTxQISM7JCIyLzoqZiE9MScpXCIhIzo3JCQiMXRYIkhRQ3J3JCEjOiQiMGM+IVskUio0OiEjODckJCIydE5yVSYpSC4jUSEjOyQiMiMpXFItTWE8XyIhIzo3JCQiMVolKil5UCJHblEhIzokIjJgLVMkcEUjPmAiISM6NyQkIjJqQlolKm9gOSNSISM7JCIyQyRvRjo/T1Y6ISM6NyQkIjEuMTdDK01xUiEjOiQiMmRkbTZeQU1iIiEjOjckJCIxRFwpcFw1KD1TISM6JCIyLjYuenNQSmMiISM6NyQkIjFoQlolSEMyMiUhIzokIjFaJEgsJGVLdDohIzk3JCQiMkRZI1wpXFtINyUhIzskIjIqUiZRLC0jSCRlIiEjOjckJCIxJ0dkOXprLzwlISM6JCIyaWJuVFVPQGYiISM6NyQkIjEpZj5SKVwjNEElISM6JCIyKHpNKEdiLDhnIiEjOjckJCIxLCsrKyNwLEYlISM6JCIyOC1xOmlEK2giISM6NyQkIjFXKFsoXC9iQlYhIzokIjIyVDEiKjRWIz47ISM6NyQkIjE0PU9zK3VwViEjOiQiMi9TJj1PRS1GOyEjOjckJCIwJGU7TDB6QVchIzkkIi8xJnAjKVFkaiIhIzc3JCQiMSNRd192Z0VaJSEjOiQiMlc2NjxNQ1BrIiEjOjckJCIxT3JVJkdcP18lISM6JCIyLkBdPyopUTlsIiEjOjckJCIyTmhBWDV4OWQlISM7JCIyXy8peiQ9ciplOyEjOjckJCIyOUdjN04hKj1pJSEjOyQiMWAlZiRwVVltOyEjOTckJCIyYz9UI1t2S3ZZISM7JCIxcip6V0cuVW4iISM5NyQkIjInR2M3RFVLQ1ohIzskIjJ1JCk+QHc+Nm8iISM6NyQkIjE8TW9PLXBzWiEjOiQiMih6MyczLiR5KG8iISM6NyQkIjF3Xi4yRlJDWyEjOiQiMjdRcV0qKkhacCIhIzo3JCQiMVQiR2NpIylmKFshIzokIjFZJVFgNyZbLDwhIzk3JCQiMC4xN2soR0FcISM5JCIwTV53aix1cSIhIzg3JCQiMSV5YzhGUnIoXCEjOiQiMjopMzA8X0I5PCEjOjckJCIxJ0hmPWRPSS0mISM6JCIydzxoVyc0Iik+PCEjOjckJCIxLTA1Pyc9bjImISM6JCIyWTdKQDhzaHMiISM6NyQkIjFuTG5NSCZ6NyYhIzokIjJhJ3oiXDEnM0s8ISM6NyQkIjFzVihbMipwdF4hIzokIjBzLicpKSpSc3QiISM4NyQkIjAmKXBSUlheQSYhIzkkIjJNa2ooNHEqR3UiISM6NyQkIjE8TG1LRiRwRiYhIzokIjFMPjdYa1dbPCEjOTckJCIwKWU8TlJaR2AhIzkkIjFjJEciZSdHUXYiISM5NyQkIjEiPlF3QWZpUCYhIzokIjJPVDhjRCdwZTwhIzo3JCQiMWEyOkkiUWRVJiEjOiQiMmk8XTtJOk93IiEjOjckJCIxWiY0Pkc0cFomISM6JCIyRGNDSj13Jm88ISM6NyQkIjF6YzhGaiJ6XyYhIzokIjI3NHduVidSdDwhIzo3JCQiMW9OclVMUSFlJiEjOiQiMmIkcGldIUcjeTwhIzo3JCQiMWxIZj1qZkVjISM6JCIydVNwJzMrUSN5IiEjOjckJCIxKWY+Uj1AJ3ljISM6JCIyUEFKOXdTcHkiISM6NyQkIjE0PE1vJ2YzdCYhIzokIjIteUdrVC05eiIhIzo3JCQiMSM+UXc3LDd5JiEjOiQiMWNVNnhKZiZ6IiEjOTckJCIxKFF4YSpmInAjZSEjOiQiMTsiKlFzJTMkKnoiISM5NyQkIjF2WyhcUnY3KWUhIzokIjFsc2wqPjxPIT0hIzk3JCQiMTxNb09cS0ZmISM6JCIyWj1GPHV3ciE9ISM6NyQkIjFFXSssOiozKWYhIzokIjI4M1BHXzo3Ij0hIzo3JCQiMXRZJHApcElHZyEjOiQiMUpPeCNRLFoiPSEjOTckJCIyJlsoXCoqR0guMychIzskIjFBKHB1LkslPT0hIzk3JCQiMTpKaUNvJylIaCEjOiQiMltTQkJWJio9Iz0hIzo3JCQiMUpoQVhTYiI9JyEjOiQiMiNlVUBHKT1hIz0hIzo3JCQiMT5QdVshPiFIaSEjOiQiMmwlXChSPHYmRz0hIzo3JCQiMSVvT3RZOy1HJyEjOiQiMkFgVFN5Jyo9JD0hIzo3JCQiMTdDWyc+J1JMaiEjOiQiMk41MlQ5ZV8kPSEjOjckJCIxeF4uMiQqb3pqISM6JCIxQStuR0I2UT0hIzk3JCQiMW5Mbk1yb0hrISM6JCIyOFl1b2hANiU9ISM6NyQkIjJEZDlIeVI4WychIzskIjJCeGI8XV9UJT0hIzo3JCQiMWEwNkE9KD1gJyEjOiQiMiQqRyVcUkcvWj0hIzo3JCQiMTQ8TW9UdyFlJyEjOiQiMmNHJHBJMnhcPSEjOjckJCIxdVgiSHldXWonISM6JCIyJlJ4S3dDc189ISM6NyQkIjE6SWc/KEhRbychIzokIjF3J2UicHNJYj0hIzk3JCQiMTE3Q1s7ImZ0JyEjOiQiMVRmPHgkKSp6Jj0hIzk3JCQiMCp5ZDpjNSR5JyEjOSQiMkJBNjw4eC4nPSEjOjckJCIxJj1QdTMsWiRvISM6JCIyZFpQKDNWIkgnPSEjOjckJCIxUndfMCVbSylvISM6JCIyTClweGhEQ2w9ISM6NyQkIjFKZz9UVipSJHAhIzokIjBSR1BmO3cnPSEjODckJCIxQFUlKW94ZyQpcCEjOiQiMiVcViI0OCEpKXA9ISM6NyQkIjFVI1snSCVbYi4oISM6JCIydUtWU1chPnM9ISM6NyQkIjEkWydIZk1kJjMoISM6JCIxTG1bVCZmVig9ISM5NyQkIjFNcFF4N3RPciEjOiQiMjxwKCo9TkFsKD0hIzo3JCQiMXVYIkhbbHU9KCEjOiQiMllKUV5JOCd5PSEjOjckJCIxPVB1W1wzTXMhIzokIjIyUidlbiUpWyEpPSEjOjckJCIyJlJ6ZTxoXihHKCEjOyQiMWBGMVRYZSMpPSEjOTckJCIvKWY+ZjBgTCghIzgkIjJ0ViVlUENUJSk9ISM6NyQkIjApZj5SImZpUSghIzkkIjFbRzlUUUonKT0hIzk3JCQiMVByVSY9R11WKCEjOiQiMnMxJCpHIiopMykpPSEjOjckJCIxJGY9UCUqeiIqWyghIzokIjFCKW9mJSk0Kyo9ISM5NyQkIjEnPVB1UXJnYCghIzokIjE9dERHPGoiKj0hIzk3JCQiMShReGFmJEgqZSghIzokIjJjKj1CKjNGTSo9ISM6NyQkIjEzOklnKSp6UHchIzokIjIvIz5ob0EtJio9ISM6NyQkIjF1WiY0NGUzcCghIzokIjIkeS51ME1zJyo9ISM6NyQkIjFzVSYzUFFtdCghIzokIjJLYlpXT2IiKSo9ISM6NyQkIjFZIkhld3coKXkoISM6JCIxLU4qNHhZKCoqPSEjOTckJCIxT3FTIikpUiJSeSEjOiQiMSIpcCpvNFg3IT4hIzk3JCQiLyp6ZjRxJSopeSEjOCQiLykqPUBfcS0+ISM3NyQkIjFwUXhhXmhSeiEjOiQiMmNvRXIpUTcvPiEjOjckJCIxWSopeWQpKXkoKXohIzokIi8tNyhvYGEhPiEjNzckJCIxbEVgMWYnKVIhKSEjOiQiMWVUdWdjJm8hPiEjOTckJCIxXy0wNSIqXCozKSEjOiQiMilvJFxCSGUiMz4hIzo3JCQiMCp6Zj5kdFQiKSEjOSQiMXZGXCF5JVw0PiEjOTckJCIxYTI6SUMtKj0pISM6JCIxJlJBdilbbjU+ISM5NyQkIjFRdVsoNCRHVCMpISM6JCIxY2lWdnElPiI+ISM5NyQkIjFfLTA1KFI4SCkhIzokIjJCVG9IJFs4OD4hIzo3JCQiMSczPE1bdTdNKSEjOiQiMmlVXHRXIUg5PiEjOjckJCIxJHBReGFKTVIpISM6JCIxeWlLd3FZOj4hIzk3JCQiMSwtLzMsWlQlKSEjOiQiMSMpZUIiKVJfOz4hIzk3JCQiMUFUI1s7bzFcKSEjOiQiMidcIj5CUyFlPD4hIzo3JCQiMUxtS2wnKSlcYSkhIzokIjJpIVEoUS08KD0+ISM6NyQkIjB0WCJIOTwlZikhIzkkIjIpWyhcJFspPig+PiEjOjckJCIwJXplPCtaVycpISM5JCIyJT11Xkwtcz8+ISM6NyQkIjFoQFYnUUFjcCkhIzokIjFMS2tAQHJAPiEjOTckJCIxJm9PdE9URXUpISM6JCIxKHouPXYsRSM+ISM5NyQkIjFoPlJ5P3cjeikhIzokIjFzZ0I+dF9CPiEjOTckJCIxJm9PdEUsRCUpKSEjOiQiMS9JI1E+QlcjPiEjOTckJCIxb01wUW5xJiopKSEjOiQiMiQ0djhjck5EPiEjOjckJCIxZDpKaSNlRSUqKSEjOiQiMUMjPjomNDtFPiEjOTckJCIxWyRwUWRJbyoqKSEjOiQiMkdtLSU+XzFGPiEjOjckJCIxOUZhM3ByWCEqISM6JCIxdmdvIVxneSM+ISM5NyQkIjFPcVMiUSgzJTQqISM6JCIyJyozRiFvJUcnRz4hIzo3JCQiMXVXKil5NjVZIiohIzokIjJiJEhYRVFWSD4hIzo3JCQiMXRYIkhRRCQpPiohIzokIjFBKWVIa0AtJD4hIzk3JCQiMShSemVuVGVDKiEjOiQiMiRHLV0keT80JD4hIzo3JCQiMHJUJG89SSdIKiEjOSQiMUQkenBGWDskPiEjOTckJCIwNkFXM1liTSohIzkkIjJkRlZeIVxMSz4hIzo3JCQiMWIzPE10IyopUiohIzokIjE8MXVPTjFMPiEjOTckJCIwI1J5Y3A2WCUqISM5JCIxOV5xKyZ5TyQ+ISM5NyQkIjAlemU8dTspXCohIzkkIjA5RmhdblYkPiEjODckJCIxJVwpcFJ3LlsmKiEjOiQiMlkmUVNfKCkqXCQ+ISM6NyQkIjFbI1wpcGhVKGYqISM6JCIyMTNQI3omM2MkPiEjOjckJCIxQ1olKikpUiZvayohIzokIjFPa3MpKlI/Tz4hIzk3JCQiMSNSeWNCbnNwKiEjOiQiMndINCE+anpPPiEjOjckJCIxPExtS1dxXSgqISM6JCIxW3B2eiEzdSQ+ISM5NyQkIjFSeGE0NnEqeiohIzokIjItSDQpPlsmeiQ+ISM6NyQkIjFHYjVAcjFbKSohIzokIjIkUUgoUmEiW1E+ISM6NyQkIjEpR2Q5ZnAoKiopKiEjOiQiMiZIYiZScEkhUj4hIzo3JCQiMV4tMDUmZjgmKiohIzokIjImNEtuelljUj4hIzo3JCQiJCsiISIiJCIydHEjKW9uYislPiEjOi0lJkNPTE9SRzYnJSRSR0JHJCIpQyllcSUhIikkIiIhISIiJCIoJz4hXCYhIiktJStfQVRUUklCVVRFRzYjLyUnc291cmNlRyUsbWF0aGRlZmF1bHRHLSUlVklFV0c2JDskIiIhISIiJCIkKyIhIiIlKERFRkFVTFRHLSYlJl9BWElTRzYjIiIiNictJStfR1JJRExJTkVTRzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSYlJl9BWElTRzYjIiIjNictJStfR1JJRExJTkVTRzYmLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUmQ09MT1JHNiYlJFJHQkckIiIhISIiJCIiISEiIiQiIiEhIiItJSpMSU5FU1RZTEVHNiMiIiEtJSpUSElDS05FU1NHNiMiIiEtJS1UUkFOU1BBUkVOQ1lHNiMkIiIhISIiLSUrQVhFU0xBQkVMU0c2JC1JI21pRzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliRzYiNjVRInQ2Ii8lJ2ZhbWlseUdRITYiLyUlc2l6ZUdRIzEwNiIvJSVib2xkR1EmZmFsc2U2Ii8lJ2l0YWxpY0dRJXRydWU2Ii8lKnVuZGVybGluZUdRJmZhbHNlNiIvJSpzdWJzY3JpcHRHUSZmYWxzZTYiLyUsc3VwZXJzY3JpcHRHUSZmYWxzZTYiLyUrZm9yZWdyb3VuZEdRKFswLDAsMF02Ii8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdNiIvJSdvcGFxdWVHUSZmYWxzZTYiLyUrZXhlY3V0YWJsZUdRJmZhbHNlNiIvJSlyZWFkb25seUdRJmZhbHNlNiIvJSljb21wb3NlZEdRJmZhbHNlNiIvJSpjb252ZXJ0ZWRHUSZmYWxzZTYiLyUraW1zZWxlY3RlZEdRJmZhbHNlNiIvJSxwbGFjZWhvbGRlckdRJmZhbHNlNiIvJTZzZWxlY3Rpb24tcGxhY2Vob2xkZXJHUSZmYWxzZTYiLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWM2IlEhNiItJSVST09URzYnLSUpQk9VTkRTX1hHNiMkIiQrJiEiIi0lKUJPVU5EU19ZRzYjJCIjXSEiIi0lLUJPVU5EU19XSURUSEc2IyQiJXFMISIiLSUuQk9VTkRTX0hFSUdIVEc2IyQiJT9RISIiLSUpQ0hJTERSRU5HNiI=

# Tanks! Eqn4A:=diff(Q(t),t)=9-(3/100)*Q(t);

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJRRzYiNiNJInRHRilGKywmIiIqIiIiRicjISIkIiQrIg==

dsolve({Eqn4A,Q(0)=0},Q(t));

Ly1JIlFHNiI2I0kidEdGJSwmIiQrJCIiIi1JJGV4cEc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkdGJTYjLCRGJyMhIiQiJCsiISQrJA==

# This looks reasonable, since Q(0)=0 and as t-> infinity, Q -> 300. Eqn4B:=diff(Q(t),t)=9-(4/(100-t))*Q(t);

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJRRzYiNiNJInRHRilGKywmIiIqIiIiKiYsJiIkKyJGLkYrISIiRjJGJ0YuISIl

dsolve({Eqn4B,Q(0)=0},Q(t));

Ly1JIlFHNiI2I0kidEdGJSwoIiQrJCIiIkYnISIkKiQsJiEkKyJGKkYnRioiIiUjRisiKCsrKyI=

# Note that t is only valid for 0<=t <=100. Eqn5:=diff(u(t),t)=k*(R-u(t)); #R is the constant environmental temp

Ly1JJWRpZmZHJSpwcm90ZWN0ZWRHNiQtSSJ1RzYiNiNJInRHRilGKyomSSJrR0YpIiIiLCZJIlJHRilGLkYnISIiRi4=

dsolve(Eqn5,u(t));

Ly1JInVHNiI2I0kidEdGJSwmSSJSR0YlIiIiKiYtSSRleHBHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2IywkKiZJImtHRiVGKkYnRiohIiJGKkkkX0MxR0YlRipGKg==