Hi All, I want to convert matlab code to julia. is there a way to access elements of symbolic expression.
function [xqq,xq,wq]=SetIntoptimal1D(n)
clc
#clear all;
syms y
# Input n: Quad pt rule
m=n-1;
P=sym('a',[1,n]);
P(1)=1;
P(2)=y;
for i=1:m
P(i+2)=((2.0*i+1)*y*P(i+1)-(i)*P(i))/(i+1.0);
# P(i)=P(i+1);
# P(i+1)=P(i+2);
end
if n==1
P(n)=P(1);
end
p=n-1;
t=p/(p+1);
#Pn=expand(P(n+1));
Pn=expand(P(n+1)-t*P(n-1));
xqq=solve(vpa(Pn,32));
xqq=sort(xqq);
#update here
xq=zeros(n);
parfor i=1:n
xq(i)=(xqq(i)+1)/2;
end
xq=sort(xq);
for k=1:1:n
P(1)=1;
P(2)=y;
m=n;
for i=1:m
P(i+2)=((2.0*i+1)*y*P(i+1)-(i)*P(i))/(i+1.0);
end
P1=((n+1)/(1-y*y))*(y*P(n+1)-P(n+2));
P2=((n-1)/(1-y*y))*(y*P(n-1)-P(n));
Pn=P(n)*(P1-t*P2);
#Pn=x*P(n+1)-P(n+2);
wqq(k)=vpa(2*(p*(1+t)+t)/(p*(p+1))/...
subs(Pn,y,xqq(k)),32);
end
parfor i=1:n
wq(i)=wqq(i)/2;
end
# disp('Nodes');
# disp(xq);
# disp('Weights')
# disp(wq)
end
Julia code converted so far
using Symbolics
@variables y
# Input n: Quad pt rule
n=10
m=n-1;
@variables P[1:n]
Symbolics.scalarize(P)