Wednesday, September 28, 2005

A good Symbolic toolbox example

Problem No 3.23, Engineering Optimization Reklaitis

clc;
close all;
syms x1 x2
a='x1*x2+(2000/(x1*x2))*(x2+x1)';
m=diff(a,x1);
n=diff(a,x2);
pretty(m)
pretty(n)
l=solve(m,x1)
l1=simplify(subs(n,x1,l(1)));
l2=simplify(subs(n,x1,l(2)));
pretty(l1)
pretty(l2)
simplify(solve(l1,x2))
simplify(solve(l2,x2))
simplify(subs(l(1),x2,solve(l1,x2)))
h=1000/(simplify(subs(l(1),x2,solve(l1,x2)))*simplify(subs(l(1),x2,solve(l1,x2))))

0 Comments:

Post a Comment

<< Home