Hi,

if I have a equation

( 15,3440 * r_1 - 0.6329 + 5.7540 * r_2  )   / 

( 21.0981 * r_1 - 3.1780 + 84.8348 *r_2 )

How do I solve for r_1 and r_2 ?

The solutions should be r_1=0.175 and r_2=0.055

Why can I only get the solution  -(-y + 2 x) ^3  (y + 2 x)  from 2) and not from 1).

I am doing the same thing in both of them but from 1) I get "Warning, solutions may have been lost" when Itry to solve for _F1

I want to go from A-B-C with commands (as in 1) ) without having to retype the expression

1)restart;

I have two differrential equations:

 

1) x ' ( t ) = x( t ) + y( t )

which can be written as

x ( t ) - x (t -1) = x ( t) + y (t )

- x (t -1) = x ( t) + y (t ) - x ( t )

- x (t -1) =  y (t )           A

 

2) y ' ( t ) = 4 * x( t )  + y (t )

which can be written as

y( t ) - y (t -1) = 4 * x( t )  + y (t )

y( t ) = 4 * x( t )  + y (t ) +  y (t -1)        B

When I want to simulate and plot the distribution from a random coin toss (normal distribution) I use the following code:

restart:
randomize():
coin:=rand(0..1):
coin_1:=proc(n) seq(coin(),i=1..n) end:
x_1:=seq([coin_1(10)],i=1..500):
x_2:=[seq(numboccur(1,x_1[i]),i=1..500)]:
with(Statistics):
Histogram(x_2,frequencyscale=relative,discrete = true );

what code should I use in order to get a coin with serial correlation. A normal distribution with fat tails? (Prefereably I want to tweak the above code)

I want to be able to do the following in Maple 1) Set up the current value hamiltonian H:=u(c(t))+lambda(t)*(f(k)-c(t)-delta*k) ; 2) Take the foc with respect to c(t) which should give me u'(c(t))=lambda(t) 3) Differentiate with respect to time which should give me u''(c(t))*diff(c(t),t)=diff(lambda(t), t) 4) Divide both sides by u'(c(t))=lambda(t) should give me u''(c(t))*diff(c(t),t) / u'(c(t)) =diff(lambda(t), t) / lambda(t) I have tried to set up the hamiltonian in maple but I find it quite hard for numerous reasons: A) I cannot differentiate a function w.r.t a function i.e