unexpected matlab operator on okada script

i tried to run Okada script from Okada
and run with the input example
[E,N] = meshgrid(linspace(-10,10,50));
[uE,uN,uZ] = okada85(E,N,2,30,70,5,3,-45,1,1,'plot');
figure, surf(E,N,uN)
and the function is
function u=uzy_ss(xi,eta,q,dip,nu)
R = sqrt(xi.^2 + eta.^2 + q.^2);
db = eta.*sin(dip) - q.*cos(dip);
yb = eta.*cos(dip) + q.*sin(dip);
u = (db.*q./R.^3).*cos(dip) ...
+ (xi.^2.*q.*A(eta,R).*cos(dip) - sin(dip)./R + yb.*q./R.^3 ...
- K2(xi,eta,q,dip,nu,R)).*sin(dip);
but then i found an error says : Error: File: okada85.m Line: 465 Column: 1 Unexpected MATLAB operator.

5 Comments

Amri - I was able to run the above code without any problems. I recommend that you replace your copy of the okada85.m file with the original and try again. Line 465 seems to correspond to
yb = eta.*cos(dip) + q.*sin(dip);
and there doesn't seem to be anything suspicious about this line.
actually line 465 is
function u=uzy_ss(xi,eta,q,dip,nu)
it is because is input 3 line of code as the input in the first line of the script
I get no error. Are you sure you have the right function without any edits? Maybe there's a function overriding matlab functions?
[E,N] = meshgrid(linspace(-10,10,50));
[uE,uN,uZ] = okada85(E,N,2,30,70,5,3,-45,1,1,'plot');
figure, surf(E,N,uN)
%It works. I get a plot.
LINE 465 of okada85.m
yb = eta.*cos(dip) + q.*sin(dip);
Amri - if I open the okada85.m (downloaded as is from the FEX) then
Note line 465. Try replacing your copy of the file with that from the FEX.
thank you ! i got this script works now. I just redownload it from the FEX

Sign in to comment.

Answers (0)

Categories

Tags

Asked:

on 5 Jul 2018

Commented:

on 6 Jul 2018

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!