Here is the code to run this in GAMS. GAMS is optimization software that is available free online:
http://www.gams.com/download/
If you download it, copy and paste this into a new program.
Variables
I total interest
A Amount invested in A
B amount invested in B
C amount invested in C
D amount invested in D
Positive variable A;
Positive variable B;
Positive variable C;
Positive variable D;
Equations
interest shows total interest
constraintA total invested equals 100000
constraintB amount invested in A less than or equal to 50000
constraintC amount invested in A must be greater than 20000
constraintD amount invested in D must be greater than C
points total points accumulated;
interest.. I =e= .04*A+.06*B+.09*C+.11*D;
constraintA.. A+B+C+D =e= 100000;
constraintB.. A =l= 50000;
constraintC.. A =g= 20000;
constraintD.. D =g= C;
points.. A/500+3*B/1000+7*C/1000+D/125 =l= 1000;
Model interest1 /all/;
solve interest1 using lp maximizing I;
display A.l, B.l, C.l, D.l, interest.l, points.l;